Added optional tolerance to some floating point asserts
Added implicit ode solvers framework

Including a Newton solver.

Currently implemented bdf1, bdf2 and midpoint solvers.