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


 

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

51361. Красноярский территориальный центр фирменного транспортного обслуживания 61.42 KB
  Сокращенное наименование Подразделения: Красноярский ТЦФТО ОАО РЖД. Подразделение для организации и осуществления своей деятельности открывает расчетный счет ОАО РЖД в банках и иных кредитных организациях в порядке установленном законодательством Российской Федерации внутренними документами ОАО РЖД Филиала. Подразделение имеет круглую печать содержащую его сокращенное наименование указание на его место нахождения а также полное фирменное наименование ОАО РЖД. Подразделение может иметь штампы и бланки со своим наименованием...
51362. Визуализация графика математической функции 1.04 MB
  Получение значения объекта step сохранение его в переменную s; получение значения объекта xmin с последующим преобразованием в формат flot и сохранением в переменную x0; получение значения объекта xmx с последующим преобразованием в формат flot и сохранением в переменную x1; получение значения объекта osx сохранение его в переменную sx; получение значения объекта osy сохранение его в переменную sy; получение значения объекта oscle сохранение его в переменную sc; создание сигнала chnge c передачей параметров...
51365. Демультиплексоры 48.24 KB
  Техническое задание Требуется спроектировать четырехразрядный демультиплексор на языке VHDL составить таблицу истинности спроектированного устройства показать логическую и техническую схемы и привести временную диаграмму с полученными результатами. Спроектировать четырехразрядный демультиплексор имеющий два входа адресный и информационный и один выход на языке VHDL. ДеМультиплексор позволяет передавать сигнал с одного входа на...
51366. Экономический анализ финансово-хозяйственной деятельности ООО «Адидас» 329 KB
  Организационная характеристика предприятия Экономический анализ финансово-хозяйственной деятельности предприятия Основные цели преддипломной практики: закрепление расширение углубление и систематизацию знаний полученных при изучении общепрофессиональных и специальных дисциплин на основе изучения деятельности предприятия отрасли; проведение системного анализа организации с целью выявления проблем управления и разработки мероприятий по их устранению; более глубоко освоить методы и приемы системного анализа; сбор и...
51368. Исследование начальной остойчивости плавучей полупогружной буровой установки 155 KB
  Ознакомление студентов с особенностями остойчивости плавучих полупогружных буровых установок (ППБУ) и их поведения на взволнованной поверхности моря, изучение основных положений теории и расчета, а также ознакомление с методикой постановки эксперимента по определению параметров начальной остойчивости плавучих технических средств для освоения шельфа.
51369. Двухфазная СМО с отказами 95.5 KB
  Для упрощения расчёта представим данную СМО как совокупность 2ух одноканальных. Т.к. в данной системе очередь не бесконечной длинны, то все расчёты будут не очень точны. Но главная цель проведения данных расчётов – это сравнение их результатов с результатами имитационной модели (программой). Для оценки соответствия результатов такой точности будет достаточно.