Adaptive Simpson method

The code computes an integral of a function over an interval using the adaptive simpson rule.

Input: function, endpoints of an interval, tolerance
Output: integral of the function over the interval

The code contains an example with the function. For more precise computation, split an interval into subintervals, implement the function on each of them, and sum the results. One way of the splitting is taking uniformly split subintervals.