使用SciPy求解常微分方程

Question

我有以下内容 ordinary differential equation 和数字参数Sigma = 0.4，x（0）= 4和dx（0）/ dt = 0
我的任务是使用ode函数得到微分方程的Cauchy问题解（初值问题解） 有人能帮我吗？我甚至不知道如何以正确的方式为SciPy编写方程式，特别是数字参数。
P.S.很抱歉没有发布图片，我刚刚注册了。

Answer 1

就像沃伦所说，scipy.integrate.odeint是解决这个问题的'SciPy'方法。

但是在你把问题带到SciPy（或者你最终使用的任何解算器）之前，你需要使用类似的东西将你的二阶ODE转换为一阶ODE：http://tutorial.math.lamar.edu/Classes/DE/SystemsDE.aspx

为了让事情进入SciPy，你需要让你的等式看起来像：

y'= f（y）

但是现在你的等式写得像：

y''= f（y，y'）

解决方案是向系统添加更多变量，但链接将更彻底地解释它。