速求C语言或者MATLAB高手!急急急!

用龙贝格方法计算SINX/X的积分,积分区间【0,1】,使截断误差不超过0.5*10的-6次方。谢谢!!!!

#include "stdio.h"
#include "math.h"
main()
{
double x,s=0,i,n=100000;
for(i=1;i<n;i++)
{
x=i/n;
s=s+sin(x)/x*(1.0/n);

}

printf("%lf\n",s);

}
用的是C语言编的,没用什么龙贝格方法,不知道能不能帮得上你
温馨提示:答案为网友推荐,仅供参考
第1个回答  2009-03-29
你有误差限制的话,只要调整龙贝格方法里面的步长就可以。

他的误差是

e=f(h),误差是步长的一个函数。

跟你你的误差,计算出步长,然后再使用Matlab数值积分即可。
第2个回答  2009-03-29
参考一下,我以前帮别人修改的龙贝格程序。

http://zhidao.baidu.com/question/79401173.html
相似回答