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

newplan

阿基米德在洗澡時(shí)發(fā)現(xiàn)浮力原理,高興得來不及穿上褲子,跑到街上大喊:Eureka(我找到了)。
posts - 39, comments - 26, trackbacks - 0, articles - 4
  C++博客 :: 首頁 :: 新隨筆 :: 聯(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>
            国产精品久久久久久久app| 久久久噜噜噜久久久| 欧美日韩精品免费观看视频完整| 91久久在线播放| 中文日韩在线| 国产日韩精品视频一区二区三区| 久久精品2019中文字幕| 亚洲电影自拍| 亚洲欧美日韩视频一区| 国产一区二区看久久| 鲁大师影院一区二区三区| 日韩天天综合| 久久综合电影一区| 一区二区三区鲁丝不卡| 国产日本欧美视频| 欧美精品成人在线| 欧美夜福利tv在线| 亚洲黄一区二区三区| 欧美一区二区三区另类 | 日韩一级欧洲| 久久免费偷拍视频| 亚洲一级二级| 亚洲成人资源网| 国产精品久久久久久妇女6080| 久久香蕉国产线看观看网| 中文国产一区| 亚洲国产精品一区二区www在线| 性欧美xxxx视频在线观看| 亚洲欧洲在线视频| 国产色视频一区| 欧美日韩免费高清| 麻豆乱码国产一区二区三区| 亚洲综合欧美| 亚洲免费黄色| 亚洲电影激情视频网站| 久久久一本精品99久久精品66| 亚洲亚洲精品在线观看| 亚洲激情在线观看视频免费| 国产农村妇女精品一二区| 欧美日韩国产一区二区| 蜜臀av国产精品久久久久| 香蕉久久夜色精品国产| 亚洲天堂av图片| 99精品热视频只有精品10| 欧美激情精品久久久久久| 久久乐国产精品| 欧美有码在线视频| 亚洲欧美春色| 亚洲一区二区三区乱码aⅴ| 亚洲精美视频| 亚洲国产激情| 亚洲成色999久久网站| 国产亚洲成av人在线观看导航| 国产精品美女久久久| 欧美三级不卡| 欧美日韩三级在线| 欧美日韩三区四区| 欧美日韩国产综合一区二区| 欧美国产在线视频| 欧美国产第一页| 欧美激情在线播放| 欧美日本免费一区二区三区| 欧美高清日韩| 欧美欧美全黄| 欧美日韩综合| 国产精品成人一区二区艾草| 国产精品爱久久久久久久| 欧美日韩一区高清| 国产精品magnet| 国产精品久久综合| 国产日韩在线看片| 黄色工厂这里只有精品| 136国产福利精品导航网址| 亚洲福利在线观看| 日韩视频免费观看| 国产精品99久久久久久白浆小说| 亚洲视频欧美在线| 亚洲欧美日韩国产一区二区三区| 午夜精品av| 久久中文字幕导航| 亚洲高清资源| 一本一本久久| 欧美亚洲自偷自偷| 久久综合激情| 欧美日韩成人免费| 国产精品永久入口久久久| 国外精品视频| 亚洲精品乱码久久久久久日本蜜臀| 日韩视频一区二区| 欧美在线播放一区| 欧美大香线蕉线伊人久久国产精品| 亚洲日本欧美日韩高观看| 一区二区三区精品在线 | 国产精品高潮呻吟久久av黑人| 国产精品亚洲视频| 亚洲福利在线视频| 亚洲图片欧美日产| 久久先锋影音av| 亚洲人精品午夜| 午夜久久影院| 欧美激情在线观看| 国产人成一区二区三区影院| 亚洲国产日韩一区| 香港成人在线视频| 欧美激情乱人伦| 亚洲欧美变态国产另类| 免费一级欧美在线大片| 国产精品美女久久久免费| 亚洲高清免费| 久久成人一区二区| 亚洲人成亚洲人成在线观看图片| 午夜精品久久久| 欧美国产日韩一区二区三区| 国产日韩精品视频一区| 一区二区三区产品免费精品久久75| 久久精品国产亚洲精品| 99精品国产一区二区青青牛奶| 久久精品国产精品亚洲| 欧美日韩一区二区在线观看视频| 在线激情影院一区| 欧美一级淫片aaaaaaa视频| 亚洲电影天堂av| 欧美亚洲在线| 国产精品美腿一区在线看| 亚洲久久在线| 蜜臀99久久精品久久久久久软件 | 欧美色区777第一页| 在线观看成人网| 欧美一区二区三区四区在线观看地址| 亚洲动漫精品| 欧美一区亚洲| 国产日韩精品一区二区浪潮av| 一区二区av| 亚洲国内自拍| 久久伊人亚洲| 好吊日精品视频| 羞羞色国产精品| 一区二区三区欧美激情| 欧美日韩福利视频| 亚洲伦理精品| 亚洲国产精品一区制服丝袜| 久久人人97超碰国产公开结果| 国产九区一区在线| 香蕉久久一区二区不卡无毒影院| 日韩午夜视频在线观看| 欧美激情成人在线视频| 亚洲精品欧美| 亚洲电影免费在线| 欧美国产第二页| 亚洲免费观看高清完整版在线观看熊| 欧美成人亚洲成人日韩成人| 久久久久.com| 亚洲国产欧洲综合997久久| 欧美jjzz| 欧美wwwwww| 日韩一级黄色片| 91久久久在线| 欧美日韩在线播放三区| 亚洲一区二区3| 一区二区精品国产| 国产精品观看| 免费观看成人鲁鲁鲁鲁鲁视频| 亚洲国产成人久久| 欧美激情在线观看| 欧美精品久久久久久| 99精品视频网| 中日韩在线视频| 国产精品一区二区你懂的| 久久精品理论片| 久久精品免费电影| 亚洲高清视频中文字幕| 亚洲黄色尤物视频| 欧美视频二区| 久久av红桃一区二区小说| 久久精品最新地址| 最新69国产成人精品视频免费 | 亚洲午夜黄色| 国产亚洲欧美激情| 美女主播一区| 欧美精品一区二区三区蜜臀 | 久久裸体艺术| 亚洲美女精品成人在线视频| 亚洲精品久久久一区二区三区| 欧美日韩在线免费观看| 欧美在线一区二区| 麻豆久久婷婷| 亚洲制服少妇| 久久久久国内| 一区二区三区免费在线观看| 午夜精品偷拍| 亚洲精品亚洲人成人网| 亚洲色无码播放| 亚洲第一精品夜夜躁人人爽| 99国产精品| 136国产福利精品导航网址应用| 亚洲精品激情| 国产婷婷97碰碰久久人人蜜臀| 欧美激情一区| 国产九九精品视频| 亚洲欧洲精品一区二区|