TOPICS
Search

Heun's Method


Heun's method, also called the explicit trapezoidal method or improved Euler method, is a two-stage second-order explicit Runge-Kutta method. For the initial value problem

 y^'=f(t,y), y(t_0)=y_0,
(1)

it advances an approximation by

k_1=f(t_n,y_n)
(2)
k_2=f(t_n+h,y_n+hk_1)
(3)
y_(n+1)=y_n+h/2(k_1+k_2).
(4)

The first stage is an Euler forward method predictor, while the update averages the slopes at the beginning of the step and at the predicted endpoint. This slope average mirrors the trapezoidal rule, while the predictor keeps Heun's method explicit. For sufficiently smooth f, its local truncation error is O(h^3) and its global error is O(h^2), where O is a Landau symbol.

For the test equation y^'=lambday, the stability function is

 R(z)=1+z+(z^2)/2, z=hlambda.
(5)

The name "Heun's method" is also used for other Runge-Kutta methods introduced by Heun, including a three-stage third-order formula (Butcher 2016).


See also

Butcher Tableau, Euler Forward Method, Initial Value Problem, Predictor-Corrector Methods, Runge-Kutta Method, Trapezoidal Rule

Explore with Wolfram|Alpha

References

Butcher, J. C. Numerical Methods for Ordinary Differential Equations, 3rd ed. Chichester, England: Wiley, 2016.Hairer, E.; Nørsett, S. P.; and Wanner, G. Solving Ordinary Differential Equations I. New York: Springer-Verlag, 1987.Heun, K. "Neue Methoden zur approximativen Integration der Differentialgleichungen einer unabhängigen Veränderlichen." Z. Math. Phys. 45, 23-38, 1900.

Cite this as:

Weisstein, Eric W. "Heun's Method." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/HeunsMethod.html

Subject classifications