курсачи / курсач 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
% двухмерный случий а не трехмерный