35304

Тема: Метод РунгеКутта вирішення задачі Коші.

Практическая работа

Информатика, кибернетика и программирование

h void min {int in; double hbdely; doublek = new double [4]; doublex=new double [n]; doubley=new double [n]; cout Vvedite bh endl; cin ; cin b; cin h; cout Vvedite y[0] endl; cin y[0]; n=b h; x[0]=; cout x y ; cout endl; cout ; cout endl; fori=0;i =n;i { k[0]=hx[i]x[i]2y[i]; k[1]=hx[i]h 2x[i]h 22y[i]k[0] 2; k[2]=hx[i]h 2x[i]h 22y[i]k[1] 2; k[3]=hx[i]hx[i]h2y[i]k[2]; dely=k[0]2k[1]2k[2]k[3] 6; x[i1]=x[i]h;...

Украинкский

2013-09-09

38 KB

0 чел.

обец А., гр. 1ПМ-06  Вариант 7

17.03.09

Лабораторна робота №26-27

Тема: Метод Рунге-Кутта вирішення задачі Коші. Складання алгоритму.

Мета: Навчитися вирішувати задачу Коші методом Рунге-Кутта; скласти алгоритм.

Устаткування: папір формату А4, ПК, програмне забезпечення Borland С++.

  1.  Індивідуальне завдання.

Вирішити задачу, використовуючи метод Рунге-Кутта для диференціального рівняння  с початковою умовою , обчислити з точністю до  рішення цього рівняння при , . Обчислення провести із двома запасними знаками.

7)

#include<iostream.h>

#include<math.h>

void main()

{int i,n;

double h,a,b,dely;

double*k = new double [4];

double*x=new double [n];

double*y=new double [n];

cout<<"Vvedite a,b,h"<<endl;

cin>>a;

cin>>b;

cin>>h;

cout<<"Vvedite y[0]"<<endl;

cin>>y[0];

n=(b-a)/h;

x[0]=a;

cout<<"    x    |   y  ";

cout<<endl;

cout<<"-------------------";

cout<<endl;

for(i=0;i<=n;i++)

{

 k[0]=h*(x[i]*x[i]+2*y[i]);

 k[1]=h*((x[i]+h/2)*(x[i]+h/2)+2*(y[i]+k[0]/2));

 k[2]=h*((x[i]+h/2)*(x[i]+h/2)+2*(y[i]+k[1]/2));

 k[3]=h*((x[i]+h)*(x[i]+h)+2*(y[i]+k[2]));

 dely=(k[0]+2*k[1]+2*k[2]+k[3])/6;

 x[i+1]=x[i]+h;

 y[i+1]=y[i]+dely;

}

cout<<"  "<<x[n]<<"    |  "<<y[n]<<"\n";

 cout<<endl;

}

  1.  Контрольні питання
  •  Поставте задачу, розв’язувану методом Рунге-Кутта.
  •  Сформулюйте алгоритм метода Рунге-Кутта.
  •  Який порядок точності методу Рунге-Кутта?
  •  Як практично досягти заданого ступеня точності?
  1.  Звіт про роботу.
  2.  Захист роботи.

PAGE  1


 

А также другие работы, которые могут Вас заинтересовать

47638. Теории вычислительных процессов. Методические указания 819.5 KB
  Третья модель предполагает проведение оптимального синтеза системы РВ основываясь на оптимизации функций Лагранжа применительно к экспоненциальным стохастическим сетям. Синтез системы разделения времени по критерию минимума стоимости при заданном времени реакции системы на запрос пользователя. Синтез системы разделения времени по критерию минимума времени реакции системы на запрос пользователя заданной стоимости. На...
47640. Правотворчество. Правотворческая деятельность 41 KB
  Правотворческая деятельность – это форма деятельности компетентных органов государства, в ходе которой устанавливаются нормы права посредством издания, изменения или отмены правовых актов. Особо следует отметить, что государство не создает право
47641. Расчёт конденсатора паровой турбины (методические указания) 1.1 MB
  Целью расчета конденсатора является: определение геометрических размеров, режимных характеристик, характеристик конденсатора на переменном режиме, гидравлический расчет, а также механические расчеты на прочность основных элементов конденсатора
47644. ДЕРЖАВНА ЕКОЛОГІЧНА ПОЛІТИКА. РОБОЧА НАВЧАЛЬНА ПРОГРАМА 325.5 KB
  Ці та інші чинники зокрема низький рівень екологічної свідомості суспільства призвели до значної деградації довкілля України надмірного забруднення поверхневих і підземних вод повітря і земель нагромадження у дуже великих кількостях шкідливих у тому числі високотоксичних відходів виробництва. Охорона навколишнього природного середовища раціональне використання природних ресурсів забезпечення екологічної безпеки життєдіяльності людини невідємна умова сталого економічного та соціального розвитку України. В основу формування державної...