仅打印最终的odeint输出

y[-1] 。例如:

import numpy as np
import scipy.integrate as si
def F(y, t):
    return [y[1], y[0]]
t = np.arange(0, 1, 0.001)
y = si.odeint(F, [1, 0], t)
print(y[-1])

收益 [ 1.54190626 1.17365875] 该系统的精确解为y(t)=[cosh(t),sinh(t)];所示数字相当接近于cosh(1)和sinh(1)。