Plot of Initial condition

> plot3d([r,theta,u(r,theta,0)],r=0..1,theta=-Pi..Pi,coords=cylindrical,axes=frame,orientation=[54,42],labels=[`x`,`y`,`u`]);

>