求解常微分方程组用：

scipy.integrate.odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, 
  ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, 
  hmin=0.0, ixpr=0, mxstep=0, mxhnil=0, mxordn=12, mxords=5, printmessg=0)

• func：梯度函数。第一个参数为y，第二个参数为t0，即计算t0时刻的梯度。其他的参数由args提供
• y0：初始的y
• t：一个时间点序列。
• args：额外提供给func的参数。
• Dfun：func的雅可比矩阵，行优先
• col_deriv：一个布尔值。如果Dfun未给出，则算法自动推导。该参数决定了自动推导的方式
• full_output：如果True，则通过字典返回更多的信息
• printmessg：布尔值。如果为True，则打印收敛信息
• 其他参数用于控制求解的细节

返回值：

• y：一个数组，形状为 (len(t),len(y0)。它给出了每个时刻的y值
• infodict：一个字典，包含额外的信息