Runge-Kutta Method

Concept

---

Theoretically, Runge-Kutta Method is a family of numerical integration methods for approximating the solutions of ordinary differential equations (ODEs) with specified initial conditions. It was developed to improve upon simpler numerical methods, particularly the Euler Method, by achieving greater accuracy through multiple evaluations of the derivative function within each integration step. The most widely used member of this family is the classical fourth-order Runge-Kutta (RK4) method, which provides an effective balance between computational efficiency and numerical accuracy. In health economics, Runge-Kutta methods are used when solving differential equation models describing disease progression, pharmacokinetics, epidemiological dynamics and other continuous-time systems that lack closed-form analytical solutions.

Mathematically, the Runge-Kutta framework approximates the solution of an initial value problem by calculating weighted averages of several estimates of the derivative over each integration interval. These intermediate evaluations substantially reduce local truncation error compared with first-order methods. The classical RK4 method has fourth-order global accuracy, meaning that the cumulative approximation error decreases proportionally to h? as the step size h becomes smaller, making it suitable for solving complex dynamic systems encountered in health economic modelling.

In practice, Runge-Kutta methods are implemented within numerical simulation software rather than calculated manually. The modeller specifies the governing differential equations, initial conditions and step size, after which successive approximations are generated throughout the simulation horizon. In health economics, Runge-Kutta methods are applied to pharmacokinetic and pharmacodynamic models, infectious disease transmission models, tumour growth models and physiological simulations that generate inputs for decision-analytic or cost-effectiveness models.

---

Purpose

---

Used to obtain accurate numerical solutions to ordinary differential equations when analytical solutions are unavailable, supporting dynamic disease modelling, pharmacokinetic simulation, epidemiological modelling and other continuous-time health economic analyses.

---

Mathematical Formulae

Primary Formula

For the classical fourth-order Runge-Kutta method:

k? = f(t?, y?)

k? = f(t? + h/2, y? + hk?/2)

k? = f(t? + h/2, y? + hk?/2)

k? = f(t? + h, y? + hk?)

y??? = y? + (h/6) ? (k? + 2k? + 2k? + k?)

Supporting Formulae

Initial value problem:

dy/dt = f(t, y)

y(t?) = y?

Global error:

Error = O(h?)

Local truncation error:

LTE = O(h?)

Related Mathematical Methods

• Euler Method
• Improved Euler Method
• Runge-Kutta-Fehlberg Method
• Adaptive Step-Size Methods
• Dormand-Prince Method
• Numerical Integration
• Initial Value Problem
• Ordinary Differential Equation

---

Example

---

A pharmacokinetic model describes drug concentration C(t) following administration using the differential equation:

dC/dt = ?0.25C

with an initial concentration of 100 mg/L.

Using the classical RK4 method with a step size of h = 1 hour, four intermediate slope estimates are calculated during each hour, producing a more accurate approximation than the Euler Method. The resulting concentration profile can then be incorporated into a pharmacodynamic model to estimate treatment effectiveness, adverse event probabilities and subsequent quality-adjusted life-years within a health economic evaluation.

---

Excel Implementation

Function	Example Formula	Health Economics Application
LET	=LET(k1,F2,k2,F3,k3,F4,k4,F5,E2+(H1/6)*(k1+2*k2+2*k3+k4))	Calculates one RK4 integration step.
LAMBDA	User-defined RK4 calculation	Creates reusable numerical solvers for differential equation models.
EXP	=EXP(-0.25*A2)	Compares RK4 approximation with analytical exponential decay when available.
ABS	=ABS(B2-C2)	Calculates approximation error.
SEQUENCE	=SEQUENCE(101,1,0,0.1)	Generates equally spaced time intervals for numerical integration.

---

VBA (Optional)

---

VBA can automate iterative Runge-Kutta integration across thousands of time steps, allowing efficient simulation of complex pharmacokinetic, epidemiological or physiological models within Excel.

---

Sources

• Butcher JC. Numerical Methods for Ordinary Differential Equations. Wiley.
• Hairer E, N?rsett SP, Wanner G. Solving Ordinary Differential Equations I: Nonstiff Problems. Springer.
• Briggs A, Claxton K, Sculpher M. Decision Modelling for Health Economic Evaluation. Oxford University Press.
• Drummond MF, Sculpher MJ, Claxton K, Stoddart GL, Torrance GW. Methods for the Economic Evaluation of Health Care Programmes. Oxford University Press.
• ISPOR-SMDM Modeling Good Research Practices Task Force reports.