// Numerical integration via the Extended Trapezoidal Formula (ETF)
#include <cmath>
#include <iostream>
using namespace std;

double f(double x) {
  double y = 1.76 + 3.21*x*x;
  return 894.*x / (y*y*y);
}

int main() {

  const int n=100;
  const double a=0.0, b=1.61;
  const double h = (b-a)/n;

  double etf = (f(a)+f(b))/2;

  for (int i=1; i<n; i++)
    etf += f(a+i*h);

  etf *= h;

  cout.precision(9);
  cout << "The integral = " << etf << endl;

}