青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

newplan

阿基米德在洗澡時(shí)發(fā)現(xiàn)浮力原理,高興得來不及穿上褲子,跑到街上大喊:Eureka(我找到了)。
posts - 39, comments - 26, trackbacks - 0, articles - 4
  C++博客 :: 首頁(yè) :: 新隨筆 :: 聯(lián)系 :: 聚合  :: 管理
/*
Romberg Algorithm 
開發(fā)者:newplan
開發(fā)日期:
08.05.07 
*/


/*=======================================*/
/*INCLUDES*/ 
#include           
<cstdlib>
#include           
<iostream>
#include           
<cmath>
/*=======================================*/
/*MACROS  USED IN THIS FILE*/
#define            MAX       
20   
#define            PRECISION 
0.000008
/*=======================================*/
/*DECLARE NAMES IN STL NAMESAPCE */
using  std::cout;
using  std::endl;

/*=======================================*/
/*CLASS FUNC (FUNCTION OBJECT): THE ORIGINAL FUNCTION WE WANT TO INTEGRAL*/
class  func{
       
public:
              func(
double x=1.0):exp(x){}
              
double operator()(const double& dnum)const{return pow(dnum,exp);}
       
private:
              
double exp;
       };
/*=======================================*/
/*CLASS ECHELONFUNC (FUNCTION OBJECT)梯形法的遞推公式*/
class  echelonFunc{
       
public:
              echelonFunc(
double begining,double ending,func & myfunc);
              
double operator()();
       
private:
               
double    h;
               
int       n;
               
double    T;
               func      myfunc;
               
double    begining;
               
double    ending ;
     };
/*=======================================*/
echelonFunc::echelonFunc(
double begining,double ending,func & myfunc)
              {
               this
->begining=begining;
               this
->ending=ending;
               this
->h=ending-begining;
               this
->n=0;
               this
->T=0;
               this
->myfunc=myfunc;//FUCNTION 
              }
/*------------------------------*/
/* INCREASE FUNCTION 遞推函數(shù)*/
double echelonFunc::operator()()
              {   
if(this->n==0)
                  {
                   this
->T=h*0.5*(myfunc(this->begining)+myfunc(this->ending));
                   this
->n=1;
                   return this
->T;
                  }
                  
double len=0.5*h;
                  
double sum=0;
                  
int k=0;
                  
for(k=0;k<this->n;k++)
                      {
                       sum
+=myfunc(len);
                       
len=len+h;
                      }
                  this
->T=0.5*this->T+0.5*h*sum;
                  this
->h/=2;
                  this
->n*=2;
                  return this
->T;
              }
/*=======================================*/
/*THE MAIN CLASS IN THIS PROGRAM*/ 
class  Romberg{
       
public:
              Romberg(
double begining,double ending,double exp);
              ~Romberg();
       
private:
              void RombergCPU();
/*THE MOST IMPORTANT FUNCTION IN THIS PROGRAM*/
              echelonFunc 
*echol; 
              
double T[MAX][MAX];/*STO THE ROMBERG TABLE*/       
       };
/*------------------------------*/
Romberg::Romberg(
double begining ,double ending ,double paraexp)
{
    func     myfunc(paraexp);
    echol 
=  new echelonFunc(begining,ending,myfunc);
    RombergCPU();
}
/*------------------------------*/
Romberg::~Romberg()
{
     delete echol;
}
/*------------------------------*/
void Romberg::RombergCPU()
{   clock_t Start; 
//TIME STARAT
    clock_t 
End;   //TIME END
    
double *p[MAX];//WE USE THIS POINTER ARRAY TO ACCELERATE ALGOTITHM
    
double **q;
    
int i,j;
    Start 
= clock();//TIME START FROM HERE
    
for(i = 0,q = p; q < p+MAX; q++,i++)
        
*q= &T[i][0];
    
double a,b,pows; 
    cout
<<"-----------------------Romberg Algorithm---------------------"<<endl;
    
*p[0]=(*echol)();
    cout
<<"  "<<*p[0]<<endl;
    p[
0]++;
    
do{
       
*p[0]=(*echol)();
        cout
<<"  "<<*p[0];
        p[
0]++;
        
for(i=1;;i++)
                  {
                    pows
=pow(4.0,double(i));
                    a
=pows/(pows-1);
                    b
=1/(pows-1);
                    
*p[i]=a*(*(p[i-1]-1))-b*(*(p[i-1]-2));//ROMBERG ALGORITHM
                    cout
<<"  "<<*p[i];
                    
if(p[i]==&T[i][0])
                      {
                      p[i]
++
                      break;
                      } 
                     p[i]
++;
                   }
                   cout
<<endl;//fabs(T[i][0]-T[i-1][0])
       }
while(fabs(T[i][0]-T[i-1][0])>PRECISION);
    
    
End = clock();//TIME END HERE
    
    cout
<<"-------------------------------------------------------------"<<
    endl
<<"  TIME SPEND:"<<(double)(End-Start)/CLOCKS_PER_SEC<<endl; 

/*=======================================*/
/*MAIN FUNCTION*/
int main(int argc, char *argv[])
{  
   Romberg  romberg(
0,1,1.5);//ROMBERG API :BEGIN(0END(1EXP(1.5)
   system(
"PAUSE");
   return EXIT_SUCCESS;
}

Feedback

# re: 龍貝格Romberg算法cpp實(shí)現(xiàn)  回復(fù)  更多評(píng)論   

2008-09-05 17:57 by matdu
good
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            亚洲一区免费在线观看| 亚洲精品久久久久久下一站 | 国产亚洲在线观看| 欧美在线观看视频一区二区三区| 中文日韩电影网站| 国产亚洲欧美日韩日本| 毛片基地黄久久久久久天堂| 久久综合影视| 亚洲婷婷综合色高清在线| 亚洲天天影视| 国产一区二区三区在线观看网站| 久久亚洲国产精品一区二区| 快射av在线播放一区| 在线视频精品一区| 性色av香蕉一区二区| 亚洲国内在线| 亚洲一区二区免费视频| 国内激情久久| 亚洲人成网站精品片在线观看| 欧美人在线观看| 久久aⅴ国产紧身牛仔裤| 久久婷婷蜜乳一本欲蜜臀| 一二美女精品欧洲| 欧美专区福利在线| 亚洲美女中文字幕| 欧美中在线观看| 一区二区三区黄色| 久久久精品免费视频| 日韩亚洲欧美综合| 欧美在线视频观看免费网站| 99国内精品久久久久久久软件| 午夜精品久久久久久久白皮肤| 亚洲国产女人aaa毛片在线| 亚洲亚洲精品三区日韩精品在线视频| 在线播放中文一区| 亚洲欧美资源在线| 亚洲视频福利| 欧美1区3d| 久久这里只有| 国产女优一区| 一本色道久久综合亚洲精品高清| 黄页网站一区| 性娇小13――14欧美| 一本大道久久a久久精二百| 老司机久久99久久精品播放免费| 亚洲欧美一级二级三级| 欧美巨乳在线观看| 欧美国产日韩精品| 国内精品美女在线观看| 午夜精品一区二区三区在线| 99伊人成综合| 欧美激情片在线观看| 欧美大片一区二区| 樱花yy私人影院亚洲| 欧美有码在线视频| 欧美一区成人| 国产区精品视频| 亚洲免费中文字幕| 午夜亚洲视频| 国产精品久久久一本精品| 亚洲精品中文字幕有码专区| 亚洲国产日韩一级| 免费成人黄色av| 亚洲国产精品成人| 99国产一区| 欧美日韩免费高清| 中文无字幕一区二区三区| 洋洋av久久久久久久一区| 欧美日韩伦理在线| 在线亚洲精品| 性一交一乱一区二区洋洋av| 国产精品久久毛片a| 亚洲欧美www| 久久久精品动漫| **性色生活片久久毛片| 久久综合99re88久久爱| 欧美激情免费观看| 亚洲精品一区中文| 欧美日韩一区高清| 亚洲欧美日韩另类精品一区二区三区| 久久国产精品电影| 亚洲大胆美女视频| 欧美精品麻豆| 亚洲欧美精品suv| 蜜桃久久av| 日韩视频在线观看国产| 欧美午夜久久久| 性色一区二区三区| 亚洲国产精品一区二区三区| 一二三区精品福利视频| 国产视频在线观看一区二区三区| 久久9热精品视频| 欧美激情久久久久| 欧美亚洲一区二区在线观看| 海角社区69精品视频| 欧美激情va永久在线播放| 亚洲视频在线观看三级| 乱人伦精品视频在线观看| 99国产精品| 国产亚洲aⅴaaaaaa毛片| 免费高清在线一区| 亚洲欧美韩国| 亚洲国产99精品国自产| 欧美在线国产| 亚洲欧洲综合另类| 国产欧美日韩另类一区| 欧美粗暴jizz性欧美20| 亚洲一区在线看| 亚洲国产婷婷香蕉久久久久久99| 亚洲欧美精品| 日韩午夜在线| 在线不卡免费欧美| 国产欧美亚洲一区| 欧美日韩国产成人在线91| 欧美在线观看天堂一区二区三区| 亚洲激情一区二区| 久久久噜噜噜久久狠狠50岁| 亚洲一区二区三区涩| 91久久嫩草影院一区二区| 国产午夜精品久久久久久久| 欧美精品一区二区在线播放| 久久久久久综合| 亚洲视频播放| 一本色道久久加勒比88综合| 欧美成人亚洲成人| 久久香蕉精品| 久久精品亚洲国产奇米99| 亚洲综合首页| 一区二区三区高清视频在线观看| 亚洲国产婷婷香蕉久久久久久99| 国产亚洲精品bt天堂精选| 国产精品麻豆成人av电影艾秋| 欧美连裤袜在线视频| 你懂的视频一区二区| 久久综合伊人77777| 久久久国产视频91| 欧美在线一级视频| 亚洲免费中文字幕| 亚洲在线免费视频| 午夜伦理片一区| 午夜精彩视频在线观看不卡| 亚洲自拍高清| 欧美亚洲综合网| 欧美在线亚洲综合一区| 欧美一乱一性一交一视频| 亚洲欧美一区二区激情| 香蕉亚洲视频| 久久精品国产清高在天天线| 久久久久久久999| 久久婷婷成人综合色| 欧美91大片| 欧美精品一区二区在线播放| 欧美精品激情在线| 欧美日韩直播| 国产欧美日韩免费| 精品88久久久久88久久久| 在线看片成人| 亚洲乱码国产乱码精品精天堂| 99综合在线| 亚洲一区免费网站| 久久久不卡网国产精品一区| 两个人的视频www国产精品| 欧美mv日韩mv国产网站| 亚洲国产欧美精品| 夜色激情一区二区| 午夜久久黄色| 女女同性女同一区二区三区91| 欧美精品一区三区| 国产精品综合av一区二区国产馆| 狠狠色狠狠色综合日日tαg| 亚洲精品乱码久久久久久日本蜜臀| 亚洲最新在线| 久久精品麻豆| 亚洲日本中文| 欧美伊久线香蕉线新在线| 欧美www视频在线观看| 欧美日韩一区二区在线观看| 国产一区二区三区黄| 最新日韩在线视频| 欧美在线免费视频| 最新热久久免费视频| 午夜久久久久久| 欧美日本三区| 在线看片第一页欧美| 亚洲欧美国产精品va在线观看 | 亚洲精品美女91| 欧美一区1区三区3区公司| 亚洲国产合集| 香蕉久久夜色精品国产| 欧美国产视频在线| 狠狠色丁香婷婷综合影院| 制服丝袜激情欧洲亚洲| 免费影视亚洲| 亚洲欧美另类中文字幕| 欧美风情在线| 亚洲国产成人不卡| 久久精品成人欧美大片古装| 艳女tv在线观看国产一区| 欧美电影免费观看大全| 狠狠狠色丁香婷婷综合激情|