Question

c语言 求定积分的通用函数

请问这里哪里错了？？？
#include<stdio.h>
#include<math.h>
void main()
{
float integral(float a,float b,float(*p)(float x) );
float f1(float x);
float f2(float x);
float f3(float x);
float f4(float x);
float f5(float x);

float x;
float a=5,b=7;
float result;
printf("请输入x的值:\n");
scanf("%f",&x);

result=integral(a,b,f1);
result=integral(a,b,f2);
result=integral(a,b,f3);
result=integral(a,b,f4);
result=integral(a,b,f5);
}
float f1(float x)
{
float s;
s=1+x;
return s;
}
float f2(float x)
{
float s;
s=2*x+3;
return s;
}
float f3(float x)
{
float s;
s=(float)exp(x)+1; //e的x次方+1
return s;
}
float f4(float x)
{
float s;
s=(1+x)*(1+x);
return s;
}
float f5(float x)
{
float s;
s=x*x*x;
return s;
}
float integral(float a,float b,float(*p)(float x) )
{
float i,result,x;
for(i=a;i<=b;i++)
{
result+=(*p)(x);
}
return result;
printf("结果是: %3.3f",result);
}

Answer 1

对于一重定积分来说其求解可以使用梯形法进行求解，计算公式如下所示：

其中，f(x)为被积函数，为横坐标的两点间的间隔，越小，则计算出的结果越精确。

对于求解此类问题可以使用C语言中的回调函数编写通用的计算函数，代码如下：

#include <stdio.h>
#include <stdlib.h>
#include<math.h>

//功能：返回f(x)在积分区间[a,b]的值
//参数：FunCallBack 指向用于计算f(x)的函数
//      a  积分区间的起始值
//      b  积分区间的结束值
//      dx 横坐标的间隔数，越小计算结果越准确
double Calculate(double (*FunCallBack)(double x),
    double a,double b,double dx)
{
    double doui;
    double total = 0;        //保存最后的计算结果

    for (doui = a; doui <= b; doui += dx)
    {
        total += FunCallBack(doui)*dx;
    }
    return total;
}

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

double f(double x)
{
    return x;
}

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

int main()
{
    double total;
    total = (Calculate(f, 2, 3, 0.000001));
    printf("total = %lf\n", total);

    total = (Calculate(f2, 2, 3, 0.000001));
    printf("total = %lf\n", total);

    total = (Calculate(f3, 2, 3, 0.000001));
    printf("total = %lf\n", total);

    return 0 ;
}
其中，函数f，f2，f3为自行编写的关于x的被积函数。
运行结果：
        total = 2.500000
        total = 6.333331
        total = 16.249991

Answer 2

float integral(float a,float b,float(*p)(float x) )
{
float i,result,x;
for(i=a;i<=b;i++)
{
result+=(*p)(x);
}
return result;
printf("结果是: %3.3f",result);
}

这里错了，return result;
printf("结果是: %3.3f",result);
换一个位置，return 就跳出函数了，不会打印出来

Answer 3

哥们儿？你这代码能运行？我复制你的怎么不行呀？

Answer 4

我也没有试，就看了一下，上面的函数调用时，参数的定义在下面，没定义先用，不服合语法吧

追问

没呢。。不过我在主函数那里声明了喔。。

追答

我的意思是，你的声明在下面，可上面调用的时候，参数还没声明啊，我是这个意思

追问

噢。。原来这样啊。。

追答

你改好了，哪里错了呢

追问

最后 return result; 那里错了，我先跳了出来所以结果输不出来，谢谢你帮忙啦。。。

追答

啊，这到没事，不过，那种声明在下面的语法居然没错？我还真没那样学过。恩，C中，return是结束语。这到是，当时只顾看上面那些，没看下面，呵呵。我也受教了。