In numerical analysis and scientific computing, the backward Euler method (or implicit Euler method) is one of the most basic numerical methods for the solution of ordinary differential equations. It is similar to the (standard) Euler method, but differs in that it is an implicit method. The backward Euler method has error of order one in time.
dy | |
dt |
=f(t,y)
y(t0)=y0.
f
t0
y0
y
t
y0,y1,y2,\ldots
yk
y(t0+kh)
h
The backward Euler method computes the approximations using
yk+1=yk+hf(tk+1,yk+1).
f(tk,yk)
f(tk+1,yk+1)
The backward Euler method is an implicit method: the new approximation
yk+1
yk+1
[0] | |
y | |
k+1 |
=yk,
[i+1] | |
y | |
k+1 |
=yk+hf(tk+1,
[i] | |
y | |
k+1 |
).
yk+1
Alternatively, one can use (some modification of) the Newton–Raphson method to solve the algebraic equation.
Integrating the differential equation
dy | |
dt |
=f(t,y)
tn
tn+1=tn+h
y(tn+1)-y(tn)=
tn+1 | |
\int | |
tn |
f(t,y(t))dt.
y(tn+1)-y(tn) ≈ hf(tn+1,y(tn+1)).
yn
y(tn)
The same reasoning leads to the (standard) Euler method if the left-hand rectangle rule is used instead of the right-hand one.
The local truncation error (defined as the error made in one step) of the backward Euler Method is
O(h2)
t
O(h2)
O(hk+1)
The region of absolute stability for the backward Euler method is the complement in the complex plane of the disk with radius 1 centered at 1, depicted in the figure. This includes the whole left half of the complex plane, making it suitable for the solution of stiff equations. In fact, the backward Euler method is even L-stable.
The region for a discrete stable system by Backward Euler Method is a circle with radius 0.5 which is located at (0.5, 0) in the z-plane.[1]
The backward Euler method is a variant of the (forward) Euler method. Other variants are the semi-implicit Euler method and the exponential Euler method.
The backward Euler method can be seen as a Runge–Kutta method with one stage, described by the Butcher tableau:
\begin{array}{c|c} 1&1\\ \hline &1\\ \end{array}
The method can also be seen as a linear multistep method with one step. It is the first method of the family of Adams–Moulton methods, and also of the family of backward differentiation formulas.