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