> with(plots):
Let us plot the cylindrical part of the surface which encloses the volume under consideration:
>
plot3d([1,theta,z],theta=0..2*Pi,z=0..exp(-1),
coords=cylindrical,labels=['x','y','z'],
axes='FRAMED',scaling='CONSTRAINED',grid=[30,5],
orientation=[45,75]);
Now let's plot the top of the surface:
>
plot3d([rho,theta,exp(-rho^2)],rho=0..1,theta=0..2*Pi,
coords=cylindrical,labels=['x','y','z'],
axes='FRAMED',scaling='CONSTRAINED',grid=[8,30],
orientation=[45,75]);
Putting the two pieces together in one picture, we obtain:
>
plot1:=plot3d([1,theta,z],theta=0..2*Pi,z=0..exp(-1),
coords=cylindrical,labels=['x','y','z'],
axes='FRAMED',scaling='CONSTRAINED',grid=[30,5],
orientation=[45,75]):
>
plot2:=plot3d([rho,theta,exp(-rho^2)],rho=0..1,theta=0..2*Pi,
coords=cylindrical,labels=['x','y','z'],
axes='FRAMED',scaling='CONSTRAINED',grid=[8,30],
orientation=[45,75]):
> display({plot1,plot2});