35257

Метод прогонки розв’язання крайової задачі. Складання алгоритму

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

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

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

Украинкский

2013-09-09

29.5 KB

1 чел.

еревко О.В.                                                                                                              2ПМС-07

                                     Лабораторна робота №29,30

Тема. Метод прогонки розв’язання крайової задачі. Складання алгоритму.

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

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

Хід роботи

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

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

16)

Программа на С++

#include<iostream.h>

#include<math.h>

double PX(double x)

{

return -3;

}

double QX(double x)

{

return 1/x;

}

double FX(double x)

{

return 1;

}

void main()

{

int n;

double h,a,b,LO,L1,A,BO,B1,B;

cout<<"Vvedite shag h:";

cin>>h;

cout<<"Vvedite a: ";

cin>>a;

cout<<"Vvedite b: ";

cin>>b;

n=(b-a)/h;

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

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

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

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

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

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

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

cout<<"Vvedite L[0]: ";

cin>>L0;

cout<<"Vvedite L[1]: ";

cin>>L1;

cout<<"Vvedite A: ";

cin>>A;

cout<<"Vvedite B[0]: ";

cin>>B0;

cout<<"Vvedite B[1]: ";

cin>>B1;

cout<<"Vvedite B: ";

cin>>B;

x[0]=a;

c[0]=L1/(L0*h-L1);

d[0]=A*h/Ll;

//Prjamoi hod

for(int i=0;i<=n;i++)

{

x[i]=x[0]+i*h;

m[i]=2*h*h*QX(x[i])-4)/(2+h*PX(x[i]));

nn[i]=(2-h*PX(x[i]))/(2+h*PX(x[i]));

f[i]=2*FX(x[i])/(2+h*PX(x[i]));

if(i!=0 {c[i]=1/(m[i]-nn[i]*c[i-1]);

d[i]=FX(x[i])*h*h-nn[i]*c[i-1]*d[i-1];}

}

//Obratnii hod

y[n]=(B*h+Bl*c[n-1]*d[n-1])/(BO*h+B1*(c[n-1]+1));

for (i=n-1;i>=0;i--) y[i]=c[i]*(d[i]-y[i+1]);

cout<<"Tablichnaja fynkchija:\n";

for(i=0;i<=n;i++) cout<<"x["<<i<<"]="<<x[i]<<"\ty["<<i<<"]="<<y[i]<<"\n";

}

Результат программы


 

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

66429. УПРАВЛІННЯ КОНКУРЕНТОСПРОМОЖНІСТЮ ПРОМИСЛОВОЇ ПРОДУКЦІЇ З УРАХУВАННЯМ ЕКОЛОГІЧНОЇ СКЛАДОВОЇ 474.5 KB
  Формування відкритої економіки, членство України в СОТ супроводжується не тільки позитивними, але і негативними явищами, серед яких слід відзначати жорсткі конкурентні умови, особливо для країн із застарілою технологією та виробництвом екологічно неякісної продукції.
66430. Внутрішні перенапруги в електричних мережах з ізольованою нейтраллю, їх обмеження та основи моделювання 4.68 MB
  Головним призначенням систем пересилання, розподілу та споживання електроенергії є ефективне забезпечення потреб національного господарства в енергії відповідної якості за мінімальних витрат і достатньому рівні надійності. Принципи побудови, розвитку сучасних і перспективних систем енергопостачання...
66431. Особенности корректировки таможенной стоимости товаров 203.5 KB
  Необходимо проанализировать изменения, связанные с образованием Таможенного союза; Выявить принципиально новые особенности Таможенного кодекса Таможенного союза; Рассмотреть что такое таможенная стоимость товаров, как она контролируется и в связи, с чем может корректироваться.
66432. Некоторые вопросы практики вексельного обращения 116.82 KB
  Наименование вексель включенное в самый текст документа и выраженное на том языке на котором этот документ составлен; простое и ничем не обусловленное предложение уплатить определенную сумму; 3 наименование того кто должен платить плательщика; указание срока платежа; указание места в котором должен быть совершен платеж; наименование того кому или приказу кого платеж должен быть совершен; указание...
66434. Акцизний збір в системі платежів в державний бюджет України 331 KB
  Необхідність збільшення податків і в той же час справедливого їх розподілу ставить державу перед складною дилемою: економічна ефективність чи справедливість? Оптимальне поєднання й урахування цих аспектів - неодмінна умова гармонійності податкової політики.
66435. Методика обліку і організація аналізу і аудиту виробничих запасів підприємства (ТзОВ „Уніплит”) 860.5 KB
  Теорикоекономічна сутність виробничих запасів та їх характеристика як обєкту обліку та аудиту Економічна сутність запасів та іх класифікація Оцінка виробничих запасів для цілей обліку Нормативноправова база з обліку та аудиту виробничих запасів Загальноекономічна характеристика ТзОВ...
66436. Облік і аудит основних засобів на ЗАТ „Лукор” 511.5 KB
  Методика обліку нематеріальних активів Визнання нематеріальних активів Відповідно до ПБО 8 нематеріальними вважаються немонетарні активи які: не мають матеріальної форми; не можуть бути ідентифіковані; утримуються підприємством з метою використання протягом періоду більше одного року...
66437. Поняття модернізму і постмодернізму в українському літературознавстві 249 KB
  Об’єктом дослідження стали два потужних явища літератури: модернізм і постмодернізм. Зважаючи на усталеність і визнаність їх у світовому літературному процесі, прагнемо “узаконити” їх в українській літературі, оскільки ще точаться суперечки щодо правомірності вживання термінів...