Algoritmos para resolução de questões em Cálculo Numérico para Eliminação Gaussiana com e sem pivotamento

Eliminação Gaussiana sem pivotamento:

- A

- b

- Lu sem

- carrega b

- y=L\b (termo independente)

- x=U\y (solução)

Eliminação Gaussiana com pivotamento:

- A

- b

- [L,U,P]=Lu (A)

- y=L\(P*b) (termo independente)

- x=U\y (solução)

- vetor de pivotamento (indica a troca de linhas que foram realizadas)

->A=[0.2641 0.1735 0.8624;0.3457 1.8423 0.4759;-0.8468 0.4759 1.2147]

 A  =

    0.2641    0.1735    0.8624 

    0.3457    1.8423    0.4759 

  - 0.8468    0.4759    1.2147 

-->[L,U,P]=lu(A)

 P  =

    0.    0.    1. 

    0.    1.    0. 

    1.    0.    0. 

 U  =

  - 0.8468    0.4759       1.2147    

    0.        2.0365827    0.9717925 

    0.        0.           1.0876289 

 L  =

    1.           0.           0. 

  - 0.4082428    1.           0. 

  - 0.3118800    0.1580705    1. 

OBSERVAÇÕES:

k= número de interações que demora para convergir à solução do método

norma_r= resíduo

vetor zero no Scilab: zeros(n,1)

matriz identidade (In): eye(n,n)

Rodar os seguintes scripts no SCILAB: sor, jacobi e lu.sem