Тогда поток через правую границу ячейки (n,m), который мы сопос-

тавляем ребру (n+1/2,m), можно вычислить как

(13)

Площадь ячейки (n+1/2,m)

(14)

Для оператора поглощения простейшей аппроксимацией будет

rim_m „m cifi

^Sn°nPn . (15)

Объединяя формулы (13,15), получаем

(16)

Индексом h подчеркиваем, что р_п - это набор чисел (рЩ, используемый для аппроксимации значений р( x, y) в узлах сетки.

При постоянных коэффициентах и шагах сетки получается классическая пятиточечная разностная аппроксимация оператора Лапласа, дополненная младшим членом при а Ф 0.

3.3.7. Способ решения дискретных уравнений диффузии

Основные вычислительные трудности при решении уравнения (6) связаны с обращением оператора диффузии с поглощением А (16). Поэтому