integral_calc.hpp 427 B

1234567891011121314151617
  1. #include <omp.h>
  2. double integral_calc(const double param, Interval_t interval){
  3. double res = 0;
  4. double increase = (interval.sup - interval.inf) / (double) interval.increase_count;
  5. double x = interval.inf + 0.5 * increase;
  6. for(int i = 0; i < interval.increase_count; ++i){
  7. //std::cout << increase << " | " << x << "\n";
  8. res += integrand(x, param);
  9. x += increase;
  10. }
  11. res *= increase;
  12. return res;
  13. }