15-Digit Accuracy Calculations of Chandrasekhar's $H$-function for Isotropic Scattering by Means of the Double Exponential Formula

Abstract

This work shows that it is possible to calculate numerical values of the Chandrasekhar $H$-function for isotropic scattering at least with 15-digit accuracy by making use of the double exponential formula (DE-formula) of Takahashi and Mori (Publ. RIMS, Kyoto Univ. Vol. 9, 721, 1974) instead of the Gauss-Legendre quadrature employed in the numerical scheme of Kawabata and Limaye (Astrophys. Space Sci. Vol. 332, 365, 2011) and simultaneously taking a precautionary measure to minimize the effects due to loss of significant digits particularly in the cases of near-conservative scattering and/or errors involved in returned values of library functions supplied by compilers in use. The results of our calculations are presented for 18 selected values of single scattering albedo ${\varpi }_0$ and 22 values of an angular variable $\mu$, the cosine of zenith angle $\theta$ specifying the direction of radiation incident on or emergent from semi-infinite media.

15 pages, 8 tables, The article is available electronically: http://link.springer.com/article/10.1007/s10509-016-2961-3

References in corpus (1)

• A Fast Iterative Method for Chandrasekhar's H-functions for General Laws of Scattering

Cited by in corpus (1)

• 15-Digit Accuracy Calculations of Ambartsumian-Chandrasekhar's $H$-Functions for Four-Term Phase Functions with the Double-Exponential Formula