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