Black–Scholes model

The Black–Scholes /ˌblæk ˈʃoʊlz/[1] or Black–Scholes–Merton model is a mathematical model for the dynamics of a financial market containing derivative investment instruments. From the partial differential equation in the model, known as the Black–Scholes equation, one can deduce the Black–Scholes formula, which gives a theoretical estimate of the price of European-style options and shows that the option has a unique price regardless of the risk of the security and its expected return (instead replacing the security's expected return with the risk-neutral rate). The formula led to a boom in options trading and provided mathematical legitimacy to the activities of the Chicago Board Options Exchange and other options markets around the world.[2] It is widely used, although often with some adjustments, by options market participants.[3]^:751

Based on works previously developed by market researchers and practitioners, such as Louis Bachelier, Sheen Kassouf and Ed Thorp among others, Fischer Black and Myron Scholes demonstrated in 1968 that a dynamic revision of a portfolio removes the expected return of the security, thus inventing the risk neutral argument.[4][5] In 1970, after they attempted to apply the formula to the markets and incurred financial losses due to lack of risk management in their trades, they decided to focus in their domain area, the academic environment.[6] After three years of efforts, the formula—named in honor of them for making it public—was finally published in 1973 in an article entitled "The Pricing of Options and Corporate Liabilities", in the Journal of Political Economy.[7][8][9] Robert C. Merton was the first to publish a paper expanding the mathematical understanding of the options pricing model, and coined the term "Black–Scholes options pricing model". Merton and Scholes received the 1997 Nobel Memorial Prize in Economic Sciences for their work, the committee citing their discovery of the risk neutral dynamic revision as a breakthrough that separates the option from the risk of the underlying security.[10] Although ineligible for the prize because of his death in 1995, Black was mentioned as a contributor by the Swedish Academy.[11]

The key idea behind the model is to hedge the option by buying and selling the underlying asset in just the right way and, as a consequence, to eliminate risk. This type of hedging is called "continuously revised delta hedging" and is the basis of more complicated hedging strategies such as those engaged in by investment banks and hedge funds.

The model's assumptions have been relaxed and generalized in many directions, leading to a plethora of models that are currently used in derivative pricing and risk management. It is the insights of the model, as exemplified in the Black–Scholes formula, that are frequently used by market participants, as distinguished from the actual prices. These insights include no-arbitrage bounds and risk-neutral pricing (thanks to continuous revision). Further, the Black–Scholes equation, a partial differential equation that governs the price of the option, enables pricing using numerical methods when an explicit formula is not possible.

The Black–Scholes formula has only one parameter that cannot be directly observed in the market: the average future volatility of the underlying asset, though it can be found from the price of other options. Since the option value (whether put or call) is increasing in this parameter, it can be inverted to produce a "volatility surface" that is then used to calibrate other models, e.g. for OTC derivatives.

Fundamental hypotheses

The Black–Scholes model assumes that the market consists of at least one risky asset, usually called the stock, and one riskless asset, usually called the money market, cash, or bond.

Now we make assumptions on the assets (which explain their names):

• (riskless rate) The rate of return on the riskless asset is constant and thus called the risk-free interest rate.
• (random walk) The instantaneous log return of stock price is an infinitesimal random walk with drift; more precisely, the stock price follows a geometric Brownian motion, and we will assume its drift and volatility are constant (if they are time-varying, we can deduce a suitably modified Black–Scholes formula quite simply, as long as the volatility is not random).
• The stock does not pay a dividend.[Notes 1]

The assumptions on the market are:

• no arbitrage opportunity (i.e., there is no way to make a riskless profit).
• ability to borrow and lend any amount, even fractional, of cash at the riskless rate.
• ability to buy and sell any amount, even fractional, of the stock (this includes short selling).
• The above transactions do not incur any fees or costs (i.e., frictionless market).

With these assumptions holding, suppose there is a derivative security also trading in this market. We specify that this security will have a certain payoff at a specified date in the future, depending on the values taken by the stock up to that date. It is a surprising fact that the derivative's price is completely determined at the current time, even though we do not know what path the stock price will take in the future. For the special case of a European call or put option, Black and Scholes showed that "it is possible to create a hedged position, consisting of a long position in the stock and a short position in the option, whose value will not depend on the price of the stock".[12] Their dynamic hedging strategy led to a partial differential equation which governed the price of the option. Its solution is given by the Black–Scholes formula.

Several of these assumptions of the original model have been removed in subsequent extensions of the model. Modern versions account for dynamic interest rates (Merton, 1976), transaction costs and taxes (Ingersoll, 1976), and dividend payout.[13]

Notation

The notation used throughout this page will be defined as follows:

${\displaystyle S(t)}$, the price of the underlying asset at time t.;

${\displaystyle V(S,t)}$, the price of the option as a function of the underlying asset S, at time t;

${\displaystyle C(S,t)}$, the price of a European call option and ${\displaystyle P(S,t)}$ the price of a European put option;

${\displaystyle K}$, the strike price of the option, also known as the exercise price;

${\displaystyle r}$, the annualized risk-free interest rate, continuously compounded Also known as the force of interest;

${\displaystyle \mu }$, the drift rate of ${\displaystyle S}$, annualized;

${\displaystyle \sigma }$, the standard deviation of the stock's returns; this is the square root of the quadratic variation of the stock's log price process;

${\displaystyle t}$, a time in years; we generally use: now ${\displaystyle =0}$, expiry ${\displaystyle =T}$;

${\displaystyle \Pi }$, the value of the portfolio.

We will use ${\displaystyle N(x)}$ to denote the standard normal cumulative distribution function,

${\displaystyle N(x)={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{x}e^{-z^{2}/2}\,dz.}$

${\displaystyle N'(x)}$ will denote the standard normal probability density function,

${\displaystyle N'(x)={\frac {1}{\sqrt {2\pi }}}e^{-x^{2}/2}.}$

Black–Scholes equation

Main article: Black–Scholes equation

Simulated geometric Brownian motions with parameters from market data

As above, the Black–Scholes equation is a partial differential equation, which describes the price of the option over time. The equation is:

${\displaystyle {\frac {\partial V}{\partial t}}+{\frac {1}{2}}\sigma ^{2}S^{2}{\frac {\partial ^{2}V}{\partial S^{2}}}+rS{\frac {\partial V}{\partial S}}-rV=0}$

The key financial insight behind the equation is that one can perfectly hedge the option by buying and selling the underlying asset and the bank account asset (cash) in just the right way and consequently "eliminate risk". This hedge, in turn, implies that there is only one right price for the option, as returned by the Black–Scholes formula (see the next section).

Black–Scholes formula

A European call valued using the Black–Scholes pricing equation for varying asset price ${\displaystyle S}$ and time-to-expiry ${\displaystyle T}$. In this particular example, the strike price is set to 1.

The Black–Scholes formula calculates the price of European put and call options. This price is consistent with the Black–Scholes equation as above; this follows since the formula can be obtained by solving the equation for the corresponding terminal and boundary conditions.

The value of a call option for a non-dividend-paying underlying stock in terms of the Black–Scholes parameters is:

${\displaystyle {\begin{aligned}C(S_{t},t)&=N(d_{1})S_{t}-N(d_{2})PV(K)\\d_{1}&={\frac {1}{\sigma {\sqrt {T-t}}}}\left[\ln \left({\frac {S_{t}}{K}}\right)+\left(r+{\frac {\sigma ^{2}}{2}}\right)(T-t)\right]\\d_{2}&=d_{1}-\sigma {\sqrt {T-t}}\\PV(K)&=Ke^{-r(T-t)}\end{aligned}}}$

The price of a corresponding put option based on put–call parity is:

${\displaystyle {\begin{aligned}P(S_{t},t)&=Ke^{-r(T-t)}-S_{t}+C(S_{t},t)\\&=N(-d_{2})Ke^{-r(T-t)}-N(-d_{1})S_{t}\end{aligned}}\,}$

For both, as above:

• ${\displaystyle N(\cdot )}$ is the cumulative distribution function of the standard normal distribution
• ${\displaystyle T-t}$ is the time to maturity (expressed in years)
• ${\displaystyle S_{t}}$ is the spot price of the underlying asset
• ${\displaystyle K}$ is the strike price
• ${\displaystyle r}$ is the risk free rate (annual rate, expressed in terms of continuous compounding)
• ${\displaystyle \sigma }$ is the volatility of returns of the underlying asset

Alternative formulation

Introducing some auxiliary variables allows the formula to be simplified and reformulated in a form that is often more convenient (this is a special case of the Black '76 formula):

${\displaystyle {\begin{aligned}C(F,\tau )&=D\left[N(d_{+})F-N(d_{-})K\right]\\d_{\pm }&={\frac {1}{\sigma {\sqrt {\tau }}}}\left[\ln \left({\frac {F}{K}}\right)\pm {\frac {1}{2}}\sigma ^{2}\tau \right]\\d_{\pm }&=d_{\mp }\pm \sigma {\sqrt {\tau }}\end{aligned}}}$

The auxiliary variables are:

• ${\displaystyle \tau =T-t}$ is the time to expiry (remaining time, backwards time)
• ${\displaystyle D=e^{-r\tau }}$ is the discount factor
• ${\displaystyle F=e^{r\tau }S={\frac {S}{D}}}$ is the forward price of the underlying asset, and ${\displaystyle S=DF}$

with d₊ = d₁ and d₋ = d₂ to clarify notation.

Given put–call parity, which is expressed in these terms as:

${\displaystyle C-P=D(F-K)=S-DK}$

the price of a put option is:

${\displaystyle P(F,\tau )=D\left[N(-d_{-})K-N(-d_{+})F\right]}$

Interpretation

The Black–Scholes formula can be interpreted fairly handily, with the main subtlety the interpretation of the ${\displaystyle N(d_{\pm })}$ (and a fortiori ${\displaystyle d_{\pm }}$) terms, particularly ${\displaystyle d_{+}}$ and why there are two different terms.[14]

The formula can be interpreted by first decomposing a call option into the difference of two binary options: an asset-or-nothing call minus a cash-or-nothing call (long an asset-or-nothing call, short a cash-or-nothing call). A call option exchanges cash for an asset at expiry, while an asset-or-nothing call just yields the asset (with no cash in exchange) and a cash-or-nothing call just yields cash (with no asset in exchange). The Black–Scholes formula is a difference of two terms, and these two terms equal the values of the binary call options. These binary options are much less frequently traded than vanilla call options, but are easier to analyze.

Thus the formula:

${\displaystyle C=D\left[N(d_{+})F-N(d_{-})K\right]}$

breaks up as:

${\displaystyle C=DN(d_{+})F-DN(d_{-})K,}$

where ${\displaystyle DN(d_{+})F}$ is the present value of an asset-or-nothing call and ${\displaystyle DN(d_{-})K}$ is the present value of a cash-or-nothing call. The D factor is for discounting, because the expiration date is in future, and removing it changes present value to future value (value at expiry). Thus ${\displaystyle N(d_{+})~F}$ is the future value of an asset-or-nothing call and ${\displaystyle N(d_{-})~K}$ is the future value of a cash-or-nothing call. In risk-neutral terms, these are the expected value of the asset and the expected value of the cash in the risk-neutral measure.

The naive, and not quite correct, interpretation of these terms is that ${\displaystyle N(d_{+})F}$ is the probability of the option expiring in the money ${\displaystyle N(d_{+})}$, times the value of the underlying at expiry F, while ${\displaystyle N(d_{-})K}$ is the probability of the option expiring in the money ${\displaystyle N(d_{-}),}$ times the value of the cash at expiry K. This is obviously incorrect, as either both binaries expire in the money or both expire out of the money (either cash is exchanged for asset or it is not), but the probabilities ${\displaystyle N(d_{+})}$ and ${\displaystyle N(d_{-})}$ are not equal. In fact, ${\displaystyle d_{\pm }}$ can be interpreted as measures of moneyness (in standard deviations) and ${\displaystyle N(d_{\pm })}$ as probabilities of expiring ITM (percent moneyness), in the respective numéraire, as discussed below. Simply put, the interpretation of the cash option, ${\displaystyle N(d_{-})K}$, is correct, as the value of the cash is independent of movements of the underlying, and thus can be interpreted as a simple product of "probability times value", while the ${\displaystyle N(d_{+})F}$ is more complicated, as the probability of expiring in the money and the value of the asset at expiry are not independent.[14] More precisely, the value of the asset at expiry is variable in terms of cash, but is constant in terms of the asset itself (a fixed quantity of the asset), and thus these quantities are independent if one changes numéraire to the asset rather than cash.

If one uses spot S instead of forward F, in ${\displaystyle d_{\pm }}$ instead of the ${\displaystyle {\frac {1}{2}}\sigma ^{2}}$ term there is ${\displaystyle \left(r\pm {\frac {1}{2}}\sigma ^{2}\right)\tau ,}$ which can be interpreted as a drift factor (in the risk-neutral measure for appropriate numéraire). The use of d₋ for moneyness rather than the standardized moneyness ${\displaystyle m={\frac {1}{\sigma {\sqrt {\tau }}}}\ln \left({\frac {F}{K}}\right)}$ – in other words, the reason for the ${\displaystyle {\frac {1}{2}}\sigma ^{2}}$ factor – is due to the difference between the median and mean of the log-normal distribution; it is the same factor as in Itō's lemma applied to geometric Brownian motion. In addition, another way to see that the naive interpretation is incorrect is that replacing N(d₊) by N(d₋) in the formula yields a negative value for out-of-the-money call options.[14]^:6

In detail, the terms ${\displaystyle N(d_{1}),N(d_{2})}$ are the probabilities of the option expiring in-the-money under the equivalent exponential martingale probability measure (numéraire=stock) and the equivalent martingale probability measure (numéraire=risk free asset), respectively.[14] The risk neutral probability density for the stock price ${\displaystyle S_{T}\in (0,\infty )}$ is

${\displaystyle p(S,T)={\frac {N^{\prime }[d_{2}(S_{T})]}{S_{T}\sigma {\sqrt {T}}}}}$

where ${\displaystyle d_{2}=d_{2}(K)}$ is defined as above.

Specifically, ${\displaystyle N(d_{2})}$ is the probability that the call will be exercised provided one assumes that the asset drift is the risk-free rate. ${\displaystyle N(d_{1})}$, however, does not lend itself to a simple probability interpretation. ${\displaystyle SN(d_{1})}$ is correctly interpreted as the present value, using the risk-free interest rate, of the expected asset price at expiration, given that the asset price at expiration is above the exercise price.[15] For related discussion – and graphical representation – see section "Interpretation" under Datar–Mathews method for real option valuation.

The equivalent martingale probability measure is also called the risk-neutral probability measure. Note that both of these are probabilities in a measure theoretic sense, and neither of these is the true probability of expiring in-the-money under the real probability measure. To calculate the probability under the real ("physical") probability measure, additional information is required—the drift term in the physical measure, or equivalently, the market price of risk.

Derivations

See also: Martingale pricing

A standard derivation for solving the Black–Scholes PDE is given in the article Black–Scholes equation.

The Feynman–Kac formula says that the solution to this type of PDE, when discounted appropriately, is actually a martingale. Thus the option price is the expected value of the discounted payoff of the option. Computing the option price via this expectation is the risk neutrality approach and can be done without knowledge of PDEs.[14] Note the expectation of the option payoff is not done under the real world probability measure, but an artificial risk-neutral measure, which differs from the real world measure. For the underlying logic see section "risk neutral valuation" under Rational pricing as well as section "Derivatives pricing: the Q world" under Mathematical finance; for detail, once again, see Hull.[16]^:307–309

The Greeks

"The Greeks" measure the sensitivity of the value of a derivative or a portfolio to changes in parameter values while holding the other parameters fixed. They are partial derivatives of the price with respect to the parameter values. One Greek, "gamma" (as well as others not listed here) is a partial derivative of another Greek, "delta" in this case.

The Greeks are important not only in the mathematical theory of finance, but also for those actively trading. Financial institutions will typically set (risk) limit values for each of the Greeks that their traders must not exceed. Delta is the most important Greek since this usually confers the largest risk. Many traders will zero their delta at the end of the day if they are not speculating on the direction of the market and following a delta-neutral hedging approach as defined by Black–Scholes.

The Greeks for Black–Scholes are given in closed form below. They can be obtained by differentiation of the Black–Scholes formula.[17]

		Calls	Puts
Delta	${\displaystyle {\frac {\partial V}{\partial S}}}$	${\displaystyle N(d_{1})\,}$	${\displaystyle -N(-d_{1})=N(d_{1})-1\,}$
Gamma	${\displaystyle {\frac {\partial ^{2}V}{\partial S^{2}}}}$	${\displaystyle {\frac {N'(d_{1})}{S\sigma {\sqrt {T-t}}}}\,}$
Vega	${\displaystyle {\frac {\partial V}{\partial \sigma }}}$	${\displaystyle SN'(d_{1}){\sqrt {T-t}}\,}$
Theta	${\displaystyle {\frac {\partial V}{\partial t}}}$	${\displaystyle -{\frac {SN'(d_{1})\sigma }{2{\sqrt {T-t}}}}-rKe^{-r(T-t)}N(d_{2})\,}$	${\displaystyle -{\frac {SN'(d_{1})\sigma }{2{\sqrt {T-t}}}}+rKe^{-r(T-t)}N(-d_{2})\,}$
Rho	${\displaystyle {\frac {\partial V}{\partial r}}}$	${\displaystyle K(T-t)e^{-r(T-t)}N(d_{2})\,}$	${\displaystyle -K(T-t)e^{-r(T-t)}N(-d_{2})\,}$

Note that from the formulae, it is clear that the gamma is the same value for calls and puts and so too is the vega the same value for calls and put options. This can be seen directly from put–call parity, since the difference of a put and a call is a forward, which is linear in S and independent of σ (so a forward has zero gamma and zero vega). N' is the standard normal probability density function.

In practice, some sensitivities are usually quoted in scaled-down terms, to match the scale of likely changes in the parameters. For example, rho is often reported divided by 10,000 (1 basis point rate change), vega by 100 (1 vol point change), and theta by 365 or 252 (1 day decay based on either calendar days or trading days per year).

(Vega is not a letter in the Greek alphabet; the name arises from reading the Greek letter ν (nu) as a V.)

The payoff, delta, and gamma for the most common option strategies as determined by the Black-Scholes models can be found on the page Option strategy.

Extensions of the model

The above model can be extended for variable (but deterministic) rates and volatilities. The model may also be used to value European options on instruments paying dividends. In this case, closed-form solutions are available if the dividend is a known proportion of the stock price. American options and options on stocks paying a known cash dividend (in the short term, more realistic than a proportional dividend) are more difficult to value, and a choice of solution techniques is available (for example lattices and grids).

Instruments paying continuous yield dividends

For options on indices, it is reasonable to make the simplifying assumption that dividends are paid continuously, and that the dividend amount is proportional to the level of the index.

The dividend payment paid over the time period ${\displaystyle [t,t+dt]}$ is then modelled as

${\displaystyle qS_{t}\,dt}$

for some constant ${\displaystyle q}$ (the dividend yield).

Under this formulation the arbitrage-free price implied by the Black–Scholes model can be shown to be

${\displaystyle C(S_{t},t)=e^{-r(T-t)}[FN(d_{1})-KN(d_{2})]\,}$

and

${\displaystyle P(S_{t},t)=e^{-r(T-t)}[KN(-d_{2})-FN(-d_{1})]\,}$

where now

${\displaystyle F=S_{t}e^{(r-q)(T-t)}\,}$

is the modified forward price that occurs in the terms ${\displaystyle d_{1},d_{2}}$:

${\displaystyle d_{1}={\frac {1}{\sigma {\sqrt {T-t}}}}\left[\ln \left({\frac {S_{t}}{K}}\right)+(r-q+{\frac {1}{2}}\sigma ^{2})(T-t)\right]}$

and

${\displaystyle d_{2}=d_{1}-\sigma {\sqrt {T-t}}={\frac {1}{\sigma {\sqrt {T-t}}}}\left[\ln \left({\frac {S_{t}}{K}}\right)+(r-q-{\frac {1}{2}}\sigma ^{2})(T-t)\right]}$.[18]

Instruments paying discrete proportional dividends

It is also possible to extend the Black–Scholes framework to options on instruments paying discrete proportional dividends. This is useful when the option is struck on a single stock.

A typical model is to assume that a proportion ${\displaystyle \delta }$ of the stock price is paid out at pre-determined times ${\displaystyle t_{1},t_{2},\ldots }$. The price of the stock is then modelled as

${\displaystyle S_{t}=S_{0}(1-\delta )^{n(t)}e^{ut+\sigma W_{t}}}$

where ${\displaystyle n(t)}$ is the number of dividends that have been paid by time ${\displaystyle t}$.

The price of a call option on such a stock is again

${\displaystyle C(S_{0},T)=e^{-rT}[FN(d_{1})-KN(d_{2})]\,}$

where now

${\displaystyle F=S_{0}(1-\delta )^{n(T)}e^{rT}\,}$

is the forward price for the dividend paying stock.

American options

The problem of finding the price of an American option is related to the optimal stopping problem of finding the time to execute the option. Since the American option can be exercised at any time before the expiration date, the Black–Scholes equation becomes a variational inequality of the form

${\displaystyle {\frac {\partial V}{\partial t}}+{\frac {1}{2}}\sigma ^{2}S^{2}{\frac {\partial ^{2}V}{\partial S^{2}}}+rS{\frac {\partial V}{\partial S}}-rV\leq 0}$[19]

together with ${\displaystyle V(S,t)\geq H(S)}$ where ${\displaystyle H(S)}$ denotes the payoff at stock price ${\displaystyle S}$ and the terminal condition: ${\displaystyle V(S,T)=H(S)}$.

In general this inequality does not have a closed form solution, though an American call with no dividends is equal to a European call and the Roll–Geske–Whaley method provides a solution for an American call with one dividend;[20][21] see also Black's approximation.

Barone-Adesi and Whaley[22] is a further approximation formula. Here, the stochastic differential equation (which is valid for the value of any derivative) is split into two components: the European option value and the early exercise premium. With some assumptions, a quadratic equation that approximates the solution for the latter is then obtained. This solution involves finding the critical value, ${\displaystyle s*}$, such that one is indifferent between early exercise and holding to maturity.[23][24]

Bjerksund and Stensland[25] provide an approximation based on an exercise strategy corresponding to a trigger price. Here, if the underlying asset price is greater than or equal to the trigger price it is optimal to exercise, and the value must equal ${\displaystyle S-X}$, otherwise the option "boils down to: (i) a European up-and-out call option… and (ii) a rebate that is received at the knock-out date if the option is knocked out prior to the maturity date". The formula is readily modified for the valuation of a put option, using put–call parity. This approximation is computationally inexpensive and the method is fast, with evidence indicating that the approximation may be more accurate in pricing long dated options than Barone-Adesi and Whaley.[26]

Perpetual put

Despite the lack of a general analytical solution for American put options, it is possible to derive such a formula for the case of a perpetual option - meaning that the option never expires (i.e., ${\displaystyle T\rightarrow \infty }$).[27] In this case, the time decay of the option is equal to zero, which leads to the Black–Scholes PDE becoming an ODE:

$${\displaystyle {1 \over {2}}\sigma ^{2}S^{2}{d^{2}V \over {dS^{2}}}+(r-q)S{dV \over {dS}}-rV=0}$$

Let ${\displaystyle S_{-}}$ denote the lower exercise boundary, below which is optimal for exercising the option. The boundary conditions are:

$${\displaystyle V(S_{-})=K-S_{-},\quad V_{S}(S_{-})=-1,\quad V(S)\leq K}$$

The solutions to the ODE are a linear combination of any two linearly independent solutions:

$${\displaystyle V(S)=A_{1}S^{\lambda _{1}}+A_{2}S^{\lambda _{2}}}$$

For ${\displaystyle S_{-}\leq S}$, substitution of this solution into the ODE for ${\displaystyle i={1,2}}$ yields:

$${\displaystyle \left[{1 \over {2}}\sigma ^{2}\lambda _{i}(\lambda _{i}-1)+(r-q)\lambda _{i}-r\right]S^{\lambda _{i}}=0}$$

Rearranging the terms in gives:

$${\displaystyle {1 \over {2}}\sigma ^{2}\lambda _{i}^{2}+\left(r-q-{1 \over {2}}\sigma ^{2}\right)\lambda _{i}-r=0}$$

Using the quadratic formula, the solutions for ${\displaystyle \lambda _{i}}$ are:

$${\displaystyle {\begin{aligned}\lambda _{1}&={-\left(r-q-{1 \over {2}}\sigma ^{2}\right)+{\sqrt {\left(r-q-{1 \over {2}}\sigma ^{2}\right)^{2}+2\sigma ^{2}r}} \over {\sigma ^{2}}}\\\lambda _{2}&={-\left(r-q-{1 \over {2}}\sigma ^{2}\right)-{\sqrt {\left(r-q-{1 \over {2}}\sigma ^{2}\right)^{2}+2\sigma ^{2}r}} \over {\sigma ^{2}}}\end{aligned}}}$$

In order to have a finite solution for the perpetual put, since the boundary conditions imply upper and lower finite bounds on the value of the put, it is necessary to set ${\displaystyle A_{1}=0}$, leading to the solution ${\displaystyle V(S)=A_{2}S^{\lambda _{2}}}$. From the first boundary condition, it is known that:

$${\displaystyle V(S_{-})=A_{2}S_{-}^{\lambda _{2}}=K-S_{-}\implies A_{2}={K-S_{-} \over {S_{-}}}}$$

Therefore, the value of the perpetual put becomes:

$${\displaystyle V(S)=(K-S_{-})\left({S \over {S_{-}}}\right)^{\lambda _{2}}}$$

The second boundary condition yields the location of the lower exercise boundary:

$${\displaystyle V_{S}(S_{-})=\lambda _{2}{K-S_{-} \over {S_{-}}}=-1\implies S_{-}={\lambda _{2}K \over {\lambda _{2}-1}}}$$

To conclude, for ${\textstyle S\geq S_{-}={\lambda _{2}K \over {\lambda _{2}-1}}}$, the perpetual American put option is worth:

$${\displaystyle V(S)={K \over {1-\lambda _{2}}}\left({\lambda _{2}-1 \over {\lambda _{2}}}\right)^{\lambda _{2}}\left({S \over {K}}\right)^{\lambda _{2}}}$$

Binary options

By solving the Black–Scholes differential equation, with for boundary condition the Heaviside function, we end up with the pricing of options that pay one unit above some predefined strike price and nothing below.[28]

In fact, the Black–Scholes formula for the price of a vanilla call option (or put option) can be interpreted by decomposing a call option into an asset-or-nothing call option minus a cash-or-nothing call option, and similarly for a put – the binary options are easier to analyze, and correspond to the two terms in the Black–Scholes formula.

Cash-or-nothing call

This pays out one unit of cash if the spot is above the strike at maturity. Its value is given by

${\displaystyle C=e^{-r(T-t)}N(d_{2}).\,}$

Cash-or-nothing put

This pays out one unit of cash if the spot is below the strike at maturity. Its value is given by

${\displaystyle P=e^{-r(T-t)}N(-d_{2}).\,}$

Asset-or-nothing call

This pays out one unit of asset if the spot is above the strike at maturity. Its value is given by

${\displaystyle C=Se^{-q(T-t)}N(d_{1}).\,}$

Asset-or-nothing put

This pays out one unit of asset if the spot is below the strike at maturity. Its value is given by

${\displaystyle P=Se^{-q(T-t)}N(-d_{1}),}$

Foreign exchange

Further information: Foreign exchange derivative

If we denote by S the FOR/DOM exchange rate (i.e., 1 unit of foreign currency is worth S units of domestic currency) we can observe that paying out 1 unit of the domestic currency if the spot at maturity is above or below the strike is exactly like a cash-or nothing call and put respectively. Similarly, paying out 1 unit of the foreign currency if the spot at maturity is above or below the strike is exactly like an asset-or nothing call and put respectively. Hence if we now take ${\displaystyle r_{FOR}}$, the foreign interest rate, ${\displaystyle r_{DOM}}$, the domestic interest rate, and the rest as above, we get the following results.

In case of a digital call (this is a call FOR/put DOM) paying out one unit of the domestic currency we get as present value,

${\displaystyle C=e^{-r_{DOM}T}N(d_{2})\,}$

In case of a digital put (this is a put FOR/call DOM) paying out one unit of the domestic currency we get as present value,

${\displaystyle P=e^{-r_{DOM}T}N(-d_{2})\,}$

While in case of a digital call (this is a call FOR/put DOM) paying out one unit of the foreign currency we get as present value,

${\displaystyle C=Se^{-r_{FOR}T}N(d_{1})\,}$

and in case of a digital put (this is a put FOR/call DOM) paying out one unit of the foreign currency we get as present value,

${\displaystyle P=Se^{-r_{FOR}T}N(-d_{1})\,}$

Skew

In the standard Black–Scholes model, one can interpret the premium of the binary option in the risk-neutral world as the expected value = probability of being in-the-money * unit, discounted to the present value. The Black–Scholes model relies on symmetry of distribution and ignores the skewness of the distribution of the asset. Market makers adjust for such skewness by, instead of using a single standard deviation for the underlying asset ${\displaystyle \sigma }$ across all strikes, incorporating a variable one ${\displaystyle \sigma (K)}$ where volatility depends on strike price, thus incorporating the volatility skew into account. The skew matters because it affects the binary considerably more than the regular options.

A binary call option is, at long expirations, similar to a tight call spread using two vanilla options. One can model the value of a binary cash-or-nothing option, C, at strike K, as an infinitesimally tight spread, where ${\displaystyle C_{v}}$ is a vanilla European call:[29][30]

${\displaystyle C=\lim _{\epsilon \to 0}{\frac {C_{v}(K-\epsilon )-C_{v}(K)}{\epsilon }}}$

Thus, the value of a binary call is the negative of the derivative of the price of a vanilla call with respect to strike price:

${\displaystyle C=-{\frac {dC_{v}}{dK}}}$

When one takes volatility skew into account, ${\displaystyle \sigma }$ is a function of ${\displaystyle K}$:

${\displaystyle C=-{\frac {dC_{v}(K,\sigma (K))}{dK}}=-{\frac {\partial C_{v}}{\partial K}}-{\frac {\partial C_{v}}{\partial \sigma }}{\frac {\partial \sigma }{\partial K}}}$

The first term is equal to the premium of the binary option ignoring skew:

${\displaystyle -{\frac {\partial C_{v}}{\partial K}}=-{\frac {\partial (SN(d_{1})-Ke^{-r(T-t)}N(d_{2}))}{\partial K}}=e^{-r(T-t)}N(d_{2})=C_{\text{no skew}}}$

${\displaystyle {\frac {\partial C_{v}}{\partial \sigma }}}$ is the Vega of the vanilla call; ${\displaystyle {\frac {\partial \sigma }{\partial K}}}$ is sometimes called the "skew slope" or just "skew". If the skew is typically negative, the value of a binary call will be higher when taking skew into account.

${\displaystyle C=C_{\text{no skew}}-{\text{Vega}}_{v}\cdot {\text{Skew}}}$

Relationship to vanilla options' Greeks

Since a binary call is a mathematical derivative of a vanilla call with respect to strike, the price of a binary call has the same shape as the delta of a vanilla call, and the delta of a binary call has the same shape as the gamma of a vanilla call.

Black–Scholes in practice

The normality assumption of the Black–Scholes model does not capture extreme movements such as stock market crashes.

The assumptions of the Black–Scholes model are not all empirically valid. The model is widely employed as a useful approximation to reality, but proper application requires understanding its limitations – blindly following the model exposes the user to unexpected risk.[31] Among the most significant limitations are:

• the underestimation of extreme moves, yielding tail risk, which can be hedged with out-of-the-money options;
• the assumption of instant, cost-less trading, yielding liquidity risk, which is difficult to hedge;
• the assumption of a stationary process, yielding volatility risk, which can be hedged with volatility hedging;
• the assumption of continuous time and continuous trading, yielding gap risk, which can be hedged with Gamma hedging.

In short, while in the Black–Scholes model one can perfectly hedge options by simply Delta hedging, in practice there are many other sources of risk.

Results using the Black–Scholes model differ from real world prices because of simplifying assumptions of the model. One significant limitation is that in reality security prices do not follow a strict stationary log-normal process, nor is the risk-free interest actually known (and is not constant over time). The variance has been observed to be non-constant leading to models such as GARCH to model volatility changes. Pricing discrepancies between empirical and the Black–Scholes model have long been observed in options that are far out-of-the-money, corresponding to extreme price changes; such events would be very rare if returns were lognormally distributed, but are observed much more often in practice.

Nevertheless, Black–Scholes pricing is widely used in practice,[3]^:751[32] because it is:

• easy to calculate
• a useful approximation, particularly when analyzing the direction in which prices move when crossing critical points
• a robust basis for more refined models
• reversible, as the model's original output, price, can be used as an input and one of the other variables solved for; the implied volatility calculated in this way is often used to quote option prices (that is, as a quoting convention).

The first point is self-evidently useful. The others can be further discussed:

Useful approximation: although volatility is not constant, results from the model are often helpful in setting up hedges in the correct proportions to minimize risk. Even when the results are not completely accurate, they serve as a first approximation to which adjustments can be made.

Basis for more refined models: The Black–Scholes model is robust in that it can be adjusted to deal with some of its failures. Rather than considering some parameters (such as volatility or interest rates) as constant, one considers them as variables, and thus added sources of risk. This is reflected in the Greeks (the change in option value for a change in these parameters, or equivalently the partial derivatives with respect to these variables), and hedging these Greeks mitigates the risk caused by the non-constant nature of these parameters. Other defects cannot be mitigated by modifying the model, however, notably tail risk and liquidity risk, and these are instead managed outside the model, chiefly by minimizing these risks and by stress testing.

Explicit modeling: this feature means that, rather than assuming a volatility a priori and computing prices from it, one can use the model to solve for volatility, which gives the implied volatility of an option at given prices, durations and exercise prices. Solving for volatility over a given set of durations and strike prices, one can construct an implied volatility surface. In this application of the Black–Scholes model, a coordinate transformation from the price domain to the volatility domain is obtained. Rather than quoting option prices in terms of dollars per unit (which are hard to compare across strikes, durations and coupon frequencies), option prices can thus be quoted in terms of implied volatility, which leads to trading of volatility in option markets.

The volatility smile

Main article: Volatility smile

One of the attractive features of the Black–Scholes model is that the parameters in the model other than the volatility (the time to maturity, the strike, the risk-free interest rate, and the current underlying price) are unequivocally observable. All other things being equal, an option's theoretical value is a monotonic increasing function of implied volatility.

By computing the implied volatility for traded options with different strikes and maturities, the Black–Scholes model can be tested. If the Black–Scholes model held, then the implied volatility for a particular stock would be the same for all strikes and maturities. In practice, the volatility surface (the 3D graph of implied volatility against strike and maturity) is not flat.

The typical shape of the implied volatility curve for a given maturity depends on the underlying instrument. Equities tend to have skewed curves: compared to at-the-money, implied volatility is substantially higher for low strikes, and slightly lower for high strikes. Currencies tend to have more symmetrical curves, with implied volatility lowest at-the-money, and higher volatilities in both wings. Commodities often have the reverse behavior to equities, with higher implied volatility for higher strikes.

Despite the existence of the volatility smile (and the violation of all the other assumptions of the Black–Scholes model), the Black–Scholes PDE and Black–Scholes formula are still used extensively in practice. A typical approach is to regard the volatility surface as a fact about the market, and use an implied volatility from it in a Black–Scholes valuation model. This has been described as using "the wrong number in the wrong formula to get the right price".[33] This approach also gives usable values for the hedge ratios (the Greeks). Even when more advanced models are used, traders prefer to think in terms of Black–Scholes implied volatility as it allows them to evaluate and compare options of different maturities, strikes, and so on. For a discussion as to the various alternative approaches developed here, see Financial economics § Challenges and criticism.

Valuing bond options

Black–Scholes cannot be applied directly to bond securities because of pull-to-par. As the bond reaches its maturity date, all of the prices involved with the bond become known, thereby decreasing its volatility, and the simple Black–Scholes model does not reflect this process. A large number of extensions to Black–Scholes, beginning with the Black model, have been used to deal with this phenomenon.[34] See Bond option: Valuation.

Interest-rate curve

In practice, interest rates are not constant – they vary by tenor (coupon frequency), giving an interest rate curve which may be interpolated to pick an appropriate rate to use in the Black–Scholes formula. Another consideration is that interest rates vary over time. This volatility may make a significant contribution to the price, especially of long-dated options. This is simply like the interest rate and bond price relationship which is inversely related.

Short stock rate

It is not free to take a short stock position. Similarly, it may be possible to lend out a long stock position for a small fee. In either case, this can be treated as a continuous dividend for the purposes of a Black–Scholes valuation, provided that there is no glaring asymmetry between the short stock borrowing cost and the long stock lending income.

Criticism and comments

Espen Gaarder Haug and Nassim Nicholas Taleb argue that the Black–Scholes model merely recasts existing widely used models in terms of practically impossible "dynamic hedging" rather than "risk", to make them more compatible with mainstream neoclassical economic theory.[35] They also assert that Boness in 1964 had already published a formula that is "actually identical" to the Black–Scholes call option pricing equation.[36] Edward Thorp also claims to have guessed the Black–Scholes formula in 1967 but kept it to himself to make money for his investors.[37] Emanuel Derman and Nassim Taleb have also criticized dynamic hedging and state that a number of researchers had put forth similar models prior to Black and Scholes.[38] In response, Paul Wilmott has defended the model.[32][39]

In his 2008 letter to the shareholders of Berkshire Hathaway, Warren Buffett wrote: "I believe the Black–Scholes formula, even though it is the standard for establishing the dollar liability for options, produces strange results when the long-term variety are being valued... The Black–Scholes formula has approached the status of holy writ in finance ... If the formula is applied to extended time periods, however, it can produce absurd results. In fairness, Black and Scholes almost certainly understood this point well. But their devoted followers may be ignoring whatever caveats the two men attached when they first unveiled the formula."[40]

British mathematician Ian Stewart FRS CMath FIMA—author of the 2012 book entitled In Pursuit of the Unknown: 17 Equations That Changed the World—[41][42] said that Black-Scholes had "underpinned massive economic growth" and the "international financial system was trading derivatives valued at one quadrillion dollars per year" by 2007. He said that the Black-Scholes equation was the "mathematical justification for the trading"—and therefore—"one ingredient in a rich stew of financial irresponsibility, political ineptitude, perverse incentives and lax regulation" that contributed to the financial crisis of 2007–08.[43] He clarified that "the equation itself wasn't the real problem", but its abuse in the financial industry.[43]

See also

• Binomial options model, a discrete numerical method for calculating option prices
• Black model, a variant of the Black–Scholes option pricing model
• Black Shoals, a financial art piece
• Brownian model of financial markets
• Financial mathematics (contains a list of related articles)
• Fuzzy pay-off method for real option valuation
• Heat equation, to which the Black–Scholes PDE can be transformed
• Jump diffusion
• Monte Carlo option model, using simulation in the valuation of options with complicated features
• Real options analysis
• Stochastic volatility

Notes

1. Although the original model assumed no dividends, trivial extensions to the model can accommodate a continuous dividend yield factor.

References

1. "Scholes on merriam-webster.com". Retrieved March 26, 2012.
2. MacKenzie, Donald (2006). An Engine, Not a Camera: How Financial Models Shape Markets. Cambridge, MA: MIT Press. ISBN 0-262-13460-8.
3. Bodie, Zvi; Alex Kane; Alan J. Marcus (2008). Investments (7th ed.). New York: McGraw-Hill/Irwin. ISBN 978-0-07-326967-2.
4. Taleb, 1997. pp. 91 and 110–111.
5. Mandelbrot & Hudson, 2006. pp. 9–10.
6. Mandelbrot & Hudson, 2006. p. 74
7. Mandelbrot & Hudson, 2006. pp. 72–75.
8. Derman, 2004. pp. 143–147.
9. Thorp, 2017. pp. 183–189.
10. https://www.nobelprize.org/nobel_prizes/economic-sciences/laureates/1997/press.html
11. "Nobel Prize Foundation, 1997" (Press release). October 14, 1997. Retrieved March 26, 2012.
12. Black, Fischer; Scholes, Myron (1973). "The Pricing of Options and Corporate Liabilities". Journal of Political Economy. 81 (3): 637–654. doi:10.1086/260062.
13. Merton, Robert (1973). "Theory of Rational Option Pricing". Bell Journal of Economics and Management Science. 4 (1): 141–183. doi:10.2307/3003143. hdl:10338.dmlcz/135817. JSTOR 3003143.
14. Nielsen, Lars Tyge (1993). "Understanding N(d₁) and N(d₂): Risk-Adjusted Probabilities in the Black–Scholes Model" (PDF). Revue Finance (Journal of the French Finance Association). 14 (1): 95–106. Retrieved Dec 8, 2012, earlier circulated as INSEAD Working Paper 92/71/FIN (1992); abstract and link to article, published article.
15. Don Chance (June 3, 2011). "Derivation and Interpretation of the Black–Scholes Model" (PDF). Retrieved March 27, 2012.
16. Hull, John C. (2008). Options, Futures and Other Derivatives (7th ed.). Prentice Hall. ISBN 978-0-13-505283-9.
17. Although with significant algebra; see, for example, Hong-Yi Chen, Cheng-Few Lee and Weikang Shih (2010). Derivations and Applications of Greek Letters: Review and Integration, Handbook of Quantitative Finance and Risk Management, III:491–503.
18. "Extending the Black Scholes formula". finance.bi.no. October 22, 2003. Retrieved July 21, 2017.
19. André Jaun. "The Black–Scholes equation for American options". Retrieved May 5, 2012.
20. Bernt Ødegaard (2003). "Extending the Black Scholes formula". Retrieved May 5, 2012.
21. Don Chance (2008). "Closed-Form American Call Option Pricing: Roll-Geske-Whaley" (PDF). Retrieved May 16, 2012.
22. Giovanni Barone-Adesi & Robert E Whaley (June 1987). "Efficient analytic approximation of American option values". Journal of Finance. 42 (2): 301–20. doi:10.2307/2328254. JSTOR 2328254.
23. Bernt Ødegaard (2003). "A quadratic approximation to American prices due to Barone-Adesi and Whaley". Retrieved June 25, 2012.
24. Don Chance (2008). "Approximation Of American Option Values: Barone-Adesi-Whaley" (PDF). Retrieved June 25, 2012.
25. Petter Bjerksund and Gunnar Stensland, 2002. Closed Form Valuation of American Options
26. American options
27. Crack, Timothy Falcon (2015). Heard on the Street: Quantitative Questions from Wall Street Job Interviews (16th ed.). Timothy Crack. pp. 159–162. ISBN 9780994118257.
28. Hull, John C. (2005). Options, Futures and Other Derivatives. Prentice Hall. ISBN 0-13-149908-4.
29. Breeden, D. T., & Litzenberger, R. H. (1978). Prices of state-contingent claims implicit in option prices. Journal of business, 621-651.
30. Gatheral, J. (2006). The volatility surface: a practitioner's guide (Vol. 357). John Wiley & Sons.
31. Yalincak, Hakan, "Criticism of the Black–Scholes Model: But Why Is It Still Used? (The Answer is Simpler than the Formula)" <<https://ssrn.com/abstract=2115141>>
32. Paul Wilmott (2008): In defence of Black Scholes and Merton Archived 2008-07-24 at the Wayback Machine, Dynamic hedging and further defence of Black–Scholes
33. Riccardo Rebonato (1999). Volatility and correlation in the pricing of equity, FX and interest-rate options. Wiley. ISBN 0-471-89998-4.
34. Kalotay, Andrew (November 1995). "The Problem with Black, Scholes et al" (PDF). Derivatives Strategy.
35. Espen Gaarder Haug and Nassim Nicholas Taleb (2011). Option Traders Use (very) Sophisticated Heuristics, Never the Black–Scholes–Merton Formula. Journal of Economic Behavior and Organization, Vol. 77, No. 2, 2011
36. Boness, A James, 1964, Elements of a theory of stock-option value, Journal of Political Economy, 72, 163–175.
37. A Perspective on Quantitative Finance: Models for Beating the Market, Quantitative Finance Review, 2003. Also see Option Theory Part 1 by Edward Thorpe
38. Emanuel Derman and Nassim Taleb (2005). The illusions of dynamic replication Archived 2008-07-03 at the Wayback Machine, Quantitative Finance, Vol. 5, No. 4, August 2005, 323–326
39. See also: Doriana Ruffinno and Jonathan Treussard (2006). Derman and Taleb's The Illusions of Dynamic Replication: A Comment, WP2006-019, Boston University - Department of Economics.
41. In Pursuit of the Unknown: 17 Equations That Changed the World. New York: Basic Books. 13 March 2012. ISBN 978-1-84668-531-6.
42. Nahin, Paul J. (2012). "In Pursuit of the Unknown: 17 Equations That Changed the World". Physics Today. Review. 65 (9): 52–53. doi:10.1063/PT.3.1720. ISSN 0031-9228.
43. Stewart, Ian (February 12, 2012). "The mathematical equation that caused the banks to crash". The Guardian. The Observer. ISSN 0029-7712. Retrieved April 29, 2020.

Primary references

• Black, Fischer; Myron Scholes (1973). "The Pricing of Options and Corporate Liabilities". Journal of Political Economy. 81 (3): 637–654. doi:10.1086/260062. (Black and Scholes' original paper.)
• Merton, Robert C. (1973). "Theory of Rational Option Pricing". Bell Journal of Economics and Management Science. The RAND Corporation. 4 (1): 141–183. doi:10.2307/3003143. hdl:10338.dmlcz/135817. JSTOR 3003143.
• Hull, John C. (1997). Options, Futures, and Other Derivatives. Prentice Hall. ISBN 0-13-601589-1.

Historical and sociological aspects

• Bernstein, Peter (1992). Capital Ideas: The Improbable Origins of Modern Wall Street. The Free Press. ISBN 0-02-903012-9.
• Derman, Emanuel. "My Life as a Quant" John Wiley & Sons, Inc. 2004. ISBN 0471394203
• MacKenzie, Donald (2003). "An Equation and its Worlds: Bricolage, Exemplars, Disunity and Performativity in Financial Economics" (PDF). Social Studies of Science. 33 (6): 831–868. doi:10.1177/0306312703336002. hdl:20.500.11820/835ab5da-2504-4152-ae5b-139da39595b8. S2CID 15524084.
• MacKenzie, Donald; Yuval Millo (2003). "Constructing a Market, Performing Theory: The Historical Sociology of a Financial Derivatives Exchange". American Journal of Sociology. 109 (1): 107–145. CiteSeerX 10.1.1.461.4099. doi:10.1086/374404.
• MacKenzie, Donald (2006). An Engine, not a Camera: How Financial Models Shape Markets. MIT Press. ISBN 0-262-13460-8.
• Mandelbrot & Hudson, "The (Mis)Behavior of Markets" Basic Books, 2006. ISBN 9780465043552
• Szpiro, George G., Pricing the Future: Finance, Physics, and the 300-Year Journey to the Black–Scholes Equation; A Story of Genius and Discovery (New York: Basic, 2011) 298 pp.
• Taleb, Nassim. "Dynamic Hedging" John Wiley & Sons, Inc. 1997. ISBN 0471152803
• Thorp, Ed. "A Man for all Markets" Random House, 2017. ISBN 9781400067961

Further reading

• Haug, E. G (2007). "Option Pricing and Hedging from Theory to Practice". Derivatives: Models on Models. Wiley. ISBN 978-0-470-01322-9. The book gives a series of historical references supporting the theory that option traders use much more robust hedging and pricing principles than the Black, Scholes and Merton model.
• Triana, Pablo (2009). Lecturing Birds on Flying: Can Mathematical Theories Destroy the Financial Markets?. Wiley. ISBN 978-0-470-40675-5. The book takes a critical look at the Black, Scholes and Merton model.

External links

Discussion of the model

• Ajay Shah. Black, Merton and Scholes: Their work and its consequences. Economic and Political Weekly, XXXII(52):3337–3342, December 1997
• The mathematical equation that caused the banks to crash by Ian Stewart in The Observer, February 12, 2012
• When You Cannot Hedge Continuously: The Corrections to Black–Scholes, Emanuel Derman
• The Skinny On Options TastyTrade Show (archives)

Derivation and solution

• Derivation of the Black–Scholes Equation for Option Value, Prof. Thayer Watkins
• Solution of the Black–Scholes Equation Using the Green's Function, Prof. Dennis Silverman
• Solution via risk neutral pricing or via the PDE approach using Fourier transforms (includes discussion of other option types), Simon Leger
• Step-by-step solution of the Black–Scholes PDE, planetmath.org.
• The Black–Scholes Equation Expository article by mathematician Terence Tao.

Computer implementations

• Black–Scholes in Multiple Languages
• Black–Scholes in Java -moving to link below-
• Black–Scholes in Java
• Chicago Option Pricing Model (Graphing Version)
• Black–Scholes–Merton Implied Volatility Surface Model (Java)
• Online Black–Scholes Calculator

Historical

• Trillion Dollar Bet—Companion Web site to a Nova episode originally broadcast on February 8, 2000. "The film tells the fascinating story of the invention of the Black–Scholes Formula, a mathematical Holy Grail that forever altered the world of finance and earned its creators the 1997 Nobel Prize in Economics."
• BBC Horizon A TV-programme on the so-called Midas formula and the bankruptcy of Long-Term Capital Management (LTCM)
• BBC News Magazine Black–Scholes: The maths formula linked to the financial crash (April 27, 2012 article)