7

与えられた初期条件でこの微分方程式を解きたい:

(3x-1)y''-(3x+2)y'+(6x-8)y=0, y(0)=2, y'(0)=3

ansはあるべきです

y=2*exp(2*x)-x*exp(-x)

ここに私のコードがあります:

def g(y,x):
    y0 = y[0]
    y1 = y[1]
    y2 = (6*x-8)*y0/(3*x-1)+(3*x+2)*y1/(3*x-1)
    return [y1,y2]

init = [2.0, 3.0]
x=np.linspace(-2,2,100)
sol=spi.odeint(g,init,x)
plt.plot(x,sol[:,0])
plt.show()

しかし、私が得るものは答えとは異なります。私は何を間違えましたか?

4

1 に答える 1