Se presenta un programa de cálculo que permite resolver las ecuaciones de equilibrio estadístico para las poblaciones de los tres primeros niveles del átomo de hidrógeno y protones. El programa se ha construido para una atmósfera de hidrógeno irradiada por una fuente de radiación externo. El método se ha comprobado para densidades entre 10¹³ y 10¹⁰ cm⁻³ y temperaturas del orden de 10³-10⁴ K.