курсачи / курсач 3 вар / kirs

sist= zeros(40,40);
ro=4;
mut=zeros(25);
h=1;
a=-4;
x0=4;
y0=4;
n=(y0-1)*5+x0;
fitrub=5;
y=1;
razrsetki=5;
%для сетки 5 на 5
%не понятно почему но матрица отображается как Y X
%составление матрицы
x=1;
for(y=1:25)
sist(y,x)=a;
if(x==n)
sist(y,n)=fitrub;
sist(n,26)=ro-a*fitrub;
sist(y,x)=0;
end
if(x-razrsetki)>0
sist(y,x-razrsetki)=1/h;
if(x-razrsetki==n)
sist(y,26)=-fitrub;
sist(y,x-razrsetki)=0;
end

end

if(x+razrsetki)<26
sist(y,x+razrsetki)=1/h;
if (x+razrsetki==n)
sist(y,26)=-fitrub;
sist(y,x+razrsetki)=0;
end
end

if(x+1)<26
sist(y,x+1)=1/h;
if(x+1==n)
sist(y,26)=-fitrub;
sist(y,x+1)=0;
end
end

if(x-1)>0
sist(y,x-1)=1/h;
if(x-1==n)
sist(y,26)=-fitrub;
sist(y,x-1)=0;
end
end
x=x+1;
end


%сортировка матрицы методом Гаусса

for(i=1:25)
for (y=i:26)
if( sist(y+1,i)==0)
while(sist(y+1,i)==0)&&(y<25)
y=y+1;

end

end
for(x=1:26)
end
for(x=1:26)

sist(y+1,x)=sist(i,i)*sist(y+1,x)/sist(y+1,i);
sist(y+1,x)=sist(y+1,x)-sist(i,x);
end
end
end
i=1;
% МАТРИЦА ПАТЕНИЦАЛОВ
% mat=zeros(25);
% mat2;
% for (y=26:-1:1)
% i=25;
% for(x=1:25-i)
% if(x~=i)
% mat2=mat2+sist(y,26)-sist(y,x);
% end
% mat2=mat2/sisr(y,x);
% a=div(y)+1;
% b= mod(y);
% if b==0
% b=razrsetki;
% end
% if y==razrsetki*razrsetki;
% mat(div(y)+1)
% end
% end
% двухмерный случий а не трехмерный