该双重积分可以用两次int（）积分函数来求解，得到其解析值。
求解过程如下：
>> lambda=632.8;z=300e6;k=2*pi/lambda;
>> syms x y u v
>> E=exp(i*k*z)/(i*lambda*z)*int(int(exp((i*k)/(2*z)*((x-u)^2+(y-v)^2)),u,-4e6,4e6),v,-4e6,4e6)
