仅打印最终的 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])

returns[ 1.54190626 1.17365875]。该系统的精确解是 y(t) = [cosh(t), sinh(t)]；显示的数字相当接近 cosh(1) 和 sinh(1)。