hide long namesshow long names
hide short namesshow short names
Integer type:  int32  int64  nag_int  show int32  show int32  show int64  show int64  show nag_int  show nag_int

PDF version (NAG web site, 64-bit version, 64-bit version)
Chapter Contents
Chapter Introduction
NAG Toolbox

NAG Toolbox: nag_specfun_opt_asian_geom_price (s30sa)

 Contents

    1  Purpose
    2  Syntax
    7  Accuracy
    9  Example

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, r, and cost of carry, b (see Kemna and Vorst (1990)). For a given strike price, X, the price of a call option with underlying price, S, and time to expiry, T, is
Pcall = S e b--r T Φ d- 1 - X e-rT Φ d- 2 ,  
and the corresponding put option price is
Pput = X e-rT Φ -d-2 - S e b--r T Φ - d-1 ,  
where
d-1 = lnS/X + b- + σ-2 / 2 T σ- T  
and
d-2 = d-1 - σ- T ,  
with
σ- = σ 3 ,  b- = 1 2 r- σ2 6 .  
Φ is the cumulative Normal distribution function,
Φx = 1 2π - x exp -y2/2 dy .  
The option price Pij=PX=Xi,T=Tj is computed for each strike price in a set Xi, i=1,2,,m, and for each expiry time in a set Tj, j=1,2,,n.

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:     calput – string (length ≥ 1)
Determines whether the option is a call or a put.
calput='C'
A call; the holder has a right to buy.
calput='P'
A put; the holder has a right to sell.
Constraint: calput='C' or 'P'.
2:     xm – double array
xi must contain Xi, the ith strike price, for i=1,2,,m.
Constraint: xiz ​ and ​ xi 1 / z , where z = x02am , the safe range parameter, for i=1,2,,m.
3:     s – double scalar
S, the price of the underlying asset.
Constraint: sz ​ and ​s1.0/z, where z=x02am, the safe range parameter.
4:     tn – double array
ti must contain Ti, the ith time, in years, to expiry, for i=1,2,,n.
Constraint: tiz, where z = x02am , the safe range parameter, for i=1,2,,n.
5:     sigma – double scalar
σ, the volatility of the underlying asset. Note that a rate of 15% should be entered as 0.15.
Constraint: sigma>0.0.
6:     r – double scalar
r, the annual risk-free interest rate, continuously compounded. Note that a rate of 5% should be entered as 0.05.
Constraint: r0.0.
7:     b – double scalar
b, the annual cost of carry rate. Note that a rate of 8% should be entered as 0.08.

Optional Input Parameters

1:     m int64int32nag_int scalar
Default: the dimension of the array x.
The number of strike prices to be used.
Constraint: m1.
2:     n int64int32nag_int scalar
Default: the dimension of the array t.
The number of times to expiry to be used.
Constraint: n1.

Output Parameters

1:     pldpn – double array
ldp=m.
pij contains Pij, the option price evaluated for the strike price xi at expiry tj for i=1,2,,m and j=1,2,,n.
2:     ifail int64int32nag_int scalar
ifail=0 unless the function detects an error (see Error Indicators and Warnings).

Error Indicators and Warnings

Errors or warnings detected by the function:
   ifail=1
On entry, calput=_ was an illegal value.
   ifail=2
Constraint: m1.
   ifail=3
Constraint: n1.
   ifail=4
Constraint: xi_ and xi_.
   ifail=5
Constraint: s_ and s_.
   ifail=6
Constraint: ti_.
   ifail=7
Constraint: sigma>0.0.
   ifail=8
Constraint: r0.0.
   ifail=11
Constraint: ldpm.
   ifail=-99
An unexpected error has been triggered by this routine. Please contact NAG.
   ifail=-399
Your licence key may have expired or may not have been installed correctly.
   ifail=-999
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 3 months, a stock price of 80 and a strike price of 85. The risk-free interest rate is 5% per year, the cost of carry is 8% and the volatility is 20% per year.
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)
Chapter Contents
Chapter Introduction
NAG Toolbox

© The Numerical Algorithms Group Ltd, Oxford, UK. 2009–2015