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