Skip to content
PW.← Selected work
02 / Quantitative financeActive

Numerical methods · Python · 2026

Monte Carlo risk, made measurable.

A numerical-finance engine that prices European calls and puts, validates simulation against Black–Scholes and estimates price and Greek uncertainty. Reproducible studies test convergence, confidence intervals and antithetic variance reduction.

View source ↗Examine results ↓
2×less noisy call estimates
600ksimulated market outcomes
95%uncertainty ranges reported

In plain English

What does this engine actually do?

An option is a contract whose value depends on what an asset might be worth in the future. Because that future price is uncertain, the engine simulates thousands of plausible market outcomes, calculates the contract's payoff in each one and averages the results into an estimated fair value.

01

Monte Carlo simulation

A way to study uncertainty by running the same model many times with different random market movements—like repeating an experiment under thousands of plausible futures.

02

Black–Scholes benchmark

A recognised analytical pricing formula. Comparing the simulation with this known answer checks that the engine is producing credible results.

03

Greeks

Risk measures that show how an option's value responds when the market changes—for example, when the asset price or expected volatility moves.

01 / The question

How closely can simulation recover analytical prices and sensitivities, and how honestly can it report the remaining uncertainty?

The project begins with risk-neutral geometric Brownian motion and the closed-form Black–Scholes solution. Those analytical prices provide a benchmark for standard and antithetic Monte Carlo estimators.

Rather than reporting a single simulated price, the engine exposes standard errors and 95% confidence intervals. Tests cover analytical benchmarks, put–call parity, finite-difference Greek validation and the statistical treatment of independent antithetic pair averages.

The latest study extends that validation to pathwise Delta and Vega and a common-random-number Gamma estimator, tracking all three from 2,000 to 600,000 simulated paths.

In practical terms, this tests more than whether the code returns a number. It checks whether that number is accurate, whether its uncertainty is measured honestly and whether the method becomes more reliable as more computation is used.

02 / Method

One model. Prices, sensitivities and uncertainty.

01

Analytical benchmark

Black–Scholes call and put prices establish the deterministic reference.

02

Standard simulation

Exact terminal GBM samples produce discounted payoff observations.

03

Antithetic estimator

Each normal draw is paired with its negative before payoffs are averaged.

04

Monte Carlo Greeks

Pathwise and common-random-number estimators measure Delta, Vega and Gamma.

03 / Results

Lower variance at the same evaluation budget.

Standard and antithetic estimators use equal numbers of terminal-price evaluations. The antithetic estimator consistently produces the lower standard error while retaining the expected N⁻¹ᐟ² scaling.

Spot price100
Strike100
Risk-free rate5%
Volatility20%
Maturity1 year
Evaluation budget600,000
Standard and antithetic Monte Carlo error comparison
OptionStandard SEAntithetic SESE reductionVariance reduction
Call0.0190170.0134601.413x1.996x
Put0.0111880.0085641.306x1.707x

What this means: both methods use the same computational budget, but pairing each simulated market move with its opposite produces a more stable answer. For calls, the remaining statistical variance is almost halved—so the engine extracts roughly twice as much information from the same number of simulated outcomes.

Call and put price estimates from standard and antithetic Monte Carlo converging towards Black-Scholes prices
Standard and antithetic estimates converge towards the analytical call and put prices as the terminal-price evaluation budget increases.
Log-log comparison showing lower standard errors for antithetic call and put estimators
Antithetic standard errors remain below their standard Monte Carlo counterparts across every tested budget.

04 / Greeks convergence

Sensitivity estimates converge with measured uncertainty.

Across a 300-fold increase in path count, standard errors fall by roughly the expected square-root rate. At 600,000 paths, the analytical benchmark lies inside the reported 95% confidence interval for each studied Greek.

Monte Carlo Greek standard-error convergence from 2,000 to 600,000 paths
GreekStarting pathsFinal pathsStarting SEFinal SESE reduction
Call Delta2,000600,0000.0128190.00074417.23x
Gamma2,000600,0000.0025040.00013918.01x
Call Vega2,000600,0001.6005930.09786316.36x
Monte Carlo call and put Delta, Gamma and Vega estimates with 95% confidence intervals converging towards their Black-Scholes benchmarks
Delta, Gamma and Vega estimates converge towards their analytical Black–Scholes values. Error bars show the Monte Carlo 95% confidence interval at each simulation count.
Log-log plots of Delta, Gamma and Vega estimator standard errors following the expected inverse square-root simulation scaling
The estimator standard errors track the N⁻¹ᐟ² reference slope, providing an empirical check of the expected Monte Carlo convergence rate.

Why this matters: a risk figure without uncertainty can look more precise than it really is. These confidence intervals show the range of values supported by the simulation, while the convergence study checks that this range narrows predictably as more paths are added. The underlying data is exported to CSV so the result can be inspected and reproduced.

Download study data ↓

05 / Impact

Why this project matters beyond the model.

Financial models influence decisions only when their outputs can be trusted. This project moves from simply producing a price to showing how accurate it is, how sensitive it is to changing markets and how efficiently that evidence can be generated.

01

More confidence per computation

Antithetic sampling nearly halves call-price variance without increasing the number of simulated market outcomes. The same computational budget therefore produces a more stable estimate.

02

Risk that can be interpreted

Delta, Gamma and Vega translate a single option price into an explanation of how exposure changes when price and volatility move—information that is closer to an actual risk decision.

03

Evidence instead of false precision

Confidence intervals, analytical benchmarks and convergence tests make uncertainty visible. A reviewer can see not only the answer, but also why it should be believed and where its limitations remain.

06 / Engineering

Designed as a reproducible numerical system.

The project demonstrates the full workflow behind trustworthy quantitative software: derive a method, implement it, compare it with an independent benchmark, quantify error and automate the checks so future changes cannot silently break the result.

A laptop displaying source code and analytical charts during a software-development workflow
Numerical research becomes useful when the analysis is implemented, tested and made reproducible.Photo: Daniil Komov / Unsplash
  • Black–Scholes calls, puts and five analytical Greeks
  • Standard and antithetic Monte Carlo estimators
  • Pathwise Delta and Vega with common-random-number Gamma
  • Monte Carlo Greek standard errors and 95% confidence intervals
  • Finite-difference validation of every Greek
  • Price and Greek convergence analysis with CSV export
  • Automated Ruff, pytest and Python 3.13 CI checks
  • Versioned releases and reproducible analysis scripts

07 / Next

From research engine to interactive proof of concept.

The next version will expose the analysis through this website, allowing a visitor to change option parameters and compare prices, confidence intervals and Greeks without installing the Python package.

Explore the repository ↗
PW.

Monte Carlo option pricing · 2026

Selected work ↑