应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 【求助】求积分程序中被积函数问题
51/1返回列表
查看: 2994  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

muqiuyan

木虫 (小有名气)

[交流] 【求助】求积分程序中被积函数问题 已有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（）表示的：


怎么区分啊
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-05-21 14:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

simthhuang

木虫 (小有名气)
★
resonant(金币+1):感谢参与：-） 2010-05-22 21:01:25
第二个被积函数就是被积函数啊，就是积分号和dx之间那部分啊。f(a) f(b) f(x) 分别是用初值a,终值b和中间迭代的数x调用被积函数。MM，我的理解对吗？
赞一下(1人) 回复此楼
一贯知足二目远眺三餐有节四季不懒五谷皆食六欲不张七分忍让八方交往九薄烟断十分坦荡
4楼2010-05-22 19:41:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

simthhuang

木虫 (小有名气)
★ ★
muqiuyan(金币+1):谢谢参与
resonant(金币+1):感谢积极参与：-） 2010-05-21 17:38:39
哪里有两个被积函数啊？明明就一个嘛！
赞一下(2人) 回复此楼
一贯知足二目远眺三餐有节四季不懒五谷皆食六欲不张七分忍让八方交往九薄烟断十分坦荡
2楼2010-05-21 17:04:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muqiuyan

木虫 (小有名气)
引用回帖:
Originally posted by simthhuang at 2010-05-21 17:04:29:
哪里有两个被积函数啊？明明就一个嘛！

主程序里有两个被积函数，我只列了两个子程序，一个是计算积分的，一个是被积函数的子程序，想问问第二个被积函数怎么弄，积分子程序里的f(a) f(b) f(x) 之类的需要调要被积函数子程序的地方要不要调整。
赞一下 回复此楼
3楼2010-05-21 21:46:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muqiuyan

木虫 (小有名气)
resonant:调用函数没有什么特别的吧？区分只要函数名不同不就可以区分么？ 2010-05-24 11:09:30
引用回帖:
Originally posted by simthhuang at 2010-05-22 19:41:04:
第二个被积函数就是被积函数啊，就是积分号和dx之间那部分啊。f(a) f(b) f(x) 分别是用初值a,终值b和中间迭代的数x调用被积函数。MM，我的理解对吗？

可能是我表达的不清楚。我要弄一个很复杂的程序，有很多公式，其中有两个计算积分的式子，所以我编了一个计算积分的子程序，就是fsimp(a,b,eps)；这其中就涉及到调用被积函数，那我的两个被积函数怎么调用呢，怎么区分开呢

PS：你是学计算机的吗，能认识下吗？
赞一下(1人) 回复此楼
5楼2010-05-23 10:37:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 2rKfVQg2jQww 2026-11-06 10/500 2026-11-08 23:01 by l48hFd8Nwae2
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 22:39 by l48hFd8Nwae2
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 21:51 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭