The plot of f4 in cylindrical coordinates

> plot3d(f4,r=1.1..10,theta=0..2*Pi,coords=cylindrical,axes=frame,orientation=[45,70],labels=[`x`,`y`,`z`]);

>