常微分方程式の解法（python）

常微分非线性方程式用scipy.integrate.odeint(func, y0, t, args=())

func：用def函数定义

y0：初始值（几个和t相关的函数就要设置几个初始值）

t：自变量

args：放定数的（自己设定的常数）

例子如下：

1.分析对t自变量有几个导数。

这整体是二阶导数。可以分为：theta对t的导数=omega（一阶），和omega对t的导数（二阶）。

因此，对t有两个相关导数如下图（也就是说要有两个初始值：theta0和omega0）：

2.编写python程序スクリプト

import numpy as np
import scipy.integrate as sciin
import matplotlib.pyplot as plt

#定义函数f。y是包括和t相关的两个函数theta和omega。t是自变量。parameters是定数。
def f(y, t, parameters):    
    theta, omega = y     #把和t相关的两个函数原始名字写上~123