b=matrix([1,2,2])
a=matrix([2,-2,1])
print a*(a.transpose()*b)         /*projection of b on a p=(a(a.transpose()*b))/(a.tranpose()*a)

Strang-3.5-3.9-U (last edited 2010-12-18 07:47:19 by manishadeepthi)