Details are provided for two subroutines with which one can compute the various characteristics of the electromagnetic radiation scattered by an absorbing, homogeneous sphere of any reasonable size. The necessary expressions for this purpose were first derived byM ie. The method of computations used is the so-called method of logarithmic derivative of one of the complex functions, introduced by Infeld. The main difference between the two subroutines is in the procedure used in computations of one of the functions. This function is computed by an upward recurrence procedure in one subroutine and by a downward recurrence procedure in the other. Sufficient results for demonstrating the reliability of these programs are presented and discussed for a sphere of 10µm radius illuminated by an unpolarized radiation of 0.4µm wavelength.
[1]
E. Wolf,et al.
Principles of Optics
,
2019
.
[2]
Wilfried Heller,et al.
Angular Scattering Functions for Spheres
,
1967
.
[3]
David M. Miller,et al.
Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55)
,
1965
.
[4]
F. J. Corbató,et al.
Generation of Spherical Bessel Functions in Digital Computers
,
1959,
JACM.
[5]
H. V. Hulst.
Light Scattering by Small Particles
,
1957
.
[6]
C. Sliepcevich,et al.
Tables of functions first and second partial derivatives of Legendre polynomials
,
1951
.