Home > atmlab > mie > besselriccati.m

besselriccati

PURPOSE ^

Computation of Riccati-Bessel Functions of Order n

SYNOPSIS ^

function result = besselriccati(n, m, x)

DESCRIPTION ^

 Computation of Riccati-Bessel Functions of Order n
 for complex argument z=m*x, used in Mie Theory. 
 input: order n, refractive index m, x value, 
 C. M�zler, August 2002

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

besselriccati.m

SOURCE CODE ^

0001 function result = besselriccati(n, m, x)
0002 
0003 % Computation of Riccati-Bessel Functions of Order n
0004 % for complex argument z=m*x, used in Mie Theory.
0005 % input: order n, refractive index m, x value,
0006 % C. M�zler, August 2002
0007 
0008 m1=real(m); m2=imag(m);
0009 nu=n+0.5;
0010 z=m.*x;
0011 sqx= sqrt(0.5*pi*z); 
0012 psz = besselj(nu, z).*sqx;
0013 chz = -bessely(nu, z).*sqx;
0014 dpic=psz-i*chz;
0015 z2=0.5*z.*z;
0016 n2=2*n;
0017 A1=4./((n2-1).*(n2+3));
0018 lez=(2*n+1)./z./(1+z2.*A1); % Low-frequency approximation of Ez
0019 ez=1./(chz.*psz);
0020 fz=psz./chz/i-1;
0021 result=[psz;chz;dpic;ez;fz;];

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005