False Position Method (Regula Falsi)

Summary

False position keeps a sign-changing bracket like bisection, but places the next sample at the intercept of the secant through (a,f(a)) and (b,f(b)) . It often reduces the residual faster than pure bisection while preserving a bracket.

Prerequisites

Problem Type

Solve f(x)=0 on a continuous bracket with opposite signs.

Method Definition

Given continuous f on [a,b] with f(a)f(b)<0 , set

c=af(b)bf(a)f(b)f(a)

and replace the endpoint whose function value has the same sign as f(c) .[1]

Assumptions / Requirements

Algorithm

  1. Evaluate f(a),f(b) with opposite signs.
  2. Compute c from the formula above.
  3. If |f(c)| or |ba| is small enough, stop.
  4. If f(a)f(c)<0 then bc , else ac .
  5. Repeat.

Formula / Iteration Rule

ck=akf(bk)bkf(ak)f(bk)f(ak)

Equivalent form:

ck=akf(ak)bkakf(bk)f(ak)

Convergence

Typically faster residual reduction than bisection when f is roughly linear near the root. Convergence is still linear in general. One endpoint can stagnate for convex/concave f ; modified regula falsi variants fix this.

Error / Accuracy

Common stops: |f(c)|<ε or |cnewcold|<ε . The interval length need not shrink as aggressively as in bisection.

Worked Example

Again use f(x)=x32x5 . Valid bracket: [2,3] because f(2)=1 and f(3)=16 .

Iteration 1:

c1=2163(1)16(1)=32+317=35172.05882 f(c1)0.3908<0new bracket [c1,3]

Iteration 2:

c2=2.05882163f(c1)16f(c1)2.08126,f(c2)0.1472

Subsequent iterates approach 2.09455 from below while the right endpoint stays at 3 for several steps (classic regula falsi stagnation pattern).

Common Failure Modes

Connections

References


  1. Burden & Faires, Numerical Analysis, regula falsi; NIST DLMF Ch. 3, https://dlmf.nist.gov/3 ↩︎