#include
double integ(double x1,double xn,double f(double x),int n)
{
double sum,h=(xn-x1)/(n-1);
int i;
for(i=2,sum=0;i sum+=f(x1+h*(i-1)); return h/2*(f(x1)+sum*2+f(xn)); } double f(double x) { return x*x+x+1; } int main() { printf("%lf\n",integ(0,1,f,1000)); return 0; }