> restart:

> r:=3+cos(t):

> x:=r*cos(s):

> y:=r*sin(s):

> z:=sin(t):

> plot3d([x,y,z],s=0..2*Pi,t=0..2*Pi,scaling=constrained,title="87121 By J.S.C.");