Estoy haciendo matrices de 3x3 con Gauss Jordan para resolver primero los tres primeros ceros de la parte posterior de la izquierda. No estoy usando numpy ni ninguna librería, todo es de forma manual. Al momento de llamar a dos variables para dividir me arroja un error de que una variable no está definida:
NameError: name 'div2' is not defined
Aquí está el problema:
numero=div2/div4
Este es mi código, les advierto que es exageradamente redundante, pero funciona para mi.
from fractions import Fraction
print('\nFila 1')
a1=int(input('Valor a: '))
b1=int(input('Valor b: '))
c1=int(input('Valor c: '))#Se definen las matrices
d1=int(input('Valor n: '))
print('\nFila 2')
a2=int(input('Valor a: '))
b2=int(input('Valor b: '))
c2=int(input('Valor c: '))#Se definen las matrices
d2=int(input('Valor n: '))
print('\nFila 3')
a3=int(input('Valor a: '))
b3=int(input('Valor b: '))
c3=int(input('Valor c: '))#Se definen las matrices
d3=int(input('Valor n: '))
print(f'\n {a1}x {b1}y {c1}z | {d1}\n', f'{a2}x {b2}y {c2}z | {d2}\n', f'{a3}x {b3}y {c3}z | {d3}')#Se imprimen de forma estética las matrices con sus incógnitas
print(f'\n {a1} {b1} {c1} | {d1}\n', f'{a2} {b2} {c2} | {d2}\n', f'{a3} {b3} {c3} | {d3}') #Se definen las matrices de forma estética sin incógnitas
if a2 and a3 > a1: #Se declara la primera condición para poner en la fila1 la variable (a1,a2 o a3) más pequeña
v1=(a1,b1,c1,d1) #Si se cumple la condición 1...
print(v1) #Se imprime la fila1 que es la más pequeña por la variable a1
if (a1 != 1) and (b1 != 1) and (c1 != 1) and (d1 != 1): #Si la fila 1 es diferente de 1....
a,b,c,d=Fraction(a1/a1),Fraction(b1/a1),Fraction(c1/a1),Fraction(d1/a1) #Se divide entre a1
print(f'\n{a}, {b}, {c}, {d}') #Se imprime la fila ya dividida
if (a2 < 0): #Terminada la condición anterior se declara una condición donde si el primer valor de la fila2 es diferente de 0...
opuesto1 = abs(a2) #Se saca su valor absoluto
multiplicacion_fila2_con_fila1=(opuesto1*a1/a1),(opuesto1*b1/a1),(opuesto1*c1/a1),(opuesto1*d1/a1) #La fila 2 se multiplica con la fila 1
print(multiplicacion_fila2_con_fila1) #Se imprime
resta_fila2_con_fila1=((opuesto1*a1/a1)+a2),((opuesto1*b1/a1)+b2),((opuesto1*c1/a1)+c2),((opuesto1*d1/a1)+d2) #Se resta el producto de la fila1 y fila2 con la fila2
print(resta_fila2_con_fila1) #Se imprime
div1=((opuesto1*b1/a1)+b2) #Se declara la variable dev1 que nos va a servir más adelant para dividir el valor 2 de la fila2 con el valor 2 de la fila 3 y multiplicar ese valor con la fila 3 para obtener el tercer 0.
else:
opuesto2 = a2 * (-1)
multiplicacion_fila2_con_fila1=(opuesto2*a1/a1),(opuesto2*b1/a1),(opuesto2*c1/a1),(opuesto2*d1/a1)
print(multiplicacion_fila2_con_fila1)
resta_fila2_con_fila1=((opuesto2*a1/a1)+a2),((opuesto2*b1/a1)+b2),((opuesto2*c1/a1)+c2),((opuesto2*d1/a1)+d2)
print(resta_fila2_con_fila1)
div2=((opuesto2*b1/a1)+b2)
if (a3 < 0):
opuesto3 = abs(a3)
multiplicacion_fila3_con_fila1=(opuesto3*a1/a1),(opuesto3*b1/a1),(opuesto3*c1/a1),(opuesto3*d1/a1)
print(multiplicacion_fila3_con_fila1)
resta_fila3_con_fila1=((opuesto3*a1/a1)+a3),((opuesto3*b1/a1)+b3),((opuesto3*c1/a1)+c3),((opuesto3*d1/a1)+d3)
print(resta_fila3_con_fila1)
div3=((opuesto3*b1/a1)+b3)
else:
opuesto4 = a3 * (-1)
multiplicacion_fila3_con_fila1=(opuesto4*a1/a1),(opuesto4*b1/a1),(opuesto4*c1/a1),(opuesto4*d1/a1)
print(multiplicacion_fila3_con_fila1)
resta_fila3_con_fila1=((opuesto4*a1/a1)+a3),((opuesto4*b1/a1)+b3),((opuesto4*c1/a1)+c3),((opuesto4*d1/a1)+d3) #Aquí condición 1 con fila3 con 2
print(resta_fila3_con_fila1)
div4=((opuesto4*b1/a1)+b3)
numero=div2/div4
print(numero)
if (numero < 0):
opuesto = abs(numero)
multiplicacion_fila3_con_fila2=(opuesto*((opuesto4*a1/a1)+a3)),(opuesto*((opuesto4*a1/a1)+b3)),(opuesto*((opuesto4*a1/a1)+c3)),(opuesto*((opuesto4*a1/a1)+d3))
print(multiplicacion_fila3_con_fila2)
resta_fila3_con_fila2=((opuesto4*a1/a1)+a3)+((opuesto2*a1/a1)+a2),((opuesto*a3)+b2),((opuesto*a3)+c2),((opuesto*a3)+d2)#Aquí
print(resta_fila3_con_fila2)
else:
opuesto = numero * (-1)
multiplicacion_fila3_con_fila2=(opuesto*((opuesto4*a1/a1)+a3)),(opuesto*((opuesto4*a1/a1)+b3)),(opuesto*((opuesto4*a1/a1)+c3)),(opuesto*((opuesto4*a1/a1)+d3))
print(multiplicacion_fila3_con_fila2)
resta_fila3_con_fila2=((opuesto4*a1/a1)+a3)+((opuesto2*a1/a1)+a2),((opuesto*a3)+b2),((opuesto*a3)+c2),((opuesto*a3)+d2)
print(resta_fila3_con_fila2)
else:
if (a2 < 0):
opuesto = abs(a2)
multiplicacion_fila2_con_fila1=(opuesto*a1/a1),(opuesto*b1/a1),(opuesto*c1/a1),(opuesto*d1/a1)
print(multiplicacion_fila2_con_fila1)
else:
opuesto = a2 * (-1)
multiplicacion_fila2_con_fila1=(opuesto*a1/a1),(opuesto*b1/a1),(opuesto*c1/a1),(opuesto*d1/a1)
print(multiplicacion_fila2_con_fila1)
elif a1 and a3 > a2:
v2=(a2,b2,c2,d2)
print(v2)
if (a2 != 1) and (b2 != 1) and (c2 != 1) and (d2 != 1):
a1,b1,c1,d1=Fraction(a2/a2),Fraction(b2/a2),Fraction(c2/a2),Fraction(d2/a2)
print(f'\n{a1},{b1},{c1},{d1}')
elif a1 and a2 > a3:
v3=(a3,b3,c3,d3)
print(v3)
else:
print('Todas las matrices son iguales')
¿Cómo puedo solucionar este problema? Por cierto, si existe alguna forma de optimizar el código estaré encantado por la ayuda. Agradezco mucho su apoyo, saludos!