redid m_rk_prim to m_ode_prim -- added different numerical solvers