naginterfaces.library.specfun.opt_heston_term¶
- naginterfaces.library.specfun.opt_heston_term(calput, x, fwd, disc, ts, t, alpha, lamda, corr, sigmat, var0)[source]¶
opt_heston_term
computes the European option price given by Heston’s stochastic volatility model with term structure.For full information please refer to the NAG Library document for s30nc
https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/s/s30ncf.html
- Parameters
- calputstr, length 1
Determines whether the option is a call or a put.
A call; the holder has a right to buy.
A put; the holder has a right to sell.
- xfloat, array-like, shape
contains the th strike price, for .
- fwdfloat
The forward price of the asset.
- discfloat
The discount factor, where the current price of the underlying asset, , is given as .
- tsfloat, array-like, shape
must contain the length of the time intervals for which the corresponding element of , , and apply. These should be ordered as they occur in time i.e., .
- tfloat
contains the time to expiry. If then the parameters associated with the last time interval are extended to the expiry time. If then the parameters specified are used up until the expiry time. The rest are ignored.
- alphafloat, array-like, shape
must contain the value of , the value of the volatility of the scaled volatility, , over time subinterval .
- lamdafloat, array-like, shape
must contain the value, , of the mean reversion parameter over the time subinterval .
- corrfloat, array-like, shape
must contain the value, , of the correlation parameter over the time subinterval .
- sigmatfloat, array-like, shape
must contain the value, , of the variance scale factor over the time subinterval .
- var0float
, the initial scaled variance.
- Returns
- pfloat, ndarray, shape
contains the computed option price at the expiry time, , corresponding to strike for the specified term structure, for .
- Raises
- NagValueError
- (errno )
On entry, .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
Quadrature has not converged to the specified accuracy. However, the result should be a reasonable approximation.
- (errno )
Solution cannot be computed accurately. Check values of input arguments.
- Notes
opt_heston_term
computes the price of a European option for Heston’s stochastic volatility model with time-dependent parameters which are piecewise constant. Starting from the stochastic volatility model given by the Stochastic Differential Equation (SDE) system defined by Heston (1993), a scaling of the variance process is introduced, together with a normalization, setting the long run variance, , equal to . This leads towhere is the drift term representing the contribution of interest rates, , and dividends, , while is the scaling parameter, is the scaled variance, is the mean reversion rate and is the volatility of the scaled volatility, . Then, , for , are two standard Brownian motions with correlation parameter . Without loss of generality, the drift term, , is eliminated by modelling the forward price, , directly, instead of the spot price, , with
If required, the spot can be expressed as, , where is the discount factor.
The option price is computed by dividing the time to expiry, , into subintervals and applying the method of characteristic functions to each subinterval, with appropriate initial conditions. Thus, a pair of ordinary differential equations (one of which is a Riccati equation) is solved on each subinterval as outlined in Elices (2008) and Mikhailov and Nögel (2003). Reversing time by taking , the characteristic function solution for the first time subinterval, starting at , is given by Heston (1993), while the solution on each following subinterval uses the solution of the preceding subinterval as initial condition to compute the value of the characteristic function.
In the case of a ‘flat’ term structure, i.e., the parameters are constant over the time of the option, , the form of the SDE system given by Heston (1993) can be recovered by setting , , and .
Conversely, given the Heston form of the SDE pair, to get the term structure form set , , and .
- References
Bain, A, 2011, Private communication
Elices, A, 2008, Models with time-dependent parameters using transform methods: application to Heston’s model, arXiv:0708.2020v2
Heston, S, 1993, A closed-form solution for options with stochastic volatility with applications to bond and currency options, Review of Financial Studies (6), 327–343
Mikhailov, S and Nögel, U, 2003, Heston’s Stochastic Volatility Model Implementation, Calibration and Some Extensions, Wilmott Magazine July/August, 74–79