Novo Link  de trabalhos academicos https://www.fenix-eng.info

no Mathlab 

>> 

%Metodo de Gauss Siel  

% entrar com a Matriz 

% Verificar qual a condição Neste caso E<10^-2

%k=0

% Apresentar a matriz

l10=[3,2,4,1];l20=[1,1,2,2];l30=[4,3,2,3];M1=[l10;l20;l30]

% a Cada iteração a partir da 2º iteração verificar a condição. 

% criação do vetor de aproximaçãovetor nulo

v1=[0];v2=0;v3=[0];V=[v1;v2;v3]

% 1º iteração para encontrar x1 - x2 - x3  

x11=(1/l10(1))*(l10(4)-l10(2)*v2-l10(3)*v3)

x21=(1/l20(2))*(l20(4)-l20(1)*x11-l20(3)*v3)

x31=(1/l30(3))*(l30(4)-l30(1)*x11-l30(2)*x21)

% apresentar valores de x1-x2-x3 da 1º iteração

x1=[x11;x21;x31]

% fim da Primeira iteração

% 2º etapa - 2º iteração

x12=(1/l10(1))*(l10(4)-l10(2)*x21-l10(3)*x31)

x22=(1/l20(2))*(l20(4)-l20(1)*x12-l20(3)*x31)

x32=(1/l30(3))*(l30(4)-l30(1)*x12-l30(2)*x22)

% apresentar valores de x1-x2-x3 da 2º iteração

x1=[x11;x21;x31]

x2=[x12;x22;x32]

% fim da 2º iteração

% iniciar Teste E=10-²

% |x1 - x2|/|x2 Maximo| 1,2 interações

Ex11=abs(x11-x12)/abs(x22)

Ex12=abs(x21-x22)/abs(x22)

Ex13=abs(x31-x32)/abs(x22)

EX1=[Ex11;Ex12;Ex13]

% não passou no teste caso 1 não passe no teste  tem que fazer mais interações 

% 3º etapa - 3º iteração

x13=(1/l10(1))*(l10(4)-l10(2)*x22-l10(3)*x32)

x23=(1/l20(2))*(l20(4)-l20(1)*x13-l20(3)*x32)

x33=(1/l30(3))*(l30(4)-l30(1)*x13-l30(2)*x23)

x1=[x11;x21;x31];x2=[x12;x22;x32];x3=[x13;x23;x33]

% iniciar Teste 2 E=10-²

% |x3 - x4|/|x4 Maximo| 2,3 interações

Ex21=abs(x13-x12)/abs(x23)

Ex22=abs(x23-x22)/abs(x23)

Ex23=abs(x33-x32)/abs(x23)

EX2=[Ex21;Ex22;Ex23]

% não passou no teste caso 1 não passe no teste  tem que fazer mais interações 

% 4º etapa - 4º iteração

x14=(1/l10(1))*(l10(4)-l10(1)*x13-l10(3)*x33)

x24=(1/l20(2))*(l20(4)-l20(1)*x14-l20(3)*x33)

x34=(1/l30(3))*(l30(4)-l30(1)*x14-l30(2)*x24)

x1=[x11;x21;x31];x2=[x12;x22;x32];x3=[x13;x23;x33];x4=[x14;x24;x34]

% iniciar Teste 2 E=10-²

% |x3 - x4|/|x4 Maximo| 2,3 interações

Ex31=abs(x14-x13)/abs(x14)

Ex32=abs(x24-x23)/abs(x14)

Ex33=abs(x34-x33)/abs(x14)

EX1=[Ex11;Ex12;Ex13]

EX2=[Ex21;Ex22;Ex23]

EX3=[Ex31;Ex32;Ex33]