Home > atmlab > mie > mie_pt.m

mie_pt

PURPOSE ^

pi_n and tau_n, -1 <= u <= 1, n1 integer from 1 to nmax

SYNOPSIS ^

function result=mie_pt(u,nmax)

DESCRIPTION ^

 pi_n and tau_n, -1 <= u <= 1, n1 integer from 1 to nmax 
 angular functions used in Mie Theory
 Bohren and Huffman (1983), p. 94 - 95

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

mie_pt.m

SOURCE CODE ^

0001 function result=mie_pt(u,nmax)
0002 % pi_n and tau_n, -1 <= u <= 1, n1 integer from 1 to nmax
0003 % angular functions used in Mie Theory
0004 % Bohren and Huffman (1983), p. 94 - 95
0005 
0006 p(1)=1; 
0007 t(1)=u;
0008 p(2)=3*u; 
0009 t(2)=3*cos(2*acos(u));
0010 for n1=3:nmax,
0011     p1=(2*n1-1)./(n1-1).*p(n1-1).*u;
0012     p2=n1./(n1-1).*p(n1-2);
0013     p(n1)=p1-p2;
0014     t1=n1*u.*p(n1);
0015     t2=(n1+1).*p(n1-1);
0016     t(n1)=t1-t2;
0017 end;
0018 
0019 result=[p;t];

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