Файловый архив студентов. Файловый архив студентов.
1302 вуза, 5124 предмета.
/ Регистрация
FAQ Обратная связь Вопросы и предложения
Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Уфимский Государственный Авиационный Технический Университет
Предмет:
Программирование
Файл:

Программирование на C / C++ / Лабораторная работа №12 / ЛабаВМ

.CPP
Скачиваний:
27
Добавлен:
02.05.2014
Размер:
7.69 Кб
Скачать
☆ 
#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
Помощь Обратная связь Вопросы и предложения Пользовательское соглашение Политика конфиденциальности