SciPy

scipy.integrate.DOP853.D

DOP853.D = array([[-8.42893828e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 5.66714954e-01, -3.06894995e+00, 2.38466766e+00, 2.11703458e+00, -8.71391584e-01, 2.24043743e+00, 6.31578779e-01, -8.89903365e-02, 1.81485055e+01, -9.19463239e+00, -4.43603639e+00], [ 1.04275086e+01, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 2.42283492e+02, 1.65200452e+02, -3.74546755e+02, -2.21136669e+01, 7.73343267e+00, -3.06740847e+01, -9.33213053e+00, 1.56972381e+01, -3.11394032e+01, -9.35292436e+00, 3.58168415e+01], [ 1.99850532e+01, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, -3.87037309e+02, -1.89178138e+02, 5.27808159e+02, -1.15739025e+01, 6.88123269e+00, -1.00060510e+00, 7.77713780e-01, -2.77820575e+00, -6.01966952e+01, 8.43204055e+01, 1.19922911e+01], [-2.56939335e+01, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, -1.54189749e+02, -2.31529379e+02, 3.57639118e+02, 9.34053242e+01, -3.74583231e+01, 1.04099650e+02, 2.98402934e+01, -4.35334566e+01, 9.63245540e+01, -3.91772617e+01, -1.49726836e+02]])