Matlab - metodo iterativo de jacobi

%programa del metodo iterativo de jacobi 2

clear all;

n=input('ingrese el orden de la matriz: \n');

A=input('ingrese la amtriz A= \n');

b=input('ingrese el vector de terminos indepemndientes: \n');

X0=input('ingrese la aproximacion inicial X0= \n');

tol=input('ingrese la tolerancia tol= \n');

N=input('ingrese el numero de iteraciones N= \n');

k=1;

while k<=N

    for i=1:n

        suma1=0;

        for j=1:n

            if i~=j

                suma1=suma1+A(i,j)*X0(j);

            end

            x(i)=(b(i)-suma1)/A(i,i);

        end

        norma=norm(x-X0);

        if norma<tol

           fprintf('%10.6f\n',x(i));

           fprintf('prosedimiento terminado exitosamente')

        end

    end

    k=k+1;

    for i=1:n

        X0(i)=x(i);

    end

end

    if k>N

        break

    end

disp('numero de iteraciones excedido');

fprintf('procedimiento terminado sin exito');

fprintf('\nTABLA:\n\t   norma   x(i)\n\n  ');