zeydel

%Исходные данные
A = [1,-0.08,0.23,-0.32;
    -0.16,1.23,-0.18,-0.16;
    -0.15,-0.12,0.68,0.18;
    -0.25,-0.21,0.16,0.97;];

B = [1.34,-2.33,0.34,0.63];


%размер системы
N = max(size(B));

%требуемая точность
e = 0.0001;


X= [zeros(1,N);
    zeros(1,N)];

root = zeros(1,N);
%Начальные приближения
for i = 1:N
    root(i) = B(1,i) / A(i,i);
end

m = max(abs(root));

%Итерации

while m>e
    for i = 1:N
        X (1,i) = root(i);
        ss = B(1,i);
        for j = 1:N
            if(i ~= j)
                ss = ss - A(i,j)*root(j);
            end;
        end;
        root(i) = ss/A(i,i);
        X(2,i) = root(i) - X(1,i);
    end;
    m = max(abs(X(2)));
end;

disp (root)