C语言习题9.6求定积分

最新推荐文章于 2022-05-06 16:35:08 发布

菜馒头

最新推荐文章于 2022-05-06 16:35:08 发布

阅读量1.5k

点赞数 1

文章标签： c语言算法

本文链接：https://blog.csdn.net/vscodenmsl/article/details/121675057

版权

按如下函数原型，采用梯形法编程实现，在积分区间[a,b]内计算函数y1=1+x²在[0,1]的定积分和函数y2= $\frac{x}{1+x^{2}}$ 在[0,3]的定积分。其中，指向函数的指针变量f用于接收被积函数的入口地址

#include<stdio.h>
float Integral(float(*f)(float),float a,float b);
float Y1(float x);
float Y2(float x);
int main(void)
{
    float a,b,ret1,ret2;
    ret1 = Integral(Y1,0,1);
    printf("函数1+x*x在积分区间[0,1]的定积分为%f\n",ret1);
    ret2 = Integral(Y2,0,3);
    printf("函数x/(1+x*x)在积分区间[0,3]的定积分为%f",ret2);
    return 0;
}
float Integral(float(*f)(float),float a,float b)
{
    float s,h;
    int n = 100,i;
    s = ((*f)(a) + (*f)(b)) / 2;
    h = (b - a) / n;
    for(i=1;i<n;i++)//1到n-1内的f(a+i*h)
    {
        s += (*f)(a + i * h);
    }
    return s * h;
}
float Y1(float x)//可以任意增加函数计算其定积分
{
    return 1 + x * x;
}
float Y2(float x)
{
    return x / (1 + x * x);
}

参考输出