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

隨筆-48 評(píng)論-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<<"請(qǐng)輸入積分的下、上限：";
cin>>a>>b;
cout<<"輸入你要的精度：";
cin>>e;
Romberg(a,b,e);
}

posted on 2007-06-08 23:28 星夢(mèng)情緣閱讀(3459) 評(píng)論(6) 編輯收藏引用所屬分類: 數(shù)據(jù)結(jié)構(gòu)的所有實(shí)現(xiàn)程序

評(píng)論:

# re: 數(shù)值分析--龍貝格算法 2007-06-11 15:30 | 農(nóng)夫

怎么用的啊？回復(fù) 更多評(píng)論

# re: 數(shù)值分析--龍貝格算法 2007-06-11 20:53 | 星夢(mèng)情緣

double f(double x)
{
return pow(x,3.0/2.0);
}

你把這里的函數(shù)改成你要的函數(shù)就可以了用了的回復(fù) 更多評(píng)論

# re: 數(shù)值分析--龍貝格算法 2007-06-20 23:39 | 曹

有關(guān)于反冪法計(jì)算矩陣特征值的程序嗎回復(fù) 更多評(píng)論

# re: 數(shù)值分析--龍貝格算法 2010-05-12 21:23 | S

對(duì) 函數(shù)的調(diào)用不明確···· 回復(fù) 更多評(píng)論

# re: 數(shù)值分析--龍貝格算法 2010-05-12 21:24 | S

k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(34) : error C2668: “pow”: 對(duì)重載函數(shù)的調(diào)用不明確
1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“l(fā)ong 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> 試圖匹配參數(shù)列表“(int, int)”時(shí)
1>k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(36) : error C2668: “pow”: 對(duì)重載函數(shù)的調(diào)用不明確
1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“l(fā)ong 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> 試圖匹配參數(shù)列表“(int, int)”時(shí)
1>k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(36) : error C2668: “pow”: 對(duì)重載函數(shù)的調(diào)用不明確
1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“l(fā)ong 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> 試圖匹配參數(shù)列表“(int, int)”時(shí)
1>k:\jisuan\longbeigea\longbeigea\longbeigea.cpp(36) : error C2668: “pow”: 對(duì)重載函數(shù)的調(diào)用不明確
1> d:\軟件程序\vc++8.0\vc\include\math.h(575): 可能是“l(fā)ong 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> 試圖匹配參數(shù)列表“(int, int)”時(shí)
1>生成日志保存在“file://k:\jisuan\longbeigea\longbeigea\Debug\BuildLog.htm”
1>longbeigea - 4 個(gè)錯(cuò)誤，0 個(gè)警告
========== 生成: 成功 0 個(gè)，失敗 1 個(gè)，最新 0 個(gè)，跳過(guò) 0 個(gè) ==========
回復(fù) 更多評(píng)論

# re: 數(shù)值分析--龍貝格算法[未登錄](méi) 2011-04-20 19:58 | C++愛(ài)好者

太謝謝你的這段龍貝格積分的程序，很贊！謝謝！回復(fù) 更多評(píng)論

刷新評(píng)論列表

只有注冊(cè)用戶登錄后才能發(fā)表評(píng)論。


相關(guān)文章: C語(yǔ)言版--圖的實(shí)現(xiàn)和各種操作數(shù)值分析--三次樣條插值公式數(shù)值分析--速度最快的牛頓插值 C語(yǔ)言版--二叉樹(shù) 數(shù)值分析--龍貝格算法數(shù)值分析--數(shù)值積分公式數(shù)值分析--解線性方程 C++版--鏈表 C++版--線性表 C語(yǔ)言版--復(fù)數(shù)的實(shí)現(xiàn)

網(wǎng)站導(dǎo)航: 博客園 IT新聞 BlogJava 博問(wèn) Chat2DB 管理

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

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

常用鏈接

留言簿(27)

隨筆分類

隨筆檔案

其他BLOG

搜索

積分與排名

最新評(píng)論

閱讀排行榜

評(píng)論排行榜