python如何求解微分方程,如何使用Python内置函数odeint求解微分方程？

I want to solve this differential equations with the given initial conditions:

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

the ans should be

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

here is my code:

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()

but what I get is different from the answer.

what have I done wrong?

解决方案

There are several things wrong here. Firstly, your equation is apparently

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

(note the sign of the term in y). For this equation, your analytical