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 , .
- Discretize time s.t. we only simulate at
- 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 .
- 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.