with(plots):

r:=1+cos(t);

x:=r*cos(t);

rho:=r*sin(t);

y:=rho*cos(s);

z:=rho*sin(s);

plot3d([x,y,z],t=0..Pi,s=0..Pi,scaling=constrained);

with(plots):

r:=1+cos(2*t);

rho:=r*cos(t);

y:=r*sin(t);

x:=rho*cos(s);

z:=rho*sin(s);      

plot3d([x,y,z],t=0..Pi/2,s=0..Pi,scaling=constrained,grid=[25,25]);

with(plots):

r:=3*cos(t)-cos(3*t);

x:=r*cos(s);

y:=r*sin(s);

z:=3*sin(t)-sin(3*t);

plot3d([x,y,z],t=-Pi/2..Pi/2,s=0..Pi,scaling=constrained);