Quant Memo
Core

Fixed-Point Iteration

A simple numerical method for solving an equation by repeatedly plugging a guess back into a rearranged version of the equation, until the output stops changing — a workhorse for equations with no clean algebraic solution.

Many equations that arise in pricing or calibration — a bond's yield to maturity given its price, an implied volatility, a nonlinear model parameter — have no formula you can solve for algebraically. Fixed-point iteration is one of the simplest ways to home in on the answer numerically: rearrange the equation into the form x=g(x)x = g(x), pick a starting guess, and repeatedly feed the output of gg back in as the next guess, hoping the sequence settles down.

Concretely, given x=g(x)x = g(x), start with an initial guess x0x_0 and compute x1=g(x0)x_1 = g(x_0), then x2=g(x1)x_2 = g(x_1), and so on. If the sequence converges — that is, if xnx_n stops changing meaningfully as nn grows — the value it settles on, called a fixed point, satisfies x=g(x)x = g(x) by construction, which is exactly the equation you were trying to solve.

As an example, solving x=cos(x)x = \cos(x) (in radians) by iteration: starting from x0=1x_0 = 1, x1=cos(1)0.540x_1 = \cos(1) \approx 0.540, x2=cos(0.540)0.858x_2 = \cos(0.540) \approx 0.858, x3=cos(0.858)0.653x_3 = \cos(0.858) \approx 0.653, and continuing this bouncing pattern converges after a couple dozen steps to approximately x0.739x \approx 0.739 — the equation's unique solution, reached with nothing more than repeated evaluation of cosine.

Whether the iteration converges at all, and how fast, depends on how "steep" gg is near the fixed point — a rearrangement where gg's slope has magnitude below 1 near the answer will converge, often slowly, while a poorly chosen rearrangement of the same underlying equation can diverge entirely, which is why picking the right form of g(x)g(x) (not just any algebraic rearrangement) matters in practice, and why faster alternatives like Newton-Raphson are often preferred when derivatives are available.

Fixed-point iteration solves x=g(x)x = g(x) by repeatedly feeding each output back in as the next guess until the sequence stops changing, but it only converges reliably when gg is chosen so its slope near the answer is shallow enough — a poor rearrangement of the same equation can fail to converge at all.

Related concepts

Practice in interviews

Further reading

  • Burden & Faires, Numerical Analysis
ShareTwitterLinkedIn