Assume the stock price is an Ito-process, satisfying Geometric Brownian Motion dynamics.

Using Ito’s lemma, the natural logarithm of satisfies Arithmetic Brownian motion dynamics:

We can solve this without Ito Integrals because the coefficients of the infinitesimals are constant. The solution is:

Implementation

When implementing this for option pricing, recall the discounted risk-neutral payout expectation formula, where the option payout at time is a function of the realized path of the stock process up to time :

We can use this closed form solution along with the discounted risk-neutral to construct a MC Algorithm

Also under , in our GBM solution (Not sure how to prove this), so we can say:

Black Scholes MC Algorithm

Given : number of simulation runs, number of time discretizations suth that for a maturity time horizon , .

  1. Discretize time s.t. we only simulate at
  2. For (number of simulation runs):
    • Sample a discretized path corresponding to .
    • Compute the payout in simulation run , say . Unlike with binomial models we don’t need to know for .
  3. Return the estimate (the discounted expected value of the average payout at maturity.)

This can be used for path-dependent (Asian) or independent options (European, rainbow?) but, as the algorithm is, the options should be exercised at maturity i.e. no american-type options (you probably just modify it a bit like with algo 2.3)

Sources of error for Monte Carlo:

  • Discretization error as we sample only at and not continuosuly (impossible on a finite resource computer)
  • Sampling error or MC error, due to only sampling paths. We can’t ever sample all uncountably many paths.