An introduction to Multilevel Monte Carlo with applications to options.

A standard problem in mathematical finance is the calculation of the price of some financial derivative such as various types of options. Since there exists analytical solutions in only a few cases it will often boil down to estimating the price with Monte Carlo simulation in conjunction with some numerical discretization scheme. The upside of using what we can call standard Monte Carlo is that it is relative straightforward to apply and can be used for a wide variety of problems. The downside is that it has a relatively slow convergence which means that the computational cost or complexity can be very large. However, this slow convergence can be improved upon by using Multilevel Monte Carlo instead of standard Monte Carlo. With this approach it is possible to reduce the computational complexity and cost of simulation considerably. The aim of this thesis is to introduce the reader to the Multilevel Monte Carlo method with applications to European and Asian call options in both the Black-Scholes-Merton (BSM) model and in the Heston model. To this end we first cover the necessary background material such as basic probability theory, estimators and some of their properties, the stochastic integral, stochastic processes and Ito’s theorem. We introduce stochastic differential equations and two numerical discretizations schemes, the Euler–Maruyama scheme and the Milstein scheme. We define strong and weak convergence and illustrate these concepts with examples. We also describe the standard Monte Carlo method and then the theory and implementation of Multilevel Monte Carlo. In the applications part we perform numerical experiments where we compare standard Monte Carlo to Multilevel Monte Carlo in conjunction with the Euler–Maruyama scheme and Milsteins scheme. In the case of a European call in the BSM model, using the Euler–Maruyama scheme, we achieved a cost O(ε-2(log ε)2) to reach the desired error in accordance with theory in comparison to the O(ε-3) cost for standard Monte Carlo. When using Milsteins scheme instead of the Euler–Maruyama scheme it was possible to reduce the cost in terms of the number of simulations needed to achieve the desired error even further. By using Milsteins scheme, a method with greater order of strong convergence than Euler–Maruyama, we achieved the O(ε-2) cost predicted by the complexity theorem compared to the standard Monte Carlo cost of order O(ε-3). In the final numerical experiment we applied the Multilevel Monte Carlo method together with the Euler–Maruyama scheme to an Asian call in the Heston model. In this case, where the coefficients of the Heston model do not satisfy a global Lipschitz condition, the study of strong or weak convergence is much harder. The numerical experiments suggested that the strong convergence was slightly slower compared to what was found in the case of a European call in the BSM model. Nevertheless, we still achieved substantial savings in computational cost compared to using standard Monte Carlo