naginterfaces.library.pde.dim1_​blackscholes_​closed

naginterfaces.library.pde.dim1_blackscholes_closed(kopt, x, s, t, tmat, tdpar, r, q, sigma)

dim1_blackscholes_closed computes an analytic solution to the Black–Scholes equation for a certain set of option types.

For full information please refer to the NAG Library document for d03nd

https://www.nag.com/numeric/nl/nagdoc_29.3/flhtml/d03/d03ndf.html

Parameters

koptint

Specifies the kind of option to be valued:

$kopt=1$

A European call option.

$kopt=2$

An American call option.

$kopt=3$

A European put option.

xfloat

The exercise price $X$.

sfloat

The stock price at which the option value and the Greeks should be evaluated.

tfloat

The time at which the option value and the Greeks should be evaluated.

tmatfloat

The maturity time of the option.

tdparbool, array-like, shape $\left(3\right)$

Specifies whether or not various arguments are time-dependent. More precisely, $r$ is time-dependent if $tdpar\left[0\right]=True$ and constant otherwise. Similarly, $tdpar\left[1\right]$ specifies whether $q$ is time-dependent and $tdpar\left[2\right]$ specifies whether $\sigma$ is time-dependent.

rfloat, array-like, shape $(:)$

Note: the required length for this argument is determined as follows: if $tdpar\left[0\right]=True$: $3$; otherwise: $1$.

If $tdpar\left[0\right]=False$ then $r\left[0\right]$ must contain the constant value of $r$. The remaining elements need not be set.

If $tdpar\left[0\right]=True$ then $r\left[0\right]$ must contain the value of $r$ at time $t$ and $r\left[1\right]$ must contain its average instantaneous value over the remaining life of the option:

$$\hat{r}={\int }_t^{tmat}r\left(\zeta \right)d\zeta \text{.}$$

The auxiliary function dim1_blackscholes_means() may be used to construct $r$ from a set of values of $r$ at discrete times.

qfloat, array-like, shape $(:)$

Note: the required length for this argument is determined as follows: if $tdpar\left[1\right]=True$: $3$; otherwise: $1$.

If $tdpar\left[1\right]=False$ then $q\left[0\right]$ must contain the constant value of $q$. The remaining elements need not be set.

If $tdpar\left[1\right]=True$ then $q\left[0\right]$ must contain the constant value of $q$ and $q\left[1\right]$ must contain its average instantaneous value over the remaining life of the option:

$$\hat{q}={\int }_t^{tmat}q\left(\zeta \right)d\zeta \text{.}$$

The auxiliary function dim1_blackscholes_means() may be used to construct $q$ from a set of values of $q$ at discrete times.

sigmafloat, array-like, shape $(:)$

Note: the required length for this argument is determined as follows: if $tdpar\left[2\right]=True$: $3$; otherwise: $1$.

If $tdpar\left[2\right]=False$ then $sigma\left[0\right]$ must contain the constant value of $\sigma$. The remaining elements need not be set.

If $tdpar\left[2\right]=True$ then $sigma\left[0\right]$ must contain the value of $\sigma$ at time $t$, $sigma\left[1\right]$ the average instantaneous value $\hat{\sigma }$, and $sigma\left[2\right]$ the second-order average $\overline{\sigma }$, where:

$$\hat{\sigma }={\int }_t^{tmat}\sigma \left(\zeta \right)d\zeta \text{,}$$

$$\overline{\sigma }={\left({\int }_t^{tmat}{\sigma }^2\left(\zeta \right)d\zeta \right)}^{1/2}\text{.}$$

The auxiliary function dim1_blackscholes_means() may be used to compute $sigma$ from a set of values at discrete times.

Returns

ffloat

The value $f$ of the option at the stock price $s$ and time $t$.

thetafloat

The values of various Greeks at the stock price $s$ and time $t$, as follows:

$$\begin{array}{c}\begin{array}{ccc}theta=\Theta =\frac{\partial f}{\partial t}\text{,}& delta=\Delta =\frac{\partial f}{\partial s}\text{,}& gamma=\Gamma =\frac{{\partial }^{2}f}{\partial s^2}\text{,}\\ & & \\ lamda=\Lambda =\frac{\partial f}{\partial \sigma }\text{,}& rho=\rho =\frac{\partial f}{\partial r}\text{.}& \end{array}\end{array}$$

deltafloat

The values of various Greeks at the stock price $s$ and time $t$, as follows:

$$\begin{array}{c}\begin{array}{ccc}theta=\Theta =\frac{\partial f}{\partial t}\text{,}& delta=\Delta =\frac{\partial f}{\partial s}\text{,}& gamma=\Gamma =\frac{{\partial }^{2}f}{\partial s^2}\text{,}\\ & & \\ lamda=\Lambda =\frac{\partial f}{\partial \sigma }\text{,}& rho=\rho =\frac{\partial f}{\partial r}\text{.}& \end{array}\end{array}$$

gammafloat

The values of various Greeks at the stock price $s$ and time $t$, as follows:

$$\begin{array}{c}\begin{array}{ccc}theta=\Theta =\frac{\partial f}{\partial t}\text{,}& delta=\Delta =\frac{\partial f}{\partial s}\text{,}& gamma=\Gamma =\frac{{\partial }^{2}f}{\partial s^2}\text{,}\\ & & \\ lamda=\Lambda =\frac{\partial f}{\partial \sigma }\text{,}& rho=\rho =\frac{\partial f}{\partial r}\text{.}& \end{array}\end{array}$$

lamdafloat

The values of various Greeks at the stock price $s$ and time $t$, as follows:

$$\begin{array}{c}\begin{array}{ccc}theta=\Theta =\frac{\partial f}{\partial t}\text{,}& delta=\Delta =\frac{\partial f}{\partial s}\text{,}& gamma=\Gamma =\frac{{\partial }^{2}f}{\partial s^2}\text{,}\\ & & \\ lamda=\Lambda =\frac{\partial f}{\partial \sigma }\text{,}& rho=\rho =\frac{\partial f}{\partial r}\text{.}& \end{array}\end{array}$$

rhofloat

The values of various Greeks at the stock price $s$ and time $t$, as follows:

$$\begin{array}{c}\begin{array}{ccc}theta=\Theta =\frac{\partial f}{\partial t}\text{,}& delta=\Delta =\frac{\partial f}{\partial s}\text{,}& gamma=\Gamma =\frac{{\partial }^{2}f}{\partial s^2}\text{,}\\ & & \\ lamda=\Lambda =\frac{\partial f}{\partial \sigma }\text{,}& rho=\rho =\frac{\partial f}{\partial r}\text{.}& \end{array}\end{array}$$

Raises

NagValueError

(errno $1$)

On entry, $sigma[⟨value⟩]=⟨value⟩$.

Constraint: $sigma[\text{i}-1]>0.0$.

(errno $1$)

On entry, $tmat=⟨value⟩$ and $t=⟨value⟩$.

Constraint: $tmat\ge t$.

(errno $1$)

On entry, $t=⟨value⟩$.

Constraint: $t\ge 0.0$.

(errno $1$)

On entry, $s=⟨value⟩$.

Constraint: $s\ge 0.0$.

(errno $1$)

On entry, $x=⟨value⟩$.

Constraint: $x\ge 0.0$.

(errno $1$)

On entry, $q\left[0\right]$ is not equal to $0.0$ with American call option. $q\left[0\right]=⟨value⟩$.

(errno $1$)

On entry, $kopt=⟨value⟩$.

Constraint: $kopt=1$, $2$ or $3$.

Notes

dim1_blackscholes_closed computes an analytic solution to the Black–Scholes equation (see Hull (1989) and Wilmott et al. (1995))

$$\frac{\partial f}{\partial t}+(r-q)S\frac{\partial f}{\partial S}+\frac{{\sigma }^2{S}^2}{2}\frac{{\partial }^{2}f}{\partial S^2}=rf$$

$$S_{min}<S<S_{max}\text{,}{t}_{min}<t<t_{max}\text{,}$$

for the value $f$ of a European put or call option, or an American call option with zero dividend $q$. In equation (1) $t$ is time, $S$ is the stock price, $X$ is the exercise price, $r$ is the risk free interest rate, $q$ is the continuous dividend, and $\sigma$ is the stock volatility. The parameter $r$, $q$ and $\sigma$ may be either constant, or functions of time. In the latter case their average instantaneous values over the remaining life of the option should be provided to dim1_blackscholes_closed. An auxiliary function dim1_blackscholes_means() is available to compute such averages from values at a set of discrete times. Equation (1) is subject to different boundary conditions depending on the type of option. For a call option the boundary condition is

$$f(S,t=t_{\text{mat}})=\text{max}(0,S-X)$$

where $t_{\text{mat}}$ is the maturity time of the option. For a put option the equation (1) is subject to

$$f(S,t=t_{\text{mat}})=\text{max}(0,X-S)\text{.}$$

dim1_blackscholes_closed also returns values of the Greeks

$$\Theta =\frac{\partial f}{\partial t}\text{,}\Delta =\frac{\partial f}{\partial x}\text{,}\Gamma =\frac{{\partial }^{2}f}{\partial x^2}\text{,}\Lambda =\frac{\partial f}{\partial \sigma }\text{,}\rho =\frac{\partial f}{\partial r}\text{.}$$

specfun.opt_bsm_greeks also computes the European option price given by the Black–Scholes–Merton formula together with a more comprehensive set of sensitivities (Greeks).

Further details of the analytic solution returned are given in Algorithmic Details.

References

Hull, J, 1989, Options, Futures and Other Derivative Securities, Prentice–Hall

Wilmott, P, Howison, S and Dewynne, J, 1995, The Mathematics of Financial Derivatives, Cambridge University Press