mie散射系数

function resault = Mie_abcd（mx）
nmax = round（2+x+4*x.^（1/3））;
n = （1:nmax）;nu = （n+0.5）;z = m.*x;m2 = m.*m;
sqx = sqrt（0.5*pi./x）;sqz = sqrt（0.5*pi./z）;
bx = besselj（nux）.*sqx;
bz = besselj（nuz）.*sqz;
yx = bessely（nux）.*sqx;
hx =bx+1i*yx;
b1x =[sin（x）/xbx（1:nmax-1）];
b1z =[sin（z）/zbz（1:nmax-1）];