# test_solve.py
from numpy import array
from numpy.linalg import solve
from numpy.linalg import qr
from numpy.linalg import eig
from numpy.linalg import svd
from numpy.matlib import rand,zeros,ones,eye
import numpy.linalg

# can initialize an array to zero
# a=array([[0.0 for j in range(n)] for i in range(m)]) # n columns, m rows 

print "test_solve.py  a x = y "

a=array([[1.0/2,1.0/3,1.0/4,1.0/5],[1.0/3,1.0/4,1.0/5,1.0/6],[1.0/4,1.0/5,1.0/6,1.0/7],[1.0/5,1.0/6,1.0/7,1.0/8]])
print "a="
print  a
y=array([1,2,3,4])
print "y="
print y

x=solve(a,y)
print "solve x="
print x
d,v=eig(a)

print "eig d="
print d
print "eig v="
print v

q1,q2=qr(a)
print "qr q1="
print q1
print "qr q2="
print q2

s1,s2,s3=svd(a)
print "svd s1="
print s1
print "svd s2="
print s2
print "svd s3="
print s3

b=rand(3,3)
print "rand(3,3)"
print b

print "a[0][0]"
print a[0][0]
print "a[1][1]"
print a[1][1]
print "a[2][2]"
print a[2][2]
print "a[3][3]"
print a[3][3]
print "a[0][3]"
print a[0][3]
a[3][3]=7.5


print "dir(numpy.linalg)"
print dir(numpy.linalg)
print "dir(numpy.matlib)"
# print dir(numpy.matlib # big