Eliminação de Gauss(python)

Prévia do material em texto

03/02/2021 Trabalho 03
localhost:8888/nbconvert/html/Trabalho 03.ipynb?download=false 1/6
1. Método exato ¶
Desenvolvendo na mão
03/02/2021 Trabalho 03
localhost:8888/nbconvert/html/Trabalho 03.ipynb?download=false 2/6
03/02/2021 Trabalho 03
localhost:8888/nbconvert/html/Trabalho 03.ipynb?download=false 3/6
03/02/2021 Trabalho 03
localhost:8888/nbconvert/html/Trabalho 03.ipynb?download=false 4/6
03/02/2021 Trabalho 03
localhost:8888/nbconvert/html/Trabalho 03.ipynb?download=false 5/6
Importar biblioteca
In [6]:
import numpy as np 
Criar função
In [7]:
def Elim_gauss (a,b): 
 n = len(b) 
 x = np.zeros(n) 
 
 for k in range (0,n-1): 
 
 for i in range ( k+1 ,n ): 
 m = a[i][k]/a[k][k] 
 a[i][k] = 0 
 b[i] = b[i] - m*b[k] 
 
 for j in range (k+1,n): 
 a[i][j] = a[i][j] - m*a[k][j] 
 
 x[n-1] = b[n-1]/a[n-1][n-1] 
 
 for i in range (n-2,-1,-1): 
 soma = 0 
 
 for j in range (i+1,n): 
 soma = soma + a[i][j]*x[j] 
 
 x[i] = (b[i]-soma)/a[i][i] 
 
 
 return (x) 
Definir os valores
In [8]:
a = [[-1,0,-70],[1,-1,-60],[0,1,-40]] 
b = [-736.46,-658.4,-477.27] 
03/02/2021 Trabalho 03
localhost:8888/nbconvert/html/Trabalho 03.ipynb?download=false 6/6
Aplicar a função
In [9]:
resultado = Elim_gauss(a,b) 
Mostrar os valores
In [10]:
print('T = %0.2f'%resultado[0]) 
print('R = %0.2f'%resultado[1]) 
print('a = %0.2f'%resultado[2]) 
In [ ]:
 
T = -34.42 
R = -36.77 
a = 11.01