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

C语言实现任意函数求定积分_第1张图片

求函数fab点之间的定积分

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

输入数据

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

    梯形法求定积分

#include
#include
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

你可能感兴趣的:(课后练习,c语言,算法,c++)