Question

15、用梯形法求定积分（数值求解算法）

用梯形法编程求函数f(x)=x2+2x+1的定积分
提示：把积分区间[a，b]等分为n等分（a,b,n的值由程序输入。），可得到若干个小梯形，积分面积就近似为这些小梯形面积之和。

Answer 1

#include<stdio.h>
void main()
{
double fun(double x);
double a,b,i,sum=0;
int n;
printf("请输入a,b,n的值:");
scanf("%lf%lf%d",&a,&b,&n);
for(i=a;i<b;i=i+(b-a)/n)//分割面积求和
{
sum=sum+(b-a)*(fun(i)+fun(i+(b-a)/n))/(2*n);
}
printf("结果为:%lf\n",sum);
}
double fun(double x)
{
double c;
c=x*x+2*x+1;
return c;
}