仅打印最终的odeint输出

Question

对不起，这可能看起来是一个非常愚蠢的问题，但我需要问一下，在解决odeint中的耦合微分方程时是否可以只打印最终输出值？实际上，我正在尝试解决随机生成的时间间隔的两个耦合微分方程，并且只获得每个区间打印的最终输出。

Answer 1

y值数组的最后一个元素是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）相当接近。