35235

Тема: Екстраполяційний метод Адамса розв’язання задачі Коші.

Практическая работа

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

h double Fdouble x double y { return cos2xy1.5xy; } void min {int n; double hb; doublek=new double [4]; doubleq=new double[n1]; doubledq1=new double[n1]; doubledq2=new double[n1]; doubledq3=new double[n1]; doublex=new double[n1]; doubley=new double[n1]; 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; for int i=0; i =2; i { k[0]=hFx[i]y[i]; k[1]=hFx[i]h 2y[i]k[0] 2;...

Украинкский

2013-09-09

42 KB

0 чел.

7.04.09  Мишура Денис, гр. 1ПМ-06

 Вариант 11

Лабораторна робота №28

Тема: Екстраполяційний метод Адамса розвязання задачі Коші.

Мета: Навчитися знаходити розв’язок диференційного рівняння екстраполяційним методом Адамса.

Устаткування: папір формату А4, ручка, калькулятор, ПЗ  С ++.

  1.  Індивідуальне завдання.

Використовуючи метод Адамса з трьома кінцевими різницями, скласти таблицю наближених  значень інтеграла диференційного рівняння , з початковими умовами на відрізку  з точністю 0,001. Початковий відрізок встановити методом Рунге-Кутта.

11)

Текст програми:

#include <iostream.h>

#include <math.h>

double F(double x, double y)

{

return (cos(2*x+y)+1.5*(x-y));

}

void main()

{int n;

double h,a,b;

double*k=new double [4];

double*q=new double[n+1];

double*dq1=new double[n+1];

double*dq2=new double[n+1];

double*dq3=new double[n+1];

double*x=new double[n+1];

double*y=new double[n+1];

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 (int i=0; i<=2; i++)

{

 k[0]=h*F(x[i],y[i]);

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

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

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

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

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

 q[i]=h*F(x[i],y[i]);

 cout<<"  "<<x[i]<<"    |  "<<y[i]<<"\n";

}

cout<<endl;

for (i=3; i<=n+1;i++)

{

q[i]=h*F(x[i],y[i]);

for (int k=0;k<i-1;k++) dq1[k]=q[k+1]-q[k];

for (k=0;k<i-2;k++) dq2[k]=dq1[k+1]-dq1[k];

for (k=0;k<i-3;k++) dq3[k]=dq2[k+1]-dq2[k];

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

y[i+1]=y[i]+q[i]+dq1[i-1]/2+5*dq2[i-2]/12+3*dq3[i-3]/8;

cout<<"  "<<x[i]<<"    |  "<<y[i]<<"\n";

}

cout<<endl;

}

  1.  Контрольні питання

  •  Для якої задачі використовується метод Адамса?
  •  Які формули методу Адамса?
  •  Як знайти початковий відрізок?
  •  Що являється результатом рішення задачі методом Адамса?
  •  Як оцінити похибку метода на практиці?
  1.  Звіт про роботу.
  2.  Захист роботи

PAGE  1


 

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

5601. Введение потребителей в заблуждение 96 KB
  Введение потребителей в заблуждение Другой формой недобросовестной конкуренции является: введение потребителей в заблуждение относительно характера, способа и места изготовления, потребительских свойств, качества товаров. Введение потребителей в з...
5602. Кинематика. Механическое движение 55.5 KB
  Кинематика Механическим движением называется изменение положения предмета относительно заданной системы отсчета. Понятие системы отсчета включает в себя тело отсчета и систему координат. Для большинства задач нашего курса достаточно ограничиться пря...
5603. Экономика предприятия. Курс лекций 1 MB
  Производственная деятельность людей осуществляется в форме предприятия. Предприятие - трудовая кооперация людей, составленная из разнородных специалистов, способных производить полезный для общества товар или услугу, обменивая его затем на...
5604. Агентура исторического духа в философии События М. М. Бахтина. Ретроактивный футуризм Ницше: генеалогия как онтология события 479 KB
  Рассмотрение аналитики историчности в значительной степени мотивировано той нередуцируемой навязчивостью, которую выявляет данный феномен как предел традиционного метафизического мышления. Для нас, данные заметки указывают, прежде вс...
5605. Контроль упруго-напряженного состояния поверхностей металла деталей машин. Расчет напряжений термической анизотропии 105 KB
  Контроль упруго-напряженного состояния поверхностей металла деталей машин. Расчет напряжений термической анизотропии, возникающих в поликристаллах на границах двух зерен. (вар.№17) Цель работы - освоение методов оценки внутренних напряжен...
5606. Графический интерфейс в MS WINDOWS 1.37 MB
  Графический интерфейс в MSWINDOWS Переход к программированию для Windows требует знаний графических возможностей этой операционной системы. Связано это с тем, что Windows является графической системой, в нем отсутствует текстовый режим...
5607. Устройства формирования и обработки сигналов 323 KB
  Предисловие Современная радиоэлектронная аппаратура и радиотехнические системы относятся к классу сложных электронных систем, для обслуживания которой требуются высококвалифицированные специалисты. Дисциплина Устройство, формирование и обрабо...
5608. Вплив держави на розміщення продуктивних сил 104 KB
  Регіональна економічна політика, її сутність та завдання. Вплив держави на розміщення продуктивних сил (РПС) в країні здійснюється за допомогою економічної регіональної політики. Отже для розгляду впливу держави на РПС необхідно розглянути регіон...
5609. Верховые болота 80.5 KB
  Верховые болота По площади и запасам торфа среди болот в России преобладают верховые болота. Хотя лучше всего они развиты в тайге, но встречаются и южнее ее границ, вплоть до степи. Почему же эти болота назвали верховыми? Это связано с особенностью...