最終的なodeint出力のみを出力する

Question

これは結構な疑問に思えるかもしれませんが、最終的な出力値だけを出力することが可能かどうかを尋ねる必要があります。実際には、無作為に生成された時間間隔の2つの結合微分方程式を解いて、各区間ごとに最終出力のみを出力しようとしています。

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）にかなり近い数値です。