a = matrix([[1, 0], [1, 1],[1, 2]])
b = matrix([[6],[0],[0]])
x1=a.transpose()*a
x2=a.transpose()*b
x=x1.inverse()*x2
print 'x cap is:'
print x
p=a*x
print 'p is :'
print p
P= a* x1.inverse()*a.transpose()
print 'P is :'
print P

Strang-4.2-3 (last edited 2010-12-18 09:30:03 by vasanthi emani)