Abstract
We present a method for the numerical solution of the kinetic equations for a gas composed by photons, electrons, atoms and ions. The gas is assumed to satisfy the statistical hypothesis. We show the integro-diferential equations that determine the distribution functions, for situations departing from thermodynarnical equilibrium as in stellar atmospheres. These functions differ from Boltzrnann' s and even from Maxwell's function. We give the equations for a one-dimensional problem and propose the use of the Newton-Raphson method to solve the equations for given boundary conditions. We also show how to compute first order devia tions from Maxwell's distribution, and, from these departures, how to compute the transport coefficients and their range of aplicability. We further suggest correction procedures for saturated fluxes.
