Estoy realizando un código para realizar matrices inversas, el código genera una matriz identidad y la adjunta a la primera matriz, así mismo para imprimirlas las separa, luego durante la ejecución desarrolla las columnas inferiores sin problema, pero al llegar a la segunda parte para desarrollar las columnas superiores las iteraciones no llegan hasta A[0][1] para reemplazar A[0]:
#Esta es la matriz con la primera parte ya desarrollada:
A = [
[1, 1/3, 1/3, 1/3, 0, 0],
[0, 1, 10, 1, -3, 0],
[0, 0, 1, 1/5, -1/10, -1/10]
]
#Función para mostrar la matriz con espacios.
def mostrar():
for i in range(len(A)):
for x in range(len(A[i])):
if x != len(A[i])-(len(A[i])/2):
print("{:.3f}".format(A[i][x]), end="\t")
else: print("| \t {:.3f}".format(A[i][x]), end="\t")
print(" ")
print(" ")
#Función para desarrollar segunda parte (columnas superiores.)
def col2(i):
verificador = i
uwu = len(A)-1
while verificador != -len(A)-1:
#print("verificador:", verificador, "len(A):", len(A))
pivote = A[verificador][uwu]
print("Cambiar columnas.")
print("")
print("range:",range(len(A[i])-1))
for x in range(len(A[i])):
print("x:",x)
print("{} - ( {} * {} ) = {}".format(A[verificador][x], pivote, A[i+1][x], A[verificador][x] - pivote * A[i+1][x]))
A[verificador][x] = A[verificador][x] - pivote * A[i+1][x]
mostrar()
verificador -= 1
uwu -= 1
for i in range(-2, -len(A)-1, -1):
col2(i)
print("")
print("RESPUESTA:")
mostrar()
El output que me dá es este:
1.000 0.333 0.000 | 0.267 0.033 0.033
0.000 1.000 0.000 | -1.000 -2.000 1.000
0.000 0.000 1.000 | 0.200 -0.100 -0.100
Pero el esperado es este:
1.000 0.000 0.000 | 0.600 0.700 -0.300
0.000 1.000 0.000 | -1.000 -2.000 1.000
0.000 0.000 1.000 | 0.200 -0.100 -0.100
La diferencia está en el A[0] pero no sé cómo solucionarlo, agradezco cualquier comentario o sugerencia.
Gracias :D.