SciPy

scipy.integrate.RK23.order

RK23.order = 3