По рекуррентным формулам (1.2.19)—(П.2.21), (1.2.23) 
—
— (1.2.25), (1П.2.29) — (1.2.33) последовательно вычисляются:
А 
А!‘‚‚г ‚ В 
В,‘_„;
9, В, ... ° 
2е .. ›
С‚ 
‚ В!‚’д;
С "Д к› 
э„ Ё’ Р, ..
Е в› Е ‚ &* Рг‚к‚...› Рл‚л-
л
Такую последовательность 
вычислений 
назовем прямой
прогонкой. В силу второго условия (П.2.11) и (1.2.28) имеем
Рп‚ К
Тогда из уравнения (1.2.28) и (1П.2.34) можно последователь-
но 
определить 
Л, 1»› Йдозко * 
Й ю йо в 
Далее, 
из
уравнений (1.2.18) и (1.2.22) 
и условий (1.2.10) и (1.2.14)
последовательно 
вычисляются 
Ир ь Инр 
в› -) Мрк Ф›
@_ к› % ‚ Эту последовательность 
вычислений назовем
обратной прогонкой.
Следовательно, такая модификация 
позволяет численно
решить систему уравнений (1.2.5) — (П.2.7) с начальными и
краевыми условиями (1.2.8) 
— (1.2.14) для любого момснта
времени т,, & = 1, 2, ... .
В этой схеме 
численного интегрирования есть один не-
достаток, 
а именно: 
при разбиении отрезка (0, 1) точками
5 в число узловых точек должна входить точка $. ПоЭ-
тому в окрестности галереи будет иметь место неравномер-
ный шаг разбиения. Этот 
недостаток метода можно устра-
нить, если вместо искомой функции Й/ ввести новую функцию:
о= й- 0*5 (-%'——аіп : ) 
(1.2.35)
0, ож ЕЕ
б 1, ЕЖЕК1.
Тогда формулы (1.2.18) 
— (1.2.21), (1.2.22) 
— (1.2.25) сох-
раняют силу. А вместо формул (1П.2.28) 
— (1.2.33) будем
иметь
= Ми вОр в 'Ё М 
(1.2.36)
б + 
25 А
‘Оі‚іг
М
’
й — 
А 
А
1+ @— Мув) б НОр М — оа, [(*—Ацнаа) 
Х
Ах
х Ак — 3] — 
тодп, [( 
— 
Сна в) Ст„ в —3]
(1.2.37)
12}