clear all;

s=10;
dx=0.1;
dy=dx;

x=([0:dx:4]-1e-3)*1e0;
szx0=size(x);
szx=szx0(2);
y=[-4:dy:4]*1e1;
szy0=size(y);
szy=szy0(2);

z=zeros(szy,szx);
for j=1:szx,
    for k=1:szy,
        z(k,j)=x(j)+y(k)*i;
    end;
end;

fz=log(z)-z;
Fz=z.^s.*exp(-z.*s);
L=abs(Fz).^2;

%=====================================
figure(1);
clf;

subplot(2,3,1);
surf(x,y,real(fz));
shading flat;
xlabel('x');
ylabel('y');
zlabel('u');

subplot(2,3,2);
surf(x,y,imag(fz));
shading flat;
xlabel('x');
ylabel('y');
zlabel('v');

%subplot(2,3,3);
%surf(x,y,(abs(fz))+1e-4);
%shading flat;
%xlabel('x');
%ylabel('y');

subplot(2,3,4);
surf(x,y,real(Fz));
shading flat;
xlabel('x');
ylabel('y');
zlabel('Re(F)');

subplot(2,3,5);
surf(x,y,imag(Fz));
shading flat;
xlabel('x');
ylabel('y');
zlabel('Im(F)');

subplot(2,3,6);
surf(x,y,log10(L));
shading flat;
xlabel('x');
ylabel('y');
zlabel('L=|F|^2');


pause;
subplot(2,3,4);
plot(x,real(Fz((szy+1)/2,:)));
xlabel('x');
ylabel('Re(F)');

pause;
subplot(2,3,5);
plot(x,imag(Fz((szy+1)/2,:)));
xlabel('x');
ylabel('Im(F)');



   