romberg integration
adaptive quadrature