Simulation Parameters

$
Current underlying asset price
$
Option exercise price
%
Enter as percentage (e.g., 5 for 5%)
%
Annualized volatility as percentage
years
Time until option expiration
%
Continuous dividend yield (0 for non-dividend stocks)
More = better accuracy, slower calculation
Steps per path (1-252). Mainly affects visualization.
GBM Formula (Hull Eq. 21.16)
S(t+Δt) = S(t) × exp((r - q - σ²/2)Δt + σ√Δt × Z)
S = Stock price | r = Risk-free rate | q = Dividend yield | σ = Volatility | Z ~ N(0,1)
Ryan O'Connell, CFA
Calculator by Ryan O'Connell, CFA

Simulation Results

Monte Carlo Price --
BSM Analytical --
Difference (MC-BSM) --
Standard Error --
95% Confidence Interval --

Simulation Paths

30 sample paths shown. Green = ITM at expiry, Red = OTM at expiry. Dashed line = Strike (K).

Price Convergence

Blue = Running MC average. Dashed green line = BSM analytical price.

Formula Breakdown

MC Price = e-rT × (1/N) × Σ max(ST - K, 0)
Discounted average of simulated payoffs
Model Assumptions
  • Geometric Brownian motion (GBM) under risk-neutral dynamics
  • European exercise only (no early exercise)
  • Continuous dividend yield (if applicable)
  • Constant volatility and risk-free rate over the option's life
  • Pseudo-random i.i.d. normal draws via Box-Muller; results vary slightly between runs

For educational purposes. Not financial advice. Market conventions simplified.

Understanding Monte Carlo Option Pricing

What is Monte Carlo Simulation?

Monte Carlo simulation is a numerical method for pricing options by generating thousands of random stock price paths using risk-neutral dynamics. Rather than solving equations analytically (like Black-Scholes), it simulates many possible future scenarios and averages the results. The approach is based on Hull Chapter 21 (Basic Numerical Procedures).

Monte Carlo Option Price
Call: C = e-rT × (1/N) × Σ max(ST(i) - K, 0)
Put: P = e-rT × (1/N) × Σ max(K - ST(i), 0)
Average discounted payoff across N simulated paths

Geometric Brownian Motion (GBM)

Each simulation path is generated using geometric Brownian motion under the risk-neutral measure. The stock price evolves as:

GBM Step (Hull Eq. 21.16)
S(t+Δt) = S(t) × exp((r - q - σ²/2) × Δt + σ × √Δt × Z)
Z ~ N(0,1) is a standard normal random variable

The -σ²/2 term is the Itô correction, ensuring the expected return under risk-neutral pricing equals the risk-free rate minus dividends.

Monte Carlo vs. Black-Scholes

Monte Carlo

Simulation-based
Flexible: handles exotic payoffs, path-dependence, multiple assets. Converges to analytical price as N increases. Tradeoff: computational cost and random estimation error (standard error).

Black-Scholes

Analytical formula
Exact and instant for European vanilla options. Limited: requires specific assumptions (constant vol, lognormal returns). Cannot price path-dependent or exotic options.

Convergence and Accuracy

The standard error (SE) of a Monte Carlo estimate decreases as 1/√N. This means:

  • Doubling accuracy requires 4× more simulations
  • 10,000 sims typically gives SE < $0.30 for ATM options with 20% volatility
  • 100,000+ sims needed for precision under $0.10
Related: Learn more about Monte Carlo Simulation in Finance.

Frequently Asked Questions

Monte Carlo simulation is a numerical method that prices options by generating thousands of random stock price paths using risk-neutral dynamics (geometric Brownian motion). Each path produces a terminal stock price, from which the option payoff is calculated. The option price is the discounted average of all simulated payoffs. This approach, described in Hull Chapter 21, is especially valuable for pricing complex derivatives where closed-form solutions don't exist, such as path-dependent or multi-asset options.

Black-Scholes-Merton (BSM) provides an exact analytical price for European vanilla options under specific assumptions (constant volatility, lognormal returns, continuous dividend yield). Monte Carlo is a numerical approximation that converges to the BSM price as the number of simulations increases. MC's advantage is flexibility: it can handle exotic payoffs, path-dependent options, stochastic volatility, and multiple underlying assets that BSM cannot. The tradeoff is computational cost and the inherent randomness of the estimate, quantified by standard error.

GBM is the standard stochastic process model for stock prices under the risk-neutral measure. It models the logarithm of stock prices as a continuous-time random walk with drift: S(t+Δt) = S(t) × exp((r - q - σ²/2)Δt + σ√Δt × Z), where r is the risk-free rate, q is the dividend yield, σ is volatility, and Z is a standard normal random variable. The -σ²/2 term is the Itô correction ensuring the expected return equals the risk-free rate minus dividends under risk-neutral pricing.

The standard error decreases as 1/√N, meaning doubling accuracy requires quadrupling simulations. For a typical at-the-money option with these default parameters (S=K=$100, σ=20%, T=1yr), 10,000 simulations typically produce a standard error under $0.30. For precision under $0.10, you generally need 100,000+ simulations. Professional applications often use variance reduction techniques (antithetic variates, control variates) to achieve the same accuracy with fewer simulations.

Standard Monte Carlo simulation cannot directly price American options because it works forward in time, making it impossible to determine the optimal early exercise boundary. The Longstaff-Schwartz least-squares Monte Carlo (LSM) method addresses this by estimating the continuation value at each time step using regression. This calculator implements European options only. For American option approximations, binomial tree methods (also covered in Hull Chapter 21) are more commonly used.

The standard error (SE) measures the precision of the Monte Carlo price estimate. It is calculated as SE = s / √N, where s is the sample standard deviation of the individual discounted payoffs and N is the number of simulations. The 95% confidence interval is MC_price ± 1.96 × SE. In repeated simulations, 95% of such intervals would contain the true option price. Standard error decreases with √N, so achieving tenfold improvement requires 100× more simulations.
Disclaimer

This calculator is for educational purposes only and simulates European option pricing under the GBM model with constant parameters. Actual option prices are affected by stochastic volatility, jumps, discrete dividends, and market microstructure. Results vary between simulation runs due to the random nature of Monte Carlo methods. This tool should not be used for trading decisions.

Related Calculators

Professional finance illustration representing Delta Hedging Calculator

Delta Hedging Calculator

Simulate dynamic delta hedging for options positions. Track P/L, hedge adjustments, and gamma exposure over...

Try Calculator →
Professional finance illustration representing Put Call Parity Calculator

Put Call Parity Calculator

Verify put-call parity relationships and identify potential arbitrage opportunities between call and put options.

Try Calculator →

Course by Ryan O'Connell, CFA, FRM

Options Mastery: From Theory to Practice

Master options trading from theory to practice. Covers fundamentals, Black-Scholes pricing, Greeks, and basic to advanced strategies with hands-on paper trading in Interactive Brokers.

  • 100 lessons with 7 hours of video
  • Black-Scholes, Binomial & Greeks deep dives
  • Basic to advanced strategies (spreads, straddles, condors)
  • Hands-on paper trading with Interactive Brokers
View Course