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


 

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

50839. Применение программы РSpice для моделирования двухтактного полумостового преобразователя постоянного напряжения 218.5 KB
  При включении питания схема запуска кратковременно замыкает ключ S1. При замыкании S1 через первичную обмотку W11 силового трансформатора TV1 начинает протекать ток, являющийся суммой токов...
50842. Изучение характеристик вакуумного диода и определение работы выхода электронов из вольфрама 737.5 KB
  Цель работы: экспериментальное изучение характеристик вакуумного диода и определение работы выхода электронов из вольфрама. Описание работы: Экспериментальная часть работы состоит из трех этапов: измерение мощностей нагрева катода Рн по значениям Uн. Сила тока накала измеряется стрелочным амперметром А1 встроенным в источник ИП а напряжение накала UН внешним цифровым прибором род работы U= подключаемым к соответствующим гнездам на панели источника ИП прибор на рис. Сила тока в лампе...