Skip to content

Euler Method

An explicit first-order ODE integrator that advances an initial state by adding the step size times the derivative evaluated at the current numerical state.

Version
v3 · 2026-09-19 · History
Domain-specific #
1794
Domain group
Formal Sciences
Origin domain
Mathematics
Subdomains
Numerical Analysis, Initial Value Methods → Mathematics
Aliases
Forward Euler method, Explicit Euler method

Core Idea

For the initial-value problem

\[ y'(t)=f(t,y(t)),\qquad y(t_0)=y_0, \]

the Euler method chooses step points \(t_{n+1}=t_n+h\) and advances by

\[ y_{n+1}=y_n+h f(t_n,y_n). \]

The derivative is evaluated only at the known left endpoint, so the update is explicit. Geometrically, the method follows the tangent line at the current numerical point for one step and then recomputes the slope. Repeating this operation produces a polygonal approximation to the unknown solution curve. MIT’s differential-equations materials present it as the simplest numerical ODE method and as the basic idea underlying more sophisticated integrators.

How would you explain it like I'm…

Follow-the-Arrow Steps

Imagine you're walking in the dark and a sign tells you which way to go, but only where you're standing. You take one small step that way, then look at the new sign and step again. Small steps keep you close to the real path. Big steps can make you drift off it.

Follow-the-Slope Stepping

The Euler method is a way to guess how something changes over time when you know where it starts and a rule for how fast it is changing. You use the rule at your current spot to find the direction, take a small straight step that way, and then repeat from your new spot. Joining all the steps makes a zig-zag path that follows the real curve. Making the steps smaller usually makes the guess more accurate. Halving the step roughly halves the error.

Tangent-Line Step Method

The Euler method approximately solves a differential equation y′ = f(t, y) when you are given a starting value y(t₀) = y₀. You choose a step size h and repeat the update y_{n+1} = y_n + h·f(t_n, y_n). The slope is calculated only at the point you already know, the left end of the step, which makes the method explicit. Geometrically, you follow the tangent line for one step and then recompute the slope, which builds a polygon-shaped approximation to the true curve. The method is 'first order', meaning the total error over a fixed interval shrinks roughly in proportion to h. This describes how the error scales and does not mean the answers have only one correct digit.

 

The Euler method applies to the initial-value problem y′ = f(t, y), y(t₀) = y₀. It takes steps t_{n+1} = t_n + h and uses the recurrence y_{n+1} = y_n + h f(t_n, y_n). Because f is evaluated only at the known left endpoint, the scheme is explicit. Each step follows the tangent line at the current numerical point, so the iterates form a polygonal approximation to the solution. A Taylor expansion, y(t_n + h) = y(t_n) + h f(t_n, y(t_n)) + O(h²), shows that one step taken from the exact value has a local defect of O(h²). Over a fixed interval, about 1/h steps accumulate. Under standard Lipschitz and smoothness assumptions, this gives a global error of O(h), which is what 'first order' means here. That term does not refer to the significant digits of the results, nor to the order of the differential equation. Stability is only conditional, so the step size must be small enough relative to the problem. The method is the simplest numerical ODE integrator and the basic idea behind more sophisticated ones.

Scope of Application

Euler Method applies to scalar and vector ordinary differential initial-value problems. A higher-order ODE can be rewritten as a first-order system and stepped componentwise. In simulation and teaching it provides a transparent baseline, a prototype for explicit Runge–Kutta methods, and a predictor within predictor–corrector constructions.

Its practical scope is limited by accuracy and stability. Achieving small error may require many steps because global error is only first order. Stiff decay modes can impose a step much smaller than the time scale of interest. Production solvers therefore usually prefer adaptive higher-order or stability-appropriate methods.

Clarity

The numerical state \(y_n\) approximates \(y(t_n)\); they are not equal by definition after the initial point. Writing \(f(t_n,y(t_n))\) inside the implemented update incorrectly assumes access to the exact solution. The algorithm uses \(f(t_n,y_n)\), so earlier error affects later slopes.

Local truncation error and global error must not be conflated. The one-step Taylor remainder from an exact starting state is proportional to \(h^2\).

Manages Complexity

Many ODEs lack usable closed-form solutions. Euler Method converts continuous evolution into a sequence of elementary evaluations and additions. It needs one derivative evaluation per step and stores only the current state unless output history is desired.

The method exposes the central computational contract shared by one-step integrators: approximate the flow map over a short interval, control local error, and propagate.

Abstract Reasoning

Applied to the linear test equation \(y'=\lambda y\), Euler gives

\[ y_{n+1}=(1+h\lambda)y_n. \]

The exact step multiplies by \(e^{h\lambda}\); Euler replaces this with the first Taylor polynomial \(1+h\lambda\). Absolute stability requires

\[ |1+h\lambda|\le1. \]

For real \(\lambda<0\), this means \(0\le h|\lambda|\le2\). A true solution can decay while Euler oscillates or grows if the step lies outside this interval.

Knowledge Transfer

The recurrence transfers literally from scalar equations to vector systems: \(y\) and \(f\) become vectors, while the update remains unchanged. Mechanical position–velocity systems, compartment models, circuit equations, and population models all use the same numerical operation after being expressed as first-order IVPs.

The method also transfers as a construction component. An Euler prediction can be corrected with a future or averaged slope; Runge–Kutta methods combine several Euler-like stages; time-stepping PDE discretizations may apply forward Euler to a large semi-discrete ODE system.

Relationships to Other Abstractions

Local relationship map for Euler MethodParents appear above the current abstraction, mutual partners to the right, and children below. Node labels state whether each abstraction is prime or domain-specific; colors identify relation types.Euler MethodDOMAINPrime abstraction: Algorithm — is a kind ofAlgorithmPRIME

Current abstraction Euler Method Domain-specific

Parents (1) — outgoing typed relations

  • Euler Method is a kind of Algorithm Prime

    prime:algorithm is the broader abstraction this entry instantiates.

Hierarchy path (1) — routes to 1 parentless root

Neighborhood in Abstraction Space

Euler Method sits in a sparse region of the domain-specific corpus (91st percentile for distinctiveness): few abstractions share its structure, so a faithful description tends to retrieve it precisely.

Family — Unclustered & Miscellaneous (2769 abstractions)

Nearest neighbors

Computed from structural-signature embeddings · 2026-10-25