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