50944

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

Лабораторная работа

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

Мета. Навчитися вирішувати задачу Коші методом Рунге-Кутта; скласти алгоритм. Устаткування: папір формату А4, ПК, програмне забезпечення Borland С++. Хід роботи Правила техніки безпеки Теоретичні дані Індивідуальне завдання.

Украинкский

2014-02-02

42 KB

1 чел.

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

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

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

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

Хід роботи

  1.  Правила техніки безпеки
  2.  Теоретичні дані
  3.  Індивідуальне завдання.

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

12

#include<iostream.h>

#include<math.h>

void main()

{

int i;

double h;

double *x = new double [6]; 

double *y = new double [6];

double *dy = new double [6];

double **k=new double*[6];

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

k[i]=new double [6];

cout<<"Input x[0] = ";

cin>>x[0];

cout<<"Input y[0] = ";

cin>>y[0];

cout<<"Input step h = ";

cin>>h;

for(i=0;x[i]<=1;i++)

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

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

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

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

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

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

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

}

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

{

  cout<<"\nx["<<i<<"] = "<<x[i]<<"        y["<<i<<"] = "<<y[i];

}

}

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


 

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

66501. Программирование и использование программных модулей 83 KB
  Разработать программный модуль (ПМ), в котором содержится не менее 4 подпрограмм (таблица 1) Задание 2 Составить Паскаль-программу, в которой используется не менее 4 подпрограмм программного модуля, разработанного в задании 1. Отчет должен содержать: -название лабораторной работы и номер варианта...
66502. Вибір та тестування оперативної пам’яті ПК 319 KB
  Мета: Набути вмінь та навиків при виборі та тестуванні оперативної пам’яті. ХІД РОБОТИ 1. Отримати від викладача материнську плату. 2. Визначити кількість роз’ємів для оперативної пам’яті. 3. Визначити тип модулів пам’яті, які встановлюються в дану материнську плату...
66503. Исследование датчика линейных ускорений 806 KB
  Датчик линейных ускорений (ДЛУ) предназначен для измерения линейных ускорений летательных аппаратов и выдачи электрического сигнала, величина которого пропорциональна линейному ускорению, действующему вдоль оси чувствительности.
66504. Электроизмерительные приборы, описание устройства, принципа действия и характеристик изучаемых приборов 5.33 MB
  Ознакомиться с классификацией назначением устройством принципом действия и характеристиками основных типов измерительных приборов. По представленным наглядным образцам и макетам уяснить принципы работы приборов и варианты их конструктивного исполнения отразив в отчете их основные технические характеристики.
66506. Тестування моніторів та відеоадаптерів 136 KB
  Крізь металеву маску або грати вони потрапляють на внутрішню поверхню скляного екрану монітора яка покрита різнокольоровими люмінофорними точками. Причини виходу монітора з ладу: Порушення втрата вакууму Часта причина відмови монітора особливо на початку експлуатації походить від того що...
66507. СИНТЕЗ СИСТЕМ УПРАВЛЕНИЯ С ОБРАТНОЙ СВЯЗЬЮ 539.92 KB
  Построим переходные и частотные характеристики непрерывной и дискретной модели: Рис. Переходная характеристика непрерывной системы Рис. Переходная характеристика дискретной системы Рис. Частотные характеристики непрерывной системы...