PDF version (NAG web site
, 64-bit version, 64-bit version)
NAG Toolbox: nag_specfun_opt_asian_geom_price (s30sa)
Purpose
nag_specfun_opt_asian_geom_price (s30sa) computes the Asian geometric continuous average-rate option price.
Syntax
[
p,
ifail] = s30sa(
calput,
x,
s,
t,
sigma,
r,
b, 'm',
m, 'n',
n)
[
p,
ifail] = nag_specfun_opt_asian_geom_price(
calput,
x,
s,
t,
sigma,
r,
b, 'm',
m, 'n',
n)
Description
nag_specfun_opt_asian_geom_price (s30sa) computes the price of an Asian geometric continuous average-rate option for constant volatility,
, risk-free rate,
, and cost of carry,
(see
Kemna and Vorst (1990)). For a given strike price,
, the price of a call option with underlying price,
, and time to expiry,
, is
and the corresponding put option price is
where
and
with
is the cumulative Normal distribution function,
The option price is computed for each strike price in a set , , and for each expiry time in a set , .
References
Kemna A and Vorst A (1990) A pricing method for options based on average asset values Journal of Banking and Finance 14 113–129
Parameters
Compulsory Input Parameters
- 1:
– string (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.
Constraint:
or .
- 2:
– double array
-
must contain
, the th strike price, for .
Constraint:
, where , the safe range parameter, for .
- 3:
– double scalar
-
, the price of the underlying asset.
Constraint:
, where , the safe range parameter.
- 4:
– double array
-
must contain
, the th time, in years, to expiry, for .
Constraint:
, where , the safe range parameter, for .
- 5:
– double scalar
-
, the volatility of the underlying asset. Note that a rate of 15% should be entered as 0.15.
Constraint:
.
- 6:
– double scalar
-
, the annual risk-free interest rate, continuously compounded. Note that a rate of 5% should be entered as 0.05.
Constraint:
.
- 7:
– double scalar
-
, the annual cost of carry rate. Note that a rate of 8% should be entered as .
Optional Input Parameters
- 1:
– int64int32nag_int scalar
-
Default:
the dimension of the array
x.
The number of strike prices to be used.
Constraint:
.
- 2:
– int64int32nag_int scalar
-
Default:
the dimension of the array
t.
The number of times to expiry to be used.
Constraint:
.
Output Parameters
- 1:
– double array
-
.
contains , the option price evaluated for the strike price at expiry for and .
- 2:
– int64int32nag_int scalar
unless the function detects an error (see
Error Indicators and Warnings).
Error Indicators and Warnings
Errors or warnings detected by the function:
-
-
On entry, was an illegal value.
-
-
Constraint: .
-
-
Constraint: .
-
-
Constraint: and .
-
-
Constraint: and .
-
-
Constraint: .
-
-
Constraint: .
-
-
Constraint: .
-
-
Constraint: .
-
An unexpected error has been triggered by this routine. Please
contact
NAG.
-
Your licence key may have expired or may not have been installed correctly.
-
Dynamic memory allocation failed.
Accuracy
The accuracy of the output is dependent on the accuracy of the cumulative Normal distribution function,
. This is evaluated using a rational Chebyshev expansion, chosen so that the maximum relative error in the expansion is of the order of the
machine precision (see
nag_specfun_cdf_normal (s15ab) and
nag_specfun_erfc_real (s15ad)). An accuracy close to
machine precision can generally be expected.
Further Comments
None.
Example
This example computes the price of an Asian geometric continuous average-rate put with a time to expiry of months, a stock price of and a strike price of . The risk-free interest rate is per year, the cost of carry is and the volatility is per year.
Open in the MATLAB editor:
s30sa_example
function s30sa_example
fprintf('s30sa example results\n\n');
put = 'P';
s = 80.0;
sigma = 0.2;
r = 0.05;
b = 0.08;
x = [85.0];
t = [0.25];
[p, ifail] = s30sa( ...
put, x, s, t, sigma, r, b);
fprintf('\nAsian Option: Geometric Continuous Average-Rate\nAsian Put :\n');
fprintf(' Spot = %9.4f\n', s);
fprintf(' Volatility = %9.4f\n', sigma);
fprintf(' Rate = %9.4f\n', r);
fprintf(' Cost of carry = %9.4f\n\n', b);
fprintf(' Strike Expiry Option Price\n');
for i=1:1
for j=1:1
fprintf('%9.4f %9.4f %9.4f\n', x(i), t(j), p(i,j));
end
end
s30sa example results
Asian Option: Geometric Continuous Average-Rate
Asian Put :
Spot = 80.0000
Volatility = 0.2000
Rate = 0.0500
Cost of carry = 0.0800
Strike Expiry Option Price
85.0000 0.2500 4.6922
PDF version (NAG web site
, 64-bit version, 64-bit version)
© The Numerical Algorithms Group Ltd, Oxford, UK. 2009–2015