This paper shows some analytic approximations of the H-function derived earlier and obtained by decomposition of the angular flux density of particles combined with the zero or der DPN method. A novel solution based on the modified DPN method has been presented too. The comparison with well-known Chandrasekhar’s results obtained by numerical procedure and with Pomraning’s analytic approximations has been performed.