restart;
f:=r->3+cos(r);
x:=r*cos(theta);
y:=r*sin(theta);
z:=f(r);
plot3d([x,y,z],r=0..2*Pi,theta=0..2*Pi);