Решение СЛАУ по методу Якоби

 function X=yakobi(A,B,X0,delta,Imax)
 
% Итерация Якоби.
% A     - невырожденная матрица коэффициентов
%         размера n x n;
 
 
                                   
%   B       -   вектор-столбец свободных членов;
%   X0      -   вектор-столбец начального приближения;
%   delta   -   допустимая относительная погрешность;
%   Imax    -   максимальное число итераций;
%   X       -   приближенное решение линейной системы AX = B.
 
n=length(B);
err=5*delta;
X=X0;
ct=0;
while err>delta
    Xp=X;
    for i=1:n
        X(i)=(B(i)-A(i,[1:i-1,i+1:n])*Xp([1:i-1,i+1:n]))/A(i,i); 
              
    end
    if max(abs(X0))==0
        error('Следует изменить начальное приближение')
        break
    end
    err=max(abs(X-Xp))/max(abs(X0));
    ct=ct+1;
    Xe(ct,:)=[X' err];
    if ct>Imax
        error('Требуемое число итераций слишком велико')
        break
    end
end
Xe
ct
и программа:
 A=[5.96 1.40 5.03; 2.99 1.25 -0.38; 1.05 -3.57 1.92];
 b=[-7.41; 4.48; -15.14];
X0=[1;1;1];
 X=yakobi1(A,B,X0,1e-2,100);