SciPy

scipy.integrate.RK45.P

RK45.P = array([[ 1. , -2.85358007, 3.07174346, -1.12701757], [ 0. , 0. , 0. , 0. ], [ 0. , 4.02313338, -6.24932157, 2.67542448], [ 0. , -3.73240196, 10.06897059, -5.68552696], [ 0. , 2.55480383, -6.39911238, 3.52193237], [ 0. , -1.37442411, 3.27265775, -1.76728126], [ 0. , 1.38246893, -3.76493786, 2.38246893]])