Блок-схема программы вычисления решения дифференциального уравнения по методу Рунге-Кутта.

Начало

Описание Описание

процедуры- процедуры-

-функции Ввод х₀, у₀, -функции

rk(x₀, y₀, h, m) h₀, E, b y′ =f(x, y)

x₀, y₀ – начальные условия

h0 – шаг таблицы

n = ] (b- x₀ )/ h₀[ + 1 E погрешность

b – конец отрезка

интегрирования [a, b]

i = 1, n n – число точек в таблице

решений

] [ - целая часть числа

h = h₀

m = 1

y = rk() Вызов процедуры-функции

y = rk(x₀, y₀, h, m)

y1 = y, h = h/2

x = x₀, y = y₀, m = 2m

Вызов процедуры-функции

y = rk() y = rk(x₀, y₀, h, m)

да |y-y1|> 15E/16

нет

Вывод

x₀, y₀, h, m

x₀ = x₀ + h₀, y₀ = y

end