last updated: May 1, 2025

Stochastic Processes and Financial Mathematics
(part two)

Chapter C Solutions to exercises (part two)

Chapter 11
  • 11.1

    • (a) We have Ct=μt+σBt. Hence,

      E[Ct]=E[μt+σBt]=μt+σE[Bt]=μtE[Ct2]=E[μ2tt+2tμσBt+σ2Bt2]=μ2t2+2tμσ(0)+σ2t=μ2t2+σ2tvar(Ct)=E[Ct2]E[Ct]2=σ2t. where we use that E[Bt]=0 and E[Bt2]=t.

    • (b) We have

      CtCu=μt+σBtμuσBu=μ(tu)+σ(BtBu)μ(tu)+σN(0,tu)

      where, in the final step, we use the definition of Brownian motion. Then, by the scaling properties normal random variables we have CtCuN(μ(tu),σ2(tu)).

    • (c) Yes. By definition, Brownian motion Bt is a continuous stochastic process, meaning that the probability that Bt is a continuous function is one. Since tμt is a continuous function, we have that μt+σBt is a continuous function with probability one; that is, Bt is a continuous stochastic process.

    • (d) We have E[Ct]=μt, but Brownian motion has expectation zero, so Ct is not a Brownian motion.

  • 11.2 We have

    cov(Bu,Bt)=E[BtBu]E[Bt]E[Bu]=E[BtBu]=E[(BtBu)Bu]+E[Bu2]=E[BtBu]E[Bu]+E[Bu2]=0×0+u=u Here, to deduce the fourth line, we use the second property in Theorem 11.2.1, which tells us that BtBu and Bu are independent.

  • 11.3 When ut the martingale property of Brownian motion (Lemma 11.4.3) implies that E[Bt|Fu]=Bu. When tu we have BtFu so by taking out what is known we have E[Bt|Fu]=Bt. Combing the two cases, for all u0 and t0 we have E[Bt|Fu]=Bmin(u,t).

  • 11.4

    • (a) We’ll use the pdf of the normal distribution to write E[Btn] as an integral. Then, integrating by parts (note that ddzez22t=ztez22t) we have

      E[Btn]=12πtznez22tdz=12πt(tzn1)(ztez22t)dz=[12πttzn1ez22t]z=+12πtt(n1)zn2tez22tdz=0+t(n1)E[Btn2].

    • (b) We use the formula we deduced in part (a). Since E[Bt0]=E[1]=1, we have E[Bt2]=t(21)(1)=t. Hence E[Bt4]=t(41)E[Bt2]=t(41)t=3t2 and therefore var(Bt2)=E[Bt4]E[Bt2]2=2t2.

    • (c) Again, we use the formula we deduced in part (a). Since E[Bt]=0 it follows (by a trivial induction) that E[Btn]=0 for all odd nN. For even nN we have E[Bt2]=t and (by induction) we obtain

      E[Btn]=tn/2(n1)(n3)(1).

    • (d) We have var(Btn)=E[Bt2n]E[Btn]2 which is finite by part (c). Hence BtnL2, which implies that BtnL1.

  • 11.5 Using the scaling properties of normal random variables, we write Z=μ+Y where YN(0,σ2). Then, eZ=eμeY and

    E[exp(σBt12σ2t)|Fu]=E[exp(σ(BtBu)+σBu12σ2t)|Fu]=exp(σBu12σ2t)E[exp(σ(BtBu))|Fu]=exp(σBu12σ2t)E[exp(σ(BtBu))]=exp(σBu12σ2t)exp(12σ2(tu))=exp(σBu12σ2u). Here, the second line follows by taking out what is known, since Bu is Fu measurable. The third line then follows by the definition of Brownian motion, in particular that BtBu is independent of Fu. The fourth line follows by (11.2), since BtBuN(0,tu) and hence σ(BtBu)N(0,σ2(tu)).

  • (b) Since Bt is adapted, Bt33tBt is also adapted. From 11.4 we have Bt3,BtL1, so also Bt3tBtL1. Using that Bu is Fu measurable, we have

    E[Bt33tBt|Fu]=E[(Bu33uBu)+Bt33tBt(Bu33uBu)|Fu]=Bu33uBu+E[Bt3Bu33tBt+3uBu|Fu] so we need only check that the second term on the right hand side is zero. To see this,

    E[Bt3Bu33tBt+3uBu|Fu]=E[Bt3Bu33tBu+3uBu3t(BtBu)|Fu]=E[Bt3Bu33tBu+3uBu|Fu]+0=E[Bt3Bu3|Fu]3Bu(tu)=E[(BtBu)3+3Bt2Bu3Bu2Bt|Fu]3Bu(tu)=E[(BtBu)3|Fu]+3BuE[Bt2|Fu]3Bu2E[Bt|Fu]3Bu(tu)=E[(BtBu)3]+3Bu(E[Bt2t|Fu]+t)3Bu33Bu(tu)=0+3Bu(Bu2u+t)3Bu33Bu(tu)=0. Here we use several applications of the fact that BtBu is independent of Fu, whilst Bu is Fu measurable. We use also that E[Z3]=0 where ZN(0,σ2), which comes from part (c) of 11.4 (or use that the normal distribution is symmetric about 0), as well as that both Bt and Bt2t are martingales (from Lemmas 11.4.3 and 11.4.4).

  • 11.7

    • (a) We have

      k=0n1(tk+1tk)=(tntn1)+(tn1tn2)++(t2t1)+(t1t0)=tnt0=t0=t. This is known as a ‘telescoping sum’. The same method shows that k=0n1(Btk+1Btk)=BtB0=0.

    • (b) We need a bit more care for this one. We have

      0k=0n1(tk+1tk)2k=0n1(tk+1tk)(maxj=0,,n1|tj+1tj|)=(maxj=0,,n1|tj+1tj|)k=0n1(tk+1tk)=(maxj=0,,n1|tj+1tj|)t. Here, the last line is deduced using part (a). Letting n0 we have maxj=0,,n1|tj+1tj|0, so the right hand side of the above tends to zero as n. Hence, using the sandwich rule, we have that k=0n1(tk+1tk)20.

    • (c) Using the properties of Brownian motion, Btk+1BtkN(0,tk+1tk) so from exercise 11.4 we have

      E[k=0n1(Btk+1Btk)2]=k=0n1E[(Btk+1Btk)2]=k=0n1(tk+1tk)=t. Here, the last line follows by part (a).

      For the last part, the properties of Brownian motion give us that each increment Btk+1Btk is independent of Ftk. In particular, the increments Btk+1Btk are independent of each other. So, using exercise 11.4,

      var(k=0n1(Btk+1Btk)2)=k=0n1var((Btk+1Btk)2)=k=0n12(tk+1tk)2 which tends to zero as n, by the same calculation as in part (b).

  • 11.8

    • (a) Let y1. Then

      P[Bty]=y12πtez22tdzy12πtyez22tdzy12πtzez22tdz=12πt[tez22t]z=y=t2πey22t. Putting y=tα where α>12 we have

      (C.1)P[Bttα]t2πe12t2α1.

      Since 2α10, the exponential term dominates the square root, and right hand side tends to zero as t.

      If α=12 then we can use an easier method. Since BtN(0,t), we have t1/2BtN(0,1) and hence

      P[Btt1/2]=P[N(0,1)1](0,1),

      which is independent of t and hence does not tend to zero as t. Note that we can’t deduce this fact using the same method as for α>12, because (C.1) only gives us an upper bound on P[Bttα].

    • (b) For this we need a different technique. For y0 we integrate by parts to note that

      y12πtez22tdz=12πty1zzey22tdz=12πt([tzez22t]z=yytz2ez22t)dz12πt([tzez22t]z=y)=12πttyey22t=t2π1yey22t Using the symmetry of normal random variables about 0, along with this inequality, we have

      P[|Bt|a]=P[Bta]+P[Bta]=2P[Bta]2t2π1aea22t. As t0, the exponential term tends to 0, which dominates the t, meaning that P[|Bt|a]0 as t0. That is, Bt0 in probability as t0.

  • Chapter 12
    • 12.1 From (12.8) we have both 0t1dBu=Bt and 0sdBu=Bs, and using (12.6) we obtain

      ut1dBu=0t1dBu0v1dBu=BtBv.

    • 12.2 Since Bt is adapted to Ft, we have that eBt is adapted to Ft. Since Bt is a continuous stochastic process and exp() is a continuous function, eBt is also continuous. From (11.2) we have that

      0tE[(eBu)2]du=0tE[e2Bu]du=0te12(22)udu=0te2udu<. Therefore, eBtH2.

    • 12.3

      • (a) We have

        E[eZ22]=12πez22ez22dz=12π1dz=.

      • (b) Using the scaling properties of normal distributions, t1/2Bt has a N(0,1) distribution for all t. Hence, if we set Ft=t1/2Bt then by part (a) we have

        0tE[e12Ft2]du=0tdu

        which is not finite. Note also that Ft is adapted to Ft. Since et, t2, t1/2 and Bt are all continuous, so is e12Ft2. Hence e12Ft2 is an example of a continuous, adapted stochastic process that is not in H2.

        Note that we can’t simply use the stochastic process Ft=Z, because we have nothing to tell us that Z is Ft measurable.

    • 12.4 We have

      E[Xt]=E[2]+E[0tt+Bu2du]+E[0tBu2dBu]=2+0tt+E[Bu2]du+0=2+t2+0tudu=2+t2+t22=2+3t22.

    • 12.5

      • (a) Xt=0+0t0du+0t0dBu so Xt is an Ito process.

      • (b) Yt=0+0t2udu+0t1dBu by (12.8), so Yt is an Ito process.

      • (c) A symmetric random walk is a process in discrete time, and is therefore not an Ito process.

    • 12.6 We have

      E[Vt]=E[ektv]+σektE[0teksdBs]=ektv. Here we use Theorem 12.2.1 to show that the expectation of the dBu integral is zero. In order to calculate var(Vt) we first calculate

      E[Vt2]=E[e2ktv2]+2σektE[0teksdBs]+σ2e2ktE[(0teksdBs)2]=e2ktv2+0+σ2e2kt0tE[(eks)2]du=e2ktv2+σ2e2kt0te2kudu=e2ktv2+σ2e2kt12k(e2kt1)=e2ktv2+σ22k(1e2kt) Again, we use Theorem 12.2.1 to calculate the final term on the first line. We obtain that

      var(Vt)=E[Vt2]E[Vt]2=σ22k(1e2kt).

    • 12.7 We have

      Xt=μt+0tσudBu.

      Since Mt=0tσudBu is a martingale (by Theorem 12.2.1), we have that Mt is adapted and in L1, and hence also Xt is adapted and in L1. For vt we have

      E[Xt|Fv]=μt+E[Mt|Fv]=μt+Mv=μt+0vσudBuμv+0vσudBu=Xv. Hence, Xt is a submartingale.

    • 12.8

      • (a) Taking Ft=0, we have E[Ft]= and 0tFsds=0, hence 0tE[Fs]ds=E[0tFsdBs]=0.

      • (b) Taking Ft=1, we have E[Ft]=1 and 0tFsds=t, hence 0tE[Fs]ds=t and E[0tFsdBs]=E[Bt]=0.

    • 12.9

      • (a) We have (|X||Y|)20, so 2|XY|X2+Y2, which by monotonicity of E means that 2E[|XY|]E[X2]+E[Y2]. Using the relationship between E and || we have

        2|E[XY]|2E[|XY|]E[X2]+E[Y2].

      • (b) Let Xt,YtH2 and let α,βR be deterministic constants. We need to show that ZtαXt+βYtH2. Since both Xt and Yt are continuous and adapted, Zt is also both continuous and adapted. It remains to show that (12.5) holds for Zt. With this in mind we note that

        Zt2=α2Xt2+2αβXtYt+β2Yt2

        and hence that

        |E[Zt2]|=|α2E[Xt2]+2αβE[XtYt]+β2E[Yt2]|α2E[Xt2]+2|αβ||E[XtYt]|+β2E[Yt2]α2E[Xt2]+|αβ|(E[Xt2]+E[Yt2])+β2E[Yt2]=(α2+|αβ|)E[Xt2]+(β2+|αβ|)|E[Yt2]. where we use part (a) to deduce the third line from the second. Hence,

        0tE[Zu2]du0t(α2+|αβ|)E[Xu2]+(β2+|αβ|)|E[Yu2]du=(α2+|αβ|)0tE[Xu2]du+(β2+|αβ|)0tE[Yu2]du< as required. The final line is < because Xt,YtH2.

    • 12.10 We have IF(t)=i=1nFti1[BtitBti1t]. We are looking to show that

      (C.2)E[IF(t)2]=0tE[Fu2]du.

      On the right hand side we have

      0tE[Fu2]du=i=1mti1ttitE[Fu2]du(C.3)=i=1m(titti1t)E[Fti12] because Ft is constant during each time interval [ti1t,tit). On the left hand side of (C.2) we have

      E[IF(t)2]=E[(i=1mFti1[BtitBti1t])2]=E[i=1mFti12[BtitBti1t]2+2i=1mj=1i1Fti1Ftj1[BtitBti1t][BtjtBtj1t]]=i=1mE[Fti12[BtitBti1t]2]+2i=1mj=1i1E[Fti1Ftj1[BtitBti1t][BtjtBtj1t]] In the first sum, using the tower rule, taking out what is known, independence, and then the fact that E[Bt2]=t, we have

      E[Fti12[BtitBti1t]2]=E[E[Fti12[BtitBti1t]2|Fti1]]=E[Fti12E[[BtitBti1t]2|Fti1]]=E[Fti12E[[BtitBti1t]2]]=E[Fti12(titti1t)] In the second sum, since j<i we have tj1ti, so using the tower rule, taking out what is known, and then the martingale property of Brownian motion, we have

      E[Fti1Ftj1[BtitBti1t][BtjtBtj1t]]=E[E[Fti1Ftj1[BtitBti1t][BtjtBtj1t]|Fti1]]=E[Ftj1[BtjtBtj1t]Fti1(E[Btit|Fti1]Bti1t)]=E[Ftj1[BtjtBtj1t]Fti1(Bti1tBti1t)]=0. Therefore,

      E[IF(t)2]=i=1mE[Fti12(titti1t)]

      which matches (C.3) and completes the proof.

    Chapter 13
    • 13.1

      • (a) Xt=X0+0t2udu+0tBudBu.

      • (b) YT=Yt+tTudu.

      By using the fundamental theorem of calculus, we obtain that Y satisfies the differential equation dYtdt=t. Using equation (13.6) from Example 13.1.2, we have that

      Xt=X0+t2+Bt22t2 which is not differentiable because Bt is not differentiable.

    • 13.2 We have Zt=f(t,Xt) where f(t,x)=t3x and dXt=αdt+βdBt. By Ito’s formula,

      n(n1)20tu(n2)/2(n3)(n5)(1)du=n(n1)2tn/2n/2(n3)(n5)(1)du=tn/2(n1)(n3)(1) as required.

    • 13.5

      • (a) By Ito’s formula, with f(t,x)=et/2cosx, we have

        dXt=(12et/2cos(Bt)+(0)(et/2sin(Bt))+12(12)(et/2cos(Bt)))dt+(1)(et/2sin(Bt))dBt=et/2sin(Bt)dBt. Hence Xt=X00teu/2sin(Bu)dBu, which is martingale by Theorem 12.2.1.

      • (b) By Ito’s formula, with f(t,x)=(x+t)ext/2 we have

        dYt={eBtt/212(Bt+t)eBtt/2+(0)(eBtt/2(Bt+t)eBtt/2)+12(12)(eBtt/2eBtt/2+(Bt+t)eBtt/2)}dt+(12)(eBtt/2(Bt+t)eBtt/2)dBt=(1tBt)eBtt/2dBt. Hence Yt=Y0+0t(1uBu)eBuu/2dBu, which is martingale by Theorem 12.2.1.

    • 13.6 We apply Ito’s formula with f(t,x)=tx to Zt=f(t,Bt)=tBt and obtain

      dZt=(Bt+(0)(t)+12(12)(0))dt+(1)(t)dBt

      so we obtain

      tBt=0B0+0tBudu+0tudBu,

      as required.

    • 13.7

      • (a) We have Xt=X0+0t2+2sds+0tBsdBs. Taking expectations, and recalling that Ito integrals have zero mean, we obtain that

        E[Xt]=X0+0t2+2sds+0=1+[2s+s2]s=0t=1+2t+t2=(1+t)2.

      • (b) From Ito’s formula,

        dYt=(0+(2+2t)(2Xt)+12(Bt)2(2))dt+Bt(2Xt)dBt=(4(1+t)Xt+(Bt)2)dt+2XtBtdBt. Writing in integral form, taking expectations, and using that Ito integrals have zero mean, we obtain

        E[Yt]=Y0+E[0t4(1+s)Xt+(Bs)2ds]+0=1+0t4(1+s)E[Xs]+E[Bs2]ds=1+0t4(1+s)3+sds=1+[(1+s)4+s22]s=0t=(1+t)4+t22 Hence, using that E[Xt2]=E[Yt],

        var(X)=E[Xt2]E[Xt]2=(1+t)4+t22(1+t)4=t22.

      • (c) If we change the dBt coefficient then we won’t change the mean, because we can see from (a) that E[Xt] depends only on the dBt coefficient. However, as we can see from part (b), the variance depends on both the dt and dBt terms, so will typically change if we alter the dBt coefficient.

    • 13.8 In integral form, we have

      Xt=X0+0tαXudu+0tσudBu.

      Taking expectations, swapping du with E, and recalling from Theorem 12.2.1 that Ito integrals are zero mean martingales, we obtain

      E[Xt]=E[X0]+0tαE[Xu]du.

      Applying the fundamental theorem of calculus, if we set xt=E[Xt], we obtain

      dxtdt=αxt

      which has solution xt=Ceαt. Putting in t=0 shows that C=E[X0]=1, hence

      E[Xt]=eαt.

    • 13.9 We have Xt=X0+0tXsdBs. Since Ito integrals are zero mean martingales, this means that E[Xt]=E[X0]=1. Writing Yt=Xt2 and using Ito’s formula,

      dYt=(0+(0)(2Xt)+12(Xt)2(2))dt+(Xt)(2Xt)dBt=(Xt2)dt+2Xt2dBt.=Ytdt+2YtdBt Writing in integral form and taking expectations, we obtain

      E[Yt]=1+0tE[Ys]ds+0. Hence, by the fundamental theorem of calculus, f(t)=E[Yt] satisfies the differential equation f(t)=f(t). The solution of this differential equation is f(t)=Aet. Since E[Y0]=1 we have A=1 and thus E[Yt]=et. Hence,

      var(Xt)=E[Xt2]E[Xt]2=E[Yt]1=et1.

    • 13.10

      Recall that cov(Xs,Xt)=E[(XtE[Xt])(XsE[Xs]). From (13.12) we have

      E[Xt]=E[eθt(Xtμ)]+E[0teθ(tu)dBu]=E[eθt(Xtμ)]+0,

      because Ito integrals are zero mean martingales. Thus XtE[Xt]=0tσeθ(ut)dBu. We can now calculate

      cov(Xs,Xt)=Cov(0sσeθ(us)dBu,0tσeθ(vt)dBv)=E[(0sσeθ(us)dBu)(0tσeθ(vt)dBv)]=σ2eθ(s+t)E[(0seθudBu)(0teθvdBv)]=σ2eθ(s+t)E[(0seθudBu+steθudBu)(0seθvdBv)]=σ2eθ(s+t)E[(0seθudBu)2]+E[(0seθvdBv)(steθudBu)]=σ2eθ(s+t)0se2θudu+E[0seθvdBv]E[steθudBu]=σ22θeθ(s+t)(e2θs1)+0. In the above, the second line and final lines follow because Ito integrals have zero mean. The penultimate line follows because 0seθvdBv and steθudBu are independent. This fact follows from the independence property in Theorem 11.2.1 (which defines Brownian motion), which implies that σ(BvBu;vus) and σ(BvBu;svu) are independent. Note that from (12.3) we know that abFudBs only depends on increments of Brownian motion between times a and b. (The same fact can be deduced from the Markov property from Section 14.2.)

      Hence var(Xt)=σ22θeθ(2t)(e2θt1)=σ22θ(1e2θt). It follows immediately that var(Xt)σ22θ as t.

    • 13.11

      • (a) Applying Ito’s formula to Zt=Bt3, with f(t,x)=x3, we have

        dZt=(0+(0)(3Bt2)+12(12)(6Bt))dt+(1)(3Bt2)dBt=3Btdt+3Bt2dt and substituting in for Z we obtain

        dZt=3Zt1/3dt+3Zt2/3dBt

        as required.

      • (b) Another solution is the (constant, deterministic) solution Xt=0.

    • 13.12 Equation (13.9) says that

      Xt=X0exp((α12σ2)t+σBt).

      Using Ito’s formula, with f(t,x)=X0exp((α12σ2)t+σx) we obtain

      dXt=((α12σ2)Xt+(0)(σXt)+12(12)(σ2Xt))dt+(1)(σXt)dBt=αXtdt+σXtdBt and thus Xt solves (13.8).

    • 13.13 Equation (13.14) says that

      Xt=X0exp(0tσudBu120tσu2du)

      We need to arrange this into a form where we can apply Ito’s formula. We write

      Xt=X0exp(Yt120tσu2du)

      where dYt=σtdBt with Y0=0. We now have Xt=f(t,Yt) where f(t,y)=X0exp(y120tσudu), so from Ito’s formula (and the fundamental theorem of calculus) we obtain

      dXt=(12σt2Xt+(0)(Xt)+12(σt2)(Xt))dt+(σt)(Xt)dBt=σtXtdBt. Hence, Xt solves (13.13).

    • 13.14

      • (a) This is essentially Example 3.3.9 but in continuous time. By definition of conditional expectation (i.e. Theorem 3.1.1) we have that MtL1 and that MtFt. It remains only to use the tower property to note that for 0ut we have

        E[Mt|Fu]=E[E[Y|Ft]Fu]=E[Y|Fu]=Mu.

      • (b)

        • (i) Note that M0=E[BT2|F0]=E[BT2]=T. We showed in Lemma 11.4.4 that Bt2t was a martingale, hence

          E[BT2|Ft]=E[BT2T+T|Ft]=Bt2t+T Using (13.6), this gives us that

          E[BT2|Ft]=20tBudBu+tt+T

          so we obtain

          Mt=T+0t2BudBu

          and we can take ht=2Bt.

        • (ii) Note that M0=E[BT3|F0]=E[BT3]=0. We showed in 11.6 that Bt33tBt was a martingale. Hence,

          E[BT3|Ft]=E[BT33TBT+3TBT|Ft]=Bt33tBt+3TBt. Using Ito’s formula on Zt=Bt3, we obtain dZt={0+(0)(3Bt2)+12(12)(6Bt)}dt+(1)(3Bt2)dBt so as

          Bt3=0+0t3Budu+0t3Bu2dBu.

          Also, from 13.6 we have

          tBt=0tBudu+0tudBu,

          so as

          E[BT3|Ft]=30tBudu+30tBu2dBu3(0tBudu+0tudBu)+3TBt=30tBu2dBu30tudBu+3T0t1dBu=0t3Bu23u+3TdBu. We can take ht=3Bt23t+3T.

        • (iii) Note that M0=E[eσBT|F0]=E[eσBT]=e12σ2T by (11.2) and the scaling properties of normal random variables. We showed in 11.6 that eσBt12σ2t was a martingale. Hence,

          E[eσBT|Ft]=E[eσBT12σ2Te12σ2T|Ft]=eσBt12σ2te12σ2T=eσBt12σ2(Tt). Applying Ito’s formula to Zt=eσBt12σ2(Tt) gives that

          dZt=(12σ2Zt+(0)(σZt)+12(12)(σ2Zt))dt+(1)(σZt)dBt=σZtdBt so we obtain that

          Zt=Z0+0tσZudBu.

          Substituting in for Zt we obtain

          E[eσBT|Ft]=e12σ2T+0tσZudBu

          so we can take ht=σZt=σeσBt12σ2(Tt).

    Chapter 14
    • 14.1 We have α(t,x)=2t, β=0 and Φ(x)=ex. By Lemma 14.1.2 the solution is given by

      F(t,x)=Et,x[eXT]

      where dXu=2udu+0dBu=du. This gives

      XT=XttT2udu=XtT2+t2. Hence,

      F(t,x)=Et,x[eXtt2]=E[exT2+t2]=exeT2+t2. Here we use that Xt=x under Et,x.

    • 14.2 We have α=β=1 and Φ(x)=x2. By Lemma 14.1.2 the solution is given by

      F(t,x)=Et,x[XT2]

      where dXu=du+dBu. This gives

      XT=Xt+tTdu+tTdBu=Xt+(Tt)+(BTBt). Hence,

      F(t,x)=Et,x[(Xt+(Tt)+(BTBt))2]=E[x2+(Tt)2+(BTBt)2+2x(Tt)+2x(BTBt)+2(Tt)(BTBt)]=x2+(Tt)2+(Tt)+2x(Tt). Here we use that Xt=x under Et,x and that BTBtBTtN(0,Tt).

    • 14.3

      • (a) We have Zt=F(t,Xt)+γ(t), where dXt=α(t,x)dt+β(t,x)dBt, so

        dZt=(Ft+γt(t)+α(t,x)Fx+12β(t,x)22Fx2)dt+β(t,x)FxdBt=β(t,x)FxdBt where as usual we have suppressed the (t,Xt) arguments of F and its partial derivatives. Note that the term in front of the dt is zero because F satisfies (14.8).

      • (b) We use the same strategy as in the proof of Lemma 14.1.2. Writing out dZt in integral form over time interval [t,T] we obtain

        ZT=Zt+tTβ(u,x)FxdBu.

        Taking expectations Et,x gives us

        Et,x[F(T,XT)+γ(T)]=Et,x[F(t,Xt)+γ(t)]

        and noting that Xt=x under Et,x, we have

        Ex,t[F(T,XT)]+γ(T)=F(t,x)+γ(t).

        Using (14.9) we have F(T,XT)=Φ(XT) and we obtain

        Ex,T[Φ(XT)]+γ(T)γ(t)=F(t,x)

        as required.

    • 14.4 Consider (for example) the stochastic process

      Mt={Bt for t1,BtBt2 for t>1.

      We now consider Mt at time 3. The intuition is that E[M3|F2] can see the value of B1, but that σ(M2)=σ(B2) and E[M3|σ(M2)] cannot see the value of B1.

      Formally: we have

      E[M3|F2]=E[B3B1|F2](C.4)=B2B1 Here, we use that Bt is a martingale. However,

      E2,M2[M3]=E[B3B1|σ(M2)]=E[B3B1|σ(B2B0)]=E[B3B1|σ(B2)]=E[B3B2|σ(B2)]+E[B2|σ(B2)]E[B1|σ(B2)]=E[B3B2]+B2E[B1|σ(B2)]=0+B2E[B1|σ(B2)](C.5)=B2E[B1|σ(B2)] If we can show that (C.4) and (C.5) are not equal, then we have that Mt is not Markov. Their difference is D=(C.4)(C.5)=E[B1|σ(B2)]B1.

      We can write B2=(B2B1)+(B1B0). By the properties of Brownian motion, B2B1 and B1B0 are independent and identically distributed. Hence, by symmetry, E[B2B1|σ(B2)]=E[B1B0|σ(B2)] and since

      B2=E[B2|σ(B2)]=E[B2B1|σ(B2)]+E[B1B0|σ(B2)] we have that E[B2B1|σ(B2)]=B22. Hence

      D=B22B1

      which is non-zero.

    Chapter 15
    • 15.1 We have

      er(Tt)EQ[Φ(ST)|Ft]=er(Tt)EQ[K]=Ker(Tt) because K is deterministic. By Theorem 15.3.1 this is the price of the contingent claim Φ(ST) at time t.

      We can hedge the contingent claim K simply by holding KerT cash at time 0, and then waiting until time T. While we wait, the cash increases in value according to (15.2) i.e. at continuous rate r. So, this answer is not surprising.

    • 15.2

      • (a) By Theorem 15.3.1 the price of the contingent claim Φ(ST)=log(ST) at time t is

        er(Tt)EQ[Φ(ST)|Ft]=er(Tt)EQ[log(ST)|Ft]. Recall that, under Q, St is a geometric Brownian motion, with S0=0, drift r and volatility σ. So from (15.20) we have ST=Ste(r12σ2)t+σBt. Hence,

        er(Tt)EQ[log(ST)|Ft]=er(Tt)EQ[log(St)+(r12σ2)(Tt)+σ(BTBt)|Ft]=er(Tt)(log(St)+(r12σ2)(Tt)+σ(EQ[BT|Ft]Bt))=er(Tt)(log(St)+(r12σ2)(Tt)). Here, we use that St,Bt are Ft measurable, and that (Bt) is a martingale.

        Putting in t=0 we obtain

        erT(log(S0)+(r12σ2)T),

      • (b) The strategy could be summarised as ‘write down Φ(ST), replace T with t=0 and hope’. The problem is that if we buy logs units of stock at time 0, then from (15.20) (with t=0) its value at time T will be (logs)exp((r12σ2)T+σBT), which is not equal to logST=logs+(r12σ2)T+σBT.

        In more formal terminology, the problem is that our excitable mathematician has assumed, incorrectly, that the log function and the ‘find the price’ function commute with each other.

    • 15.3

      • (a) We have (in the risk neutral-world Q) that dSt=rStdt+σStdBt. Hence, by Ito’s formula,

        dYt=((0)+rSt(βStβ1)+12σ2St2(β(β1)Stβ2))dt+σSt(βStβ1)dBt=(rβ+12σ2β(β1))Ytdt+(σβ)YtdBt. So Yt is a geometric Brownian motion with drift rβ+12σ2β(β1) and volatility σβ.

      • (b) Applying (15.20) and replacing the drift and volatility with those from part (a), we have that

        YT=Ytexp((rβ+12σ2β(β1)12σ2β2)(Tt)+σβ(BTBt))=Ytexp((rβ12σ2β)(Tt)+σβ(BTBt)). By Theorem 15.3.1, the arbitrage free price of the contingent claim Yt=Φ(ST) at time t is

        er(Tt)EQ[YT|Ft]=er(Tt)EQ[Stβexp((rβ12σ2β)(Tt)+σβ(BTBt))|Ft]=er(Tt)Stβe(rβ12σ2β)(Tt)EQ[eσβ(BTBt)|Ft]=er(Tt)Stβe(rβ12σ2β)(Tt)EQ[eσβ(BTBt)]=er(Tt)Stβe(rβ12σ2β)(Tt)e12σ2β2(Tt)=Stβer(Tt)(1β)12σ2β(Tt)(1β). Here, we use that St is Ft measurable. We then use (11.2) along with the properties of Brownian motion to tell us that σβ(BTBt) is independent of Ft with distribution N(0,σ2β2(Tt)).

    • 15.4 From Theorem 15.3.1 the price at time t of the binary option is

      er(Tt)EQ[Φ(ST)|Ft]=Ker(Tt)EQ[𝟙{ST[α,β]}|Ft] We have

      ST=Ste(r12σ2)(Tt)+σ(BTBt)=Ste(r12σ2)(Tt)+σTtZ where ZN(0,1) is independent of Ft. Here we use that BTBtN(0,Tt)TtN(0,1) is independent of Ft. Therefore,

      EQ[𝟙{ST[α,β]}|Ft]=EQ[𝟙{Ste(r12σ2)(Tt)+σTtZ[α,β]}|Ft]=EQ[𝟙{log(αSt)(r+12σ2)(Tt)σTtZlog(βSt)(r+12σ2)tσTt}|Ft]=EQ[𝟙{log(αSt)(r+12σ2)(Tt)σTtZlog(βSt)(r+12σ2)tσTt}]=N(e2)N(e1) where

      e1=log(αSt)(r+12σ2)(Tt)σTte2=log(βSt)(r+12σ2)(Tt)σTt. Hence the price at time t is given by Ker(Tt)[N(e1)N(e2)].

    • 15.5 If we take Ft=eμt then (in world P) we have

      StFt=S0exp((μ12σ2)t+σBt)exp(μt)=S0exp(σBt12σ2t),

      which we showed was a martingale in 11.6.

    • 15.6 We have

      ft=0,fs=c,2fx2=0

      which, put into (15.10), gives (0)+rs(c)+12s2σ2(0)r(cs)=0. Similarly,

      gt=rcert,gs=0,2gx2=0

      which, put into (15.10), gives (rcert)+rs(0)+12s2σ2(0)r(cert)=0.

      The solution f corresponds to a portfolio containing c units of stock, with value F(t,St)=cSt at time t. The solution g corresponds to a portfolio that starts at time 0 with c units of cash, with value F(t,St)=cert at time t.

    • 15.7 We have the risk neutral pricing formula Πt(Φ)=er(Tt)Et,StQ[Φ(ST)] for any contingent claim Φ. Hence, using linearity of EQ,

      Πt(αΦ1+βΦ2)=er(Tt)Et,StQ[αΦ1(ST)+βΦ2(ST)]=αer(Tt)Et,StQ[Φ1(ST)]+βer(Tt)Et,StQ[Φ2(ST)]=αΠt(Φ1)+βΠt(Φ2) as required.

    • 15.8 In the ‘new’ model we could indeed buy a unit of Yt=Stβ and hold onto it for as long as we liked. But in the ‘old’ model we can only buy (linear multiples of) the stock St; we can’t buy a commodity whose price at time t is Stβ. This means that the hedging strategy suggested within the ‘new’ model doesn’t work within the ‘old’ model. Consequently, there is no reason to expect that prices in the two models will be equal. In general they will not be.

    Chapter 16
    • 16.1

      • (a) The functions Φcash(ST)=1, Φstock(ST)=ST, Φcall(ST)=max(STK,0) and Φput(ST)=max(KST,0) look like:

        (A plot of the payoff functions)

      • (b) In terms of functions, the put-call parity relation states that

        max(KST,0)=max(STK)+KST.

        To check that this holds we consider two cases.

        • If KST then put-call parity states that 0=STK+KST, which is true.

        • If KST then put-call parity states that KST=0+KST, which is true.

    • 16.2 The put-call parity relation (16.1) says that

      Φput(ST)=Φcall(ST)+KΦcash(ST)Φstock(ST).

      Hence,

      Πtput=Πtcall+KΠtcashΠtstock.

      From here, using that Πtcash=er(Tt) (corresponding to Φ(ST)=1) and Πtstock=St (corresponding to Φ(ST)=ST), as well as the Black-Scholes formula (15.23), we have

      Πtput=StN[d1]Ker(Tt)N[d2]+Ker(Tt)St=St(N[d1]1)Ker(Tt)(N[d2]1)=StN[d1]+Ker(Tt)N[d2]. In the last line we use that N[x]+N[x]=1, which follows from the fact that the N(0,1) distribution is symmetric about 0 (i.e. P[Xx]+P[Xx]=P[Xx]+P[Xx]=P[Xx]+P[Xx]=1). The formula stated in the question follows from setting t=0.

    • 16.3 Write Φcall,K(ST)=max(STK,0) for the contingent claim of a European call option, and Φput,K(ST) for the contingent claim of a European put options, both with strike price K and exercise date T. Then

      Φ(ST)=Φcall,1(ST)+Φput,1(ST)

      so we can hedge Φ(ST) by holding a single call option with strike price 1 and a single put option with strike price 1.

    • 16.4

      • (a) A sketch of Φ(ST) for general A looks like

        (A plot of the payoff function)

        Let Φput,K(ST)=max(KST,0) denote the contingent claim corresponding to a European put option with strike price K and exercise date T. Then, we claim that

        (C.6)Φ(ST)=Φput,K+A(ST)Φput,A(ST).

        The right-hand side of this equation corresponds to a portfolio of one put option with strike price K+A and minus one put option with strike price A. We can see that the relation (C.6) holds by using a diagram

        (A plot of the payoff functions)

        in which the purple line Φput,A(ST) is subtracted from the red line Φput,K+A(ST) to obtain the green line Φ(ST). Alternatively, we can check it by considering three cases:

        • If STA then we have K=(K+AST)(AST) which is true.

        • If ASTK+A then we have K+AST=(K+AST)(0) which is true.

        • If STK+A then we have 0=(0)(0) which is true.

        Therefore, the portfolio of one put option with strike price K+A and minus one put option with strike price A is a replicating portfolio for Φ(ST).

      • (b) We use put-call parity to replicate our portfolio of put options with a portfolio of cash, stock and call options. This tells us that Φput,K(ST)=Φcall,K(ST)+KΦcash(ST)+Φstock(ST). Therefore, replacing each of our put options with the equivalent amounts of cash and stock, we have

        Φput,K+A(ST)Φput,A(ST)=Φcall,K+A(ST)Φcall,A(ST)+(K+AA)Φcash(ST)(11)Φstock(ST)=Φcall,K+A(ST)Φcall,A(ST)+KΦcash(ST). Hence, Φ(ST) can be replicated with a portfolio containing one call option with strike price K+A, minus one call option with strike price K, and K units of cash.

      • (c) Corollary 15.3.3 refers to portfolios consisting only of stocks and cash. It does not apply to portfolios that are also allowed to contain derivatives.

    • 16.5 We can write

      Φbull(ST)=A+max(STA,0)max(STB,0).

      It may help to draw a picture, in the style of 16.4. If we write Φcall,K(ST)=max(STK,0) and recall that Φcash(ST)=1 then we have

      Φbull(ST)=AΦcash(ST)+Φcall,A(ST)Φcall,B(ST).

      So, in our (constant) hedging portfolio at time 0 we will need AerT in cash, one call option with strike price A, and minus one call option with strike price B.

    • 16.9 The price computed for the contingent claim Φ(ST)=Stβ in 15.3 is

      F(t,St)=Stβexp(r(Tt)(1β)12σ2β(Tt)(1β)).

      (Recall that Theorem 15.3.1 also tells us that there is a hedging strategy ht=(xt,yt) for the contingent claim with value F(t,St).)

      We have F(t,s)=sβexp(r(Tt)(1β)12σ2β(Tt)(1β)) and we calculate

      ΔF=Fs(t,St)=βStβ1exp(r(Tt)(1β)12σ2β(Tt)(1β))ΓF=2Fs2(t,St)=β(β1)Stβ2exp(r(Tt)(1β)12σ2β(Tt)(1β))ΘF=Ft(t,St)=Stβ(1β)(r+12βσ2)exp(r(Tt)(1β)12σ2β(Tt)(1β))ρF=Fr(t,St)=Stβ(Tt)(1β)exp(r(Tt)(1β)12σ2β(Tt)(1β))VF=Fσ(t,St)=Stβσβ(Tt)(1β)exp(r(Tt)(1β)12σ2β(Tt)(1β)).

    • 16.10

      • (a) The underlying stock St has ΔS=1 and ΓS=0. If we add 2 stock into our original portfolio with value F, then its new value is V(t,St)=F(t,St)2St, which satisfies ΔV=ΔF2=0.

        The cost of adding 2 units of stock into the portfolio is 2St.

      • (b) After including an amount wt of stock and an amount dt of D, we have

        V(t,St)=F(t,St)+wtSt+dtD(t,St).

        Hence, we require that

        0=ΔF+wt+dtΔD=2+wt+dt,0=ΓF+dtΓD=3+2dt. The solution is dt=32 and wt=12.

        The cost of the extra stock and units of D that we have had to include is 32D(t,St)12St.

    • 16.11

      • (a) It doesn’t work because adding in an amount zt of Z in the second step will (typically i.e. if ΔZ0) destroy the delta neutrality that we gained from the first step.

      • (b) Since Z(t,St)=St we have

        V(t,St)=F(t,St)+wtW(t,St)+ztSt,

        so we require that

        0=Vs=ΔF+wtΔW+zt0=2Vs2=ΓF+wtΓW. The solution is easily seen to be

        wt=ΓFΓWzt=ΔWΓFΓWΔF

      • (c) This idea works. As we can see from part (b), if we use stock as our financial derivative Z(t,St)=St, then ΓZ=0. Hence, adding in a suitable amount of stock in the second step can achieve delta neutrality without destroying the gamma neutrality obtained in the first step.

    • 16.12 Omitted (good luck).

    Chapter 19
    • 19.1 The degrees of the nodes, in alphabetical order, are (1,1),(0,2),(2,2),(2,0),(1,1),(1,2),(1,0). This gives a degree distribution

      P[DG=(a,b)]={17 for (a,b){(0,2),(2,2),(2,0),(1,2),(1,0)},27 for (a,b)=(1,1).

      Sampling a uniformly random and moving along it means that the chance of ending up at a node A is proportional to degin(A). Since the graph has 8 edges, we obtain

      P[O=n]={18+28 for n=0 (nodes D and G)18+18 for n=1 (nodes A and E)0+28+18 for n=2 (nodes B,C and F)={38 for n=014 for n=138 for n=2

    • 19.2 Node Y can fail only if the cascade of defaults includes XDY. Given that X fails, the probability that X fails is η3=13, because D has 3 in-edges. Similarly, given that D fails, the probability that Y also fails is 12. Hence, the probability that Y fails, given that X fails, is 16.

    • 19.3 Node Y can fail if the cascade of defaults includes XBDY or XBCDY. Given that X fails, node B is certain to fail as well, since B has only one in-edge. Given that B fails, C is certain to fail for the same reason. Therefore, the edges (B,D) and (C,D) are both certain to default. For each of these edges, independently, there is a chance 13 that their own default causes D to fail. The probability that D fails is therefore

      13+(113)×13=59.

      Here, we condition first on if the link BD causes D to fail (which it does with probability 13) and then, if it doesn’t (which has probability 113) we ask if the link BC, which fails automatically and causes failure of CD, causes D to fail.

      Given that D fails, Y is certain to fail. Hence, the probability that Y fails, given that X fails, is 59.

    • 19.4 For any node of the graph (except for the root node), if its single incoming loan defaults, then its own probability of default if α, independently of all else. Hence, each newly defaulted loan leads to two further defaulted loans with probability α, and leads to no further defaulted loans with probability 1α. Hence, the defaulted loans form a Galton-Watson process (Zn) with off-spring distribution G, given by P[G=2]=α and P[G=0]=1α, with initial state Z0=2 (representing the two loans which initially default when V0 defaults).

      The total number of defaulted edges is given by

      S=n=0Zn.

      Combining Lemmas 7.4.7, 7.4.6 and 7.4.8, we know that either:

      • If E[G]>1 then there is positive probability that Zn as n; in this case for all large enough n we have Zn1, and hence S=.

      • If E[G]1 then, almost surely, for all large enough n we have Zn=0, which means that S<.

      We have E[G]=2α. It is clear that the number of defaulted banks is infinite if and only if the number of defaulted loans in infinite, which has positive probability if E[G]>1. So, we conclude that there is positive probability of a catastrophic default if and only if α>12.

    • 19.5

      • (a) The probability that both B and C fail is 1212=14. Given this event, the probability that both D and E fail is 12(1343434)=37128; here the first term 12 is the probability of D failing (via its link to B) and the second term is one minus the probability of E not failing despite all its inbound loans defaulting. Given that D and E both fail, the probability that F fails is 12323=59.

        Hence, the probability that every node fails is

        143712859=1854608.

      • (b) Our strategy comes in three stages: first work out the probabilities of all the possible outcomes relating to B and C; secondly, do the same for D and E; finally, do the same for F. Each stage relies on the information obtained in the previous stage.

        Stage 1: The probability that B fails and C does not is 1212=14. This is also the probability that C fails and B does not. The probability that both B and C fail is also 14.

        Stage 2: Hence, the probability that E fails and D does not is

        14(1412)+14(14)+14(12)(14+3414)=19128

        The three terms in the above correspond respectively to the three cases considered in the first paragraph. Similarly, the probability that D fails and E does not is

        14(123434)+14(0)+14(12)(343434)=63512

        and the probability that both D and E fail is

        14(12)(14+3414)+14(0)+14(12)(1343434)=65512

        .

        Stage 3: Finally, considering these three cases in turn, the probability that F fails is

        19128(13)+63512(13)+65512(12323)=3712304.