PQ := 3*cos(t)*i+2*j; F := 2*j; T := crossprod(PQ,F); t := arcsin(2/3); T := T[3]*k;