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},Shading->False]
