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


 

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

10513. Суперглобальный массив $_SERVER 54 KB
  Суперглобальный массив SERVER Одним из важнейших предопределенных массивов является массив SERVER в него PHPинтерпретатор помещает переменные полученные от сервера. Без данных переменных сложно организовать полноценную поддержку Webприложений. Ниже приводится описан...
10514. Ответственность за неисполнение денежных обязательств 36.5 KB
  Ответственность за неисполнение денежных обязательств. За их неисполнение начисляется процент. Основание ответственности. Для возникновения ответственности необходимы четыре элемента: Противоправное поведение: это неправильное поведение нарушающее дого
10515. Схема коммерческих отношений в трамповом судоходстве 52.5 KB
  Схема коммерческих отношений в трамповом судоходстве. Общая схема: Стивидорная компания чаще называется оператор терминала. Договор со стивидорной компанией заключает продавец лучше знает свойства груза. Договор с агентом заключается судовладе...
10516. Особенности работы агента в трамповом судоходстве 51 KB
  На капитана возложен контроль за качеством агентского обслуживания судна. Во всех спорных случаях право окончательного решения принадлежит капитану. Судовладелец вправе требовать от агента услуг и выполнения работ только в рамках законов правил и обычаев страны пр
10517. Коммерческое значение порта 45 KB
  Коммерческое значение порта. Порт это либо понятие географическое либо административно-управленческое. Соотношение портовых властей и коммерческих предприятий: такие как город и мэрия. За рубежом порт управляется общественным советом. Сюда входят представите...
10518. Буксирная помощь 36 KB
  Буксирная помощь. Как правило частные буксирные компании. С буксирными компаниями можно заключить отдельный долгосрочный договор тогда судовладелец обязан обращаться только к этой компании. Буксировка: внутрипортовая: содействие в маневрировании судна в пор
10519. Способы оплаты в морском бизнесе 42.5 KB
  Способы оплаты в морском бизнесе. оплата наличными: это только оплата. Предусматривает что получатель платежа может беспрепятственно и немедленно получить платеж в банке в день согласованного платежа. Так оплачиваются: фрахт, демеридж, мертвый фрахт...
10520. Согласование фрахтовой ставки 38.5 KB
  Согласование фрахтовой ставки. Надо согласовать 11 пунктов. ставка фрахта: арифметическая величина; валюта фрахта: доллары валюта в которой выражена фрахтовая ставка; фрахтовая единица: физическая единица на которой установлена ставка фрахта. Единица: ...
10521. От чего зависит оборачиваемость коносамента 38 KB
  Оборотные документы имеют два свойства: Освобождение правомерного держателя от каких бы то ни было возражений в платеже со стороны плательщика. Коносамент является носителем обязательного права. Что делать кредитору если не платят по векселю Можно обратиться к ка