Matlab Metode Jacobi

disp(' ');
disp('=====================================================')
disp('- Program Menentukan Solusi Sistem Persamaan Linear -')
disp('-             dengan Metode Jacobi                  -')
disp('-----------------------------------------------------')

function [X,g,H]= jacobi(A,b,X0,T,N)
H = X0';
n = length(b);
X1 = X0;
for k=1:N,
    for i = 1:n,
        S = b(i)-A(i,[1:i-1,i+1:n])*X0([1:i-1,i+1:n]);
        X1(i)=S/A(i,i);
    end
    g = abs(X1-X0)
    err = norm(g);
    relerr = err/(norm(X1)+ eps);
    X0 = X1
    H = [H;X0']
if (err<T)|(relerr<T),break,end
end
disp(' ')
disp('===================Terima Kasih======================')