v := i+j; w := i+j+k; len := vec -> sqrt(dotprod(vec,vec)): sint := len(crossprod(v,w))/(len(v)*len(w)); convert(angle(v,w),degrees); evalf(%);