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

no Mathlab 

>> %Metodo de Gauss Jacob - 4 iterações

% entrar com a Matriz 

>> %k=0

>> % Apresentar a matriz

>> l10=[10,5,2,1,18],l20=[5,10,-3,2,14],l30=[2,3,12,-5,6],l40=[1,2,-5,12,10];M1=[l10;l20;l30;l40]

>> % criação do vetor de aproximação

>> v10=[l10(5)/l10(1);l20(5)/l20(2);l30(5)/l30(3);l40(5)/l40(4)]

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

>> x11=(1/l10(1))*(l10(5)-l10(2)*v10(2)-l10(3)*v10(3)-l10(4)*v10(4))

>> x21=(1/l20(2))*(l20(5)-l20(1)*v10(1)-l20(3)*v10(3)-l20(4)*v10(4))

>> x31=(1/l30(3))*(l30(5)-l30(1)*v10(1)-l30(2)*v10(2)-l30(4)*v10(4))

>> x41=(1/l40(4))*(l40(5)-l40(1)*v10(1)-l40(2)*v10(2)-l40(3)*v10(3))

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

>> X1= [x11;x21;x31;x41]

% fim da Primeira iteração

% 2º etapa - 2º iteração

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

x22=1/l20(2)*(l20(5)-l20(1)*x11-l20(3)*x31-l20(4)*x41)

x32=1/l30(3)*(l30(5)-l30(1)*x11-l30(2)*x21-l30(4)*x41)

x42=1/l40(4)*(l40(5)-l40(1)*x11-l40(2)*x21-l40(3)*x31)

X1=[x11;x21;x31;x41],X2=[x12;x22;x32;x42]

% fim da segunda iteração

% 3º etapa - 3º iteração

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

x23=1/l20(2)*(l20(5)-l20(1)*x12-l20(3)*x32-l20(4)*x42)

x33=1/l30(3)*(l30(5)-l30(1)*x12-l30(2)*x22-l30(4)*x42)

x43=1/l40(4)*(l40(5)-l40(1)*x12-l40(2)*x22-l40(3)*x32)

X1=[x11;x21;x31;x41],X2=[x12;x22;x32;x42],X3=[x13;x23;x33;x43]

% fim da terceira iteração

% 4º etapa - 4º iteração

x14=1/l10(1)*(l10(5)-l10(2)*x23-l10(3)*x33-l10(4)*x43)

x24=1/l20(2)*(l20(5)-l20(1)*x13-l20(3)*x33-l20(4)*x43)

x34=1/l30(3)*(l30(5)-l30(1)*x13-l30(2)*x23-l30(4)*x43)

x44=1/l40(4)*(l40(5)-l40(1)*x13-l40(2)*x23-l40(3)*x33)

X= [x11;x21;x31;x41],X2=[x12;x22;x32;x42],X3=[x13;x23;x33;x43],X4=[x14;x24;x34;x44]

% Estimativa de precisão "E" 3 interação 

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