C语言.VC6.0软件写用变步长梯形求积法计算定积分f(x)=e^(-x^2)在[1,2]上积分.C语言C语言

代码如下：
#include
#include
#include
#define e 2.7182818
double f(double x)
{
return pow(e,-x*x); // 函数e^(-x^2)
}
double Calc(double a,double b,double esp)//变步长梯形求[a,b]定积分,esp是精度
{
int done(0);
int n=1;
double h,Tn,T2n,k,temp,x;
h=b-a;
Tn=h*(f(a)+f(b))/2.0;
while(!done)
{
temp = 0;
for(k=0;k