代码如下:#include #include #include #define e 2.7182818double 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;
w... |