SciPy

scipy.integrate.RK23.P

RK23.P = array([[ 1. , -1.33333333, 0.55555556], [ 0. , 1. , -0.66666667], [ 0. , 1.33333333, -0.88888889], [ 0. , -1. , 1. ]])