1.mws

>    restart;

>    with(plots):with(plottools); 

Warning, the name changecoords has been redefined



Warning, the assigned name arrow now has a global binding

[arc, arrow, circle, cone, cuboid, curve, cutin, cutout, cylinder, disk, dodecahedron, ellipse, ellipticArc, hemisphere, hexahedron, homothety, hyperbola, icosahedron, line, octahedron, parallelepiped,...
[arc, arrow, circle, cone, cuboid, curve, cutin, cutout, cylinder, disk, dodecahedron, ellipse, ellipticArc, hemisphere, hexahedron, homothety, hyperbola, icosahedron, line, octahedron, parallelepiped,...
[arc, arrow, circle, cone, cuboid, curve, cutin, cutout, cylinder, disk, dodecahedron, ellipse, ellipticArc, hemisphere, hexahedron, homothety, hyperbola, icosahedron, line, octahedron, parallelepiped,...
[arc, arrow, circle, cone, cuboid, curve, cutin, cutout, cylinder, disk, dodecahedron, ellipse, ellipticArc, hemisphere, hexahedron, homothety, hyperbola, icosahedron, line, octahedron, parallelepiped,...

>    a:=2;b:=5;c:=sqrt(a^3+b^3);

a := 2

b := 5

c := 133^(1/2)

>    z:=a*sin(t);r:=c+a*cos(t);

z := 2*sin(t)

r := 133^(1/2)+2*cos(t)

>    x:=r*cos(s);y:=r*sin(s)*cos(s);

x := (133^(1/2)+2*cos(t))*cos(s)

y := (133^(1/2)+2*cos(t))*sin(s)*cos(s)

>    w:=solve(z=a/b*x,s);

w := arccos(5*sin(t)/(133^(1/2)+2*cos(t)))

>    plot3d([x,y,z],s=-w..w,t=0..2*Pi,scaling=constrained,grid=[40,40]);

[Maple Plot]

>   

>