两次使用Matlab的“ Int”：
符号x1 y1 x y z lambda
f = exp（ - （x1^2+x2^2）/a^2）*exp（1i*pi/（lambda*z）*（（x-x1）^2+（y-y-y1）^2））;
e = int（int（f，x1，-3*a，3*a），x2，-3*a，3*a）
或者：
符号x1 y1 x y z lambda
f1 = exp（-x1^2/a^2+1i*pi/（lambda*z）*（x-x1）^2）;
f2 = exp（-y1^2/a^2+1i*pi/（lambda*z）*（y-y-y1）^2）;
e = int（f1，x1，-3*a，3*a）*int（f2，x2，-3*a，3*a）
最好的祝愿
托斯滕。
