UC知道c语言用梯形法求积分
来源：百度知道 编辑：UC知道 时间：2024/05/24 12:57:35
#include <math.h>
float integral(float(*fun)(float x),float a,float b,int,n)
{float s,h,y;
int i;
s=(fun(a)+fun(b))/2;
h=(b-a)/n;
for(i=1;i<n;i++)
s=s+fun(a+i*h);
y=s*h;
return y;
高手可以帮我解释下每一行语句的意思吗？
float integral(float(*fun)(float x),float a,float b,int,n)
{float s,h,y;
int i;
s=(fun(a)+fun(b))/2;
h=(b-a)/n;
for(i=1;i<n;i++)
s=s+fun(a+i*h);
y=s*h;
return y;
高手可以帮我解释下每一行语句的意思吗？
a区间起点，b积分区间中点， n是这个区间划分数。每一个划分的长度为h=(b-a)/n。每一个划分利用梯形算法来求积分，即：上底＋下底乘以高再除以2。第一个小梯形面积为 (f(a) + f(a+h))*h/2，第二个为 (f(a+h)+f(a+2h))*h/2，依次类推。将这些加和加起来有：
h*(f(a)/2 + f(a+h) + ... +f(a+(n-1)h) + f(b))/2)。预先将f(a)/2 和f(b)/2加起来，其它数据即f(a+h) + ... f(a+(n-1)n)用循环处理，然后将所有结果乘以h即可。