#svd version to solve Ax=b as in gauss.py


import numpy as np
import numpy.linalg as la





A = np.loadtxt("matrix")
b = np.loadtxt("vector")
b = b.reshape(len(b),1)


U,s,V = la.svd(A)

print A
print b
print
print
print U
print s
print V

print "Condition = ", s[0]/s[len(s)-1]

x=np.dot(np.transpose(U),b)
for i in range(len(s)):
     x[i,0]=x[i,0]/s[i]
x=np.dot(np.transpose(V),x)
print x