и т. д. Этот алгоритм повторяется до тех пор, пока не вы-
полнится условие
шах',д(р) _ (9—1›|<е
где е `>> () — малое.
Если последовательность { 9”,} начинает 
„колебаться“
(осциллировать), начиная с некоторого р, то 
целесообразно
применить операцию демпфирования —
() 
(р) 
(2-1).
Ора ° 
„ + (1 — ›) ор 
;
здесь ©, „ — значение искомой функции, 
полученное опи-
санным. выще алгоритмом.
Для оценки точности численного решения. применяется
метод балансовых уравнений, введенных в $ 2 данной главы.
Интегрируя уравнение (1.4.5) по Е в интервале (0, 1), по-
лучеем.
а0. 
аис
аір = — & ©* + 9„ (5) — 1 д.ср — 4 (*),
ИЛИ
‘0ср ('С) — — атиср (Т) —- ЕГ е*‘с + уаисп (’Ё› Сі"ъ' —
0
— 3 9 (*) @* + о, (0) + стар (0),
где 
'
1
оер (т)=3` & (& <) ах,
1
йср () = || # (г <) ак,
0
кр
7
1 д'ш(ё‚ — т)
да (=) = [Е — — —- 
&*
0
1 
п
1 —
на () = ® |+(1 - т 
) &
0
Для иллюстрации рассмотрим задачу со следующими дан-
ными: В, ==500 м; т, = 40 м; т = 100 м; т 
= 10 м;
Но = 39 м; л = 0, 1, 2, 3; д9 = 0; 0,0036; 0,036;
55