... 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