gap=0.01; %gap between data points
x=-1:gap:1;
y=1:-gap:-1;

[xs,ys]=meshgrid(x,y);
%ys=y'*(zeros(1,numel(x))+1);
%xs=(zeros(1,numel(y))+1)'*x;
zs=xs+1i*ys; %sets up complex matrix with entries like complex grid

%functions to plot
f1=@(z) (1 + 25*z.^2).^(-1); 
f2=@(z) exp(-1./z);

z1=f1(zs);
z2=f2(zs);

figure;
surf(x,y,abs(z1));
axis([-1 1 -1 1 0 10]);
xlabel('Re');
ylabel('Im');
print -depsc 1a.eps;

figure;
surf(x,y,abs(z2));
axis([-1 1 -1 1 0 10]);
xlabel('Re');
ylabel('Im');
print -depsc 1b.eps;

figure;
surf(xs, ys, 0*xs, z2rgb_kawski(z1)); view(2); shading interp;
xlabel('Re');
ylabel('Im');
print -depsc 1c.eps;

figure;
surf(xs, ys, 0*xs, z2rgb_kawski(z2)); view(2); shading interp;
xlabel('Re');
ylabel('Im');
print -depsc 1d.eps;
