About half a year ago, I learned of Gaussian numerical integration. It is a simple method to accurately evaluate a bounded integral. The idea is that there is a polynomial that approximates the function being integrated. The integral can be approximated by the sum of weights times the polynomial evaluated at different points. The more points used the more precise the answer is.
Link to a more comprehensive explanation.