function [x]= demoJacobi(A,b)
//Se determina la dimensión de la matriz A
[n m] = size (A); 
if (n<>m) then
    error('La matriz debe ser cuadrada');
end
// Se inicializa el vector  inicial con x=[0, 0, 0]
x = zeros(n,1);
//Se inicializa la cantidad maxima de tolerancia, para cuando el método no converga
iteraciones=1000; 
//Se inicializa la tolerancia al error
tolerancia=0.0001;
//Se inicializa contador de iteraciones
it=0;
//Se inicializa condicion de salida
salir=0;
while (salir==0)
    //Se inicia el método de Jacobi
    for i=1:n
        s=0
        for j=1:n
            if i<>j then
                s=s+A(i,j)*x(j);
            end
        end
        y(i)=(1/A(i,i))*(b(i)-s);
    end
    error = abs(x(1)-y(1));
    for i=2:n
        if error<abs(x(i)-y(i)) then
            error=abs(x(i)-y(i));
        end
    end
    it=it+1;
    if((it==iteraciones)|(error<tolerancia)) then
       if (it==iteraciones) then
         disp('El método no converge');
       end
       salir=1;
    end
    x=y;
end
endfunction