Добавил:
Studfiles2
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Программирование на C / C++ / Лабораторная работа №12 / ЛабаВМ
.CPP#include <iostream.h>
#include <iomanip.h>
#include <math.h>
#include <conio.h>
#define N 65536
void main()
{
double a,b,a1,b1,a2,b2,h,hx,sum,sumx,c,x,xx,U,fa,fb,fa1,fb1,d,fa2,fb2;
double y1[20][6],y2[20][6],u1[20][6],u2[20][6],ux1[20][6],ux2[20][6];
long int i,n,j,v;
const int m=4,
k1=1, //¬Ґв®¤ «Ґўле Їап¬.
k2=2; //¬Ґв®¤ ва ЇҐжЁ©
//ўлзЁб«ҐЁҐ Ї®бв®п®© бy ¤«п дгЄжЁЁ y
a=0;
b=2;
U=6.4; //в®з®Ґ § 票Ґ ЁвҐЈа « дгЄжЁЁ y
clrscr();
c=32; //Ї®бв®п п ¤«п дгЄжЁЁ y
cout<<"c="<<c<<endl;
fa=pow(a,m);
fb=pow(b,m);
//ўлзЁб«ҐЁҐ ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа « ¤«п дгЄжЁЁ y
i=0;
for (n=2;n<=N;n=2*n)
{
sum=0;
h=(b-a)/n;
y1[i][5]=c*(pow(h,k1));//ўлзЁб«ҐЁҐ ⥮аЁв. Ї®Ја-вЁ ¬Ґв®¤®¬ «Ґўле Їап¬.
y2[i][5]=c*(pow(h,k2));//ўлзЁб«ҐЁҐ ⥮аЁв. Ї®Ја-вЁ ¬Ґв®¤®¬ ва ЇҐжЁ©
///////////////////////
for (j=1;j<=n;j++)
{
x=a+double(j)*h;
sum=sum+h*(pow(x,m));
}
y1[i][0]=n;
y1[i][1]=sum;//¬ ббЁў ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа « ¬Ґв®¤®¬ «Ґўле Їап¬.
///////////////////////
sum=0;
for (j=1;j<=n-1;j++)
{
x=a+double(j)*h;
sum=sum+h*(pow(x,m));
}
y2[i][0]=n;
y2[i][1]=h*(fa+fb)/2+sum;//¬ ббЁў ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа «
//¬Ґв®¤®¬ ва ЇҐжЁ©
i++;
}
v=i;//Є®«ЁзҐбвў® § 票© ЁвҐЈа «
cout<<"v="<<v<<endl;
getch();
////////////////////////////////////////////////
a1=0;
b1=1.5;
a2=0.001;
b2=1.5;
fa1=(pow(a1,1.0/m))*(sqrt(1+a1*a1));
fb1=(pow(b1,1.0/m))*(sqrt(1+b1*b1));
fa2=(pow(a2,1.0/m))*(sqrt(1+a2*a2));
fb2=fb1;
//ўлзЁб«ҐЁҐ ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа « ¤«п дгЄжЁЁ u
i=0;
for (n=2;n<=N;n=2*n)
{
sum=0;
sumx=0;
h=(b1-a1)/n;
hx=(b2-a2)/n;
for (j=1;j<=n;j++)
{
x=a1+double(j)*h;
sum=sum+h*(pow(x,1.0/m))*(sqrt(1+x*x));
xx=a2+double(j)*hx;
sumx=sumx+hx*(pow(xx,1.0/m))*(sqrt(1+xx*xx));
}
u1[i][0]=n;
u1[i][1]=sum;//¬ ббЁў ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа «
//¬Ґв®¤®¬ «Ґўле Їап¬. [a1,b1]
ux1[i][0]=n;
ux1[i][1]=sumx;//¬ ббЁў ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа «
//¬Ґв®¤®¬ «Ґўле Їап¬. [a2,b2]
//////////////////
sum=0;
sumx=0;
for (j=1;j<=n-1;j++)
{
x=a1+double(j)*h;
sum=sum+h*(pow(x,1.0/m))*(sqrt(1+x*x));
xx=a2+double(j)*hx;
sumx=sumx+hx*(pow(xx,1.0/m))*(sqrt(1+xx*xx));
}
u2[i][0]=n;
u2[i][1]=h*(fa1+fb1)/2+sum;//¬ ббЁў ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа « ¬Ґв®¤®¬
//ва ЇҐжЁ© [a1,b1]
ux2[i][0]=n;
ux2[i][1]=hx*(fa2+fb2)/2+sumx;//¬ ббЁў ЇаЁЎ«Ё¦Ґле § 票© ЁвҐЈа « ¬Ґв®¤®¬
i++; //ва ЇҐжЁ© [a2,b2]
}
v=i;//Є®«ЁзҐбвў® § 票© ЁвҐЈа «
//ўлзЁб«ҐЁҐ Ї®ЈаҐи®бвЁ Љ
for(i=2;i<v;i++)
{
//¬Ґв®¤ «Ґўле Їап¬.¤«п дгЄжЁЁ y
y1[i][2]=(y1[i-1][1]-y1[i-2][1])/(y1[i][1]-y1[i-1][1]);
//¬Ґв®¤ ва ЇҐжЁ© ¤«п дгЄжЁЁ y
y2[i][2]=(y2[i-1][1]-y2[i-2][1])/(y2[i][1]-y2[i-1][1]);
//¬Ґв®¤ «Ґўле Їап¬.¤«п дгЄжЁЁ u [a1,b1]
u1[i][2]=(u1[i-1][1]-u1[i-2][1])/(u1[i][1]-u1[i-1][1]);
//¬Ґв®¤ ва ЇҐжЁ© ¤«п дгЄжЁЁ u [a1,b1]
u2[i][2]=(u2[i-1][1]-u2[i-2][1])/(u2[i][1]-u2[i-1][1]);
//¬Ґв®¤ «Ґўле Їап¬. ¤«п дгЄжЁЁ u [a2,b2]
ux1[i][2]=(ux1[i-1][1]-ux1[i-2][1])/(ux1[i][1]-ux1[i-1][1]);
//¬Ґв®¤ ва ЇҐжЁ© ¤«п дгЄжЁЁ u [a2,b2]
ux2[i][2]=(ux2[i-1][1]-ux2[i-2][1])/(ux2[i][1]-ux2[i-1][1]);
}
y1[0][2]=0; y1[1][2]=0; y2[0][2]=0; y2[1][2]=0;
u1[0][2]=0; u1[1][2]=0; u2[0][2]=0; u2[1][2]=0;
ux1[0][2]=0; ux1[1][2]=0; ux2[0][2]=0; ux2[1][2]=0;
//ўлзЁб«ҐЁҐ Ї®ЈаҐи®бвЁ в®з®© ¤«п дгЄжЁЁ y
for(i=0;i<v;i++)
{
y1[i][3]=U-y1[i][1];//¬Ґв®¤ «Ґўле Їап¬.¤«п дгЄжЁЁ y
y2[i][3]=U-y2[i][1];//¬Ґв®¤ ва ЇҐжЁ© ¤«п дгЄжЁЁ y
}
//ўлзЁб«ҐЁҐ Ї®ЈаҐи®бвЁ ђгЈҐ
for(i=1;i<v;i++)
{
//¬Ґв®¤ «Ґўле Їап¬.¤«п дгЄжЁЁ y
y1[i][4]=(y1[i][1]-y1[i-1][1])/(pow(2,k1)-1);
//¬Ґв®¤ ва ЇҐжЁ© ¤«п дгЄжЁЁ y
y2[i][4]=(y2[i][1]-y2[i-1][1])/(pow(2,k2)-1);
//¬Ґв®¤ «Ґўле Їап¬.¤«п дгЄжЁЁ u [a1,b1]
u1[i][3]=(u1[i][1]-u1[i-1][1])/(pow(2,k1)-1);
//¬Ґв®¤ ва ЇҐжЁ© ¤«п дгЄжЁЁ u [a1,b1]
u2[i][3]=(u2[i][1]-u2[i-1][1])/(pow(2,k2)-1);
//¬Ґв®¤ «Ґўле Їап¬. ¤«п дгЄжЁЁ u [a2,b2]
ux1[i][3]=(ux1[i][1]-ux1[i-1][1])/(pow(2,k1)-1);
//¬Ґв®¤ ва ЇҐжЁ© ¤«п дгЄжЁЁ u [a2,b2]
ux2[i][3]=(ux2[i][1]-ux2[i-1][1])/(pow(2,k2)-1);
}
y1[0][4]=0; y2[0][4]=0;
u1[0][3]=0; u2[0][3]=0;
ux1[0][3]=0; ux2[0][3]=0;
//ўлў®¤ нЄа в Ў«Ёжл дгЄжЁЁ y ¬Ґв®¤®¬ «Ґўле Їап¬.
cout<<" ‡ зҐЁп Ї®ЈаҐи®б⥩ ¤«п дгЄжЁЁ y ¬Ґв®¤®¬ «Ґўле Їаאַ㣮«мЁЄ®ў"<<endl;
cout<<setw(11)<<setiosflags(ios::left)<<"n"<<
setw(11)<<setiosflags(ios::left)<<"€вҐЈа «"<<
setw(11)<<setiosflags(ios::left)<<"Љ"<<
setw(11)<<setiosflags(ios::left)<<"’®з."<<
setw(11)<<setiosflags(ios::left)<<"ђгЈҐ"<<
setw(11)<<setiosflags(ios::left)<<"’Ґ®а."<<endl;
for(i=0;i<v;i++)
{
for(j=0;j<6;j++)
cout<<setw(10)<<setiosflags(ios::left)<<y1[i][j]<<" ";
cout<<endl;
}
getch();
cout<<endl;
//ўлў®¤ нЄа в Ў«Ёжл дгЄжЁЁ y ¬Ґв®¤®¬ ва ЇҐжЁ©
cout<<" ‡ зҐЁп Ї®ЈаҐи®б⥩ ¤«п дгЄжЁЁ y ¬Ґв®¤®¬ ва ЇҐжЁ©"<<endl;
cout<<setw(10)<<setiosflags(ios::left)<<"n"<<
setw(10)<<setiosflags(ios::left)<<"€вҐЈа «"<<
setw(10)<<setiosflags(ios::left)<<"Љ"<<
setw(15)<<setiosflags(ios::left)<<"’®з."<<
setw(15)<<setiosflags(ios::left)<<"ђгЈҐ"<<
setw(15)<<setiosflags(ios::left)<<"’Ґ®а."<<endl;
for(i=0;i<v;i++)
{
for(j=0;j<6;j++)
{
if (j<3) d=10; else d=15;
cout<<setw(d)<<setiosflags(ios::left)<<y2[i][j];
}
cout<<endl;
}
getch();
cout<<endl;
//ўлў®¤ нЄа в Ў«Ёжл дгЄжЁЁ u ¬Ґв®¤®¬ «Ґўле Їап¬. [a1,b1]
cout<<" ‡ зҐЁп Ї®ЈаҐи®б⥩ ¤«п дгЄжЁЁ u ¬Ґв®¤®¬ «Ґўле Їаאַ㣮«мЁЄ®ў [a1,b1]"<<endl;
cout<<setw(14)<<setiosflags(ios::left)<<"n"<<
setw(14)<<setiosflags(ios::left)<<"€вҐЈа «"<<
setw(14)<<setiosflags(ios::left)<<"Љ"<<
setw(14)<<setiosflags(ios::left)<<"ђгЈҐ"<<endl;
for(i=0;i<v;i++)
{
for(j=0;j<4;j++)
cout<<setw(13)<<setiosflags(ios::left)<<u1[i][j]<<" ";
cout<<endl;
}
getch();
cout<<endl;
//ўлў®¤ нЄа в Ў«Ёжл дгЄжЁЁ u ¬Ґв®¤®¬ ва ЇҐжЁ© [a1,b1]
cout<<" ‡ зҐЁп Ї®ЈаҐи®б⥩ ¤«п дгЄжЁЁ u ¬Ґв®¤®¬ ва ЇҐжЁ© [a1,b1]"<<endl;
cout<<setw(14)<<setiosflags(ios::left)<<"n"<<
setw(14)<<setiosflags(ios::left)<<"€вҐЈа «"<<
setw(14)<<setiosflags(ios::left)<<"Љ"<<
setw(14)<<setiosflags(ios::left)<<"ђгЈҐ"<<endl;
for(i=0;i<v;i++)
{
for(j=0;j<4;j++)
cout<<setw(13)<<setiosflags(ios::left)<<u2[i][j]<<" ";
cout<<endl;
}
getch();
cout<<endl;
//ўлў®¤ нЄа в Ў«Ёжл дгЄжЁЁ u ¬Ґв®¤®¬ «Ґўле Їап¬. [a2,b2]
cout<<" ‡ зҐЁп Ї®ЈаҐи®б⥩ ¤«п дгЄжЁЁ u ¬Ґв®¤®¬ «Ґўле Їаאַ㣮«мЁЄ®ў [a2,b2]"<<endl;
cout<<setw(14)<<setiosflags(ios::left)<<"n"<<
setw(14)<<setiosflags(ios::left)<<"€вҐЈа «"<<
setw(14)<<setiosflags(ios::left)<<"Љ"<<
setw(14)<<setiosflags(ios::left)<<"ђгЈҐ"<<endl;
for(i=0;i<v;i++)
{
for(j=0;j<4;j++)
cout<<setw(13)<<setiosflags(ios::left)<<ux1[i][j]<<" ";
cout<<endl;
}
getch();
cout<<endl;
//ўлў®¤ нЄа в Ў«Ёжл дгЄжЁЁ u ¬Ґв®¤®¬ ва ЇҐжЁ© [a2,b2]
cout<<" ‡ зҐЁп Ї®ЈаҐи®б⥩ ¤«п дгЄжЁЁ u ¬Ґв®¤®¬ ва ЇҐжЁ© [a2,b2]"<<endl;
cout<<setw(14)<<setiosflags(ios::left)<<"n"<<
setw(14)<<setiosflags(ios::left)<<"€вҐЈа «"<<
setw(14)<<setiosflags(ios::left)<<"Љ"<<
setw(14)<<setiosflags(ios::left)<<"ђгЈҐ"<<endl;
for(i=0;i<v;i++)
{
for(j=0;j<4;j++)
cout<<setw(13)<<setiosflags(ios::left)<<ux2[i][j]<<" ";
cout<<endl;
}
getch();
}
Соседние файлы в папке Лабораторная работа №12