... työstettävää matskua ...
addpath c:\usr\heikki\opemate\complex x=[-1 0 1]; y=linspace(-pi/6,3*pi/4,20); m=length(y); n=length(x); [X,Y]=meshgrid(x,y); Z=X+i*Y; W=exp(Z); W=W+i*eps; tol=0.2 ztaso=[min(x)-tol,max(x)+tol,min(y)-tol,max(y)+tol]; Z=Z+i*eps; % Tällä niksillä varmistetaan, että alkioissa myös imag. osa. % plot ei nimittäin toimi samalla tavoin reaalisille alkioille. clf subplot(1,2,1) axis(ztaso) axis square hold on plot(Z(:,1),'b') plot(Z(:,2),'g') plot(Z(:,3),'c') for k=1:m plot(Z(k,:),'r') end title('z-taso') shg %%% z-taso valmis, w-taso alkaa %%%%%%%%%%%%%%%%%% subplot(1,2,2) axis square hold on plot(W(:,1),'b') plot(W(:,2),'g') plot(W(:,3),'c') for k=1:m plot(W(k,:),'r') end title('w-taso') shg