“…The analytical consideration of the NE distribution function formation in the case of non-stationary, non-localized, mismatched on intensity sources and sinks is very problematic, thus numerical treatment must be used. For the numerical investigation of the problem, we develop the original conservative schemes whose efficiency is proved in different problems (see, for instance, [8,9]).…”