Quantitative Finance Interviews

Chapter 6

Finance

137

It used to be common for candidates with no finance knowledge to get hired into quantitative finance positions. Although this still happens for candidates with specialized knowledge that is in high demand, it’s more likely that you are required, or at least expected, to have a basic grasp of topics in finance. So you should expect to answer some finance questions and be judged on your answers.

Besides classic textbooks,1 there are a few interview books in the market to help you prepare for finance interviews.2 If you want to get prepared for general finance problems, you may want to read a finance interview book to get a feel for what types of questions are asked. The focus of this chapter is more on the intuitions and mathematics behind derivative pricing instead of basic finance knowledge. Derivative problems are popular choices in quantitative interviews—even for divisions that are not directly related to derivative markets—because these problems are complex enough to test your understanding of quantitative finance.

6.1. Option Pricing

Let’s begin with some notations that we will use in the following sections.

TT: maturity date; tt: the current time; τ=T−t\tau=T-t: time to maturity; SS: stock price at time tt; rr: continuous risk-free interest rate; yy: continuous dividend yield; σ\sigma: annualized asset volatility; cc: price of a European call; pp: price of a European put; CC: price of an American call; PP: price of an American put; DD: present value, at tt, of future dividends; KK: strike price; PVPV: present value at tt.

Price direction of options

How do vanilla European/American option prices change when SS, KK, τ\tau, σ\sigma, rr, or DD changes?

Solution: The payoff of a call is max⁡(S−K,0)\max(S-K,0) and the payoff of a put is max⁡(K−S,0)\max(K-S,0). A European option can only be exercised at the expiration time, while an American option can be exercised at any time before maturity. Intuitively we can figure out that the price of a European/American call should decrease when the strike price increases138 since a call with a higher strike has no higher—and sometimes lower—payoff than a call with a lower strike. Using similar analyses, we summarize the effect of changing market conditions on an option’s value in Table 6.1.

The impact of time to maturity on the price of a European call/put is uncertain. If there is a large dividend payoff between two different maturity dates, a European call with shorter maturity that expires before the ex-dividend date may be worth more than a call with longer maturity. For deep in-the-money European puts, the one with shorter maturity is worth more since it can be exercised earlier (time value of the money).

Table 6.1 Impact of SS, KK, τ\tau, σ\sigma, rr, and DD on option prices
VariableEuropean callEuropean putAmerican callAmerican Put
Stock price ↑↑↓↑↓
Strike price ↑↓↑↓↑
Time to maturity ↑??↑↑
Volatility ↑↑↑↑↑
Risk-free rate ↑↑↓↑↓
Dividends ↑↓↑↓↑

↑: increase; ↓: decrease; ?: increase or decrease

It is also worth noting that Table 6.1 assumes that only one factor changes value while all others stay the same, which in practice may not be realistic since some of the factors are related. For example, a large decrease in interest rate often triggers a stock market rally and increases the stock price, which has an opposite effect on option value.

Put-call parity

Put-call parity: c+K−rτ=p+S−Dc+K^{-r\tau}=p+S-D, where the European call option and the European put option have the same underlying security, the same maturity TT and the same strike price KK. Since p≥0p\geq0, we can also derive boundaries for cc, S−D−Ke−rτ≤c≤SS-D-Ke^{-r\tau}\leq c\leq S, from the put-call parity.

For American options, the equality no longer holds and it becomes two inequalities: S−D−K≤C−P≤S−K−rτS-D-K\leq C-P\leq S-K^{-r\tau}.

Can you write down the put-call parity for European options on non-dividend paying stocks and prove it?

139

Solution: The put-call parity for European options on non-dividend paying stocks is c+K−rτ=p+Sc+K^{-r\tau}=p+S. We can treat the left side of the equation as portfolio AA—a call and a zero-coupon bond with face value KK—and the right side as portfolio BB—a put and the underlying stock, which is a protective put. Portfolio AA has payoff max⁡(ST−K,0)+K=max⁡(ST,K)\max(S_T-K,0)+K=\max(S_T,K) at maturity TT; portfolio BB has payoff max⁡(K−ST,0)+ST=max⁡(ST,K)\max(K-S_T,0)+S_T=\max(S_T,K) at TT. Since both portfolios have the same payoff at TT and no payoff between tt and TT, the no-arbitrage argument3 dictates that they must have the same value at tt. Hence, c+K−rτ=p+Sc+K^{-r\tau}=p+S.

If we rearrange the put-call parity equation into c−p=S−K−rτc-p=S-K^{-r\tau}, it will give us different insight. The portfolio on the left side of the equation—long a call and short a put—has the payoff max⁡(ST−K,0)−max⁡(K−ST,0)=ST−K\max(S_T-K,0)-\max(K-S_T,0)=S_T-K, which is the payoff of a forward with delivery price KK. A forward with delivery price KK has present value S−K−rτS-K^{-r\tau}. So we again have the put-call parity c−p=S−K−rτc-p=S-K^{-r\tau}. This expression shows that when the strike price K=SrτK=S^{r\tau} (forward price), a call has the same value as put; when K<SrτK<S^{r\tau}, a call has higher value; and when K>SrτK>S^{r\tau}, a put has higher value.

American v.s. European options

A. Since American options can be exercised at any time before maturity, they are often more valuable than European options with the same characteristics. But when the stock pays no dividend, the theoretical price for an American call and European call should be the same since it is never optimal to exercise the American call. Why should you never exercise an American call on a non-dividend paying stock before maturity?

Solution: There are a number of solutions to this popular problem. We present three arguments for the conclusion.

Argument 1. If you exercises the call option, you will only get the intrinsic value of the call S−KS-K. The price of the American/European call also includes time value, which is positive for a call on a non-dividend paying stock. So the investor is better off selling the option than exercising it before maturity.

In fact, if we rearrange the put-call parity for European options, we have c=S−K−rτ+p=(S−K)+(K−K−rτ)+pc=S-K^{-r\tau}+p=(S-K)+(K-K^{-r\tau})+p. The value of a European call on a non-dividend paying stock includes three components: the first component is the intrinsic value S−KS-K; the second component is the time value of the strike (if you exercise now,140 you pay KK now instead of KK at the maturity date, which is lower in present value); and the third component is the value of the put, which is often considered to be a protection against falling stock price. Clearly the second and the third components are both positive. So the European call should be worth more than its intrinsic value. Considering that the corresponding American call is worth at least as much as the European call, it is worth more than its intrinsic value as well. As a result, it is not optimal to exercise the American call before maturity.

Argument 2. Let’s compare two different strategies. In strategy 1, we exercise the call option4 at time tt (t<Tt<T) and receive cash S−KS-K. Alternatively, we can keep the call, short the underlying stock and lend KK dollars with interest rate rr (the cash proceedings from the short sale, SS, is larger than KK). At the maturity date TT, we exercise the call if it’s in the money, close the short position and close the lending. Table 6.2 shows the cash flow of such a strategy:

It clearly shows that at time tt, we have the same cash flow as exercising the call, S−KS-K. But at time TT, we always have positive cash flow as well. So this strategy is clearly better than exercising the call at time tt. By keeping the call alive, the extra benefit can be realized at maturity.

Table 6.2 Payoff of an alternative strategy without exercising the call
Cash flowttTT
ST≤KS_T\leq KST>KS_T>K
Call option000ST−KS_T-K
Short StockSS−ST-S_T−ST-S_T
Lend K at tt−K-KKerτKe^{r\tau}KerτKe^{r\tau}
TotalS−KS-KKerτ−ST>0Ke^{r\tau}-S_T>0Kerτ−K>0Ke^{r\tau}-K>0

Argument 3. Let’s use a mathematical argument relying on risk-neutral pricing and Jensen’s inequality—if f(X)f(X) is a convex function,5 then E[f(X)]≥f(E[X])E[f(X)]\geq f(E[X]). From Figure 6.1, it’s obvious that the payoff (if exercised when S>KS>K) of a call option C(S)=(S−K)+C(S)=(S-K)^+ is a convex function of stock price with property

C(λS1+(1−λ)S2)≤λC(S1)+(1−λ)C(S2),0<λ<1.C(\lambda S_1+(1-\lambda)S_2)\leq\lambda C(S_1)+(1-\lambda)C(S_2),\quad 0<\lambda<1. 141

Let S1=SS_1=S and S2=0S_2=0, then C(λS)≤λC(S)+(1−λ)C(0)=λC(S)C(\lambda S)\leq\lambda C(S)+(1-\lambda)C(0)=\lambda C(S) since C(0)=0C(0)=0.

Call payoff convexity diagram: horizontal stock-price axis S and vertical call-payoff axis C; the chord through S₂ and S₁ lies above C at λS₁+(1−λ)S₂. Dashed levels mark λC(S₁)+(1−λ)C(S₂) and C(λS₁+(1−λ)S₂).
Figure 6.1 Payoff of a European call option

If the option is exercised at time tt, the payoff at tt is C(St−K)C(S_t-K). If it is not exercised until maturity, the discounted expected payoff (to tt) is E~[e−rτC(ST)]\widetilde E[e^{-r\tau}C(S_T)] under risk-neutral measure. Under risk-neutral probabilities, we also have E~[ST]=Sterτ\widetilde E[S_T]=S_te^{r\tau}.

So E~[e−rτC(ST)]=e−rτE~[C(ST)]≥e−rτC(E~[ST])=e−rτC(erτSt)\widetilde E[e^{-r\tau}C(S_T)]=e^{-r\tau}\widetilde E[C(S_T)]\geq e^{-r\tau}C(\widetilde E[S_T])=e^{-r\tau}C(e^{r\tau}S_t), where the inequality is from Jensen’s inequality.

Let S=erτStS=e^{r\tau}S_t and λ=e−rτ\lambda=e^{-r\tau}, we have C(λS)=C(St)≤e−rτC(erτSt)≤e−rτE~[C(ST)]C(\lambda S)=C(S_t)\leq e^{-r\tau}C(e^{r\tau}S_t)\leq e^{-r\tau}\widetilde E[C(S_T)].

Since the discounted payoff e−rτE~[C(ST)]e^{-r\tau}\widetilde E[C(S_T)] is no less than C(St)C(S_t) for any t≤Tt\leq T under the risk neutral measure, it is never optimal to exercise the option before expiration.

I should point out that the payoff of a put is also a convex function of the stock price. But it is often optimal to exercise an American put on a non-dividend paying stock. The difference is that P(0)=KP(0)=K, so it does not have the property that P(λS)≤λP(S)P(\lambda S)\leq\lambda P(S). In fact, P(λS)≥λP(S)P(\lambda S)\geq\lambda P(S). So the argument for American calls does not apply to American puts.

Similar analysis can also show that early exercise of an American call option for dividend-paying stocks is never optimal except possibly for the time right before an ex-dividend date.

B. A European put option on a non-dividend paying stock with strike price $80 is currently priced at $8 and a put option on the same stock with strike price $90 is priced at $9. Is there an arbitrage opportunity existing in these two options?

142

Solution: In the last problem, we mentioned that the payoff of a put is a convex function in stock price. The price of a put option as a function of the strike price is a convex function as well. Since a put option with strike 0 is worthless, we always have P(0)+λP(K)=λP(K)>P(λK)P(0)+\lambda P(K)=\lambda P(K)>P(\lambda K).

For this specific problem, we should have 8/9×P(90)=8/9×9=8>P(80)8/9\times P(90)=8/9\times9=8>P(80). Since the put option with strike price $80 is currently price at 8, it is overpriced and we should short it. The overall arbitrage portfolio is to short 9 units of put with K=$80K=\$80 and long 8 units of put with K=90K=90. At time 0, the initial cash flow is 0. At the maturity date, we have three possible scenarios:

ST≥90,payoff=0S_T\geq90,\quad\text{payoff}=0 (No put is exercised.)

90>ST≥80,payoff=8×(90−ST)>090>S_T\geq80,\quad\text{payoff}=8\times(90-S_T)>0 (Puts with K=90K=90 are exercised.)

ST<80,payoff=8×(90−ST)−9×(80−ST)=ST>0S_T<80,\quad\text{payoff}=8\times(90-S_T)-9\times(80-S_T)=S_T>0 (All puts are exercised.)

The final payoff ≥0\geq0 with positive probability that payoff >0>0. So it is clearly an arbitrage opportunity.

Black-Scholes-Merton differential equation

Can you write down the Black-Scholes-Merton differential equation and briefly explain how to derive it?

Solution: If the evolution of the stock price is a geometric Brownian motion, dS=μSdt+σSdW(t)dS=\mu Sdt+\sigma SdW(t), and the derivative V=V(S,t)V=V(S,t) is a function of SS and tt, then applying Ito’s lemma yields:

dV=(∂V∂t+μS∂V∂S+12σ2S2∂2V∂S2)dt+σS∂V∂SdW(t),dV=\left(\frac{\partial V}{\partial t}+\mu S\frac{\partial V}{\partial S}+\frac12\sigma^2S^2\frac{\partial^2V}{\partial S^2}\right)dt+\sigma S\frac{\partial V}{\partial S}dW(t),

where W(t)W(t) is a Brownian motion.

The Black-Scholes-Merton differential equation is a partial differential equation that should be satisfied by VV:

∂V∂t+rS∂V∂S+12σ2S2∂2V∂S2=rV.\frac{\partial V}{\partial t}+rS\frac{\partial V}{\partial S}+\frac12\sigma^2S^2\frac{\partial^2V}{\partial S^2}=rV.

To derive the Black-Scholes-Merton differential equation, we build a portfolio with two components: long one unit of the derivative and short ∂V∂S\dfrac{\partial V}{\partial S} unit of the underlying stock. Then the portfolio has value Π=V−∂V∂SS\Pi=V-\dfrac{\partial V}{\partial S}S and the change of Π\Pi follows equation

143 dΠ=dV−∂V∂SdS=(∂V∂t+μS∂V∂S+12σ2S2∂2V∂S2)dt+σS∂V∂SdW(t)−∂V∂S(μSdt+σSdW(t))=(∂V∂t+12σ2S2∂2V∂S2)dt.\begin{aligned} d\Pi&=dV-\frac{\partial V}{\partial S}dS\\ &=\left(\frac{\partial V}{\partial t}+\mu S\frac{\partial V}{\partial S}+\frac12\sigma^2S^2\frac{\partial^2V}{\partial S^2}\right)dt+\sigma S\frac{\partial V}{\partial S}dW(t)-\frac{\partial V}{\partial S}(\mu Sdt+\sigma SdW(t))\\ &=\left(\frac{\partial V}{\partial t}+\frac12\sigma^2S^2\frac{\partial^2V}{\partial S^2}\right)dt. \end{aligned}

It is apparent that this portfolio is risk-free since it has no diffusion term. It should have risk-free rate of return as well: dΠ=r(V−∂V∂SS)dtd\Pi=r(V-\dfrac{\partial V}{\partial S}S)dt. Combining these results we have

(∂V∂t+12σ2S2∂2V∂S2)dt=r(V−∂V∂SS)dt⇒∂V∂t+rS∂V∂S+12σ2S2∂2V∂S2=rV,\left(\frac{\partial V}{\partial t}+\frac12\sigma^2S^2\frac{\partial^2V}{\partial S^2}\right)dt=r\left(V-\frac{\partial V}{\partial S}S\right)dt\Rightarrow\frac{\partial V}{\partial t}+rS\frac{\partial V}{\partial S}+\frac12\sigma^2S^2\frac{\partial^2V}{\partial S^2}=rV,

which is the Black-Scholes-Merton differential equation.

The Black-Scholes-Merton differential equation is a special case of the discounted Feynman-Kac theorem. The discounted Feynman-Kac theorem builds the bridge between stochastic differential equations and partial differential equations and applies to all Ito processes in general:

Let XX be an Ito process given by equation dX(t)=β(t,X)dt+γ(t,X)dW(t)dX(t)=\beta(t,X)dt+\gamma(t,X)dW(t) and f(x)f(x) be a function of XX. Define function V(t,x)=E[e−r(T−t)f(XT)∣Xt=x]V(t,x)=E[e^{-r(T-t)}f(X_T)\mid X_t=x], then V(t,x)V(t,x) is a martingale process that satisfies the partial differential equation

∂V∂t+β(t,x)∂V∂x+12γ2(t,x)∂2V∂x2=rV(t,x)\frac{\partial V}{\partial t}+\beta(t,x)\frac{\partial V}{\partial x}+\frac12\gamma^2(t,x)\frac{\partial^2V}{\partial x^2}=rV(t,x)

and boundary condition V(T,x)=f(x)V(T,x)=f(x) for all xx.

Under risk-neutral measure, dS=rSdt+σSdW(t)dS=rSdt+\sigma SdW(t). Let S=XS=X, β(t,X)=rS\beta(t,X)=rS and γ(t,X)=σS\gamma(t,X)=\sigma S, then the discounted Feynman-Kac equation becomes the Black-Scholes-Merton differential equation ∂V∂t+rS∂V∂S+12σ2S2∂2V∂S2=rV\dfrac{\partial V}{\partial t}+rS\dfrac{\partial V}{\partial S}+\dfrac12\sigma^2S^2\dfrac{\partial^2V}{\partial S^2}=rV.

Black-Scholes formula

The Black-Scholes formula for European calls and puts with continuous dividend yield yy is:

c=Se−yτN(d1)−Ke−rτN(d2)andp=Ke−rτN(−d2)−Se−yτN(−d1),c=Se^{-y\tau}N(d_1)-Ke^{-r\tau}N(d_2)\quad\text{and}\quad p=Ke^{-r\tau}N(-d_2)-Se^{-y\tau}N(-d_1), 144

where

d1=ln⁡(Se−yτ/K)+(r+σ2/2)τστ=ln⁡(S/K)+(r−y+σ2/2)τστd2=ln⁡(S/K)+(r−y−σ2/2)τστ=d1−στ.\begin{aligned} d_1&=\frac{\ln(Se^{-y\tau}/K)+(r+\sigma^2/2)\tau}{\sigma\sqrt\tau}=\frac{\ln(S/K)+(r-y+\sigma^2/2)\tau}{\sigma\sqrt\tau}\\ d_2&=\frac{\ln(S/K)+(r-y-\sigma^2/2)\tau}{\sigma\sqrt\tau}=d_1-\sigma\sqrt\tau. \end{aligned}

N(x)N(x) is the cdf of the standard normal distribution and N′(x)N'(x) is the pdf of the standard normal distribution: N(x)=∫−∞x12πe−y2/2dyN(x)=\displaystyle\int_{-\infty}^x\frac1{\sqrt{2\pi}}e^{-y^2/2}dy and N′(x)=12πe−x2/2N'(x)=\dfrac1{\sqrt{2\pi}}e^{-x^2/2}.

If the underlying asset is a futures contract, then yield y=ry=r. If the underlying asset is a foreign currency, then yield y=rfy=r_f, where rfr_f is the foreign risk-free interest rate.

A. What are the assumptions behind the Black-Scholes formula?

Solution: The original Black-Scholes formula for European calls and puts consists of the equations c=SN(d1)−Ke−rτN(d2)c=SN(d_1)-Ke^{-r\tau}N(d_2) and p=Ke−rτN(−d2)−SN(−d1)p=Ke^{-r\tau}N(-d_2)-SN(-d_1), which require the following assumptions:

  1. The stock pays no dividends.
  2. The risk-free interest rate is constant and known.
  3. The stock price follows a geometric Brownian motion with constant drift μ\mu and volatility σ\sigma: dS=μSdt+σSdW(t)dS=\mu Sdt+\sigma SdW(t).
  4. There are no transaction costs or taxes; the proceeds of short selling can be fully invested.
  5. All securities are perfectly divisible.
  6. There are no risk-free arbitrage opportunities.

B. How can you derive the Black-Scholes formula for a European call on a non-dividend paying stock using risk-neutral probability measure?

Solution: The Black-Scholes formula for a European call on a non-dividend paying stock is

c=SN(d1)−Ke−rτN(d2),where d1=ln⁡(S/K)+(r+σ2/2)τστandd2=d1−στ.c=SN(d_1)-Ke^{-r\tau}N(d_2),\quad\text{where }d_1=\frac{\ln(S/K)+(r+\sigma^2/2)\tau}{\sigma\sqrt\tau}\quad\text{and}\quad d_2=d_1-\sigma\sqrt\tau. 145

Under the risk-neutral probability measure, the drift of stock price becomes the risk-free interest rate r(t)r(t): dS=r(t)Sdt+σSdW(t)dS=r(t)Sdt+\sigma SdW(t). Risk-neutral measure allows the option to be priced as the discounted value of its expected payoff with the risk-free interest rate:

V(t)=E[e−∫tTr(u)duV(T)∣S(t)],0≤t≤T,V(t)=E\left[e^{-\int_t^T r(u)du}V(T)\mathrel{\Big|}S(t)\right],\quad0\leq t\leq T,

where V(T)V(T) is the payoff at maturity TT.

When rr is constant, the formula can be further simplified as V(t)=e−rτE[V(T)∣S(t)]V(t)=e^{-r\tau}E[V(T)\mid S(t)]. Under risk-neutral probabilities, dS=rSdt+σSdW(t)dS=rSdt+\sigma SdW(t). Applying Ito’s lemma, we get

d(ln⁡(S))=(r−σ2/2)dt+σdW(t)⇒ln⁡ST∼N(ln⁡S+(r−σ2/2)τ,σ2τ).d(\ln(S))=(r-\sigma^2/2)dt+\sigma dW(t)\Rightarrow\ln S_T\sim N\left(\ln S+(r-\sigma^2/2)\tau,\sigma^2\tau\right).

So ST=Se(r−σ2/2)τ+στεS_T=Se^{(r-\sigma^2/2)\tau+\sigma\sqrt\tau\varepsilon}, where ε∼N(0,1)\varepsilon\sim N(0,1). For a European option, we have

V(T)={Se(r−σ2/2)τ+στε−K,if Se(r−σ2/2)τ+στε>K0,otherwiseV(T)=\begin{cases}Se^{(r-\sigma^2/2)\tau+\sigma\sqrt\tau\varepsilon}-K,&\text{if }Se^{(r-\sigma^2/2)\tau+\sigma\sqrt\tau\varepsilon}>K\\0,&\text{otherwise}\end{cases} Se(r−σ2/2)τ+στε>K⇒ε>ln⁡(K/S)−(r−σ2/2)τστ=−d2andSe^{(r-\sigma^2/2)\tau+\sigma\sqrt\tau\varepsilon}>K\Rightarrow\varepsilon>\frac{\ln(K/S)-(r-\sigma^2/2)\tau}{\sigma\sqrt\tau}=-d_2\quad\text{and} E[V(T)∣S]=E[max⁡(ST−K,0)∣S]=∫−d2∞(Se(r−σ2/2)τ+στε−K)12πe−ε2/2dε=Serτ∫−d2∞12πe−(ε−τσ)2/2dε−K∫−d2∞12πe−ε2/2dε.\begin{aligned} E[V(T)\mid S]&=E[\max(S_T-K,0)\mid S]=\int_{-d_2}^\infty\left(Se^{(r-\sigma^2/2)\tau+\sigma\sqrt\tau\varepsilon}-K\right)\frac1{\sqrt{2\pi}}e^{-\varepsilon^2/2}d\varepsilon\\ &=Se^{r\tau}\int_{-d_2}^\infty\frac1{\sqrt{2\pi}}e^{-(\varepsilon-\sqrt\tau\sigma)^2/2}d\varepsilon-K\int_{-d_2}^\infty\frac1{\sqrt{2\pi}}e^{-\varepsilon^2/2}d\varepsilon. \end{aligned}

Let ε~=ε−στ\widetilde\varepsilon=\varepsilon-\sigma\sqrt\tau, then dε=dε~d\varepsilon=d\widetilde\varepsilon, ε=−d2⇒ε~=−d2−στ=−d1\varepsilon=-d_2\Rightarrow\widetilde\varepsilon=-d_2-\sigma\sqrt\tau=-d_1 and we have

Serτ∫−d2∞12πe−(ε−τσ)2/2dε=Serτ∫−d1∞12πe−ε~2/2dε~=SerτN(d1),Se^{r\tau}\int_{-d_2}^\infty\frac1{\sqrt{2\pi}}e^{-(\varepsilon-\sqrt\tau\sigma)^2/2}d\varepsilon=Se^{r\tau}\int_{-d_1}^\infty\frac1{\sqrt{2\pi}}e^{-\widetilde\varepsilon^2/2}d\widetilde\varepsilon=Se^{r\tau}N(d_1), K∫−d2∞12πe−ε2/2dε=K(1−N(−d2))=KN(d2)K\int_{-d_2}^\infty\frac1{\sqrt{2\pi}}e^{-\varepsilon^2/2}d\varepsilon=K(1-N(-d_2))=KN(d_2) ∴E[V(T)]=SerτN(d1)−KN(d2)andV(t)=e−rτE[V(T)]=SN(d1)−Ke−rτN(d2).\therefore E[V(T)]=Se^{r\tau}N(d_1)-KN(d_2)\quad\text{and}\quad V(t)=e^{-r\tau}E[V(T)]=SN(d_1)-Ke^{-r\tau}N(d_2).

From the derivation process, it is also obvious that 1−N(−d2)=N(d2)1-N(-d_2)=N(d_2) is the risk-neutral probability that the call option finishes in the money.

C. How do you derive the Black-Scholes formula for a European call option on a non-dividend paying stock by solving the Black-Scholes-Merton differential equation?

146

Solution: You can skip this problem if you don’t have background in partial differential equations (PDE). One approach to solving the problem is to convert the Black-Scholes-Merton differential equation to a heat equation and then apply the boundary conditions to the heat equation to derive the Black-Scholes formula.

Let y=ln⁡Sy=\ln S (S=eyS=e^y) and τ~=T−t\widetilde\tau=T-t, then ∂V∂t=−∂V∂τ~\dfrac{\partial V}{\partial t}=-\dfrac{\partial V}{\partial\widetilde\tau}, ∂V∂S=∂V∂ydydS=1S∂V∂y\dfrac{\partial V}{\partial S}=\dfrac{\partial V}{\partial y}\dfrac{dy}{dS}=\dfrac1S\dfrac{\partial V}{\partial y} and

∂2V∂S2=∂V∂S(∂V∂S)=∂V∂S(1S∂V∂y)=−1S2∂V∂y+1S∂V∂S(∂V∂y)=−1S2∂V∂y+1S2∂2V∂y2.\frac{\partial^2V}{\partial S^2}=\frac{\partial V}{\partial S}\left(\frac{\partial V}{\partial S}\right)=\frac{\partial V}{\partial S}\left(\frac1S\frac{\partial V}{\partial y}\right)=\frac{-1}{S^2}\frac{\partial V}{\partial y}+\frac1S\frac{\partial V}{\partial S}\left(\frac{\partial V}{\partial y}\right)=\frac{-1}{S^2}\frac{\partial V}{\partial y}+\frac1{S^2}\frac{\partial^2V}{\partial y^2}.

6

The Black-Scholes-Merton differential equation ∂V∂t+rS∂V∂S+12σ2S2∂2V∂S2−rV=0\dfrac{\partial V}{\partial t}+rS\dfrac{\partial V}{\partial S}+\dfrac12\sigma^2S^2\dfrac{\partial^2V}{\partial S^2}-rV=0 can be converted to −∂V∂τ~+(r−12σ2)∂V∂y+12σ2∂2V∂y2−rV=0-\dfrac{\partial V}{\partial\widetilde\tau}+\left(r-\dfrac12\sigma^2\right)\dfrac{\partial V}{\partial y}+\dfrac12\sigma^2\dfrac{\partial^2V}{\partial y^2}-rV=0.

Let u=erτ~Vu=e^{r\widetilde\tau}V, the equation becomes −∂u∂τ~+(r−12σ2)∂u∂y+12σ2∂2u∂y2=0-\dfrac{\partial u}{\partial\widetilde\tau}+\left(r-\dfrac12\sigma^2\right)\dfrac{\partial u}{\partial y}+\dfrac12\sigma^2\dfrac{\partial^2u}{\partial y^2}=0.

Finally, let x=y+(r−12σ2)τ~=ln⁡S+(r−12σ2)τ~x=y+\left(r-\dfrac12\sigma^2\right)\widetilde\tau=\ln S+\left(r-\dfrac12\sigma^2\right)\widetilde\tau and τ=τ~\tau=\widetilde\tau, then ∂u∂y=∂u∂x\dfrac{\partial u}{\partial y}=\dfrac{\partial u}{\partial x} and ∂u∂τ~=∂u∂τ+(r−12σ2)∂u∂x\dfrac{\partial u}{\partial\widetilde\tau}=\dfrac{\partial u}{\partial\tau}+\left(r-\dfrac12\sigma^2\right)\dfrac{\partial u}{\partial x}, which transforms the equation to

−∂u∂τ−(r−12σ2)∂u∂x+(r−12σ2)∂u∂x+12σ2∂2u∂x2=0⇒∂u∂τ=12σ2∂2u∂x2.-\frac{\partial u}{\partial\tau}-\left(r-\frac12\sigma^2\right)\frac{\partial u}{\partial x}+\left(r-\frac12\sigma^2\right)\frac{\partial u}{\partial x}+\frac12\sigma^2\frac{\partial^2u}{\partial x^2}=0\Rightarrow\frac{\partial u}{\partial\tau}=\frac12\sigma^2\frac{\partial^2u}{\partial x^2}.

So the original equation becomes a heat/diffusion equation ∂u∂τ=12σ2∂2u∂x2\dfrac{\partial u}{\partial\tau}=\dfrac12\sigma^2\dfrac{\partial^2u}{\partial x^2}. For heat equation ∂u∂τ=12σ2∂2u∂x2\dfrac{\partial u}{\partial\tau}=\dfrac12\sigma^2\dfrac{\partial^2u}{\partial x^2}, where u=u(x,τ)u=u(x,\tau) is a function of time τ\tau and space variable xx, with boundary condition u(x,0)=u0(x)u(x,0)=u_0(x), the solution is

u(x,τ)=12πτσ∫−∞∞u0(ψ)exp⁡(−(x−ψ)22σ2τ)dψ.u(x,\tau)=\frac1{\sqrt{2\pi\tau}\sigma}\int_{-\infty}^\infty u_0(\psi)\exp\left(-\frac{(x-\psi)^2}{2\sigma^2\tau}\right)d\psi.

7

147

For European calls, the boundary condition is u0(ST)=max⁡(ST−K,0)u_0(S_T)=\max(S_T-K,0).

S=exp⁡(x−(r−0.5σ2)τ)S=\exp\left(x-(r-0.5\sigma^2)\tau\right). When x=ψx=\psi and τ=0\tau=0, ST=eψS_T=e^\psi.

u(S,τ)=u(x,τ)=12πτσ∫−∞∞max⁡(eψ−K,0)exp⁡(−(x−ψ)22σ2τ)dψ=12πτσ∫ln⁡K∞(eψ−K)exp⁡(−(x−ψ)22σ2τ)dψ.\begin{aligned} u(S,\tau)=u(x,\tau)&=\frac1{\sqrt{2\pi\tau}\sigma}\int_{-\infty}^\infty\max(e^\psi-K,0)\exp\left(-\frac{(x-\psi)^2}{2\sigma^2\tau}\right)d\psi\\ &=\frac1{\sqrt{2\pi\tau}\sigma}\int_{\ln K}^\infty(e^\psi-K)\exp\left(-\frac{(x-\psi)^2}{2\sigma^2\tau}\right)d\psi. \end{aligned}

Let ε=ψ−xστ\varepsilon=\dfrac{\psi-x}{\sigma\sqrt\tau}, then dε=dψστd\varepsilon=\dfrac{d\psi}{\sigma\sqrt\tau}, eψ=ex+εστe^\psi=e^{x+\varepsilon\sigma\sqrt\tau}, exp⁡(−(x−ψ)22σ2τ)=e−ε2/2\exp\left(-\dfrac{(x-\psi)^2}{2\sigma^2\tau}\right)=e^{-\varepsilon^2/2} and when ψ=ln⁡K\psi=\ln K, ε=ln⁡(K/S)−(r−σ2/2)τστ=−d2\varepsilon=\dfrac{\ln(K/S)-(r-\sigma^2/2)\tau}{\sigma\sqrt\tau}=-d_2.

∴u(S,τ)=∫−d2∞(Se(r−σ2/2)τ+στε−K)12πe−ε2/2dε.\therefore u(S,\tau)=\int_{-d_2}^\infty\left(Se^{(r-\sigma^2/2)\tau+\sigma\sqrt\tau\varepsilon}-K\right)\frac1{\sqrt{2\pi}}e^{-\varepsilon^2/2}d\varepsilon.

Now, it’s clear that the equation for u(S,τ)u(S,\tau) is exactly the same as the equation for E[V(T)∣S]E[V(T)\mid S] in question BB. Hence, we have V(S,t)=e−rτu(S,τ)=SN(d1)−Ke−rτN(d2)V(S,t)=e^{-r\tau}u(S,\tau)=SN(d_1)-Ke^{-r\tau}N(d_2) as well.

D. Assume zero interest rate and a stock with current price at $1 that pays no dividend. When the price hits level $HH (H>1H>1) for the first time you can exercise the option and receive $1. What is this option worth to you today?

Solution: First let’s use a brute-force approach to solve the problem by assuming that the stock price follows a geometric Brownian motion under risk-neutral measure: dS=rSdt+σSdW(t)dS=rSdt+\sigma SdW(t). Since r=0r=0, dS=σSdW(t)⇒d(ln⁡S)=−12σ2dt+σdW(t)dS=\sigma SdW(t)\Rightarrow d(\ln S)=-\frac12\sigma^2dt+\sigma dW(t). When t=0t=0, we have S0=1⇒ln⁡(S0)=0S_0=1\Rightarrow\ln(S_0)=0.

148

Hence, ln⁡S=−12σ2t+σW(t)⇒ln⁡S+12σ2tσ=W(t)\ln S=-\frac12\sigma^2t+\sigma W(t)\Rightarrow\dfrac{\ln S+\frac12\sigma^2t}{\sigma}=W(t) is a Brownian motion.

Whenever SS reaches $HH, the payoff is $1. Because the interest rate is 0, the discounted payoff is also $1 under risk-neutral measure. So the value of the option is the probability that SS ever reaches $HH, which is equivalent to the probability that ln⁡S\ln S ever reaches ln⁡H\ln H. Again we can apply the exponential martingale Z(t)=exp⁡{λW(t)−12λ2t}Z(t)=\exp\{\lambda W(t)-\frac12\lambda^2t\} as we did in Chapter 5:

E[Z(t)]=E[exp⁡{λln⁡S+12σ2tσ−12λ2t}]=1.E[Z(t)]=E\left[\exp\left\{\lambda\frac{\ln S+\frac12\sigma^2t}{\sigma}-\frac12\lambda^2t\right\}\right]=1.

To remove the terms including time tt, we can set λ=σ\lambda=\sigma and the equation becomes E[exp⁡(ln⁡S)]=1E[\exp(\ln S)]=1. The Let PP be the probability that ln⁡S\ln S ever reaches ln⁡H\ln H (using −∞-\infty as the negative boundary for stopping time), we have

Pexp⁡(ln⁡H)+(1−P)exp⁡(−∞)=P×H=1⇒P=1/H.P\exp(\ln H)+(1-P)\exp(-\infty)=P\times H=1\Rightarrow P=1/H.

So the probability that SS ever reaches $HH is 1/H1/H and the price of the option should be $1/H1/H. Notice that SS is a martingale under the risk-neutral measure;8 but ln⁡S\ln S has a negative drift. The reason is that ln⁡S\ln S follows a (symmetrical) normal distribution, but SS itself follows a lognormal distribution, which is positively skewed. As T→∞T\to\infty, although the expected value of STS_T is 1, the probability that ST≥1S_T\geq1 actually approaches 0.

It is simpler to use a no-arbitrage argument to derive the price. In order to pay $1 when the stock price hits $HH, we need to buy 1/H1/H shares of the stock (at $1/H1/H). So the option should be worth no more than $1/H1/H. Yet if the option price CC is less than $1/H1/H (C<1/H⇒CH<1C<1/H\Rightarrow CH<1), we can buy an option by borrowing CC shares of the stock. The initial investment is 0. Once the stock price hits $HH, we will excise the option and return the stock by buying CC shares at price $HH, which gives payoff 1−CH>01-CH>0. That means we have no initial investment, yet we have possible positive future payoff, which is contradictory to the no arbitrage argument. So the price cannot be less than $1/H1/H. Hence, the price is exactly $1/H1/H.

E. Assume a non-dividend paying stock follows a geometric Brownian motion. What is the value of a contract that at maturity TT pays the inverse of the stock price observed at the maturity?

149

Solution: Under risk-neutral measure dS=rSdt+σSdW(t)dS=rSdt+\sigma SdW(t). Apply Ito’s lemma to V=1SV=\dfrac1S:

dV=(∂V∂SrS+∂V∂t+12∂2V∂S2σ2S2)dt+∂V∂SσSdW(t)=(−1S2rS+0+122S3σ2S2)dt−1S2σSdW(t)=(−r+σ2)Vdt−σVdW(t).\begin{aligned} dV&=\left(\frac{\partial V}{\partial S}rS+\frac{\partial V}{\partial t}+\frac12\frac{\partial^2V}{\partial S^2}\sigma^2S^2\right)dt+\frac{\partial V}{\partial S}\sigma SdW(t)\\ &=\left(-\frac1{S^2}rS+0+\frac12\frac2{S^3}\sigma^2S^2\right)dt-\frac1{S^2}\sigma SdW(t)=(-r+\sigma^2)Vdt-\sigma VdW(t). \end{aligned}

So VV follows a geometric Brownian motion as well and we can apply Ito’s lemma to ln⁡V\ln V:

d(ln⁡V)=(VV(−r+σ2)+0−12V2V2σ2)dt+VVσdW(t)=(−r+12σ2)dt−σdW(t).d(\ln V)=\left(\frac VV(-r+\sigma^2)+0-\frac12\frac{V^2}{V^2}\sigma^2\right)dt+\frac VV\sigma dW(t)=\left(-r+\frac12\sigma^2\right)dt-\sigma dW(t).

Hence, ln⁡(VT)∼ln⁡(Vt)+N((−r+12σ2)τ,σ2τ)\ln(V_T)\sim\ln(V_t)+N\left((-r+\frac12\sigma^2)\tau,\sigma^2\tau\right) and E[VT]=E[eln⁡VT]=1Ste−rτ+σ2τE[V_T]=E[e^{\ln V_T}]=\dfrac1{S_t}e^{-r\tau+\sigma^2\tau}.

Discounting the payoff by e−rτe^{-r\tau}, we have V=e−rτE[VT]=1Ste−2rτ+σ2τV=e^{-r\tau}E[V_T]=\dfrac1{S_t}e^{-2r\tau+\sigma^2\tau}.