u := vector([1,0]); v := vector([cos(t),sin(t)]); angle(u,u+v); ang := simplify(",trig); ang := subs(t=2*s,ang); ang := arccos(expand(cos(ang))); simplify(");