a=vector([2,-2,1])
b=vector([1,2,2])
arccos ((a.dot_product(b))/norm(a)*norm(b))

Strang-3.5-3.17-U (last edited 2010-12-17 12:52:52 by manishadeepthi)