IMSL MATH/LIBRARY, Fortran Subroutines for Mathematical Applications, Version 2.0; IMS: Houston, TX, 1991;DIVPAG, pp. 755–771 to solve an initial-value problem for ordinary differential equations using Adams-Moulton’s method and DRUN, pp. 1319–1320 to generate random numbers from a uniform (0,1) distribution
DRFM: a New Package for the Evaluations of Gas-Phase-Transport Properties