When using an Euler discretisation to simulate a mean-reverting square root process, one runs into the problem that while the process itself is guaranteed to be nonnegative, the discretisation is not. Although an exact and efficient simulation algorithm exists for this process, at present this...

This paper presents a tailor-made discrete-time simulation model for valuing path-dependent options, such as lookback option, barrier option and Asian option. In the context of a real-life application that is interest to many students, we illustrate the option pricing by using Quasi Monte Carlo...

We present a robust method for simulating an increment of a Levy process, based on decomposing the jump part of the process into the sum of its positive and negative jump components. The characteristic exponent of a spectrally one-sided Levy process has excellent analytic properties, which we...

In this article we suggest a new method for solutions of stochastic integrals where the dynamics of the variables in integrand are given by some stochastic differential equation. We also propose numerical simulation of stochastic differential equations which is based on iterated integrals method...

Unlike tranches of synthetic CDOs, that depend only on the defaults of the underlying securities, tranches of cashflow CDOs also depend on the interest cash flows from the coupons of the securities. Whilst fast, accurate, (semi-)analytic methods exist for pricing synthetic CDO tranches (Hull and...

In this article we propose an efficient Monte Carlo scheme for simulating the stochastic volatility model of Heston (1993) enhanced by a non-parametric local volatility component. This hybrid model combines the main advantages of the Heston model and the local volatility model introduced by...

The derivation of Asian option value has posed a challenge to financial mathematicians for the last two decades. Fu, Madan and Wang (1999) made a comparison between the Laplace transform approach and the Monte Carlo approach, and found that the numerical inversion method encountered severe...

Monte Carlo simulations of diffusion processes often introduce bias in the final result, due to time discretization. Using an auxiliary Poisson process, it is possible to run simulations which are unbiased. In this article, we propose such a Monte Carlo scheme which converges to the exact value....

Chen and Shen (2003) argue that it is possible to improve the Least Squares Monte Carlo Method (LSMC) of Longstaff and Schwartz (2001) to value American options by removing the least squares regression module. This would make not only faster but also more accurate. We demonstrate, using a large...

Exact simulation schemes under the Heston stochastic volatility model (e.g., Broadie-Kaya and Glasserman-Kim) suffer from computationally expensive Bessel function evaluations. We propose a new exact simulation scheme without the Bessel function, based on the observation that the conditional...