# 1D complex function plot -- djm -- 08 sept 2006
#
restart; with(plots):ang := Pi/3;# define POLYNOMIAL here;
P := (r*exp(I*ang))^3 -2*(r*exp(I*ang))^2 -2*(r*exp(I*ang)) - 3;# real-valued contours;
plotRe := plot(Re(P),r=-5..5,color=black):# imag-valued contours;
plotIm := plot(Im(P),r=-5..5,color=blue):display( [ plotRe, plotIm ]);