We should test different numerical methods for integration. See: https://en.wikipedia.org/wiki/Simpson%27s_rule