SciPy

scipy.integrate.DOP853.order

DOP853.order = 8