PDF version (NAG web site
, 64-bit version, 64-bit version)
NAG Toolbox: nag_specfun_opt_binary_aon_price (s30cc)
Purpose
nag_specfun_opt_binary_aon_price (s30cc) computes the price of a binary or digital asset-or-nothing option.
Syntax
[
p,
ifail] = s30cc(
calput,
x,
s,
t,
sigma,
r,
q, 'm',
m, 'n',
n)
[
p,
ifail] = nag_specfun_opt_binary_aon_price(
calput,
x,
s,
t,
sigma,
r,
q, 'm',
m, 'n',
n)
Description
nag_specfun_opt_binary_aon_price (s30cc) computes the price of a binary or digital asset-or-nothing option which pays the underlying asset itself,
, at expiration if the option is in-the-money (see
Option Pricing Routines in the S Chapter Introduction). For a strike price,
, underlying asset price,
, and time to expiry,
, the payoff is therefore
, if
for a call or
for a put. Nothing is paid out when this condition is not met.
The price of a call with volatility,
, risk-free interest rate,
, and annualised dividend yield,
, is
and for a put,
where
is the cumulative Normal distribution function,
and
The option price is computed for each strike price in a set , , and for each expiry time in a set , .
References
Reiner E and Rubinstein M (1991) Unscrambling the binary code Risk 4
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 continuous yield rate. Note that a rate of 8% should be entered as 0.08.
Constraint:
.
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: .
-
-
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 asset-or-nothing put with a time to expiry of years, a stock price of and a strike price of . The risk-free interest rate is per year, there is an annual dividend return of and the volatility is per year.
Open in the MATLAB editor:
s30cc_example
function s30cc_example
fprintf('s30cc example results\n\n');
put = 'P';
s = 70;
sigma = 0.27;
r = 0.07;
q = 0.05;
x = [65.0];
t = [0.5];
[p, ifail] = s30cc( ...
put, x, s, t, sigma, r, q);
fprintf('\nBinary (Digital): Asset-or-Nothing\n European Put :\n');
fprintf(' Spot = %9.4f\n', s);
fprintf(' Volatility = %9.4f\n', sigma);
fprintf(' Rate = %9.4f\n', r);
fprintf(' Dividend = %9.4f\n\n', q);
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
s30cc example results
Binary (Digital): Asset-or-Nothing
European Put :
Spot = 70.0000
Volatility = 0.2700
Rate = 0.0700
Dividend = 0.0500
Strike Expiry Option Price
65.0000 0.5000 20.2069
PDF version (NAG web site
, 64-bit version, 64-bit version)
© The Numerical Algorithms Group Ltd, Oxford, UK. 2009–2015