FINM3405 Revision

lecture 4-8 Options

  1. Definitions and Concepts

  2. Moneyness (ITM, ATM, OTM).

  3. Payoffs vs profits (including the premium).

  4. Main contracts and markets, contract specifications.

  5. Pricing bounds and put-call-parity.

  6. American call premium = European call premium when no dividends.

  7. Time value and intrinsic value.

  8. European option premiums/prices via the Black-Scholes model:

  1. Black-Scholes assumption of geometric Brownian motion and consequent log-normal distribution of the underlying asset’s price, or normal distribution of the asset’s returns.

  2. Simulating geometric Brownian motion in preparation for the Monte Carlo approach to derivative security pricing.

  3. Heuristic (non-rigorous) discussion of risk-neutral pricing.

  4. Heuristic (hand-waving) interpretation of the factors affecting option prices, and quantification of this via the Black-Scholes model Greeks.

  5. Using delta ∆ and gamma Γ to predict small changes in option prices due to small changes in the price of the underlying asset (in preparation for delta and delta-gamma hedging).

  6. More detailed discussion of theta θ and associated concept of time decay, and how it relates to moneyness.

  7. Scholes Hedging

  1. Implied volatility
  1. Trading strategies and their payoff and profit/loss diagrams.

  2. Binomial and Monte Carlo numerical option pricing methods:

  1. 1-period binomial model for European option price and delta.

Aside - Notation

Black Scholes


1. Definitions and Concepts

Recall that there is two types of plain vanilla European options:

2. Moneyness

At time tt an option is:

3 Payoffs and Profits

3.1 Payoffs

At any time tt, an option’s intrinsic or exercise value (IV) is

call IVt=max{0,StK}IV_{t} = \max \{ 0, S_{t} − K \}

and put IVt=max{0,KSt}IV_{t} = \max \{ 0, K − S_{t} \}

3.2 Profits

and the writer’s (short position) profits are the negative of these:

4. Main contracts and markets, contract specifications.

5. Pricing bounds and put-call parity

5.1 put-call parity

Developing mathematical models for pricing options (calculating the fair value of the option premium), such as the Black-Scholes option pricing model, is a major part of the theory and practice of options.

Form a portfolio long 1 European call and short 1 European put, both with the same underlying, strike and expiry.

This portfolio’s current price is C − P and its payoff at expiry is

portfolio payoff=max{0,STK}long callmax{0,KST}short put\text{portfolio payoff} = \underbrace{\max\{0, S_T - K\}}_{\text{long call}} - \underbrace{\max\{0, K - S_T\}}_{\text{short put}}

5.2 Pricing bounds

  1. First note that American options are worth at least as much as European options over the same underlying asset and with the same strike and expiry:

0CEuCAm0 \le C_{Eu} \le C_{Am}

and

0PEuPAm0 \le P_{Eu} \le P_{Am}

American options can be exercised any time up to and including expiry

5.2.1 Option Pricing Bounds - American

Lower Pricing bounds - American 2. American options are worth at least their intrinsic (exercise) value:

max{0,SK}CAm\max \{0, S − K \} ≤ C_{Am}

and

max{0,KS}PAm\max \{0, K − S \} ≤ P_{Am}

because they can be exercised immediately

Upper Pricing bounds - American 3. American calls can never be worth more than the underlying spot price, and American puts can never be worth more than the strike:

CAmSC_{Am} ≤ S

and

PAmKP_{Am} ≤ K

Combined

Combining the lower and upper bounds from the previous two slides leads to the important pricing bounds for American options:

max{0,SK}CAmS \max \{ 0, S − K \} ≤ C^{Am} ≤ S

And

max{0,KS}PAmK \max \{ 0, K − S \} ≤ P^{Am} ≤ K

5.2.2 Option Pricing Bounds - European

Upper Pricing bounds - European

0PEuKerT0 \le P^{Eu} \le Ke^{-rT}

because a European put can only be exercised at expiry, where it is known with certainty that their maximum payoff at expiry is KK.

From this, deep in-the-money European puts can have negative time value (their premium can be less than their intrinsic value).

Combined

European calls

max{0,eqTSerTK}CEuS\max \{ 0, e^{-qT}S-e^{-rT}K \} \le C^{Eu} \le S

European puts

max(0,erTKeqTS)PEuKerT\max(0, e^{-rT}K - e^{-qT}S) \le P^{Eu} \le Ke^{-rT}

6. No early exercise of American calls

early exercise is never optimal for an American call option on a non-dividend-paying stock

max(0,SerTK)CEuCAm\max (0, S − e^{−rT}K) ≤ C^{Eu} ≤ C^{Am}

max(0,SK)<max(0,SerTK)CrTCAm\max(0, S − K) < \max (0, S − e^{−rT}K) ≤ C^{-rT} ≤ C^{Am}

But max(0,SK)\max(0, S − K) is just the call’s intrinsic value. So if there’s no dividends, call options will always have strictly positive time value.

7. Time value and intrinsic value.

From above, we noticed that call options on non-dividend-paying stocks have a premium that is strictly larger than the option’s intrinsic value.

time value = premium − intrinsic value\text{time value = premium − intrinsic value}

8. European option premiums/prices via the Black-Scholes model

The Black-Scholes European option pricing model is a mathematical equation for pricing plain vanilla European call and put options.

8.1 European option premiums on non-dividend paying assets

Black-Scholes European call option pricing model is

C=SN(d1)KerTN(d2)C = S\mathcal{N} (d_{1}) − Ke^{−rT}\mathcal{N} (d_{2})

Where:

d1=logSK+(r+12σ2)TσTd_{1} = \frac{\log \frac{S}{K} + (r+\frac{1}{2}\sigma ^{2})T}{\sigma \sqrt{T}}

And

d2=d1σTd_{2} = d_{1} - \sigma \sqrt{T}

=logSK+(r12σ2)TσT= \frac{\log \frac{S}{K} + (r -\frac{1}{2}\sigma ^{2})T}{\sigma \sqrt{T}}

8.2 European option premiums incorporating dividends

pricing equations for European call CC and put #PP options become

C=SeqTN(d1)KerTN(d2)C = Se^{-qT}\mathcal{N} (d_{1}) − Ke^{−rT}\mathcal{N} (d_{2})

P=KerTN(d2)SeqTN(d1)P = Ke^{−rT}\mathcal{N} (-d_{2}) - Se^{-qT}\mathcal{N} (-d_{1})

Where:

d1=logSK+(rq+12σ2)TσTd_{1} = \frac{\log \frac{S}{K} + (r - q + \frac{1}{2}\sigma ^{2})T}{\sigma \sqrt{T}}

And

d2=d1σTd_{2} = d_{1} - \sigma \sqrt{T}

d2=logSK+(rq12σ2)TσTd_{2} = \frac{\log \frac{S}{K} + (r - q - \frac{1}{2}\sigma ^{2})T}{\sigma \sqrt{T}}

8.3. European option premiums for currencies

We view the foreign risk-free rate rfr_f as the “dividend” on the underlying asset (the foreign currency):

Cf:d=Sf:derfTN(d1)Kf:derdTN(d2)C_{f:d} = S_{f:d}e^{-r_fT}\mathcal{N} (d_{1}) − K_{f:d}e^{−r_dT}\mathcal{N} (d_{2})

Pf:d=Kf:derdTN(d2)Sf:derfTN(d1)P_{f:d} = K_{f:d}e^{−r_dT}\mathcal{N} (-d_{2}) - S_{f:d} e^{-r_fT}\mathcal{N} (-d_{1})

Where:

d1=logSf:dKf:d+(rdrf+12σ2)TσTd_{1} = \frac{\log \frac{S_{f:d}}{K_{f:d}} + (r_d - r_f + \frac{1}{2}\sigma ^{2})T}{\sigma \sqrt{T}}

and

d2=d1σTd_{2} = d_{1} - \sigma \sqrt{T}

=logSf:dKf:d+(rdrf12σ2)TσT = \frac{\log \frac{S_{f:d}}{K_{f:d}} + (r_d - r_f - \frac{1}{2}\sigma ^{2})T}{\sigma \sqrt{T}}

9. key features/Assumptions of Black Scholes

  1. We’re using historical volatility for σ but in practice options traders tend to use the Black-Scholes model with “rules of thumb” and their own estimates of σ.
  2. Our dividend yield may not be accurate: We assumed a continuously compounded yield. It should be a forecast. And should we incorporate dividend imputation/franking?

Geometric Brownian Motion

10. Simulating geometric Brownian

So in the risk-neutral pricing approach the value of a European option is

12. Factors affecting options prices

How each of the input variables KK, SS, rr, TT, σσ and qq impact option premiums

We’re interested in the sensitivity of option prices to the other variables SS, rr, TT and σσ, and we give these sensitivities special Greek names:

12.1 Delta and Gamma

We know that as SS increases, calls premiums rise and put premiums fall.

alt text

We quantify this mathematically with the delta ∆, given by

C=N(d1)andP=N(d1)1∆_{C} = N(d_{1}) \: \: \text{and} \: \: ∆_{P} = N(d_{1}) − 1

Importantly, note the following

0<C<1and1<P<0 0 < ∆_{C} < 1 \: \: \text{and} \: \: -1 < ∆_{P} < 0

We interpret ∆ as the change in the premium due to a change in S, giving us the approximations

dCΔcdSanddPΔpdSdC \approx \Delta _{c} dS \: \: \text{and} \: \: dP \approx \Delta _{p} dS

where dCdC and dPdP are a change in the premium and dSdS a change in SS.

We calculate the “new” premiums due to a change in S as

CnewC+dCandPnewP+dPC_{new} \approx C + dC \: \: \text{and} \: \: P_{new} \approx P+dP

12.2 Rho

Rho ρ is the change in the premium from a change in rr, and is given by

ρc=KTerTN(d2)\rho_c = KTe^{-rT}\mathcal{N}(d_2)

and

ρp=KTerT[N(d2)1]\rho_p = KTe^{-rT} \left[ \mathcal{N}(d_2) -1 \right]

Where

0<ρc0 < \rho_c

and

ρp<0\rho_p < 0

Example for call prices:
alt text
(x axis is underlying stock price)

As rr increases, call premiums increase but put premiums decrease.

12.3 Vega

Vega ν is the change in the premium from a change in σ, and is given by:

ν=Sf(d1)T\nu = Sf(d_{1})\sqrt{T}

with f(x)f(x) the PDF of a standard normal random variable. Note that

0<ν0 < \nu

As σ increases, option premiums increase.

12.4 Theta

Theta θ is a bit ambiguous. It gives the negative of the change in the premium from a change in TT. And the equations are more complex:

θc=Sf(d1)σ2TrKerTN(d2)\theta_c = -\frac{Sf(d_1)\sigma}{2\sqrt{T}} -rKe^{-rT}\mathcal{N}(d_2)

and

θp=θc+rKerT\theta _p = \theta _c + rKe^{-rT}

θ is the negative of the partial derivative of the premium with respect to T, telling us the impact of approaching expiry.

θc<0\theta _c < 0

but we may have

θp0\theta _p \le 0

and

0θp 0 \le \theta _p

12.4.1 Time Decay

Premiums usually fall as expiry approaches

Theta Call Theta Put
alt text alt text
alt text alt text

13. Scholes Hedging

13.1 Static Delta Hedging

Recall the basic long call and put option payoff and premium plots:

alt text

Delta hedging involves taking a position in the asset to hedge an option position against small movements in the asset’s price.

An option’s delta is an approximation of the change in the premium due to a change dS in the price S of the underlying asset:

dCCdS    and    dPPdS.dC ≈ ∆_C dS \; \; \text{and} \; \; dP ≈ ∆_P dS.

Static delta hedging with h options

The value of a portfolio of hh calls and QQ units in the underlying asset is

V=QS+hCV = QS + hC

The value of a portfolio of hh puts and QQ assets is V=QS+hPV = QS + hP. Its change is dV(Q+hP)dSdV ≈ (Q + h∆_P )dS so we set

Q=hPQ = −h∆_P

or

Q=hCQ = −h∆_C

for hedging options

13.2 Delta-Gamma hedging

We can use delta-gamma hedging to improve the delta hedging of an option position by taking a position in the asset and a different option.

For delta-gamma hedging we calculate how many units QQ in the asset and kk in another option we need in order to hedge against small changes ddS in the price SS of the underlying asset.

The value of a portfolio of QQ units in the asset, hh units of one option, and kk units of another, different option is given by

V=QS+hV1+kV2V = QS + hV_1 + kV_2

where V1V_1 and V2V_2 are the respective option premiums.

We can show that the solution is

Q=hΔ1Γ2Δ2Γ1Γ2Q = -h \frac{\Delta _1 \Gamma _2 - \Delta _2 \Gamma _1}{\Gamma _2}

And

k=QΓ1Δ1Γ2Δ2Γ1k = \frac{Q \Gamma _1}{\Delta _1 \Gamma _2 - \Delta _2 \Gamma _1}

IMPORTANT:

13.3 Dynamic delta hedging

14. Implied Volatility

In the Black-Scholes European option the only unobservable parameter is the volatility σ.

Given observed option prices CC or PP and observable market variables SS, KK, rr and TT, an option’s implied volatility is the volatility parameter σ that yields Black-Scholes prices equal to CC or PP.

We have to code numerical (iterative) techniques on computers.

In markets, Black-Scholes implied vols are not constant but display:

Volatility smirk/smiles
alt text

Remark We see this from the volatility smile and term structure:

Consequently traders use:

14.1 VIX index

It “is a calculation designed to produce a measure of constant, 30-day expected volatility of the U.S. stock market, derived from real-time, mid-quote prices of S&P 500 Index (SPX) call and put options.”

VIX index example
alt text

A value of 17.55 means that average implied vols of 30 day S&P 500 Index options, thus the market’s 30 day volatility expectations, is 17.55%.

15 Trading Strategies

  1. Directional strategies: Speculate on the direction of the price of the underlying asset, similar to taking calls and puts.
  2. Volatility strategies: Speculate on high or low asset volatility, or changes in implied vols, often incorporating delta neutrality.
  3. Time: Strategies that seek to take advantage of time decay, typically assuming low asset volatility and relatively constant implied vols.

15.1 Directional strategies

Directional strategies seek to speculate on a directional change in the price of the asset, but at a lower upfront cost to the basic strategies of taking calls and puts, while still maintaining limited downside risk.

15.1.1 Bull Spread

Suppose we want to speculate on an increase in the asset price:

We could take an ATM call with strike K=50K = 50, costing C1=4.13C_1 = 4.13. But also write an OTM call with strike K2=55K_2 = 55.

alt text

A bull spread has lower downside risk and “profits sooner” compared to the basic long call, but its upside profits are capped.

15.1.2 Bear Spread

Speculate on a decrease in the asset price:

We can lower the cost of taking a put with strike K=50K = 50 by also writing a put with strike K2=45K_2 = 45.

alt text

15.2 Volatility Strategies

15.2.1 Long Straddle

Suppose we think that the asset price will be volatile:

We could take ATM calls and puts both with strike K = 50, a strategy called a long straddle with payoff and profit:

alt text

15.2.2 Short Butterfly

As shown above

15.2.3 Short Straddle

Suppose that we expect a period of low volatility

alt text

15.2.4 Long Butterfly

As shown above

16. Binomial and Monte Carlo numerical option pricing methods

16.1 & 16.3 Rationale

  1. When applying the Black-Scholes framework to derive pricing models for more complex, exotic derivatives (even just for American options), we immediately run into the scenario of being unable to derive closed-form, analytical solution equations.
  2. We may also want to apply more complex derivative security pricing frameworks than the Black-Scholes framework (say the Heston stochastic volatility, or the Merton jump-diffusion, frameworks), but again we immediately run into the same problem.

Stylised features of financial returns:

16.2 Methods

  1. Lattice type methods like the binomial and trinomial models.
  2. Monte Carlo methods:
  1. Partial differential equation numerical solution methods:

17. Binomial model for European option price and delta

17.1 1-Period

The call option’s payoffs in each possible outcome of the price of the underlying asset are:

Cu=max(0,SuK)C_u = \max(0, S_u − K)

Cd=max(0,SdK)C_d = \max(0, S_d − K)

Call option

C=erT[qCu+(1q)Cd]C = e^{-rT}[qC_u + (1-q)C_d]

Where

q=erTdudq = \frac{e^{rT} - d}{u-d}

17.1.1. Up and down factors

  1. The Jarrow-Rudd (JR) scheme derives

u=e(rσ2/2)T+σT      and      d=e(rσ2/2)TσTu=e^{(r - \sigma^2 /2)T + \sigma \sqrt{T}} \;\;\; \text{and} \;\;\; d = e^{(r - \sigma ^2 / 2) T - \sigma \sqrt{T}}

17.2 Multi-period

For the multi-period binomial model we simply:

  1. Discretise the interval [0, T] into N + 1 equally-spaced dates t0,t1,...,tN{t_0 , t_1 , . . . , t_N } with t0=0,tN=Tt_0 = 0, t_N = T and spacing dt=T/Ndt = T/N.
  2. Build an asset price tree starting at SS as per the next slide.
  3. Calculate the option premium by recursively stepping backwards in time in the tree and using the 1-period pricing formula each step.

Build asset price tree of S

To build the asset price tree, we note that the asset price can go up by a factor of uu or down by a factor of dd at each time step.

alt text

Asset price SiN_{iN} took ii up steps and thus NiN − i down steps.

Calculate Asset Price tree of C/P:

We calculate call prices as follows:

CiN=max{0,SiNK}C_{iN} = \max \{0, S_{iN} − K \}

for i=0,1,...,Ni = 0, 1, . . . , N.

alt text

Cij=erdt[qCi+1,j+1+(1q)Ci,j+1]C_{ij} = e^{-rdt}\left[ qC_{i+1, j+1} + (1-q)C_{i, j+1}\right]

17.3 Binal price deltas

Δ=CuCdS(ud)\Delta = \frac{C_u - C_d}{S(u-d)}

Δ=PuPdS(ud)\Delta = \frac{P_u - P_d}{S(u-d)}

17.3 Rationale of Risk-neutral approach

-- to complete

18. Monte Carlo Pricing of European Options

Option prices are given by

C=erTE[max{0,STK}]C = e^{-rT} \mathbb{E} \left[ \max \{ 0, S_T - K \} \right]

and

P=erTE[max{0,KST}]P = e^{-rT} \mathbb{E} \left[ \max \{ 0, K - S_T \} \right]

Where STS_T is log-normally distributed as per geometric Brownian motion.

To price options via the Monte Carlo method we:

19. Pricing American and path dependent options

19.1 American Options (using Binomial Tree)

An American option gives the holder the right (but not the obligation) to exercise it at any point up to and including the expiry date TT.

American options are worth at least as much as European options:

0CEuCAm0 ≤ C_{Eu} ≤ C_{Am}

and

0PEuPAm0 ≤ P_{Eu} ≤ P_{Am}

American options are worth at least as much as their intrinsic value:

max0,SKCAm\max{0, S − K} ≤ C_{Am}

and

max0,KSPAm\max{0, K − S} ≤ P_{Am}

The adjustment to price American options is:

A each node, set the American option price equal to the maximum of the “1-step European price” and the option’s intrinsic value:

The (American and European) option payoffs at expiry are

CiNAm=max{0,SiNK}C^{Am}_{iN} = \max \{ 0, S_{iN} - K \}

and

PiNAm=max{0,KSiN}P^{Am}_{iN} = \max \{ 0, K - S_{iN} \}

The “1-step European” option prices at node ij on the tree are

CijEu=erdt[qCi+1,j+1Am+(1q)Ci,j+1Am]C^{Eu}_{ij} = e^{-rdt}\left[ qC^{Am}_{i+1, j+1} + (1-q)C_{i, j+1}^{Am} \right]

PijEu=erdt[qPi+1,j+1Am+(1q)Pi,j+1Am]P^{Eu}_{ij} = e^{-rdt}\left[ qP^{Am}_{i+1, j+1} + (1-q)P_{i, j+1}^{Am} \right]

Hence, the American option prices at node ijij are

Cijiv=max{CijEu,Cijiv}C^{iv}_{ij} = \max \{C^{Eu}_{ij}, C^{iv}_{ij} \}

Pijiv=max{PijEu,Pijiv}P^{iv}_{ij} = \max \{P^{Eu}_{ij}, P^{iv}_{ij} \}

19.2 Path Dependent Options

A European option is path dependent if its payoff is calculated from the underlying asset’s path travelled over the option’s life.

19.3 European Chooser Option (binomial method)

A European chooser option is similar to a plain vanilla European option except that it allows the holder to choose at some date tt choose (the choice date) over the option’s life if the option is a call or a put

Modification to original Binomial Tree

19.4 Lookback options (Monte-Carlo)

A lookback option’s payoff at expiry depends on the maximum or minimum prices of the underlying asset over the life of the option

19.4.1 Fixed-Strike lookback Options

European fixed-strike lookback options are very similar to plain vanilla European options except the “final price” of the underlying asset used in calculating the payoff is not the asset price STS_{T} but instead the maximum SmaxS_{max} or minimum SminS_{min} asset price reached over the option’s life:

The maximum and minimum asset prices of path ii are given by

Si,max=maxj=0,...,MSijS_{i, max} = \underbrace{\max}_{j=0,...,M}S_{ij}

Si,min=minj=0,...,MSijS_{i, min} = \underbrace{\min}_{j=0,...,M}S_{ij}

The Monte Carlo Option prices are then simply

C=erT1Ni=1Nmax{0,Si,maxK}C = e^{-rT} \frac{1}{N} \sum ^N _{i=1} \max \{ 0, S_{i, \max } - K \}

P=erT1Ni=1Nmax{0,KSi,min}P = e^{-rT} \frac{1}{N} \sum ^N _{i=1} \max \{ 0, K - S_{i, \min } \}

19.4.2. (European) Floating-Strike Lookback Options

European floating-strike lookback options differ from fixed-strike options by instead setting the strike price KK to be the maximum SmaxS_{max} or minimum SminS_{min} asset prices over the life of the option. Their payoffs are

call payoff=max{0,STSmin}\text{call payoff} = \max \{ 0, S_T - S_{\min } \}

put payoff=max{0,SmaxST}\text{put payoff} = \max \{ 0, S_{\max } - S_T\}

So, after calculating the N asset price paths {Si0,Si1,...,SiM}\{ S_{i0} , S_{i1} , . . . , S_{iM} \} for i = 1, . . . , N by simulating geometric Brownian motion, and calculating the minimum S i,min and maximum S i,max prices for each path, the Monte Carlo floating-strike lookback option prices are then simply

C=erT1Ni=1Nmax{0,SiMSi,min}C = e^{-rT} \frac{1}{N} \sum ^N _{i=1} \max \{ 0, S_{iM} - S_{i, \min } \}

P=erT1Ni=1Nmax{0,Si,maxSiM}P = e^{-rT} \frac{1}{N} \sum ^N _{i=1} \max \{ 0, S_{i, \max } - S_{iM} \}

19.5 (European) Barrier Options

They are effectively plain vanilla European options whose payoff's “knock in” or “knock-out” if the price of the underlying asset hits a barrier BB at some point over the life of the option:

19.5.1 Knock-in Barrier Options

European knock-in options are activated if the price of the underlying asset hits the barrier B at some point of the option’s life.

The payoffs are then the usual

If the barrier is never hit, the option is never activated and expires worthless

There’s two kinds of knock-in options depending on the relation between BB and SS

  1. Up-and-in options set B>SB > S, noting that typically KSK ≈ S.
  2. Down-and-in options set B<SB < S, noting that typically KSK ≈ S.

The Monte Carlo pricing of knock-in options is a simple modification to the above, simply to check each path to see if the barrier was hit.

19.5.2. Knock-out Barrier Options

European knock-out options are deactivated if the price of the underlying asset hits the barrier B at some point of the option’s life.

If the barrier is never hit, the option stays alive and the payoff's are the usual max{0,STK}\max \{ 0, S_T − K \} for a call and max{0,KST}\max \{ 0, K − S_T \} for a put. There’s two kinds of knock-outs depending on the relation between B and S:

  1. Up-and-out options set B>SB > S, noting that typically KSK ≈ S.
  2. Down-and-out options set B<SB < S, noting that typically KSK ≈ S.

The Monte Carlo pricing of knock-out options is a simple modification to the above for knock-in options, namely a reversal of the logical conditions on the payoffs if the barrier was hit.

19.6. (European) Asian Options

A European Asian or average option’s payoffs depend on the average underlying asset price Sˉ\bar{S} over the life of the option.

There’s two general types of Asian options:

  1. Fixed-strike Asian options.
  2. Average-strike Asian options.

19.6.1. Fixed-Strike Asian Options

European fixed-strike Asian options are very similar to plain vanilla European options except the “final price” of the underlying asset used in calculating an option’s payoff is not the asset price itself but the average asset price Sˉ\bar{S} over the life of the option. Their payoffs are

It is simple to use Monte Carlo simulation to price them:

So, after calculating the NN asset price paths {Si0,Si1,...,SiM}\{ S_{i0} , S_{i1} , . . . , S_{iM} \} for i=1,...,Ni = 1, . . . , N by simulating geometric Brownian motion, the average price of path ii is

Sˉi=1Mj=0MSij\bar{S}_{i} = \frac{1}{M} \sum^{M}_{j=0}S_{ij}

(arithmetic average). The payoffs for path ii are

Put and call prices are then:

C=erT1Ni=1Nmax{0,SˉiK}C = e^{-rT} \frac{1}{N} \sum ^N _{i=1} \max \{ 0, \bar{S}_{i} − K \}

P=erT1Ni=1Nmax{0,KSˉi}P = e^{-rT} \frac{1}{N} \sum ^N _{i=1} \max \{ 0, K - \bar{S}_{i} \}

19.6.2. Average-strike Asian options

European average-strike Asian options differ from fixed-strike Asian options by instead setting the strike price K=SˉK = \bar{S} to be the average asset price Sˉ\bar{S} over the life of the option. So, their payoff s are simply