Jacobi Method

%Metodo de Jacobi

clear;
clc;

format long

disp('Método de Jacobi')
disp('')

A=input('Introduza a matriz A= ');

[m,n]=size(A);

b=input('Introduza o vector b(vector linha 1*n)=');
x0=input('Introduza a aproximação inicial x0 (vector linha 1*n)= ');
erro=input('Qual o erro max pretendido?  ');
iter=input('Qual o numero max de iteraçoes? ');
disp('')
it=1;
e=erro+1;
x=x0;
while (e>erro && it<=iter)

    for i=1:m;
        soma=0;
        for j=1:n
            if j~=i
                soma=soma+A(i,j)*x(j);
            end
        end
        x(i)=(1/A(i,i))*(b(i)-soma);
    end

    e=norm(x-x0, inf);

    disp([num2str(it),'ª iteração: Solução aproximada = ', mat2str(x,5), 'com estimativa de erro = ', num2str(e)])

    it=it+1;
    x0=x;

end