Jacobi

function[x] = jacobi(A,b)

eps = 10^(-6)

n=size(A,1);
L=zeros(n);
D=zeros(n);
U=zeros(n);
x=zeros(n,1);

for i=n-1:n-1
    for j=2:n
        L(i,j) = A(i,j)
    end
end

for i=1:n
    for j=1:n
        D(i,j) = A(i,j)
    end
end

for i=2:n
    for j=1:i-1
        U(i,j) = A(i,j)
    end
end

MJ= -D^(-1) * (L+U)

p=abs(eigs(MJ,1,'largestabs'))
if p < 1
    disp('Zbiezne')
end

x0=x;
k=0;
do=0;

while do==0
    x= -D^(-1)*(L+U)*x0+D^(-1)*b
    if (norm(x-x0) < eps)
        do=do+1;
    else
        k = k + 1;
    end
end