This routine solves Burgers' equations for partly ionized elements for arbitrary
compositons. A modified version of this routine which allows to use plasma-parameter
dependent resistance coefficients (e.g., from package below) is available here.
This routine computes plasma-parameter dependent resistance coefficients needed for Burgers' scheme using fitting
formulae to the collision integrals of Paquette et al., ApJS 61, 177 (1986). In addition, quantum
corrections of Hahn et al., Phys. Fluids 14, 278 (1971) are included utilizing again appropriate fitting
formulae.