私はこの方程式を解くために必要パイソン

Question

にodeint用いた二次ODEを解く：

Y '' + Y」 - B Y、A及びBは、の関数である Y = 0

それは意味がありませんプロットを返す

import numpy as np 
import matplotlib.pyplot as pyplot 
from scipy.integrate import odeint 


x0 = 0.,0.1 # initial conditions 
oe = 0.0 # cosmological constant density parameter 
om = 1. # matter density parameter 
h = 2.3E-18 # Hubble constant 
w = -1. 

a = np.linspace(0.1,1,10) # scale factor 


def H(a): # Hubble rate equation 
    return h*np.sqrt((om/(a**3.))+(oe/(a**(3.*(1.+w))))) 

def A(a): # differential equation term 
    return (-3./(2.*a)) 

def B(a): # differential equation term 
    return (3.*om*(h**2.))/((2.*(a**5.))*(H(a)**2.)) 



def system(X,a): # differential equation system 
    X0 = X[0] 
    X1 = X[1] 
    X2 = -A(a)*X1 + B(a)*X0 

    return X1,X2 



x = odeint(system,x0,a) 



pyplot.semilogx(a,x, linestyle='-', c="k", linewidth="2") 

：同じ変数「」

は、私は、次のコードを試してみました。私はa = 1で "x"の最大値が1である1つのプロットを取得する必要があります。私が手

プロット：しかし、私は次のプロットを取得

を、期待される結果が次のグラフの連続線のようなものです：

任意の提案？

Answer 1

機能システム（X、a）でB（a）* X0の前に符号を変更してみてください。それはあなたの最初の方程式に従って負でなければなりません。