| 查看: 254 | 回复: 1 | ||
zhangjin_7铜虫 (小有名气)
|
[求助]
求助romberg调试
|
|
#include #include #define jifen(x) pow(3,x)*pow(x,1.4)*(5*x+7)*sin(pow(x,2))//定义待求积分函数 #define E 0.00001//定义控制迭代误差大小 main() { double newarray[20]={0},oldarray[20]={0};//定义两个数列分别表示迭代是用到的两组数据 double fistarray(double a,double b,int j);//声明一个函数计算新数组的第一个元素 int i=1,j,k;//定义i为新数组中元素个数,j控制将新数组赋值给旧数组，k控制计算新数组 oldarray[0]=jifen(1)+jifen(3);//初始化旧数组中第一个元素 newarray[0]=(3-1)/2*(jifen(1)+jifen(3)-2*jifen(2));//初始化新数组中第一个元素 newarray[1]=(4*newarray[0]-oldarray[0])/(4-1);//初始化新数组中第二个元素 while(fabs(newarray-oldarray[i-1])>=E) { for(j=0;j<=i;j++)//将新数组数据赋给旧数组 oldarray[j]=newarray[j]; i++;//新数组长度增加一位 newarray[0]=0.5*(oldarray[0]+fistarray(1,3,i)); for(k=1;k<=i;k++)//给新数组第一个元素以后的赋值 { newarray[k]=(pow(4,k)*newarray[k-1]-oldarray[k-1])/(pow(4,k)-1); } } printf("%f",newarray); } double fistarray(double a,double b,int j)//帮助计算新数组第一个元素 { int i; double s=0,ss=0; for(i=1;i<=pow(2,j-1);i++) { s+=jifen(a+(2*i-1)*((b-a)/pow(2,j))); } ss=(b-a)/pow(2,j-1)*s; return(ss); } |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Romberg.c
2012-11-29 17:05:35, 1.19 K
» 猜你喜欢
filecode
已经有10人回复
-
听说今天filecode变了
已经有3人回复
-
filecode
已经有12人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有5人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有11人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有30人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab程序调试
已经有5人回复
- 求助一下程序调试过程出现的错误
已经有4人回复
- 求助：污水处理现场调试工程师未来发展前景
已经有17人回复
- 求助调试ACS的latex投稿模板。
已经有5人回复
- 焦化废水的调试问题求助
已经有8人回复
- 【求助】厌氧调试？
已经有8人回复
- 【求助】VS2005数值计算程序调试【急求！！！！加金币】
已经有9人回复
- 【求助】VC++程序算法集调试报错
已经有8人回复
- C数值算法程序大全
已经有7人回复
zhangjin_7
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 207.3
- 散金: 1155
- 红花: 1
- 帖子: 257
- 在线: 99.2小时
- 虫号: 2087044
- 注册: 2012-10-26
- 专业: 化学计量学与化学信息学
2楼2012-11-29 17:06:29
50