24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

muqiuyan

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2962
散金: 355
红花: 1
帖子: 162
在线: 64小时
虫号: 829827
注册: 2009-08-18
性别: MM
专业: 能源化工

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

怎么区分啊

回复此楼

» 猜你喜欢

Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有6人回复
2025冷门绝学什么时候出结果已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
论文投稿，期刊推荐已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复

1楼2010-05-21 14:44:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tjyl

金虫 (正式写手)

程序强帖: 2
应助: 0 (幼儿园)
金币: 3218.1
红花: 2
帖子: 576
在线: 156.6小时
虫号: 765184
注册: 2009-05-07
专业: 生物无机化学
管辖: 程序语言

★ ★ ★
resonant(金币+3):专家辛苦了；—） 2010-06-04 19:06:53

我也不知道是我没有弄明白你的意图还是你没有理解我建议用函数指针的用意。
直接写个很简单的范例，不知道对你有启发没有
#include
double plus(double (*f) (double), double a, double b);
double F1(double);
double F2(double);
int main()
{
double x1, x2;
double y1, y2;
x1 = 0.5;
x2 = 1.5;
y1 = plus(F1, x1, x2);
y2 = plus(F2, x1, x2);
printf("plus of F1=%0.6f\n", y1);
printf("plus of F2=%0.6f\n", y2);
return 0;
}

double plus(double (*f) (double), double x1, double x2)
{
return (f(x1) + f(x2));
}

double F1(double x)
{
return x * x;
}

double F2(double x)
{
return x;
}
只要你的2个积分函数的参数个数是一样的就行了
如果参数个数不一样的话建议用C++来写。C的可变参数估计你一时半会搞不清楚。

引用回帖:

Originally posted by muqiuyan at 2010-06-03 10:56:34:

万分感谢，因为要编程，才开始学，很多都不懂。
麻烦在帮忙看看这个，我的程序里有两个积分，每个的积分变量还不同，怎么在程序中体现区分呢。

double fsimp(double (*f)(double ),double a, d ...

[ Last edited by tjyl on 2010-6-3 at 12:39 ]