Euler（欧拉）方法 | 显示Euler格式、隐式Euler格式、两步Euler格式、改进的Euler格式

微分方程初值问题的数值解法，就是将一个连续的微分方程初值问题转化为一个离散的差分方程初值问题，然后通过解差分方程而获得其数值解。即：对自变量的一系列离散节点$x_k(k=0,1,2,\cdots ,n\cdots )$，计算出去精确解$y(x_k)$的近似值$y_k(k=1,2,\cdots ,n,\cdots )$，（其中$y_0$是准确值）从而得到函数$y=y(x)$的近似数值解$y_k(k=1,2,\cdots ,n,\cdots )$

关于微分方程解的存在性与惟一性，有下面的定理：

定理1：对于一阶常微分方程初值问题：
$$\{\begin{array}{l}{y}^{\prime }(x)=f(x,y)x\in [a,b]\\ y(x_0)=y_0\end{array}\text{\hspace{1em}(1)}$$
如果$f(x,y)$对x连续，且关于y满足Lipschitz（李普希兹）条件，即存在常数$L>0$，使得
$$|f(x,y_1)-f(x,y_2)|\le L|y_1-y_2|$$
对所有$x\in [a,b]$以及任何实数$y_1,y_2$均成立，则初值问题在$[a,b]$上有惟一解。

Euler（欧拉）方法

将微分方程（1）转化成差分方程，有几何方法、数值微分、数值积分和Taylor（泰勒）展开法等数值方法，这些方法将在推导微分方程的各种数值解法时得到具体运用。

1. 显示Euler格式

从导数的定义可知，当步长$h=x_{n+1}-x_n\to 0$时，有：
$$y^{\prime }(x)=li{m}_{\Delta x\to 0}\frac{\Delta y}{\Delta x}=li{m}_{\Delta x\to 0}\frac{f(x+\Delta x)-f(x)}{\Delta x}\cong \frac{y(x_{n+1})-y(x_n)}{x_{n+1}-x_n}\text{\hspace{1em}(2)}$$
（2）式右边是导数的近似差商表示。因此，如果向前差商、向后差商和中心差商分别代替微分方程（1）中的导数项$y^{\prime }(x)$，则可相应地导出Euler方法的各种格式。

1. 显示Euler格式

若用在$x_n$点的向前差商$\frac{y(x_{n+1})-y(x_n)}{x_{n+1}-x_n}$代替微分方程（1）的导数项$y^{\prime }(x_n)$，当$h=x_{n+1}-x_n$很小时，则有：
$$\frac{y(x_{n+1})-y(x_n)}{h}=f(x_n,y(x_n))$$
若将$y(x_n)$的近似值$y_n$代入上式右端，并将计算结果记为$y_{n+1}$，则有：
$$y_{n+1}=y_n+hf(x_n,y_n)(n=0,1,2,\cdots )\text{\hspace{1em}(3)}$$
这样就得到了$y(x_{n+1})$的近似值$y_{n+1}$的递推格式，即所谓的显示Euler格式。

若给出初始条件$y(x_0)=y_0$，则按此（3）式可逐步算出$y_1,y_2,\cdots ,$即：
$$\begin{gathered} y_1=y_0+h{y}_0^{\prime }=y_0+hf(x_0,y_0) \\ {y}_2=y_1+hf( \end{gathered}$$