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