As the title suggests, I'm looking for a very accurate algorithm to find the area under an oscillatory function.

For instance, I would like to integrate the following graphenter image description here

I'm currently using simpson's method to integrate this, but I was wondering if there are more accurate methods.

  • 3
    $\begingroup$ Your test function looks smooth (infinitely often differentiable), to me. Did you mean to write oscillatory instead? $\endgroup$ Jun 4, 2016 at 19:11
  • 1
    $\begingroup$ In what form is the function given? Do you have an analytic or partially analytic form, or is black box or defined by interpolation? $\endgroup$ Jun 6, 2016 at 4:15

1 Answer 1


You can look up various quadratures. One method that should fair better is Gauss Quadrature.

I would also recommend looking into any adaptive quadrature schemes. There are many of them out there, so searching up algorithms for adaptive integration would help you.

This wikipedia link should give you some insight into adaptive integration.

  • $\begingroup$ Thanks for your reply, I'm going to look into it and come back to this answer $\endgroup$
    – Hunter
    Jun 4, 2016 at 20:28

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.