SciPy

scipy.integrate.RK45.order

RK45.order = 4