Jacobi

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

//programa para aplicação do método de Jacobi a um sistema linear de ordem 3: AX=B
clear 
//Aqui a(i,j) denota o elemento na matriz A da i-ésima linha e j-ésima coluna
a(1,1)=10;
a(1,2)=2;
a(1,3)=1;
a(2,1)=1;
a(2,2)=5;
a(2,3)=1;
a(3,1)=2;
a(3,2)=3;
a(3,3)=10;
// Aqui b(j) denota o j-ésimo elemento do vetor B de termos independentes no sistema AX=B.
b(1)=7;
b(2)=-8;
b(3)=6;
i=1 // i é um índice usado nas iterações a seguir
dk=1; // dk está relacionado com o erro
// Abaixo x(i,j) representa o valor de xi na (j-1)-ésima iteração (Aqui não podemos usar j=0)
x(1,1)=0;
x(2,1)=0;
x(3,1)=0;
while dk >0.1
 x(1,i+1)=(b(1)-a(1,2)*x(2,i)-a(1,3)*x(3,i))/a(1,1);
 x(2,i+1)=(b(2)-a(2,1)*x(1,i)-a(2,3)*x(3,i))/a(2,2);
 x(3,i+1)=(b(3)-a(3,1)*x(1,i)-a(3,2)*x(2,i))/a(3,3);
 d(1)=abs(x(1,i+1)-x(1,i));
 d(2)=abs(x(2,i+1)-x(2,i));
 d(3)=abs(x(3,i+1)-x(3,i));
 dk=max(d);
 i=i+1;
end
disp('número de iterações'), disp(i-1);
disp('Solução aproximada')
disp(x(:,i))
disp('dk='); disp(dk);