When simulating discrete time approximations of solutions of stochastic differential equations (SDEs), numerical stability is clearly more important than numerical efficiency or some higher order of convergence. Discrete time approximations of solutions of SDEs are widely used in simulations in finance and other areas of application. The stability criterion presented is designed to handle both scenario simulation and Monte Carlo simulation, that is, strong and weak simulation methods. The symmetric predictor-corrector Euler method is shown to have the potential to overcome some of the numerical instabilities that may be experienced when using the explicit Euler method. This is of particular importance in finance, where martingale dynamics arise for solutions of SDEs and diffusion coefficients are often of multiplicative type. Stability regions for a range of schemes are visualized and discussed. For Monte Carlo simulation it turns out that schemes, which have implicitness in both the drift and the diffusion terms, exhibit the largest stability regions. It will be shown that refining the time step size in a Monte Carlo simulation can lead to numerical instabilities.
Download Info
To download:
If you experience problems downloading a file, check if you have the
proper application to
view it first. Information about this may be contained
in the File-Format links below. In case of further problems read
the IDEAS help
page. Note that these files are not on the IDEAS
site. Please be patient as the files may be large.
Publisher Info
Paper provided by Quantitative Finance Research Centre, University of Technology, Sydney in its series Research Paper Series with number
234.
References listed on IDEAS Please report citation or reference errors to , or , if you are the registered author of the cited work, log in to your RePEc Author Service profile, click on "citations" and make appropriate adjustments.: