a = matrix([[1],[2],[2]])
t=a.transpose()
x=a*t
y=t*a
z=y.inverse()
x=x.with_rescaled_row(0,1/9)
x=x.with_rescaled_row(1,1/9)
P=x.with_rescaled_row(2,1/9)
print 'Projection matrix P is'
print P

Strang-4.2-2 (last edited 2010-12-18 09:43:58 by Siva Jyoshna Poreddy)