create new Heston process
HestonProcess(x0, v0, mu, theta, kappa, sigma, rho)
algebraic expression; initial value of the state variable
algebraic expression; initial value of the variance
algebraic expression; risk-neutral drift
algebraic expression; long-run mean of the variance process
algebraic expression; speed of mean reversion of the variance process
algebraic expression; volatility of the variance process
algebraic expression; instantaneous correlation between the return process and the volatility process
The HestonProcess command creates a new stochastic process S⁡t governed by the following stochastic differential equation (SDE)
μ is the drift parameter
θ is the long-run mean
κ is the speed of mean reversion
σ is the volatility of the variance process
W⁡t is the two-dimensional Wiener process with instantaneous correlation ρ.
This is a stochastic volatility process that was introduced by Heston in A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options.
The parameter X0 defines the initial value of the underlying stochastic process.
The parameter mu is the drift parameter of state process.
The volatility of this process evolves as a SquareRootDiffusion.
The parameter kappa is the speed of mean-reversion of the variance process. The parameter theta is the long-term running mean of the variance process. The parameter sigma is the volatility of the variance process. In general, kappa, theta, and sigma can be any algebraic expressions. However, if the process is to be simulated, these parameters must be assigned numeric values.
The parameter rho is the instantaneous correlation between the state process and the volatility process.
S := HestonProcess(S0, V0, mu, theta, kappa, sigma, rho):
You can now simulate the Heston process.
S0 := 700;
V0 := 0.1;
kappa := 1.0;
theta := 0.1;
mu := 0.05;
sigma := 0.1;
rho := 0.5;
A := SamplePath(S(t), t = 0..3, timesteps = 100, replications = 5);
These are sample paths for the state variables.
PathPlot(A, timegrid = 0..3, 1, thickness = 3, color = red..blue, axes = BOXED, gridlines = true);
And these are the corresponding sample paths for the volatility.
PathPlot(A, timegrid = 0..3, 2, thickness = 3, color = red..blue, axes = BOXED, gridlines = true);
Brigo, D., Mercurio, F., Interest Rate Models: Theory and Practice. New York: Springer-Verlag, 2001.
Gatheral, J., The Volatility Surface: A Practitioner's Guide, (with foreword by Nassim Taleb), Wiley, 2006.
Glasserman, P., Monte Carlo Methods in Financial Engineering. New York: Springer-Verlag, 2004.
Heston, Steven L., A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options, The review of Financial Studies, Volume 6, Issue 2, 327-343, 1993.
Hull, J., Options, Futures, and Other Derivatives, 5th. edition. Upper Saddle River, New Jersey: Prentice Hall, 2003.
The Finance[HestonProcess] command was introduced in Maple 15.
For more information on Maple 15 changes, see Updates in Maple 15.
Download Help Document
What kind of issue would you like to report? (Optional)