Главная страница Случайная страница
КАТЕГОРИИ:
АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника
Блок-схема программы вычисления решения дифференциального уравнения по методу Рунге-Кутта.
|
|
 |
Начало
Описание Описание
процедуры- процедуры-
-функции Ввод х0, у0, -функции
rk(x0, y0, h, m) h0, E, b y′ =f(x, y)
x0, y0 – начальные условия
h0 – шаг таблицы
n = ] (b- x0 )/ h0[ + 1 E погрешность
b – конец отрезка
интегрирования [a, b]
 |
i = 1, n n – число точек в таблице
решений
] [ - целая часть числа
h = h0
m = 1
y = rk() Вызов процедуры-функции
y = rk(x0, y0, h, m)
 |
y1 = y, h = h/2
x = x0, y = y0, m = 2m
 |
Вызов процедуры-функции
y = rk() y = rk(x0, y0, h, m)
 |
да |y-y1|> 15E/16
нет
 |
Вывод
x0, y0, h, m
|
x0 = x0 + h0, y0 = y
end
Блок-схема процедуры-функции метода Рунге-Кутта.
 |
rk(x, y, h, m) x, y – начальная точка
h – шаг интегрирования
m – число дроблений шага
| |||
 |
j = 1, m
Вычисление у
по формулам
Рунге-Кутта
x = x + h
rk = y Значение решения ОДУ
в очередной точке.