> restart;
 

> f:=z->sin(z);
 

proc (z) options operator, arrow; sin(z) end proc (1)
 

> with(plots):
 

> Re(1+I);
 

1 (2)
 

> P1:=contourplot(Re(f(x+I*y)),x=-Pi..Pi,y=-Pi..Pi,    contours=[0.1,-0.1,.5,-.5,.9,-.9,1.5,-1.5,2,-2,3,-3],    numpoints=1000,color=blue);
 

PLOT(CURVES([[.263033270396252, -3.14159265358979], [.283242623819873, -3.03213619571094]], [[.321506108199784, -2.94524311274043], [.283242623819873, -3.03213619571094]], [[.321506108199784, -2.94524... (3)
 

> P2:=contourplot(Im(f(x+I*y)),x=-Pi..Pi,y=-Pi..Pi,     contours=[0.1,-0.1,.5,-.5,.9,-.9,1.5,-1.5,2,-2,3,-3],     numpoints=1000,color=red);
 

PLOT(CURVES([[-3.14159265358979, -1.81496042220010], [-3.08827677241409, -1.82046174881995]], [[-2.94524311274043, -1.83270742790713], [-3.08827677241409, -1.82046174881995]], [[-2.94524311274043, -1.... (4)
 

> display([P1,P2]);
 

Plot_2d
 

>