в обратном случае, когда

На практике установлено, что выгоднее (численное решение быстрее сходится к точному) чередоватьчерез шаг по времени. В большинстве моделируемых случаев поглощение малозначительно, поэтому для определенияиспользуется приближение. Значения т для нечетного и четного шагов по времени шагов

(10)

где N - число делений разностной сетки.

3.3.6. Дискретная модель для диффузии и поглощения

Горизонтальные и вертикальные линии сетки нумеруем целыми индексами n и m соответственно. Соответственно узлы сетки нумеруем парой чисел (n,m), горизонтальные отрезки - (n+1/2,m), вертикальные отрезки (n,m+1/2).

Определим также два типа ячеек сетки. Ячейка (n+1/2,m+1/2) ограничена ребрами сетки, а ячейка с целыми индексами (n,m), заштрихованная на рис. 9, имеет в качестве узлов центры полуцелых ячеек.

Построим оператор An, аппроксимирующий оператор A в уравнении (4), в проинтегрированном по ячейке (n,m) виде. Считаем функциюопределенной в узлах сетки. Соответствующие значения обозначаемПоскольку сетка прямоугольная, расстояние между узлами (n+1,m) и (n,m), то есть шаг сетки по x, не зависит от y, и шаг сетки по y не зависит от x. Поэтому шаги можно пронумеровать одним полуцелым индексоми