ソースコード:

clear;
clc;
clf;
figure(1);
hold on;
lambda=632.8*10^-3;%um
L=100*lambda;
k=2*pi/lambda;
a=25*lambda;
C=sqrt(1i/(lambda*L)*exp(-1i*k*L));
for x=-200:0.1:200 
    u=1;
    fun=@(x0) exp(-1i.*k.*((x-x0).^2./(2.*L)));
    for j=1:1
        u=integral(@(x0) C*fun(x0)*u,-a,a);
    end
    scatter(x,abs(u))
end

計算方法:

計算结果: