clear all; close all; [x,y] = meshgrid(-3:0.1:3) z=x+i*y; fz=1+1./z; mesh(x,y,abs(fz)) zlim([0,2]) surf(x,y,abs(fz),'FaceColor','interp',... 'EdgeColor','none',... 'FaceLighting','phong') %daspect([2 2 1]) %axis tight %view(-50,30) colormap hsv camlight left zlim([0,2]) hold on omega=0:pi/20:2*pi plot3(cos(omega),sin(omega),0*omega) plot3(cos(omega),sin(omega),abs(1+exp(-i*omega)),'LineWidth',4) figure plot(omega,abs(1+exp(-i*omega)),'LineWidth',4) grid on