Введение в численные методы. Гришанова Л.И - 46 стр.

UptoLike

Составители:

Рубрика:

согласно рис. 3.5, б), а также шаг H. Выходными являются парамет$

ры XN, YN в конце интервала.

Первое обращение к подпрограмме RK при входных параметрах

X, Y, N дает выходные параметры X2, U (см. рис. 3.5, б). Второе

обращение к RK, отличающееся от первого значением шага H/2, по$

рождает промежуточные выходные параметры X1, V в середине ин$

тервала. Третье обращение с выходными параметрами X1, V, H/2

дает выходные параметры X2, W (конец интервала – см. рис. 3.5, б).

Затем вычисляем модуль разности значений U (шаг H) и W (шаг

H/2), т. е. D = |U–W| (блок 4). Если D>EPS (блок 5), то точность еще

не достигнута, и, присвоив U = V, уменьшаем шаг вдвое H = H/2

(блок 8). Снова повторяем операции блоков 1–4 при H = H/2. Если

D £ EPS, то вычисляем окончательно Y = W–(U–W)/15 и присваива$

Рис. 3.8. Фрагмент схемы алгоритма метода Рунге – Кутты