Ang := (v,n) -> convert(Pi/2 - angle(v,n),degrees); Ang([1,0,0],[-1,1,0]); Ang([1,-3,0],[1,1,-1]); evalf(%);