Crank-Nicolson Pricing
Also known as: CN Method, Implicit Finite Difference
The Crank-Nicolson method is a widely-used implicit finite difference scheme for solving PDEs in option pricing. It provides second-order accuracy in both space and time, unconditional stability, and can efficiently price derivatives with early exercise features (American options) or complex boundary conditions.
Key highlights
- Unconditional stability: large time steps do not cause oscillations or divergence
- Second-order accuracy: O(Delta_t^2 + Delta_x^2) is better than explicit or simple implicit
- Fast solver: tridiagonal matrix solves in O(N) time via Thomas algorithm
- Flexible boundaries: easily incorporates American option early exercise and exotic boundary conditions
Intuition
This section is available to Pro members. Upgrade to Pro
How it works
This section is available to Pro members. Upgrade to Pro
When to use it
Use Crank-Nicolson for American and European option pricing when you need high accuracy and stability. Ideal for single-factor PDEs (spot price or interest rate). For multi-dimensional PDEs (two factors), use alternating direction implicit (ADI) methods. Crank-Nicolson is slow for high-dimensional problems; use Monte Carlo instead.
Strengths & limitations
- Unconditional stability: large time steps do not cause oscillations or divergence
- Second-order accuracy: O(Delta_t^2 + Delta_x^2) is better than explicit or simple implicit
- Fast solver: tridiagonal matrix solves in O(N) time via Thomas algorithm
- Flexible boundaries: easily incorporates American option early exercise and exotic boundary conditions
- One-dimensional: extending to 2+ dimensions requires ADI or splitting methods (slower)
- Grid construction: choosing appropriate spatial/temporal grids requires domain knowledge
- Boundary condition sensitivity: errors in boundary placement or shape propagate backward and degrade interior accuracy
- Not naturally parallel: implicit nature makes parallelization harder than explicit methods
Common pitfalls
This section is available to Pro members. Upgrade to Pro
Applications
This section is available to Pro members. Upgrade to Pro
Frequently asked
How do I choose spatial and temporal grids?
Spatial grid: typically 101-501 points. Start with 101, double until price converges. Temporal grid: match the number of decision points (reset dates, dividends). Stability ratio r = sigma^2 * Delta_t / (2 * Delta_x^2) should be ~0.5 for Crank-Nicolson (not critical for stability, but affects accuracy).
What boundary conditions should I use?
Far OTM: option value ≈ 0. Far ITM: option value ≈ spot - strike (American), or discounted payoff (European). At spot = 0: option value = 0 (stock cannot recover). At spot = infinity: option value ≈ spot - strike for calls (ITM premium dominates).
How do I handle dividends?
Continuous dividend yield q: modify drift term in PDE. Discrete dividends: requires jump condition (option value drops by dividend amount). Alternatively, use an implicit dividend method that updates the solution at dividend dates.
Can I use Crank-Nicolson for two-factor models?
In principle yes, but computational cost is O(M^3) where M is grid size per dimension. ADI methods split the 2D PDE into two 1D problems and reduce cost to O(M^2). Alternating Direction Implicit (ADI) is the standard for 2-factor problems.
Sources
- 1.Crank, J., & Nicolson, P. (1947). A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type. Mathematical Proceedings of the Cambridge Philosophical Society, 43(1), 50-67.
- 2.Fornberg, B. (1996). A Practical Guide to Pseudospectral Methods. Cambridge University Press.
You have read it. What now?
Cite this page
ScholarGate. (2026, June 3). Crank-Nicolson Pricing. ScholarGate. https://scholargate.app/quantitative-finance/crank-nicolson-pricing