Método de Jacobi em SCILAB (Cálculo Numérico)

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

A=[225 0 -25 0; 0 175 0 -125; -225 0 275 -50; 0 -25 -250 275;]
b=[1400; 100; 2000; 0]
x0=[a; b; c; d]
l=-%inf
n=-%inf
function[x,cont]=jacobi(A, b, x0, L, n)
 [l,c]=size(A)
 erro=1
 x=x0
 cont=0
 while erro>L & cont<n
 cont=cont+1
 xa=x
 for i=1:l
 soma=0
 for j=1:l
 if j~=i
 soma = soma + A(i,j)*xa(j)
 end
 end
 x(i)=(b(i)-soma)/A(i,i)
 end
 erro=max(abs(x-xa))/max(abs(x))
 end
endfunction
r=jacobi(A, b, x0, l, n)
disp(r)