c语言计算定积分

一、问题提出。

 二、设计思路。

假设有这样一个函数，我们要求它的面积（梯形法），我们把它分成n个小梯形
公式是：（上底+下底）*h/2；
a0和a1组成一个梯形先一个一个的求出面积，最后在相加起来。
先根据x=（an-a0）/n求出平均分了以后每一小块的梯形的高，再利用a0+x求出每一小块的位置，根据位置可以利用函数公式y=x+b求出底是多少,然后再与相邻一块的y结合求出面积，最后在全部面积一起相加。

三、代码实现。

#include<stdio.h>
float Fun1(float x);
float Fun2(float x);
float Integral(float(*f)(float), float a, float b);
int main()
{
float y1, y2;
y1 = Integral(Fun1, 0.0, 1.0);
y2 = Integral(Fun2, 0.0, 3.0);
printf("y1=%f\ny2=%f\n", y1, y2);
return 0;
}
float Fun1(float x)
{
return 1 + x*x;
}
float Fun2(float x)
{
return x / (1 + x*x);
}
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++)
{
s += ((*f)(a + i*h));
}
return s*h;
}