получим
д 
= бы(М 
-) а). в 
= @ы () 
Вы) + 4)).
{7 
,, 
„) — 
К,
( 4; Вн, ) + 
)) () СЛ )+
С 
— + ( @; С1-1,; + Л, ) )(”Ё.і _“Ё‚і АЁ—Ы)
7 ° 
[\›і.] 
’
а’ 
— С, ; ( & ; —) Ан,7) . В 
— ©)( 4) ВнЕ 
4)) .
і’] 
Кі.] 
› 
і'] Н 
К‘‚і 
›
( 4; В-) 
+ )) ( 4,7 С,7 + Л, 
;) +
С° = +( @; С,7 Л у) ( ©) — 
7 Аа ))
в) — 
В/ 
›
К, ‚ = ( 6) — @, ‚ Ан. 
)) (ді‚і — @,, Аы.і) т
+ (а В и) + @;) (а„ В„ +4,)). — (1.5.21)
При этом значения А,); В,;; Соу; А,у; Во 
;; С)о,у
трудно получать из краевых условий.
Вычисления производятся следующим образом.
По рекуррентным соотношениям (1.5.17) находим 4, 
;
В, (1=1, р); затем по формулам (1.5.21) для всех # и / —
,
А,у; В, 
Сиу: А„ув В„уз С,. По формулам (1.5.20),
(1.5.16) определяем все 2 ›; /1„ и Н, ле Аналогичным 0б-
разом вычислительный процесс строится для а — 2.
Предложенный метод обеспечивает совместное решение ли-
нейных конечно-разностных систем, не прибегая к итерацион-
ному процессу, как это предлагается в работе [112]. Согласно
[112] каждое из уравнений конечно-разностной системы в обо-
их координатных направлениях решается независимо друг от
друга. При этом неизвестные значения переменных задаются
из предыдущей итерации. Такое построение итерационного
процесса помимо нерациональной затраты машинного вре-
мени приводит к тому, что несмотря на применение абсолют-
но устойчивых схем для аппроксимации каждого дифферен-
циального уравнения системы, вычислительный процесс для
системы в целом сходится только при специальных значениях
шага по времени.
Разработанный алгоритм реализован в виде стандартной
программы для ЭВМ типа М-20, имеющих два куба МОЗУ.
Предлагаемый алгоритм и соответствующая ему програм-
ма опробированы 
нами 
на конкретном 
примере, 
т. е. для
не-
5—242 
65