Many phenomena in our world can be described as differential equations in high dimensions. However, they are notoriously challenging to solve numerically due to the exponential growth in computational cost with increasing dimensions. This thesis explores an algorithm, known as deep BSDE, for solving high-dimensional partial differential equations and applies it to finance, namely European option pricing. In addition, an implementation of the method is provided that seemingly shortens the runtime by a factor of two, compared with the results in previous studies. From the results, we can conclude that the deep BSDE method does handle high-dimensional problems well. Lastly, the thesis gives the relevant prerequisites required to be able to digest the theory from an undergraduate level.

This thesis investigated Stochastic Runge–Kutta Lawson (SRKL) schemes and their application to the Heston model. Two distinct SRKL discretization methods were used to simulate a single asset’s dynamics under the Heston model, notably the Euler–Maruyama and Midpoint schemes. Additionally, standard Monte Carlo and variance reduction techniques were implemented. European and Asian option prices were estimated and compared with a benchmark value regarding accuracy, effectiveness, and computational complexity. Findings showed that the SRKL Euler–Maruyama schemes exhibited promise in enhancing the price for simple and path-dependent options. Consequently, integrating SRKL numerical methods into option valuation provides notable advantages by addressing challenges posed by the Heston model’s SDEs. Given the limited scope of this research topic, it is imperative to conduct further studies to understand the use of SRKL schemes within other models.