Главная

Популярная публикация

Научная публикация

Случайная публикация

Обратная связь

ТОР 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:

 

 

x y
   
   
  -140
  -245
  -346
  -446
  -546
  -646
  -746
  -846
  -946
  -46

Таблица - 2

 


 

График программы 2:

 






Не нашли, что искали? Воспользуйтесь поиском:

vikidalka.ru - 2015-2024 год. Все права принадлежат их авторам! Нарушение авторских прав | Нарушение персональных данных