欧拉方法

考虑计算这样的一个未知曲线的形状：它具有给定的起点并且满足一个给定的微分方程。 这里，所谓“微分方程”可以看作能够通过曲线上任意点的位置而计算出这一点的切线斜率的公式。

思路是，一开始只知道曲线的起点（假设为${\displaystyle A_{0}}$ ），曲线其他部分是未知的，不过通过微分方程，${\displaystyle A_{0}}$ 的斜率可以被计算出来，也就得到了切线。

顺着切线向前走一小步到点${\displaystyle A_{1}}$ 。如果假设${\displaystyle A_{1}}$ 是曲线上的一点（实际上通常不是），那么同样的道理就可以确定下一条切线，依此类推。在经过几步之后，一条折线（英语：Polygonal chain）${\displaystyle A_0{A}_1{A}_2{A}_3\dots }${\displaystyle A_{0}A_{1}A_{2}A_{3}\dots }  就被计算出来了。一般情况下，这条折线与原先的未知曲线偏离不远，并且任意小的误差都可以通过减少步长来得到。

这是一个以${\displaystyle \mathbf {z} (t)}$ 为变量的一阶系统，因此可以用欧拉法求解，也可以使用其他的一阶数值方法。[1]

设微分方程为 ${\displaystyle y'=y}$  ，初始值为 ${\displaystyle y(0)=1}$ ，试用欧拉方法求 ${\displaystyle y_{3}}$ 的近似值，步长为 ${\displaystyle h=1}$ 。

欧拉法为：

${\displaystyle y_{n+1}=y_{n}+hf(t_{n},y_{n}).}$ 

首先求${\displaystyle f(t_{0},y_{0})}$ （当${\displaystyle n=0}$ ），${\displaystyle f}$ 的定义为${\displaystyle f(t,y)=y}$ ，因此有

${\displaystyle f(t_{0},y_{0})=f(0,1)=1.}$ 

透过以上步骤，求得解曲线在点${\displaystyle (0,1)}$ 的切线斜率。回顾直线斜率的定义：${\displaystyle y}$ 变化量和${\displaystyle t}$ 变化量的比值，亦记作${\displaystyle \Delta y/\Delta t}$ 。

接着是

${\displaystyle h\cdot f(y_{0})=1\cdot 1=1.}$ 

${\displaystyle y_{0}+hf(y_{0})=y_{1}=1+1\cdot 1=2.}$ 

重复以上步骤求出${\displaystyle y_{2}}$  和 ${\displaystyle y_{3}}$ 的值。

${\displaystyle y_{2}=y_{1}+hf(y_{1})=2+1\cdot 2=4}$ 

${\displaystyle y_{3}=y_{2}+hf(y_{2})=4+1\cdot 4=8}$ 

由于欧拉法属于递归算法，把运算整理成表格也许有助于避免计算错误。

欧拉法的局部截尾误差（Local truncation error, LTE）是指在实施一次欧拉法所产生的误差，是指经过一步的数值解${\displaystyle y_{1}}$ 与在${\displaystyle t_{1}=t_{0}+h}$ 时精确解的误差。数值解${\displaystyle y_{1}}$ 由以下给出：

${\displaystyle y_{1}=y_{0}+hf(t_{0},y_{0}).\quad }$ 

对于精确解，使用泰勒级数展开给出：

${\displaystyle y(t_{0}+h)=y(t_{0})+hy'(t_{0})+{\frac {1}{2}}h^{2}y''(t_{0})+O(h^{3}).}$ 

欧拉法的局部截尾误差为：

${\displaystyle \mathrm {LTE} =y(t_{0}+h)-y_{1}={\frac {1}{2}}h^{2}y''(t_{0})+O(h^{3}).}$ 

当${\displaystyle y}$ 拥有三阶有界导数时，这个结果是成立的。^[2]

结果显示：当步长${\displaystyle h}$ 很小时，局部截尾误差近似与 ${\displaystyle h^{2}}$  成比例。也就是说，欧拉法没有其他的高阶方法如龙格－库塔法和线性多步法（英语：Linear multistep method）精确，这些方法的局部截尾误差与${\displaystyle h^{p}}$ （p>2）成比例。

全局截尾误差（Global truncation error, GTE）是指在一个固定时间${\displaystyle t}$ 时的误差，但是很多步之后该方法需要以从初始时间到达该时间来计算。全局截尾误差可以看做是一个每一步的局部截尾误差的累积效应。^[3] 经过的步骤数为${\displaystyle (t-t_{0})/h}$ ，而每步的误差则正比于${\displaystyle h^{2}}$ 。因此，可以预期全局截尾误差是正比于${\displaystyle h}$ 的。^[4]

这个直观的推测可以被严谨地证明。如果解${\displaystyle y}$ 存在二阶有界导数，并且${\displaystyle f}$ 关于${\displaystyle y}$ 是利普希茨连续的，那么全局截尾误差是有界的：

${\displaystyle |{\text{GTE}}|\leq {\frac {hM}{2L}}(e^{L(t-t_{0})}-1)\qquad \qquad }$ 

其中 ${\displaystyle M}$  是在给定区间内${\displaystyle y}$ 的二阶导数的上界，${\displaystyle L}$  是${\displaystyle f}$ 的利普希茨常数。^[5]

这种精确的形式其实是没有什么意义的，通常情况下这个上界都会严重高估了欧拉法所造成的实际误差。^[6]重要的是，这显示了全局截尾误差是近似正比于${\displaystyle h}$ 的，所以欧拉法被称为是一阶的。^[7]

• 克兰克－尼科尔森方法
• 梯度下降法，也是用有限步进行，计算函数的最小值。
• 龙格－库塔法列表（英语：List of Runge-Kutta methods）
• 线性多步法（英语：Linear multistep method）
• 数值积分（计算定积分）
• 常微分方程数值方法（英语：Numerical methods for ordinary differential equations）

• Atkinson, Kendall A., An Introduction to Numerical Analysis 2nd, New York: John Wiley & Sons, 1989, ISBN 978-0-471-50023-0 .
• Ascher, Uri M.; Petzold, Linda R., Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations, Philadelphia: Society for Industrial and Applied Mathematics, 1998, ISBN 978-0-89871-412-8 .
• Butcher, John C., Numerical Methods for Ordinary Differential Equations, New York: John Wiley & Sons, 2003, ISBN 978-0-471-96758-3 .
• Hairer, Ernst; Nørsett, Syvert Paul; Wanner, Gerhard, Solving ordinary differential equations I: Nonstiff problems, Berlin, New York: Springer-Verlag, 1993, ISBN 978-3-540-56670-0 .
• Iserles, Arieh, A First Course in the Numerical Analysis of Differential Equations, Cambridge University Press, 1996, ISBN 978-0-521-55655-2 .
• Lakoba, Taras I., Simple Euler method and its modifications (PDF) (Lecture notes for MATH334, University of Vermont), 2012 [2016-01-02], （原始内容存档 (PDF)于2012-07-12） .