lmie.exp: Calculates the Lorentz-Mie coefficients a_n and b_n.

In wendellopes/rvswf: RVSWF: Expansion in Vector Spherical Wave Functions written in R

Description

Calculates the Lorentz-Mie coefficients a_n and b_n.

Usage

lmie.exp(m, x, by = "LD", ...)

Arguments

m	The ratio between the refractive indices.
x	The form factor value.
by	LD using logarithmic derivative; RB using the ratio between Ricatti-Bessel functions.

Details

Parser to the two ways to calculate the Lorentz-Mie coefficents.

Value

The coefficients a_n=C_nT_a and b_n=C_nT_b and also the coefficients C_n, T_a and T_b.

Examples

# Table 4.1 - Bohren and Hoffman Book, pg 114
# The values showed bellow are conjugated. Probably
# there is tipography error in the book.
m<-1.33+1e-8i
x<-3
n<-c(1,2,3,4,5)
an<-c(5.1631e-1+4.9973e-1i, 3.4192e-1+4.7435e-1i,4.8467e-2+2.1475e-1i,
      1.0346e-3+3.2148e-2i, 9.0375e-6+3.0062e-3i)
bn<-c(7.3767e-1+4.3990e-1i, 4.0079e-1+4.9006e-1i,9.3553e-3+9.6269e-2i,
      6.8810e-5+8.2949e-3i, 2.8390e-7+5.3204e-4i)
bh<-data.frame(n,an,bn)
a<-lmie.exp(m,x,by="LD",NMAX=5)
b<-lmie.exp(m,x,by="RB",NMAX=5)
print(a)
print(b)
print(bh)

wendellopes/rvswf documentation built on May 4, 2019, 4:19 a.m.