Bug List

Global math_simpson (double(*f)(double), double a, double b, double epsilon)

Probably does not work with SIMD