ТОР 5 статей: Методические подходы к анализу финансового состояния предприятия Проблема периодизации русской литературы ХХ века. Краткая характеристика второй половины ХХ века Характеристика шлифовальных кругов и ее маркировка Служебные части речи. Предлог. Союз. Частицы КАТЕГОРИИ:
|
Метод Рунге-Кутта четвертого порядка
В вычислительной практике наиболее часто используется метод Рунге-Кутта четвертого порядка. В этом методе величины yi+1 вычисляются по следующим формулам:
yi+1 = yi + Dyi Dyi=h×(k1 + 2k2 + 2k3 + k4)/6, i = 0, 1,... k1 = f(xi, yi), k2 = f(xi+h/2, yi+h×k1/2), k3 = f(xi+h/2, yi+h×k2/2), k4 = f(xi+h, yi+h×k3).
Листинг программы метода Рунге – Кутта четвертого порядка:
Программа 5 var k1,k2,k3,k4,x,y,a,b,h,d:real; function f(x,y:real):real; begin f:= 2 – x – y; end; begin writeln('Введите y, a, b, h'); readln(y,a,b,h); x:=a; repeat writeln(x:0:3,' ',y:0:3); k1:=f(x,y); k2:=f(x+h/2,y+h*k1/2); k3:=f(x+h/2,y+h*k2/2); k4:=f(x+h,y+h*k3); d:=(k1+2*k2+2*k3+k4)/6; y:=y+d; x:=x+h; until x>b+0.1; readln; end.
Блок схема программы 2:
Демонстрация рабочей программы 2:
Таблица - 2
График программы 2:
Не нашли, что искали? Воспользуйтесь поиском:
|