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

这里有几个问题.首先,你的等式显然是

(3x-1)y” – (3x 2)y’ – (6x-8)y = 0; y(0)= 2,y'(0)= 3

(注意y中术语的符号).对于此等式,您的分析解决方案和y2的定义是正确的.

其次,正如@Warren Weckesser所说,你必须将2个参数作为y传递给g：y [0](y),y [1](y’)并返回它们的导数y’和y”.

第三,你的初始条件是x = 0,但你要整合的x网格从-2开始.从odeint的文档中,这个参数,在他们的调用签名描述中：

odeint(func,y0,t,args =(),…)：

t : array

A sequence of time points for which to solve for y. The initial

value point should be the first element of this sequence.

因此,您必须从0开始积分或提供从-2开始的初始条件.

最后,您的整合范围涵盖了x = 1/3处的奇点. odeint可能在这里度过了不愉快的时光(但显然没有).

这是一种似乎有效的方法：

import numpy as np

import scipy as sp

from scipy.integrate import odeint

import matplotlib.pyplot as plt

def g(y, x):

y0 = y[0]

y1 = y[1]

y2 = ((3*x+2)*y1 + (6*x-8)*y0)/(3*x-1)

return y1, y2

# Ini