#include <stdio.h>
double func(double x) {
return x * x * x + 2 * x * x - 1;
}
double integral(int n, double left, double right) {
double h = (right-left)/n;
double sum = 0;
int cnt = 0;
for (double xi = left; xi < right; xi = xi + h) {
sum += func(xi);
++cnt;
}
printf("count: %d\n", cnt);
return h * sum;
}
int main()
{
double res10 = integral(10,1,4);
printf("%lf\n",res10);
double res100 = integral(100,1,4);
printf("%lf\n",res100);
return 0;
}