Metodo_Gauss_Seidel

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

clear all
close all
clc
%%METODO GAUSS SEIDEL
%%matriz
A=[ 5 1 1 ; 3 4 1 ; 3 3 6]
b=[5 6 0]
xanterior=[0 ; 0 ; 0]
e=0.05;
n=size (A,1);
 for i=1:n
 for j=1:n 
 %Matriz C - nova matriz a com 0 na diagonal
 C(i,j)= (-A(i,j))/A(i,i);
 C(i,i)=0;
 
 %Novo b
 B(i)=b(i)/A(i,i);
 
 end
 end
for k=1:20
 
 h = xanterior;
 
 for i=1:n
 
 %Novos x
 xposterior(i) = C(i,:)*xanterior + B(i);
 
 xanterior(i) = xposterior(i)
 
 end
 
 if abs(xposterior - h')<e
 break
 
 end
 
 end