• <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++愛好者
            
				
            太謝謝你的這段龍貝格積分的程序，很贊！謝謝！  回復  更多評論
              
            
			
            
			
			
		
            
	








            






            

				

            






    
    







            
	   
	



AV无码久久久久不卡蜜桃|
久久精品不卡|
狠狠狠色丁香婷婷综合久久五月|
国产精品女同久久久久电影院|
色成年激情久久综合|
人人妻久久人人澡人人爽人人精品|
久久无码人妻一区二区三区|
国产99久久久久久免费看|
久久精品国产精品亚洲精品|
一本久久a久久精品综合夜夜
|
久久九九久精品国产免费直播|
亚洲va久久久久|
国产精品免费久久久久久久久|
国内精品久久久久影院薰衣草
|
久久久SS麻豆欧美国产日韩|
久久人妻少妇嫩草AV无码专区|
久久久久亚洲av成人无码电影|
久久综合噜噜激激的五月天|
中文字幕无码久久精品青草
|
国产成人久久精品二区三区|
亚洲va国产va天堂va久久|
亚洲国产成人精品女人久久久
|
精品无码久久久久久久久久
|
亚洲精品WWW久久久久久|
色成年激情久久综合|
久久精品aⅴ无码中文字字幕不卡|
少妇久久久久久被弄到高潮|
国产成人久久精品区一区二区|
国产69精品久久久久观看软件|
欧美精品一区二区久久|
婷婷久久综合九色综合98|
久久97精品久久久久久久不卡|
亚洲va久久久噜噜噜久久男同|
久久精品一区二区三区AV|
久久久久久久综合狠狠综合|
亚洲精品tv久久久久久久久久|
亚洲国产精品成人久久蜜臀|
日韩十八禁一区二区久久
|
亚洲欧洲精品成人久久曰影片|
久久精品国产精品亜洲毛片|
精品国产乱码久久久久久浪潮|