import numpy as np
import os
os.system("clear") 

# creo datos como una lista anidada
a=[[214.0,2.0,-4.0,8.5],
    [1.0,15.0,-2.0,2.0],
    [-10.0,-2.0,5.0,-1.0],
    [82.0,13.0,-17.0,1.0]] 
# trasnformo la lista en array NP
aa=np.array(a)
print "a"
print aa
# idem con un vector
# son los mismos del ejemplo de Gauss
#b=[-7.0,18.0,-23.0,12.0] 
b=[[-7.0],[18.0],[-23.0],[12.0]]
bb=np.array(b)
print "b"
print bb
print
#resolucion de sistema de ecuaciones
#verificar que da el mismo resultado que Gauss programado
x=np.linalg.solve(aa,bb)
print "Solucion (x):"
print x
print
# utilizando multiplicacion matricial del modulo numpy
print "Verificacion (a.x);"
print aa.dot(x)
print
# multiplico  matrices, elemento a elemento
print "multip elemento a elemento (a*x)"
print aa*x
print