ID: 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";

}

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