b = matrix([[1],[1],[1]])
a = matrix([[1],[2],[2]])
t=a.transpose()
x1=t*b
x2=t*a
x=x1[0]/x2[0]
P=x*a 
print 'Projection matrix P is'
print P

Strang-4.2-1 (last edited 2010-12-18 09:11:58 by vasanthi emani)