Question

西工大poj梯形积分那题结果正确为啥是PE？查了几遍了 没错啊 求AC了的大神帮忙

Answer 1

用c和梯形法会无限PE，但用c++和牛顿-莱布尼兹公式写就AC了，附代码：

#include <iostream>
#include <math.h>
using namespace std;

int main()
{
    double a,b,sum;
    cin>>a>>b;
    sum=atan(b)-atan(a);
    cout << sum << endl;
    return 0;
}

Answer 2

#include<stdio.h>
#include<math.h>
double integral(double a,double b,double(*f)(double x))
{
int n=1000,i;
double h,x,s=0.0;
h=(b-a)/n;
for(i=1;i<=n;i++)
{
x=a+(i-1)*h;
s=s+(f(x)+f(x+h))*h/2;
}
return s;
}
double f(double x)
{
return 1/(1+x*x);
}
int main()
{
double a,b,t;
scanf("%lf %lf",&a,&b);
t=integral(a,b,f);
printf("%lf\n",t);
return 0;
}

追问

仍是PE