[求助]为什么结果总是2倍呢？
//自选步长梯形求积分#include<iostream>
#include<cmath>
using namespace std;
#define N 100//设置最大对分次数
//设置函数
double f(double x)
{
return 1/x;//cos(x)+sin(x)*sin(x);
}
//main()函数
int main()
{
double a,b,h,s=0,xk,T[N];
int i,k,m;
cout<<"请输入积分区间:[a,b]"<<endl;
cin>>a>>b;
for(m=0;m<N;m++)
{
i=pow(2,m);
h=(b-a)/i;
for(k=1;k<i;k++)
{
xk=a+k*h;
s+=f(xk);
}
T[m]=h*(f(a)+2*s+f(b))/2;
if(fabs(T[m]-T[m-1])>=1E-005)
{
cout.precision(16);
cout<<"对分"<<m<<"次时的值为: "<<T[m]<<endl;
}
else
break;
}
return 0;
}
数值计算中的问题，为什么计算结果是实际值的2倍呢？
公式都没有错误，请个位大虾指点迷津，不胜感激~~
