| 查看: 2226 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】求积分程序中被积函数问题 已有5人参与
|
|||
|
我用Simpson方法计算定积分。在程序中有两个地方用到积分，而被积分函数作为子程序调用，在计算积分的程序中该怎么定义被积函数呢？ 这是 辛普森算法 的程序： double fsimp(a,b,eps) /*辛普森算法*/ double a,b,eps; /*a为积分下限，b为积分上限，eps是希望达到的精度*/ { int n,k; double h,t1,t2,s1,s2,ep,p,x; n=1; h=b-a; t1=h*(fsimpf(a)+fsimpf(b))/2.0; /*用梯形公式求出一个大概的估值*/ s1=t1; printf("s1=%e\n",s1); ep=eps+1.0; while (ep>=eps) { p=0.0; for (k=0;k<=n-1;k++) { x=a+(k+0.5)*h; p=p+fsimpf(x); } t2=(t1+h*p)/2.0; s2=(4.0*t2-t1)/3.0; printf("s2=%e\n",s2); ep=fabs(s2-s1); t1=t2; s1=s2; n=n+n; h=h/2.0; } return(s2); } 这是被积函数的程序： double fsimpf(r) double r; /* 被积分的函数*/ { double y; y=被积函数f(x)的表达式； return(y); } 两个被积函数不能都用fsimpf（）表示吧，可不用它，计算积分时调用被积函数又是用fsimpf（）表示的： 怎么区分啊 |
回复此楼» 猜你喜欢
26年申博自荐-计算机视觉
已经有4人回复
-
考博
已经有4人回复
-
药化及相关博士的申请
已经有3人回复
-
一篇MDPI论文改变了学习工作和生活
已经有4人回复
-
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇
已经有4人回复
-
太白金星有点烦
已经有3人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有3人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有3人回复
-
有没有快的中文核心比较快录用的，纳米材料光催化
已经有4人回复
15楼2010-06-04 08:54:10
simthhuang
木虫 (小有名气)
- 程序强帖: 1
- 应助: 0 (幼儿园)
- 金币: 2231.2
- 散金: 24
- 红花: 3
- 帖子: 252
- 在线: 69.3小时
- 虫号: 843987
- 注册: 2009-09-09
- 性别: GG
- 专业: 认知科学及智能信息处理
2楼2010-05-21 17:04:29
3楼2010-05-21 21:46:17
simthhuang
木虫 (小有名气)
- 程序强帖: 1
- 应助: 0 (幼儿园)
- 金币: 2231.2
- 散金: 24
- 红花: 3
- 帖子: 252
- 在线: 69.3小时
- 虫号: 843987
- 注册: 2009-09-09
- 性别: GG
- 专业: 认知科学及智能信息处理
4楼2010-05-22 19:41:04