使用scipy.optimize模块的root和fsolve函数进行数值求解线性及非线性方程，下面直接贴上代码，代码很简单 from scipy.integrate import odeint import numpy as np import matplotlib.pyplot as plt from scipy.optimize import root,fsolve #plt.rc('text', usetex=True) #使用latex ## 使用scipy.optimize模块的ro