C语言求定积分_c语言求定积分的代码

题目要求：
求下面函数的定积分：

思路：
求一个函数的定积分，其实就是求它的面积，如对函数

求积分。
也就是要求出绿色部分的面积，如下：

我们可以通过矩形的方法来无限逼近定积分的求解，如下：

因为被分成n等分，就可以认为每一等分是一个矩形，那么每一矩形的面积为：
每一个矩形面积为：***Sn=f(x)(b-a)/n
总面积为：****S=S1+S2+…+Sn
通过这样的一个思路就可以完成，定积分的求解。
这样这三个定积分的求解代码，如下：

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

void main( )
{
    float f1(float x);
    float f2(float x);
    float f3(float x);
    float s1 = 0;
    float s2 = 0;
    float s3 = 0;
    int n;  //  被分成的份数
    float a,b;  //   积分范围
    int i;      // 循环变量

    printf("请输入极限下限：");
    scanf("%f",&a);
    printf("请输入极限上限：");
    scanf("%f",&b);
    printf("请输入被分成的份数：");
    scanf("%d",&n);
    for(i=1; i<=n;i++)
    {
        s1 = s1 + f1(a+((b-a)/n)*i)*((b-a)/n);
        s2 = s2 + f2(a+((b-a)/n)*i)*((b-a)/n);
        s3 = s3 + f3(a+((b-a)/n)*i)*((b-a)/n);
    }
    printf("积分的结果为：%0.4f\n",s1);
    printf("积分的结果为：%0.4f\n",s2);
    printf("积分的结果为：%0.4f\n",s3);
}

/* y=x */
float f1(float x)
{
    float y=0;
    y = x;

    return y;

}

/* y = x*x */
float f2(float x)
{
    float y=0;
    y = x*x;
    return y;
}

/* y = sin(x) */
float f3(float x)
{
    float y=0;
    y = (float)sin(x);
    return y;
}

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59

结果为：

这样就把定积分求完了，当然还有其他的方法。