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


 

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

41853. Изучение и анализ конструкций мостов транспортных автомобилей 86.06 KB
  1 запорное кольцо подшипника; 16 полуосевая шестерня; 2 тормозная колодка; 17 болты крепления к балке заднего моста; 3 тормозной барабан; 18 подшипники ведущей шестерни; 4 шпилька крепления колеса; 19 сальник ведущей шестерни; 5 колпак колеса; 20 фланец; 6 тормозной цилиндр; 21 гайка ведущей шестерни; 7 тормозной щит; 22 кольцо грязеотражательное; 8 подшипник полуоси; 23 распорная втулка; 9 сальник полуоси; 24 регулировочное кольцо; 10 опорная чашка пружины; 25 ведущая шестерня; 11 ...
41854. Использование систем проверки орфографии и грамматики. Форматирование текста 201.77 KB
  Форматирование текста Цель: научиться использовать системы проверки орфографии и грамматики форматировать текст. Обратите внимание что в раскладке продуктов левый край ровный но текст отодвинут от левого края. Задание: Набрать следующий текст: Тесто рассыпчатое 400 г муки 200 г масла 05 стакана воды Растереть масло добавить муку воду всыпать 05 чайной ложки соли и замесить тесто. Порядок выполнения задания №2: Заголовок выровнять по центру с помощью элемента вкладки Главная шрифт полужирный вкладки Главная разрядка 3 пт Команда:...
41855. Ознакомление с устройством и функционированием счётчиков и испытание синхронного суммирующего, реверсивного и десятичного счётчиков 576.67 KB
  Между собой ячейки счётчика соединяют таким образом чтобы каждому числу импульсов соответствовали состояния 1 или 0 определенных ячеек. Каждый разряд счётчика может находиться в двух состояниях. Максимальное число N которое может быть записано в счётчике равно 2п 1 где п число разрядов счётчика.1 Условное изображение трехразрядного суммирующего счётчика показано на рис.
41856. Ознакомление с принципом работы и испытание интегрального цифроаналогового преобразователя 354.81 KB
  При построении устройств связывающих цифровое устройство с объектами использующими информацию в непрерывно изменяющейся форме требуется преобразование информации из аналоговой формы в цифровую и из цифровой в аналоговую. называют цифро-аналоговым преобразователем ЦАП. Сменяющиеся входные цифровые коды обуславливают сменяющееся ступенчатое напряжение на выходе L идеальная передаточная характеристика ЦАП. ЦАП с весовыми двоичновзвешенными сопротивлениями рис.
41857. АНАЛОГО-ЦИФРОВОЙ ПРЕОБРАЗОВАТЕЛЬ 234.35 KB
  Входным сигналом АЦП в течение некоторого промежутка времени t является постоянное напряжение равное отсчёту uвхkt входной аналоговой функции uвх. За это время на выходе АЦП формируется цифровой обычно двоичный код соответствующий дискретному отсчёту напряжения uвхkt. Количественная связь для любого момента времени определяется соотношением где u шаг квантования входного аналогового напряжения uвх; i погрешность преобразования напряжения uвхkt на данном шаге. Процесс квантования по уровню дискретизированной функции uвхkt...
41858. Изучение и анализ конструкций рамы 95.3 KB
  Ознакомились с устройством рамы различных автомобилей, научились анализировать их конструктивные особенности.
41860. Окислительно-восстановительное титрование. Иодометрическое определение пероксида водорода. Иодометрическое определение растворённого в воде кислорода 65.63 KB
  Сформировать умения по стандартизации раствора тиосульфата натрия; выполнению иодометрического определения пероксида водорода; иодометрического определения растворенного в воде кислорода. При этом к определяемому веществу добавляют взятое в заведомом избытке точное количество стандартного раствора иода. Какую среду сильнокислую слабокислую должен иметь раствор после добавления серной кислоты Почему при добавлении крахмала амилозы к раствору иода появляется синее окрашивание Какие ещё вещества могут взаимодействовать с иодом...
41861. Определение удельной теплоты плавления олова 286.55 KB
  Температура при которой вещество плавится называется температурой плавления вещества. Температура плавления для данного вещества при одинаковых условиях одинакова. Однако это не значит что в процессе плавления к телу не надо подводить энергию.