Scipy 的 integrate 模块的 odeint 函数可以用来以数值积分法求解常微分方程。...import pyplot as plt
# 上篇的向量场绘图函数
def plot_directtion_field(x, y_x, f_xy, x_lim=(-1,1), y_lim=(-1,1...('x')
y = sympy.Function('y')
f = y(x)**2 + x
f_np = sympy.lambdify((y(x), x), f)
x0...轴负方向延伸
xp = np.linspace(x0, x0+2, 100) # 初值处向x轴正方向延伸
yn = integrate.odeint(f_np, y0, xn) # 数值积分法求解常微分方程...,负方向积分
yp = integrate.odeint(f_np, y0, xp) # 数值积分法求解常微分方程,正方向积分
fig, ax = plt.subplots(1,