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）];