1.一阶微分方程

1.数值方法

1.1 欧拉方法

  一阶微分方程$$dy/dx=f(x,y)\text{\hspace{1em}(1)}$$

以及初值$y\left(x_0\right)=y_0$给定了函数$y(x)$在初值点$(x_0,y_0)$的切线斜率$f(x_0,y_0)$。以该初值为起点，以切线为方向，按照某步长$\Delta x=x_1-x_0$变化到新点$（x_1,y_1）$。这就是欧拉方法，用公式表述为式(1):
$$y_1=y_0+\Delta xf\left(x_0,y_0\right)\text{\hspace{1em}(2)}$$

1.2 龙格库塔方法

  欧拉方法也别称为一阶龙格库塔法。
  二阶龙格库塔法的通式为：
$$k_1=\Delta xf\left(x_n,y_n\right),k_2=\Delta xf\left(x_n+\alpha \Delta x,y_n+\beta k_1\right),y_{n+1}=y_n+a{k}_1+b{k}_2\text{\hspace{1em}(3)}$$

并满足

$$a+b=1，\alpha b=\beta b=1/2\text{\hspace{1em}(4)}$$

2. 可分离变量的一阶方程

  可分离变量的一阶方程可写为
$$g(y)\frac{dy}{dx}=f(x),y\left(x_0\right)=y_0\text{\hspace{1em}(5)}$$

两边从$x_0$到$x_1$进行积分，$${\int }_{x_0}^{x}g(y(x))y^{\prime }(x)dx={\int }_{x_0}^{x}f(x)dx\text{\hspace{1em}(6)}$$

进行变量替换$u=y(x),du=y^{\prime }(x)dx$，并改变积分上下限为$y\left(x_0\right)=y_0$和$y(x)=y$。因此有$${\int }_{y_0}^{y}g(u)du={\int }_{x_0}^{x}f(x)dx\text{\hspace{1em}(7)}$$

3. 一阶线性微分方程

  形如$$\frac{dy}{dx}+p(x)y=g(x),y\left(x_0\right)=y_0\text{\hspace{1em}(8)}$$

的方程。微分方程的求解借助于名为积分因子的函数$\mu (x)$的求解。积分因子是这样一个函数，它乘以式(8)后得到的式子为
$$\mu (x)\left[\frac{dy}{dx}+p(x)y\right]=\mu (x)g(x)\text{\hspace{1em}(8)}$$

而该式左边等于一个积分
$$\mu (x)\left[\frac{dy}{dx}+p(x)y\right]=\frac{d}{dx}[\mu (x)y]\text{\hspace{1em}(9)}$$

由$\frac{d}{dx}[\mu (x)y]=\mu (x)g(x)$，利用初值$y\left(x_0\right)=y_0$并令$\mu \left(x_0\right)=1$，可得$\mu (x)y-y_0={\int }_{x_0}^x\mu (x)g(x)dx$，求得：
$$y(x)=\frac{1}{\mu (x)}\left(y_0+{\int }_{x_0}^x\mu (x)g(x)dx\right)\text{\hspace{1em}(10)}$$

为求积分因子，展开式(9)，得
$$\frac{d\mu }{dx}=p(x)\mu ,\mu \left(x_0\right)=1\text{\hspace{1em}(11)}$$

这是可分离变量微分方程，得
$$\mu (x)=\exp \left({\int }_{x_0}^{x}p(x)dx\right)\text{\hspace{1em}(12)}$$

注解：积分因子是构造出来的。