Stabilized scheme for calculating radiation transfer in the $P_1$ approximation

We consider interpolation-characteristic scheme approximating approximation to the radiative transfer equations corresponding to the $P_1$ model. The explicit finite-difference scheme is corrected by special term adjusting the rate of radiation energy transfer. Such correction can reduce the influence of non-physical effects when calculating radiative heat transfer in a medium with non-uniform opacity.