xmin:=0
xmax:=2*Pi
ymin:=0
ymax:=2*Pi
xxx[x_,y_]:=(3+Cos[x]+Cos[3*y])*Cos[y]
yyy[x_,y_]:=(3+Cos[x]+Cos[5*y])*Sin[y]
zzz[x_,y_]:=Sin[x]
ParametricPlot3D[{yyy[x,y],-xxx[x,y],zzz[x,y]},{x,xmin,xmax},{y,ymin,ymax},PlotPoints->{15,51},Shading->False]
