C语言实现任意函数求定积分

求函数f在a和b点之间的定积分

要求具有通用性，可求任意函数的定积分（以下代码中my_fun可改为想求的函数）

输入数据

•积分上限a、积分下限b、分割段数n

    梯形法求定积分

#include<stdio.h>
#include<math.h>
double integrate(double (*fun)(double),double a,double b,int n);

double my_fun(double x){
    return x*x;
}
int main(){
    double a,b;
    int n;
    printf("积分上限：a="); scanf("%lf",&a);
    printf("积分下限：b="); scanf("%lf",&b);
    printf("分割段数：n="); scanf("%d",&n);

    printf("sin函数积分值：%f\n",integrate(sin,a,b,n));
    printf("cos函数积分值：%f\n",integrate(cos,a,b,n));
    printf("x^2函数积分值：%f\n",integrate(my_fun,a,b,n));

    return 0;
}
double integrate(double (*fun)(double),double a,double b,int n){
    int i;
    double h=(b-a)/n;
    double sum=(fun(a)+fun(b))/2;

    for(i = 1;i <= n-1;i++)
        sum += fun(a+i*h);
    sum = h*sum;

    return sum;
}

 PS：在Linux系统下使用gcc编译通常会出现以下错误

wuyy@ubuntu:~$ gcc jifen.c -o jifen
/tmp/ccl4v1pg.o: In function `main':
jifen.c:(.text+0xb7): undefined reference to `sin'
jifen.c:(.text+0xeb): undefined reference to `cos'
collect2: error: ld returned 1 exit status

这通常是缺少数学库而导致

这时只需要在编译时手动加入gcc libm.so库即可，如下：

wuyy@ubuntu:~$ gcc jifen.c -lm -o jifen
wuyy@ubuntu:~$ ./jifen