Рассмотрим методику численного решения интегро-дифференциального уравнения (14.2), описывающего процесс агрегации включений. В левой части данного уравнения находятся две производные, описывающие изменение функции распределения включений по размерам f во времени и пространстве; правая часть, содержащая два интеграла, характеризует изменение функции f за счёт процесса агрегации включений*. Таким образом, в данном случае функция f является функцией трёх независимых переменных - t, x и r. Также в правой части уравнения (14.2) имеется вероятностная функция k, значения которой зависят от размера включений.
     Таким образом, для записи разностной схемы, аппроксимирующей уравнение (14.2), требуется ввести трёхмерную разностную сетку (по осям которой отложены независимые переменные - время t, координата по длине реактора x и размер включений r), а также обозначения, аналогичные обозначениям (14.8).
     Для аппроксимации производной по времени будем использовать, как и обычно, правую конечную разность. Выбор конечной разности для аппроксимации производной по координате определяется (согласно известному правилу) знаком величины v, характеризующей скорость перемещения включений. Для аппроксимации интегралов будем использовать выражение (14.4) и принцип замороженных коэффициентов. Учитывая всё сказанное, запишем схему расщепления, аппроксимирующую уравнение (14.2) для случая v > 0:

(14.10)