s1Eva_IIT2011_T2_AN-Sistema-de-ecuaciones

Propuesta de solución

matriz aumentada: 
[[ -2   5   9   1]
 [  7   1   1   6]
 [ -3   7  -1 -26]]
diagonal dominante: 
[[  7   1   1   6]
 [ -3   7  -1 -26]
 [ -2   5   9   1]]
nuevas matrices: 
A: 
[[ 7  1  1]
 [-3  7 -1]
 [-2  5  9]]
B: 
[  6 -26   1]
solucion X: 
[ 1.00000245 -2.99999785  1.99999935]
verifica A.X=B
[  6.00001868 -25.99999164   1.        ]
>>> 

Desarrolle manualmente el ejercicio para verificar sus respuestas.
el archivo de matg1013.py con las funciones desarrolladas en clase se encuentra en la seccion de Recursos de estudio.

# 1ra Evaluación II Término 2011
# Tema 2. Sistema de ecuaciones
import numpy as np
import matg1013 as fcnm

A = np.array([[-2, 5, 9],
              [ 7, 1, 1],
              [-3, 7,-1]])
B = np.array([1,6,-26])
tolera = 1e-4
X = np.zeros(len(B), dtype=float)
iteramax = 100

# PROCEDIMIENTO

# Diagonal dominante:
B = np.transpose([B]) # en columna
aumentada = np.concatenate((A,B),axis=1)
dominante = fcnm.pivoteafila(aumentada)
# Separa A y B de dominante
A = dominante[:,:-1]
B = dominante[:,-1]

respuesta = fcnm.gauss_seidel(A,B,X,tolera, iteramax)
verifica = np.dot(A,respuesta)

# SALIDA
print('matriz aumentada: ')
print(aumentada)
print('diagonal dominante: ')
print(dominante)
print('nuevas matrices: ')
print('A: ')
print(A)
print('B: ')
print(B)
print('solucion X: ')
print(respuesta)
print('verifica A.X=B')
print(verifica)