clear;
dl=0.002;?
lamuda=0.02;
l=（lamuda*3）/2;
a=0.0001;k?=2*pi/lamuda;
n=l/dl;
z=（-l/2+dl/2）:dl:（l/2-dl/2）;
Z=zeros（n）;
zzz=120*pi;
for?h=1:n
????for?l=1:n
????????xx1=exp（-j*k*sqrt（a^2+（abs（h-l）*dl）^2））;
????????xx2=sqrt（a^2+（abs（h-l）*dl）^2）;
????????qqq=（sqrt（dl^2+4*a^2）+dl）/（sqrt（dl^2+4*a^2）-dl）;
????????znn1=（1/（4*pi））*log（qqq）;
????????znn2=（1/（4*pi