фильтрации, / — максимальная 
в плане 
протяженность 0б-
ласти С (звездочки у безразмерных переменных для удоб-
ства опущены).
Дифференциальные операторы системы (1.5.12) представ-
ляем 
конечно-разностными аналогами 
на сетке 
‚ При
этом — дифференциальные 
операторы Г) аппроксимируем
двухслойной шеститочечной схемой с весом о()
Г) 
—- А 
[ °) цуа Р (1 — ©) Яе
Г) &) 
7 А 
[ °) “лвуьн Р (1 — 9) 2а
| = Т,+о,5; /Рач 
°7 (Т:н0,5; /РО П1-0,5; /) н,7 
РО Т1-0,5; 
/ Ын,
() — 
62
Аг [7#
и Т. Д.
Таким образом, 
получим систему 
— конёчно-разностных
уравнений
ё То 
‚ З@) в+1 — бе в — д [ 
+ 
т
/гОНОрг(„) в+1 
Ах 
(@) 
Х а) в +-1
213^ 
Й
+(1__оа)2‹а›й]___.___._2_ё_ 
| — — © *+5
Но 
® (а) в +1
21° 
21°
т 
2 -/.іг 1° 
5 Ю
/гОНО 
+ 
/гОНО 
+
й 
— й
(а) 
& +1 
(2)® _ _
- 
Ах 
— А(ч) [ 
(а)&+1 + (1 ) (а)’г] Т
+121‚/г 
1 
йа) & +1 
&оЁ? 
х
То!о 
® (а) в +-1 
2ТоНо $
х (З (а) 
р &+1 — 4'1[7'(0‹›17---—1‚іг+1 + 
(а) р——-2‚іг+1)
т? Т _ 
Н а) в +1 — /Т(а)в
аП Р‚[‘З 
Ах
‚° (1.5.13)
Н(а)!—-].. Ё+1 ° 2Н(Ц) , Ё1 + Н(а) 1+-1, & +1
‹62
1
р
Так 
как эта система 
нелинейна, 
то 
для 
нее 
применим
итерационный процесс, основанный на методе квазилинеари-
зации нелинейных членов [42]. Согласно этому методу не-
линейные члены системы $ (Л, &) представляются в виде
н 2—о(3, )+ (в- 2ЕВ + (—а Ь,
(1.5.14)
63