锘??xml version="1.0" encoding="utf-8" standalone="yes"?>国产精品无人区,亚洲最新合集,亚洲福利专区http://www.shnenglu.com/newplan/archive/2008/05/07/49127.html灞辨硥寮歡灞辨硥寮歡Wed, 07 May 2008 07:05:00 GMThttp://www.shnenglu.com/newplan/archive/2008/05/07/49127.htmlhttp://www.shnenglu.com/newplan/comments/49127.htmlhttp://www.shnenglu.com/newplan/archive/2008/05/07/49127.html#Feedback0http://www.shnenglu.com/newplan/comments/commentRss/49127.htmlhttp://www.shnenglu.com/newplan/services/trackbacks/49127.html/*
Romberg Algorithm 
寮鍙戣?newplan
寮鍙戞棩鏈?
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 閫掓帹鍑芥暟*/
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;
}


灞辨硥寮歡 2008-05-07 15:05 鍙戣〃璇勮
]]>
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            久久久久久日产精品| 黄色一区二区在线| 欧美诱惑福利视频| 一区二区三区四区五区在线| 国产在线观看91精品一区| 国产精品成人在线| 久久久久久久999精品视频| 亚洲校园激情| 亚洲午夜激情网站| 一区二区三区四区国产| 亚洲精选国产| 一区二区三区精品视频| 中文有码久久| 亚洲欧美激情四射在线日| 一区二区电影免费在线观看| 亚洲精品一区在线观看| 亚洲每日在线| 亚洲精选大片| 亚洲网站在线| 午夜精品福利一区二区三区av| 在线亚洲国产精品网站| 亚洲视频欧美在线| 亚洲精品欧洲| 一区二区三区日韩精品视频| 一区二区三区精品在线| 亚洲亚洲精品在线观看| 午夜欧美视频| 久久久久国内| 欧美日韩国产限制| 国产一区二区三区奇米久涩| 国产精品99久久久久久www| 久久只有精品| 亚洲一区二区高清| 欧美www视频在线观看| 国产精品视频专区| 夜夜嗨av一区二区三区中文字幕| 久久久久久久一区二区三区| 日韩一区二区精品| 欧美jizzhd精品欧美巨大免费| 国产精品天美传媒入口| 中国女人久久久| 欧美激情精品久久久久久黑人 | 极品日韩久久| 久久精品人人做人人爽| 亚洲午夜一二三区视频| 欧美女同视频| 亚洲精品在线一区二区| 欧美岛国在线观看| 久久午夜影视| 亚洲福利免费| 女生裸体视频一区二区三区| 久久久www成人免费毛片麻豆| 国产区在线观看成人精品| 性欧美办公室18xxxxhd| 亚洲先锋成人| 国产精品久久99| 亚洲影院在线观看| 亚洲校园激情| 国产亚洲欧美aaaa| 久久噜噜亚洲综合| 久久午夜精品一区二区| 亚洲第一精品福利| 欧美国产视频在线观看| 久久综合色8888| 亚洲人成在线免费观看| 亚洲人成人一区二区在线观看| 久久蜜桃精品| 亚洲精选在线观看| 日韩午夜电影av| 国产精品网站视频| 久久蜜桃av一区精品变态类天堂| 久久精品国产96久久久香蕉| 在线观看视频一区二区| 欧美激情一区二区在线| 欧美福利网址| 亚洲欧美一区二区三区久久 | 中国日韩欧美久久久久久久久| 国产精品毛片a∨一区二区三区| 欧美伊人久久| 久久影院午夜论| 99国产精品99久久久久久| 亚洲视频免费在线观看| 激情一区二区三区| 亚洲日本成人在线观看| 国产精品夜夜夜| 欧美阿v一级看视频| 欧美日韩综合在线| 久久久精品国产一区二区三区 | 欧美在线观看www| 亚洲日本国产| 亚洲欧美日韩成人| 最新亚洲电影| 亚洲综合精品自拍| 亚洲国内在线| 亚洲综合第一| 亚洲三级毛片| 欧美一区二区三区免费观看| 日韩亚洲在线| 久久精品日产第一区二区| 亚洲美女诱惑| 性18欧美另类| 中日韩美女免费视频网址在线观看| 亚洲欧美日韩在线播放| 亚洲精品视频免费| 午夜精品久久久久久久久| 日韩亚洲一区在线播放| 香蕉精品999视频一区二区| 亚洲欧洲视频在线| 亚洲欧美国产精品va在线观看| 国产精品成人v| 亚洲国产天堂久久综合| 国产一区二区三区四区三区四| 亚洲福利视频专区| 红桃视频国产精品| 亚洲综合国产激情另类一区| 亚洲人人精品| 久久午夜精品一区二区| 久久久久久一区| 国产日韩欧美在线观看| 在线视频你懂得一区二区三区| 亚洲国产另类久久精品| 欧美在线免费一级片| 午夜欧美不卡精品aaaaa| 欧美日韩一区二区三| 亚洲黄色高清| 亚洲精品美女在线| 欧美激情1区2区| 亚洲福利视频一区| 亚洲精品资源美女情侣酒店| 免费91麻豆精品国产自产在线观看| 久久免费视频网站| 精品91久久久久| 久久久爽爽爽美女图片| 老司机67194精品线观看| 激情五月婷婷综合| 麻豆成人综合网| 欧美激情小视频| 日韩网站在线| 欧美日韩中字| 亚洲综合视频网| 久久精品99国产精品| 国产自产精品| 美乳少妇欧美精品| 亚洲久色影视| 亚洲一区欧美| 国产午夜精品一区二区三区视频| 性欧美xxxx视频在线观看| 久久久久欧美精品| 亚洲高清视频在线观看| 欧美成人高清视频| 日韩亚洲在线观看| 久久精品电影| 亚洲第一福利在线观看| 欧美激情一区二区三区在线视频观看| 亚洲国产精品精华液2区45| 99精品视频一区二区三区| 欧美色网一区二区| 欧美一区二区视频免费观看| 美国十次成人| 欧美激情第9页| 翔田千里一区二区| 亚洲欧美日本日韩| 久久久国产一区二区| 久久综合伊人77777蜜臀| 老司机精品视频一区二区三区| 欧美一级一区| 女生裸体视频一区二区三区| 欧美日韩成人激情| 国产精品卡一卡二卡三| 加勒比av一区二区| 99re6热只有精品免费观看| 一区电影在线观看| 久久精品视频在线观看| 欧美大色视频| 亚洲综合视频1区| 男女av一区三区二区色多| 国产精品jvid在线观看蜜臀| 黄色综合网站| 亚洲永久网站| 欧美国产日韩精品| 亚洲欧美日韩国产中文在线| 美国十次了思思久久精品导航| 欧美日韩免费高清一区色橹橹| 国产午夜精品美女毛片视频| 亚洲人午夜精品免费| 欧美在线在线| 亚洲精品一区在线观看香蕉| 久久精品99无色码中文字幕| 欧美日韩免费一区| 亚洲国产高清视频| 久久国产天堂福利天堂| 日韩一级黄色大片| 美乳少妇欧美精品| 伊人久久久大香线蕉综合直播 | 国产亚洲网站| 先锋影音久久| 99精品热视频只有精品10| 免费h精品视频在线播放| 国内精品久久久久影院色| 香蕉成人啪国产精品视频综合网|