a=vector([3,4])
b=vector([4,3])
p=a.dot_product(b)/(a.norm() * b.norm())
print p
float(arccos(p))

Strang-3.3-U (last edited 2010-12-18 10:14:08 by pratyureepu)