restart;
p:=sqrt(1-x^2);
f:=p*sin(x)*cos(x);
g:=x*sinh(p)*cosh(p);



plot([f,g],x=-1..1);

sinh(1.0);
sin(1.0);


plot3d([r,theta,abs(sin(r*exp(I*theta)))], r=0..1, theta=0..2*Pi, coords=cylindrical);

plot3d([r,theta,abs( arcsin(r/2*exp(I*theta)))], r=0..1, theta=0..2*Pi, coords=cylindrical);

h:=z->z^2+2*z+I;

plot3d([r,theta,abs( h(r*exp(I*theta)))], r=0..1, theta=0..2*Pi, coords=cylindrical);