# 
restart;
f:=z->sin(z);
with(plots):
Re(1+I);
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);
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);
display([P1,P2]);

#