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


 

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

63067. Проектування скребкового конвеєра 414.14 KB
  Скребкові конвеєри використовують для транспортування сипких, зернистих матеріалів і коренеплодів у горизонтальному напрямі і під кутом близько 0-65° до горизонту, їх широко застосовують на зернових токах і елеваторах, у кормоцехах тваринницьких ферм...
63070. Создание транспортной накладной 202.6 KB
  Для поведите курсор к листу Лист 1 вызовите контекстное меню выберите команду Переименовать. 6 Выделите ячейки А5 и В5 на вкладке Главная выполните команду Объединить и поместить в центре нажав на кнопу.
63071. Детский массаж 1.1 MB
  При небольших покраснениях на коже ребенка в результате диатеза массаж нужно проводить осторожно, избегая участки с высыпаниями. Если же высыпания стали значительными, массаж в это время производить не следует.