Crate lsode

Functions

dlsode_⚠	Call DLSODE subroutine from ODEPACK
fake_jacobian	A dummy function to pass to dlsode_ in case the user does not want to specify a Jacobian.
solve_ode	Solves system of ODEs for times in t_dense. First time in t_dense has to be the initial time.