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=[5,0,0,-3,-1,2];l20=[-1,4,0,0,-1,3];l30=[0,0,2,-1,0,-1];l40=[-1,0,0,4,-2,0];l50=[0,0,0,-1,2,-1]; M1=[l10;l20;l30;l40;l50]

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

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

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

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

x11=(1/l10(1))*(l10(6)-l10(2)*v2-l10(3)*v3-l10(4)*v4-l10(5)*v5)

x21=(1/l20(2))*(l20(6)-l20(1)*x11-l20(3)*v3-l20(4)*v4-l20(5)*v5)

x31=(1/l30(3))*(l30(6)-l30(1)*x11-l30(2)*x21-l30(4)*v4-l30(5)*v5)

x41=(1/l40(4))*(l40(6)-l40(1)*x11-l40(2)*x21-l40(3)*x31-l40(5)*v5)

x51=(1/l50(5))*(l50(6)-l50(1)*x11-l50(2)*x21-l50(3)*x31-l50(4)*x41)

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

x1=[x11;x21;x31;x41;x51]

% fim da Primeira iteração

% 2º etapa - 2º iteração

x12=(1/l10(1))*(l10(6)-l10(2)*x21-l10(3)*x31-l10(4)*x41-l10(5)*x51)

x22=(1/l20(2))*(l20(6)-l20(1)*x12-l20(3)*x31-l20(4)*x41-l20(5)*x51)

x32=(1/l30(3))*(l30(6)-l30(1)*x12-l30(2)*x22-l30(4)*x41-l20(5)*x51)

x42=(1/l40(4))*(l40(6)-l40(1)*x12-l40(2)*x22-l40(3)*x32-l20(5)*x51)

x52=(1/l50(5))*(l50(6)-l50(1)*x12-l50(2)*x22-l50(3)*x32-l50(4)*x42)

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

x1=[x11;x21;x31;x41;x51]

x2=[x12;x22;x32;x42;x52]

% 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)

Ex14=abs(x41-x42)/abs(x22)

Ex15=abs(x51-x52)/abs(x22)

EX1=[Ex11;Ex12;Ex13;Ex14;Ex15]

% ate aqui*******************

% 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]