function tiaoxing(bochang,N,m,L,handles)
O=jLZ2os global n
|0[Buh[_:c global point
<{#_;7h" k=2*pi/bochang;
\OW:- a=m*1e+006;
3X gJZ
deltax=2*a/point;%取單位長(zhǎng)度
x0# Bc7y xx=-a:deltax:a;
xSudDhRP L=L*1e+006;
6l\5J6x nn=length(xx);
Lap?L/NS u=ones(1,nn);
&l+Qn'N jj=1:nn;
U<'N=#A
J for n=1:N %進(jìn)行迭代
UyRy>:n c5^HGIe1 for b=1:nn
Jj=qC{] uu(b)=sqrt(i*exp(-i*k*L)/(bochang*L)).*...
6 - 3?&+ sum(exp(-i*k*((b-jj)*deltax).^2./(2*L)).*u)*deltax;%
HTL6;87w+] &qbEF3p^@ end
it}h8:^< %Xh}{