Allgemeine Drehflächen
> DR:=evalm(Drehmatrix&*vector([1,0,u,5*sin(u)]));
> plot3d([DR[2],DR[3],DR[4]],t=0..2*Pi,u=0..4*Pi,axes=FRAMED,numpoints=2000,scaling=constrained,orientation=[30,60]);
>