c语言计算定积分

发布时间 2023-04-19 23:43:01作者: liuxuechao

一、问题提出。

 二、设计思路。

假设有这样一个函数,我们要求它的面积(梯形法),我们把它分成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;
}