Home > atmlab > mie > besselplot3.m

besselplot3

PURPOSE ^

Computation and plot of Inverse Products of Riccati-Bessel Functions

SYNOPSIS ^

function result = besselplot3(n, m, xmin, dx, nx)

DESCRIPTION ^

 Computation and plot of Inverse Products of Riccati-Bessel Functions 
 of Order n for complex argument z=m*x, used in Mie Theory, 
 input: order n, refractive index m, minimum x value xmin, 
 x interval dx, number of x values nx.
 C. M�zler, August 2002

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

besselplot3.m

SOURCE CODE ^

0001 function result = besselplot3(n, m, xmin, dx, nx)
0002 
0003 % Computation and plot of Inverse Products of Riccati-Bessel Functions
0004 % of Order n for complex argument z=m*x, used in Mie Theory,
0005 % input: order n, refractive index m, minimum x value xmin,
0006 % x interval dx, number of x values nx.
0007 % C. M�zler, August 2002
0008 
0009 
0010 m1=real(m);m2=imag(m);
0011 nn=(1:nx);
0012 x=xmin+dx*nn;
0013 nu=n+0.5;
0014 z=m.*x;
0015 sqz= sqrt(0.5*pi*z); 
0016 pz = besselj(nu, z).*sqz;       % Psi_n Function
0017 chz = -bessely(nu, z).*sqz;     % Chi_n Function
0018 p1z= besselj(nu-1, z).*sqz;     % Psi_n-1 Function
0019 ch1z= -bessely(nu-1, z).*sqz;;  % Chi_n-1 Function
0020 pc=1./(pz.*chz);
0021 dz1=real(pc);
0022 dz2=imag(pc);
0023 r=[dz1;dz2];
0024 
0025 plot(x,r(1:2,:))
0026 legend('real(1/(Psi_n(mx)*Chi_n(mx)))','imag(1/(Psi_n(mx)*Chi_n(mx)))')
0027 title(sprintf('Inverse product of Riccati-Bessel Functions of Order n=%g, for m=%g+%gi',n,m1,m2))
0028 xlabel('x')
0029 result=[x;r];

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