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}]