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


 

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

53979. Розвиток мовлення. Переказ 57.5 KB
  Мета: вчити учня розповідати в усній і писемній формі про те, що бачив, про що читав, що знає, добирати для розповіді найбільш вдалі слова, удосконалювати навички використання іменників різних відмінків у зв’язному мовленні.
53981. The Numerals (числівники). Time. Cardinal numerals (кількісні числівники) 117.5 KB
  Let me introduce myself! Where are you from? – I am from Ukraine. I came from Rivne region. How old are you? – I am 25. What is your occupation (profession)? – I am businessman. Are you married? – Yes, I am. No, I am not. Do you have a business here? – Yes, I do. No, I don’t. I came here to have a rest. Glad to meet you! – Glad to meet you too.
53982. Prepositions of time, place, move 724 KB
  Look at these time expressions. Some use in, some use on or at. Put them in the correct column.
53983. To be. There is/are. To have. Pronouns 140.5 KB
  Whether the weather be fine, or whether the weather be not. Whether the weather be cold, or whether the weather be hot. We'll weather the weather whether we like it or not.
53984. Present Simple 65.5 KB
  We use the simple present tense when: the action is general the action happens all the time, or habitually, in the past, present and future the action is not only happening now the statement is always true