• <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>
            隨筆-48  評論-259  文章-1  trackbacks-0

            #include"iostream.h"
            #include"math.h"

            double f(double x)
            {
               return pow(x,3.0/2.0);
            }
            double T(double n,double a,double b)
            {
              double sum=0;
              for(int i=0;2*i+1<=n;i++)
              {
                sum+=f(a+double(2*i+1)/n);
              }
              if(n==1) return 0.5*(f(a)+f(b));
              return 0.5*T(n/2,a,b)+sum*(1.0/n);
            }
            void Romberg(double a,double b,double e)
            {
              double Old[10]={0.0},New[10]={0.0};
              double h=b-a;
              New[0]=h*(f(a)+f(b))/2.0;
              int j;
              for(int i=1;;i++)
              {
                for(j=0;j<i;j++) Old[j]=New[j];
             New[0]=T(pow(2,i),a,b);
             for(j=1;j<=i;j++)
              New[j]=(pow(4,j)*New[j-1]/(pow(4,j)-1)-Old[j-1]/(pow(4,j)-1));
             if(fabs(New[i]-Old[i-1])<e) break;
              }
              cout<<New[i];
            }

            void main()
            {
              double a,b; 
              double e;
              cout<<"請輸入積分的下、上限:";
              cin>>a>>b;
              cout<<"輸入你要的精度:";
              cin>>e;
              Romberg(a,b,e);
            }

            posted on 2007-06-08 23:28 星夢情緣 閱讀(3444) 評論(6)  編輯 收藏 引用 所屬分類: 數據結構的所有實現程序

            評論:
            # re: 數值分析--龍貝格算法 2007-06-11 15:30 | 農夫
            怎么用的啊?  回復  更多評論
              
            # re: 數值分析--龍貝格算法 2007-06-11 20:53 | 星夢情緣
            double f(double x)
            {
            return pow(x,3.0/2.0);
            }

            你把這里的函數改成你要的函數就可以了用了的  回復  更多評論
              
            # re: 數值分析--龍貝格算法 2007-06-20 23:39 |
            有關于反冪法計算矩陣特征值的程序嗎  回復  更多評論
              
            # re: 數值分析--龍貝格算法 2010-05-12 21:23 | S
            對 函數的 調用不明確····  回復  更多評論
              
            # re: 數值分析--龍貝格算法 2010-05-12 21:24 | S
            k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(34) : error C2668: “pow”: 對重載函數的調用不明確
            1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“long double pow(long double,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(527): 或 “float pow(float,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(489): 或 “double pow(double,int)”
            1> 試圖匹配參數列表“(int, int)”時
            1>k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(36) : error C2668: “pow”: 對重載函數的調用不明確
            1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“long double pow(long double,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(527): 或 “float pow(float,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(489): 或 “double pow(double,int)”
            1> 試圖匹配參數列表“(int, int)”時
            1>k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(36) : error C2668: “pow”: 對重載函數的調用不明確
            1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“long double pow(long double,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(527): 或 “float pow(float,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(489): 或 “double pow(double,int)”
            1> 試圖匹配參數列表“(int, int)”時
            1>k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(36) : error C2668: “pow”: 對重載函數的調用不明確
            1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“long double pow(long double,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(527): 或 “float pow(float,int)”
            1> d:\軟件程序\vc++8.0\vc\include\math.h(489): 或 “double pow(double,int)”
            1> 試圖匹配參數列表“(int, int)”時
            1>生成日志保存在“file://k:\jisuan\longbeigea\longbeigea\Debug\BuildLog.htm”
            1>longbeigea - 4 個錯誤,0 個警告
            ========== 生成: 成功 0 個,失敗 1 個,最新 0 個,跳過 0 個 ==========
              回復  更多評論
              
            # re: 數值分析--龍貝格算法[未登錄] 2011-04-20 19:58 | C++愛好者
            太謝謝你的這段龍貝格積分的程序,很贊!謝謝!  回復  更多評論
              
            久久精品九九亚洲精品天堂| 情人伊人久久综合亚洲| 伊人色综合久久| 久久久久久无码国产精品中文字幕 | 少妇久久久久久被弄到高潮| 国产成人精品三上悠亚久久| 成人综合伊人五月婷久久| 精品人妻伦九区久久AAA片69| 亚洲人成网亚洲欧洲无码久久| 亚洲国产高清精品线久久| 国内精品久久久久| 久久亚洲精品国产精品婷婷| 久久综合久久自在自线精品自| 亚洲伊人久久综合中文成人网| 国产成人久久精品激情| 久久午夜福利无码1000合集| 国产精自产拍久久久久久蜜| 久久不见久久见免费视频7| 久久久久久久久66精品片| 久久综合中文字幕| 久久久久亚洲AV成人片 | 国产精品丝袜久久久久久不卡| 久久久久久综合网天天| 久久人人爽人人人人爽AV| 国产精品免费久久久久影院| 久久国产精品国产自线拍免费| 性高湖久久久久久久久| 久久人人爽人人人人爽AV| 人妻丰满?V无码久久不卡| 久久精品国产亚洲一区二区三区| 国产精品99久久99久久久| avtt天堂网久久精品| 99久久免费国产特黄| 国产精品18久久久久久vr| 久久国产高潮流白浆免费观看| 久久人人爽人人爽人人片av麻烦 | 亚洲精品午夜国产VA久久成人| 久久午夜福利电影| 午夜欧美精品久久久久久久| 国产欧美久久久精品影院| 日产精品久久久一区二区|