数值分析中梯形积分公式的C语言程序
/*#include<stdio.h>
void main()
{printf("hello");
}*/
#include<stdio.h>
#include<math.h>
double f(double x)
{return sin(x) ;//此处可以更改积分函数
}
void main()
{
double h,T=0.0,s;
int k,n=100;double a=0,b=3.141593;//变量改变行,积分上下边界a,b可改，n是等分数 h=(b-a)/n;
for(k=1;k<n;k++)
{T+=f(a+k*h);}
s=(h/2.0)*(f(a)+2*T+f(b));
printf("结果是:%lf\n",s);
}