锘??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>
            国产真实精品久久二三区| 国产精品免费网站| 亚洲激情视频网| 欧美韩日高清| 欧美午夜在线观看| 欧美一级日韩一级| 久久九九免费| 亚洲精品资源| 亚洲自拍偷拍麻豆| 亚洲高清视频在线观看| 亚洲欧洲一区二区在线观看| 欧美日韩性生活视频| 久久精品国产在热久久| 免费观看久久久4p| 亚洲小少妇裸体bbw| 欧美一级片久久久久久久| 伊人久久亚洲热| 亚洲精品在线一区二区| 国产亚洲福利一区| 亚洲国产精品综合| 国产精品一区二区在线观看| 欧美电影打屁股sp| 国产精品嫩草影院一区二区| 久久综合给合久久狠狠色 | 欧美一区二区三区婷婷月色| 久久久一二三| 午夜在线精品| 欧美精品在线观看| 久久亚洲精品网站| 国产精品久久777777毛茸茸| 欧美国产日本韩| 国产亚洲精品资源在线26u| 91久久久亚洲精品| 在线日本高清免费不卡| 亚洲综合丁香| 亚洲视频日本| 你懂的视频欧美| 老司机一区二区三区| 国产精品美女久久福利网站| 亚洲黄色成人网| 影音先锋久久久| 亚洲欧美一区二区激情| 一区二区三区久久精品| 免费久久久一本精品久久区| 久久免费视频在线| 国产午夜亚洲精品不卡| 中文网丁香综合网| 一区二区三区日韩精品| 欧美电影专区| 欧美高清在线| 亚洲国产一二三| 麻豆精品在线视频| 欧美成人激情视频| 在线观看日韩一区| 久久xxxx精品视频| 久久久91精品| 国模精品一区二区三区| 午夜激情久久久| 久久久久久有精品国产| 国产丝袜一区二区三区| 午夜精品一区二区三区在线| 午夜在线精品| 国产精品一区二区女厕厕| 亚洲一级黄色av| 欧美一区二区三区在| 国产乱码精品一区二区三区不卡| 亚洲午夜女主播在线直播| 亚洲一区二区精品在线| 国产精品一区二区三区免费观看| 亚洲一区在线直播| 久久午夜av| 最新亚洲一区| 欧美日精品一区视频| 亚洲一区精品视频| 久久国产夜色精品鲁鲁99| 激情综合亚洲| 欧美国产视频日韩| 一区二区激情小说| 久久精品国语| 亚洲黄色免费网站| 欧美网站大全在线观看| 午夜精品久久久久久久蜜桃app| 久久久久国产免费免费| 亚洲国产裸拍裸体视频在线观看乱了中文 | 国产精品有限公司| 久久久久久电影| 亚洲精品免费一二三区| 篠田优中文在线播放第一区| 国产一区二区中文| 欧美成人在线影院| 亚洲香蕉成视频在线观看 | 亚洲一二区在线| 国外成人在线视频| 欧美精品一区二区三| 亚洲伊人一本大道中文字幕| 美女视频网站黄色亚洲| 一本一本大道香蕉久在线精品| 国产精品一区二区久久国产| 久久亚洲精选| 亚洲一区视频在线| 欧美激情无毛| 欧美一区二区三区在线看| 亚洲国产精品免费| 国产精品高潮呻吟久久| 狼狼综合久久久久综合网| 亚洲色图制服丝袜| 亚洲成色777777女色窝| 欧美伊人久久久久久午夜久久久久| 亚洲国产成人久久综合| 国产精品无码永久免费888| 欧美wwwwww| 久久国产综合精品| 亚洲无限av看| 日韩视频一区二区三区在线播放免费观看| 久久激情五月激情| 亚洲女同精品视频| 日韩视频在线一区二区| 亚洲大胆在线| 狠色狠色综合久久| 国产乱人伦精品一区二区| 欧美日韩黄色一区二区| 欧美成在线观看| 久久久中精品2020中文| 性8sex亚洲区入口| 亚洲一区二区三区国产| 夜夜嗨av一区二区三区四区| 亚洲国产日韩欧美在线动漫| 老鸭窝毛片一区二区三区| 久久av老司机精品网站导航| 亚洲无限av看| 亚洲色图制服丝袜| 一区二区三区你懂的| 亚洲免费观看高清在线观看| 亚洲日本理论电影| 亚洲精品日韩综合观看成人91| 伊人久久综合97精品| 激情av一区二区| 激情久久久久久久| 一区二区三区在线免费观看| 韩日午夜在线资源一区二区| 国产一区二区中文字幕免费看| 国产三级精品在线不卡| 国产欧美精品在线| 国产偷国产偷精品高清尤物| 国产综合精品一区| 一区三区视频| 亚洲国产精彩中文乱码av在线播放| 精品动漫av| 91久久精品网| 99精品国产一区二区青青牛奶 | 合欧美一区二区三区| 精品成人国产| 亚洲精品韩国| 亚洲视频高清| 久久国产一区二区| 麻豆精品视频在线观看| 亚洲高清免费| 夜夜嗨av一区二区三区网站四季av | 亚洲伦理在线观看| 一本久道久久久| 亚洲影院一区| 久久五月天婷婷| 欧美精品久久一区| 国产精品中文字幕在线观看| 韩国女主播一区二区三区| 亚洲黄一区二区三区| 亚洲视频在线观看视频| 久久精品午夜| 亚洲三级性片| 欧美在线|欧美| 欧美久久精品午夜青青大伊人| 国产精品久久久久久久久久尿 | 国产一区二区三区免费在线观看| 永久555www成人免费| 一区二区三区国产在线| 久久精品国产精品亚洲综合| 亚洲高清自拍| 午夜精品一区二区三区在线视| 免费在线亚洲欧美| 国产精品无人区| 亚洲免费av电影| 午夜精品久久久久久久99热浪潮| 另类欧美日韩国产在线| 亚洲无亚洲人成网站77777 | 亚洲乱亚洲高清| 久久国产乱子精品免费女| 亚洲国产一区在线| 久久成人国产精品| 国产精品国产三级欧美二区 | 欧美特黄一区| 亚洲电影视频在线| 久久精品91| 亚洲天堂激情| 欧美精品一区二区在线观看| 影音先锋久久久| 久久久久久免费| 亚洲欧美国产日韩天堂区| 欧美美女bb生活片| 亚洲日本免费电影| 欧美国产高清|