艾里涡旋………

%% 子程序调用
function y = AI( params, r, phi )
r0=params(1);
rk1=params(2);
rk2 =params(3);
b= params(4);
w0=params(5);
m= params(6);
l= params(7);
a= params(8);
y= airy((r0-r)./(b.*w0)).exp(a.(r0-r)./(b.*w0)).exp(-(r0-r).^2./w0.2).((r.*exp(i.*phi)-rk1)./w0).^{m.*((r.*exp(i.*phi)+rk2)./w0).}l;
end