Metodo de gauss jacobi

Prévia do material em texto

A = [10 2 1;1 5 1;2 3 10] 
Ao = [0 2 1;1 0 1;2 3 0] 
B = [7;-8;6] 
 
[ni,nj] = size(A) 
m = 1 
j = 0 
h = 0 
s = 0 
i = 0 
d = 0 
k = 0 
alfa = 0 
sa = 0 
AUX = 0 
 //verificando a convergencia criterio das linhas 
 
 for i=1:ni 
 
 for j=1:nj 
 
 sa = sa + Ao(i,j) 
 //disp(sa) 
 end 
 
 alfa(i) = sa/A(i,i) 
 sa = 0 
 
 end 
 
 if max(alfa) < 1 then 
 
 disp("Esta função converge :)") 
 disp("") 
 end 
 
 //Interações 
 e = input("Digite a tolerancia-->") 
 d = 1+e 
 for h = h+1:ni 
 
 xo(h,1) = input("Insira o chute-->") 
 
 end 
 
 while d > e, 
 
 k = k+1 
 
 for i = 1:ni 
 for j = 1:nj 
 s = s + (Ao(i,j)/A(i,i))*xo(j) 
 end 
 
 x(i,1) = (B(i,1)/A(i,i)) - s 
 sub(i,1) = abs(x(i,1) - xo(i,1)) 
 //Não esta indo xo(i,1) = x(i,1) 
 s = 0 
 end 
 
 xo = x 
 
 d = (max(abs(sub))/max(abs(x))) 
 disp(d) 
 end 
 
 disp("Resultado do sistema-->") 
 disp(x)