Thanks @oscardssmith !
I see that FLINT and RKC does not explicitly require Jacobian, but they are not particularly suitable for stiff problem.
Uhm in your experience, how do you switch between non-stiff and stiff ode solver?
Because in real case, even if we use the same model of ODE, if some parameters of the ODE changes, the problem may become a stiff problem. Therefore being able to efficiently switch between non-stiff and stiff can be very useful.