Truncation Error#

The accuracy of the algebraic approximation to \(du/dt\) (3) can be determined with the help of the Taylor series expansion of \(u(t)\):

(4)#\[ u(t+\Delta t)=u(t)+\Delta t\frac{du}{dt} + \displaystyle\sum_{p=2}^{\infty} \frac{\Delta t^p}{p!}\frac{d^pu}{dt^p} \]

Rearranging for \(du/dt\), we obtain, using the discretized \(t^n\):

(5)#\[ \frac{du}{dt} = \frac{u^{n+1}-u^n}{\Delta t} - \displaystyle\sum_{p=2}^{\infty} \frac{\Delta t^{p-1}}{p!}\frac{d^pu}{dt^p}. \]

The difference between the exact derivative (1) and our algebraic approximation (3) is the second term of the right hand side of (5). Expanding this term up to \(p=3\), we see its general form:

(6)#\[ \frac{\Delta t}{2!}\frac{d^2u}{dt^2}+\frac{\Delta t^2}{3!}\frac{d^3u}{dt^3}+O(\Delta t^3), \]

where \(O(\Delta t^3)\) represents terms proportional to \(\Delta t^3\) and to higher powers of \(\Delta t\).

Expression (6) is the truncation error, which is the error we incurr when we approximate the exact derivative (1) with the algebraic expression (3).

The algebraic approximation (3) can thus be written as:

(7)#\[ \frac{du}{dt} = \frac{u^{n+1}-u^n}{\Delta t} + O(\Delta t), \]

where \(O(\Delta t)\) represents the truncation error of the approximation, whose first term, the leading order term is proportional to \(\Delta t\). Since the leading order term is proportional to the first power of \(\Delta t\), we say that this approximation is a first order approximation to the first derivative.