ALNP2EJ1

#EJERCICIO 1
#APARTADO 1: CAMBIAR REALFIELD POR LO QUE DIGA LA MANTISA Y EL CORCHETE COPYPASTEAR MATRIZ, REESCRIBIR LOS MENOS, NADA MAS
print '\n Apartado 1:'
print'\nMatriz A:'
A1 = matrix(RealField(7), 3, [0.872,0.845,-2.966,-1.866,-0.661,-1.118,2.616,-1.289,-0.205]).transpose(); show(A1)
print '\nVector b:'
b1 = vector(RealField(7), [4.238,9.078,-30.573]); show(b1)

#--------------------------------------A PARTIR DE AQUI NO TOCAR---------------------------------------

print'\nMatriz ampliada:'
M1 = A1.augment(matrix(b1).transpose()); show(M1)
M11 = block_matrix([[A1, b1.column()]])
show(M11)
F11, U11 = forma_escalonada(M11)
print'\nMatriz ampliada escalonada:'
show(U11)
print '\n===========================La solución al apartado 1 es:==========================='
show(sustitucion_regresiva(U11))
F12, U12 = forma_escalonada(M11, False, False, algoritmo = 'pivoteo parcial escalado')
print '\n===========================La solución al apartado 2 es:==========================='
show(sustitucion_regresiva(U12))