SciPy

scipy.integrate.RK23.A

RK23.A = array([[0. , 0. , 0. ], [0.5 , 0. , 0. ], [0. , 0.75, 0. ]])