programa jacobi

el codigo del programa es:

function jacobi(a,b,x0,Es)
D=diag(diag(a));
L=(tril(a)-D)*(-1);
U=(triu(a)-D)*(-1);
m=inv(D)*(L+U);
c=inv(D)*b;
E=99999999999;
cont=0;
n=length(m);
while E>=Es
x=m*x0+c;
cont=cont+1;
E=max(abs(x-x0));
x0=x;
end
disp('numero de iteraciones')
disp(cont)
disp('valor del error')
disp(E)
disp('vector x')
disp(x0)


para ejecuytar se tiene los siguiente
A es la matriz, B es el vector, xo es el vector nulo, Es es el error deseado