数值分析中梯形积分公式的C语言程序

/*#include

void main()

{printf("hello");

}*/

#include

#include

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

{T+=f(a+k*h);}

s=(h/2.0)*(f(a)+2*T+f(b));

printf("结果是:%lf\n",s);

}