In the present work, a compact, fast, and yet accurate algorithm is developed to calculate numerical values of Ambartsumian-Chandrasekhar’s \(H\)-function for isotropic scattering its moments on basis double exponential (DE) formula Takahashi Mori (RIMS, Kyoto Univ., 9:721, 1974). The main improvement made in new method an elimination iterative procedure automatic adjustment step-size integrati...