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

隨筆-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 星夢情緣 閱讀(3459) 評論(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++愛好者
太謝謝你的這段龍貝格積分的程序,很贊!謝謝!  回復  更多評論
  
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            美女主播一区| 亚洲福利视频三区| 亚洲欧美久久| 亚洲欧美成人网| 久久精品午夜| 老司机免费视频一区二区三区| 久久成人一区| 一本一本a久久| 影音先锋亚洲精品| 久久av资源网| 亚洲免费观看| 欧美一区亚洲二区| 欧美不卡三区| 亚洲精品一区在线观看| 精东粉嫩av免费一区二区三区| 亚洲天堂av图片| 久久动漫亚洲| 久久久精品国产免费观看同学| 日韩视频在线播放| 日韩视频精品在线观看| 日韩视频在线观看国产| 国产精品高清在线观看| 亚洲一区二区三区四区五区黄| 亚洲国产高清自拍| 亚洲人精品午夜| 日韩天堂在线视频| 午夜精品福利在线| 欧美成人一区二区三区在线观看| 亚洲女人天堂av| 亚洲特级片在线| 久久久久久九九九九| 亚洲一区二区在线免费观看视频| 9人人澡人人爽人人精品| 亚洲精品日韩激情在线电影| a4yy欧美一区二区三区| 欧美性猛交xxxx免费看久久久| 久久久精品网| 欧美日韩精品免费看| 国产精品国产馆在线真实露脸| 激情久久久久| 午夜视频在线观看一区二区三区| 久久一区欧美| 欧美在线视频不卡| 在线观看一区二区精品视频| 国产日韩欧美综合| 亚洲视频一区二区| 久久精品夜色噜噜亚洲a∨ | 欧美韩日一区二区三区| 亚洲精品中文字| 欧美高清视频在线播放| 亚洲国产裸拍裸体视频在线观看乱了中文 | 韩国女主播一区二区三区| 久久成人免费电影| 久久一区视频| 欧美一级片久久久久久久| 亚洲综合日本| 韩国三级在线一区| 欧美激情中文字幕一区二区 | 亚洲黄色性网站| 亚洲一区二区在线观看视频| 久久久99免费视频| 黄色国产精品| 亚洲伦理在线观看| 亚洲第一天堂无码专区| 久久九九久精品国产免费直播| 亚洲免费av片| 欧美激情视频一区二区三区在线播放| 国内精品国产成人| 国产精品国产成人国产三级| 在线亚洲精品| 日韩午夜电影av| 欧美性大战久久久久久久蜜臀| 欧美成人一区二区三区片免费| 亚洲欧美一区二区原创| 日韩网站免费观看| 一区二区三区久久久| 亚洲人成人99网站| 久久国产66| 小嫩嫩精品导航| 国产精品成人v| 国产精品永久| 亚洲激情专区| 日韩一区二区精品视频| 麻豆精品在线观看| 久久久999国产| 国产精品自拍网站| 欧美夜福利tv在线| 亚洲在线第一页| 国产精品xxxav免费视频| 亚洲宅男天堂在线观看无病毒| 久久国产一区二区三区| 日韩午夜免费视频| 老牛嫩草一区二区三区日本 | 韩日欧美一区二区三区| 亚洲男人第一av网站| 久久精品视频免费观看| 欧美激情综合五月色丁香| 在线免费观看欧美| 亚洲欧美卡通另类91av | 亚洲香蕉伊综合在人在线视看| 亚洲国产一区二区精品专区| 老司机午夜免费精品视频 | 欧美一区激情| 亚洲国产婷婷香蕉久久久久久99 | 亚洲国产一区二区a毛片| 欧美三级第一页| 日韩亚洲欧美在线观看| 蜜臀久久99精品久久久画质超高清| 暖暖成人免费视频| 在线综合视频| 欧美亚洲一级| 欧美不卡激情三级在线观看| 日韩视频一区| 欧美在线观看视频一区二区| 欧美一区二区高清在线观看| 国产精品免费看| 亚洲国产精品va| 欧美亚洲系列| 亚洲网友自拍| 欧美美女bb生活片| 亚洲国产另类久久精品| 亚洲免费福利视频| 欧美第一黄网免费网站| 久久久久久电影| 午夜久久久久| 欧美区亚洲区| 一区二区在线视频| 亚洲美女在线观看| 国产精品丝袜久久久久久app| 欧美亚洲自偷自偷| 欧美成人午夜激情视频| 国产精品区一区二区三区| 影音欧美亚洲| 久久黄色影院| 亚洲小视频在线观看| 欧美成人一区二区在线| 国产一区二区激情| 欧美在线亚洲| 亚洲主播在线| 国产日韩欧美中文| 亚洲深夜福利| 亚洲国产黄色| 亚洲永久视频| 欧美成人一区二区| 久久午夜羞羞影院免费观看| 国产欧美日韩亚洲精品| 久久精品久久99精品久久| 欧美一区二区精品在线| 永久久久久久| 日韩午夜电影| 欧美一进一出视频| 欧美大片免费久久精品三p | 久久精品30| 亚洲私人影吧| 国产精品永久免费观看| 亚洲一区国产| 亚洲一区二区视频| 国产自产女人91一区在线观看| 欧美亚洲视频| 亚洲一区二区三区中文字幕 | 欧美中文字幕不卡| 国产精品视频你懂的| 久久久久久电影| 欧美福利视频| 亚洲欧美精品伊人久久| 先锋影音国产精品| 亚洲国产欧美另类丝袜| 亚洲精品在线观看免费| 欧美性jizz18性欧美| 亚洲伊人一本大道中文字幕| 欧美亚洲色图校园春色| 亚洲高清网站| 午夜精品久久久久| 亚洲欧洲在线一区| 亚洲一级一区| 亚洲全部视频| 亚洲欧美日韩国产中文| 亚洲欧洲一区| 久久不射2019中文字幕| 亚洲天堂免费观看| 麻豆精品精华液| 久久精品综合| 国产精品一区二区久久| 亚洲激情另类| 一色屋精品视频免费看| 亚洲欧美欧美一区二区三区| 亚洲看片网站| 欧美国产第一页| 欧美激情国产高清| 亚洲国产精品va在线观看黑人| 亚洲一二三区在线| 欧美亚洲综合另类| 在线亚洲一区二区| 亚洲欧洲精品一区二区| 噜噜噜在线观看免费视频日韩| 男人天堂欧美日韩| 国产一区二区三区在线观看网站| 亚洲视频在线观看| 欧美影院成人| 极品尤物久久久av免费看|