00001 00029 #ifndef INTEGRATION_H 00030 #define INTEGRATION_H 00031 00032 #include <limits> 00033 00034 00035 namespace itpp 00036 { 00037 00044 00081 double quad(double(*f)(double), double a, double b, 00082 double tol = std::numeric_limits<double>::epsilon()); 00083 00120 double quadl(double(*f)(double), double a, double b, 00121 double tol = std::numeric_limits<double>::epsilon()); 00122 00124 00125 } // namespace itpp 00126 00127 #endif // #ifndef INTEGRATION_H