1

Pythonで対応する初期条件で次の頌歌をセットアップする方法は?

x'(t) =x(t) - y(t) - e^t

y'(t) =x(t) + y(t) + 2e^t

とx(0)= -1とy(0)= -1と0 <= t <= 4

以下は私がこれまでに持っているものです：

def f(u, t):
    x, y = u
    return [x+y-e**t, x+y+2*e**t]

x0, y0 = [-1.0,-1.0]
t = numpy.linspace( 0,4,50 )
4

1 に答える 1

1

odeintでそれらを解決しようとしていると思います。まず、スクリプトでこのプレリュードを使用すると仮定しています。

import numpy as np
from scipy.integrate import odeint

あなたの方程式は次のとおりです。

def equation(X, t):
    x, y = X
    return [ x+y-np.exp(t), x+y+2*np.exp(t) ]

そして、あなたはそれらを解決することができます

init = [ -1.0, -1.0 ]
t = np.linpsace(0, 4, 50)
X = odeint(equation, init, t)

x(t) と y(t) を抽出できます

x = X[:, 0]
y = X[:, 1]
于 2013-05-25T08:17:43.927 に答える