Make some updates to LMM ODE solving (fix simpsons rule), and add an
example of gaussian quadrature to calculate an integral