arxiv: v5 [cs.ce] 3 Dec 2008

Transcription

1 PRICING AMERICAN OPTIONS FOR JUMP DIFFUSIONS BY ITERATING OPTIMAL STOPPING PROBLEMS FOR DIFFUSIONS ERHAN BAYRAKTAR AND HAO XING arxiv: v5 [cs.ce] 3 Dec 2008 Abstract. We approximate the price of the America put for jump diffusios by a sequece of fuctios, which are computed iteratively. This sequece coverges to the price fuctio uiformly ad expoetially fast. Each elemet of the approximatig sequece solves a optimal stoppig problem for geometric Browia motio, ad ca be umerically computed usig the classical fiite differece methods. We prove the covergece of this umerical scheme ad preset examples to illustrate its performace. 1. Itroductio Jump diffusio models are heavily used i modelig stock prices sice they ca capture the excess kurtosis ad skewess of the stock price returs, ad they ca produce the smile i the implied volatility curve (see [10]). Two well-kow examples of these models are i) the model of [18], i which the jump sizes are log-ormally distributed, ad ii) the model of [17], i which the logarithm of jump sizes have the so called double expoetial distributio. Based o the results of [6] we propose a umerical algorithm to calculate the America optio prices for jump diffusio models ad aalyze the covergece behavior of this algorithm. As observed by [6], we ca costruct a icreasig sequece of fuctios, which are value fuctios of optimal stoppig problems (see (2.8) ad also (2.11)), that coverge to the price fuctio of the America put optio uiformly ad expoetially fast. Because each elemet of this sequece solves a optimal stoppig problem it shares the same regularity properties, such as covexity ad smoothess, with the origial price fuctio. Eve the correspodig free boudaries have the same smoothess properties (whe they have a discotiuity, which ca oly happe at maturity, the magitude of the discotiuity is the same; see (2.13)). Therefore, the elemets i this approximatig sequece provide a good imitatio to the value fuctio besides beig close to it umerically (see Remark 2.1). O the other had, each of these fuctios ca be represeted as classical solutios of free boudary problems (see (2.9)) for geometric Browia motio, ad therefore ca be implemeted usig classical fiite differece methods. We build a iterative umerical algorithm based o discretizig these free boudary problems (see (3.10)). Whe the mesh sizes are fixed, we show that the iterative sequece we costructed is mootoous ad coverges uiformly ad expoetially fast (see Propositio 3.3). We also show, i a rather direct way, that whe the mesh sizes go to zero our algorithm coverges to the true price fuctio (see Propositio 3.4). The pricig i the cotext of jump models is difficult sice the prices of optios satisfy itegro-partial differetial equatios (itegro-pdes), i.e. they have o-local itegral terms, ad the usual fiite-differece methods are ot directly applicable because the itegral term leads to full matrices. Recetly there has bee a lot of iterest i developig umerical algorithms for pricig i jump models, see e.g. [1], [2], [4], [9], [11], [13], [14], [16], [17], [19], [23], amog them [2], [9], [13] ad [14] treated specific or geeral jump models with ifiite activity jumps. These algorithms have bee extesively discussed i Chapter 12 of [10]. I this paper, relyig o the results of [6] as Key words ad phrases. Pricig derivatives, America optios, jump diffusios, barrier optios, fiite differece methods. Departmet of Mathematics, Uiversity of Michiga, A Arbor, MI 48109, USA; {erha,haoxig}@umich.edu. This research is supported i part by the Natioal Sciece Foudatio. 1

2 2 ERHAN BAYRAKTAR AND HAO XING described above, we give a efficiet umerical algorithm (ad aalyze its error versus accuracy characteristics) to efficietly compute America optio prices for jump diffusio models with fiite activity. Oe ca hadle ifiite activity models by icreasig the volatility coefficiet appropriately as suggested o p. 417 of [10]. A ideal umerical algorithm, which is most ofte a iterative scheme, *should mootoically coverge to the true price uiformly (across time ad space) ad expoetially fast*, that is, the error bouds should be very tight. This is the oly way oe ca be sure that the price output of the algorithm is close to the true price after a reasoable amout of rutime ad without havig to compare the price obtaied from the algorithm to other algorithms output. It is also desirable to obtai a scheme **that does ot deviate from the umerical pricig schemes, such as fiite differece methods, that were developed for models that do ot accout for jumps**. Fiacial egieers workig i the idustry are already familiar with fiite differece schemes such as projected successive over relaxatio, PSOR, (see e.g. [20]) ad Brea-Schwartz algorithm (see [8] ad [15]) to solve the partial differetial equatios associated with free boudary problems, but may ot be familiar with the itricacies ivolved i solvig itegro-partial differetial equatios developed i the literature. It would be ideal for them if they could use what they already kow with oly a slight modificatio to solve for the prices i a jump diffusio model. I this paper, we develop a algorithm which establishes both * ad **. I Sectio 4, we will ame this algorithm, depedig o which classical method we use to solve the sparse liear systems i (3.10), as either Iterated PSOR or Iterated Brea-Schwartz. I the ext sectio we itroduce a sequece of optimal stoppig problems that approximate the price fuctio of the America optios, ad discuss their properties. I Sectio 3, we itroduce a umerical algorithm ad aalyze its covergece properties. I the last sectio we give umerical examples to illustrate the competitiveess of our algorithm ad price America, Barrier ad Europea optios for the models of [17] ad [18]. 2. A Sequece of Optimal Stoppig Problems for Geometric Browia Motio Approximatig the America Optio Price For Jump Diffusios We will cosider a jump diffusio model for the stock price S t with S 0 = S, ad assume that retur process X t := log(s t /S), uder the risk eutral measure, is give by dx t = (µ 12 ) N t σ2 dt+σdw t + Z i, X 0 = 0. (2.1) I (2.1), µ = r + λ λξ, r is the risk-free rate, W t is a Browia motio, N t is a Poisso process with rate λ idepedet of the Browia motio, Z i are idepedet ad idetically distributed, ad come from a commo distributio F o R, that satisfies ξ := R ez F(dz) <. The last coditio guaratees that the stock prices have fiite expectatio. We will assume that the volatility σ is strictly positive. The price fuctio of the America put with strike price K is i=1 V(S,t) := sup τ S t,t E{e r(τ t) (K S τ ) + St = S}, (2.2) i which S t,t is the set of stoppig times of the filtratio geerated by X that belog to the iterval [t,t] (t it the curret time, T is the maturity of the optio). Istead of workig with the pricig fuctio V directly, which is the uique classical solutio of the followig itegro-differetial free boudary problem (see Theorem 3.1 of [6]) t V(S,t)+AV(S,t)+λ V(S,t) = K S, S s(t), V(S,T) = (K S) +, R V(e z S,t)F(dz) (r +λ) V(S,t) = 0 S > s(t), (2.3)

3 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 3 i which, A is the differetial operator A := 1 2 σ2 S 2 d2 ds 2 +µs d ds, (2.4) ad t s(t), t [0,T], is the exercise boudary that eeds to be determied alog with the pricig fuctio V; we will costruct a sequece of pricig problems for the geometric Browia motio dst 0 = µs0 t dt+σs0 t dw t, S0 0 = S. (2.5) To this ed, let us itroduce a fuctioal operator J, whose actio o a test fuctio f : R + [0,T] R + is the solutio of the followig pricig problem for the geometric Browia motio: (St) 0 t 0 { τ } Jf(S,t) = sup E e (r+λ)(u t) λ Pf(Su 0,u)du+e (r+λ)(τ t) (K Sτ 0 )+ S 0 t = S, (2.6) τ S t,t t i which Pf(S,u) = R f(e z S,u)F(dz) = E[f(e Z S,u)], S 0, (2.7) for a radom variable Z whose distributio is F, ad S t,t is the set of stoppig times of the filtratio geerated by W that take values i [t,t]. Let us defie a sequece of pricig fuctios by v 0 (S,t) = (K S) +, v +1 (S,t) = Jv (S,t), 0, for all (S,t) R + [0,T]. (2.8) For each 1, the pricig fuctio v is the uique solutio of the classical free-boudary problem (istead of a free boudary problem with a itegro-diffetial equatio) t v (S,t)+Av (S,t) (r+λ) v (S,t) = λ (Pv 1 )(S,t), S > s (t), v (S,t) = K S, S s (t), v (S,T) = (K S) +, i which t s (t) is the free-boudary (the optimal exercise boudary) which eeds to be determied (see Lemma 3.5 of [6]). Now startig from v 0, we ca calculate {v } 0 sequetially. For v, the solutio of (2.9) ca be determied usig a classical fiite differece method (we use the Crak-Nicolso discretizatio alog with Bera-Schwartz algorithm or PSOR i the the followig sectios) give that the fuctio v 1 is available. The term o the right-had-side of (2.9) ca be computed either usig Mote-Carlo or a umerical itegrator (we use the umerical itegratio with the Fast Fourier Trasformatio (FFT) i our examples). Iteratig the solutio for (2.9) a few times we are able to obtai the America optio price V accurately sice the sequece of fuctios {v } 0 coverges to V uiformly ad expoetially fast: ( v (S,t) V(S,t) v (S,t)+K 1 e (r+λ)(t t)) ( ) λ, S R +, t (0,T), (2.10) λ+r see Remark 3.3 of [6]. Note that the usual values of T for the traded optios is 0.25, 0.5, 0.75, 1 year. Remark 2.1. The approximatig sequece {v } 0 goes beyod approximatig the value fuctio V. Each v ad its correspodig free boudary have the same regularity properties which V ad its correspodig free boudary have. I a sese, for large eough, v provides a good imitatio of V. Below we list these properties: 1) The fuctio v ca be writte as the value fuctio of a optimal stoppig problem: (2.9) v (S,t) := sup τ S t,t E{e r(τ σ t) (K S τ σ ) + St = S}, (2.11) i which σ is the -th jump time of the Poisso process N t.

4 4 ERHAN BAYRAKTAR AND HAO XING 2) Each v is a covex fuctio i the S-variable, which is a property that is also shared by V. Moreover, the sequece {v } 0 is a mootoe icreasig sequece covergig to the value fuctio V (see (2.10)). 3) The free boudaries s(t) ad s (t) have the same regularity properties (see [7] ad its refereces): a) They are strictly decreasig. b) They may exhibit discotiuity at T: If the parameters satisfy r < λ (e z 1)F(dz), (2.12) R + we have lim s(t) = lim s (t) = S < K, 1, (2.13) t T t T where S is the uique solutio of the followig itegral equatio [ ] rk +λ (K Se z ) + (K Se z ) F(dz) = 0. (2.14) R We will see such a example i Sectio 4, where the equatio (2.14) ca be solved aalytically for some jump distributio F. c) Both s(t) ad s (t) are cotiuously differetiable o [0,T). 3. A Numerical Algorithm ad its covergece aalysis 3.1. The umerical algorithm. I this sectio, we will discretize the algorithm itroduced i the last sectio ad give more details. For the coveiece of the umerical calculatio, we will first chage the variable: x logs, x(t) logs(t) ad u(x,t) V(S,t). u satisfies the followig itegro-differetial free boudary problem t u+ 1 ( 2 σ2 2 x 2u+ µ 1 ) 2 σ2 u (r +λ)u+λ (Iu)(x,t) = 0, x > x(t) x u(x,t) = K e x, x x(t) (3.1) i which u(x,t) = (K e x ) +, (Iu)(x,t) = R u(x + z, t)ρ(z)dz, (3.2) with ρ(z) as the desity of the distributio F. Similarly, u (x,t) v (S,t) satisfies the similar free boudary problem where u i (3.1) is replaced by u i differetial parts ad by u 1 i the itegral part. I additio, it was show i Theorem 4.2 of [21] that the free boudary problem (3.1) is equivalet to the followig variatioal iequality L D u(x,t)+λ (Iu)(x,t) 0 u(x,t) g(x) [L D u(x,t)+λ (Iu)(x,t)] [u(x,t) g(x)] = 0, (x,t) R [0,T], (3.3) i which L D u t u+ 1 ( 2 σ2 2 x 2u+ µ 1 ) 2 σ2 u (r +λ)u x g(x) = (K e x ) +.

5 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 5 Sice the secod spacial derivative of u does ot exist alog the free boudary x(t), the variatioal iequality (3.3) does ot have a classical solutio. However, Theorem 3.2 of [21] showed that u is the solutio of (3.3) i the Sobolev sese. I the same sese, u (x,t) satisfies a similar variatioal iequality L D u (x,t)+λ (Iu 1 )(x,t) 0 u (x,t) g(x) [L D u (x,t)+λ (Iu 1 )(x,t)] [u (x,t) g(x)] = 0, (x,t) R [0,T]. Let us discretize (3.3) usig Crak-Nicolso scheme. For fixed t, x, x mi ad x max, let M t = T ad L x = x max x mi. Let us deote x l = x mi + l x, l = 0,,L. By ũ l,m we will deote the solutio of the followig differece equatio θp ũ l 1,m +(1+θp 0 )ũ l,m θp + ũ l+1,m b l,m 0 ũ l,m g l [ θp ũ l 1,m +(1+θp 0 )ũ l,m θp + ũ l+1,m b l,m] [ũ l,m g l] = 0, for m = M 1,,0, l = 0,,L, satisfyig the termial coditio ũ l,m = g l = (K e x l ) + ad Dirichlet boudary coditios. θ is the weight factor. Whe θ = 1, the scheme (3.5) is the completely implicit Euler scheme; whe θ = 1/2, it is the classical Crak-Nicolso scheme. The coefficiets p, p + ad p 0 are give by p = 1 t 2 σ2 ( x) 2 1 ( µ 1 ) t 2 2 σ2 x, The term b is defied by p + = 1 2 σ2 t ( x) p 0 = p +p + +(r+λ) t. ( µ 1 2 σ2 ) t x, bl,m = (1 θ)p ũ l 1,m+1 +(1 (1 θ)p 0 )ũ l,m+1 +(1 θ)p + ũ l+1,m+1 +λ t (3.4) (3.5) (3.6) [ (1 θ)(ĩũ)l,m+1 +θ(ĩũ)l,m]. (3.7) Ĩ i (3.7) is the discrete versio of the covolutio operator I i (3.2). It will be coveiet to approximate this covolutio itegral usig Fast Fourier Trasformatio (FFT). Discretizig a sufficietly large iterval [z mi,z max ] ito J sub-itervals. For the coveiece of the FFT, we will choose these J sub-itervals equally spaced, such that J z = z max z mi. We also choose x = α z, where α is a positive iteger, so that the umerical itegral may have fier grid tha the grid i x. Let z j = z mi +j z, j = 0,,J. Ĩ is defied by ) l,m J 1 (Ĩũ = ũ iterp (x l +z j,m t)ρ(z j ) z, (3.8) j=0 i which the value of ũ iterp is determied by the liear iterpolatio ũ. That is if there is some l satisfyig the x l x l +z j x l +1, ũ iterp (x l +z j,m t) = (1 w)ũ l,m +wũ l +1,m, for some w [0,1]. O the other had, if x l + z j is outside the iterval [x mi,x max ], the value of ũ iterp is determied by the boudary coditios. Moreover, i (3.8) we also assume ρ(z j ) 0, for all j, ad J 1 ρ(z j ) 1. (3.9) j=0

6 6 ERHAN BAYRAKTAR AND HAO XING Now (3.8) ca be calculated usig FFT. See Sectio 6.1 i [2] for implemetatio details. Note that umerically solvig the system (3.5) is difficult due to the cotributio of the itegral term Ĩũ. Therefore, followig the results i Sectio 2, we will discretize (3.4) recursively (usig the Crak-Nicoslo scheme) to obtai the sequece {ũ } 0 recursively. Let ũ l,m 0 = g l. For 1, ũ is defied recursively by θp ũ l 1,m +(1+θp 0 )ũ l,m θp +ũ l+1,m b l,m 0 ũ l,m g l [ θp ũ l 1,m +(1+θp 0 )ũ l,m θp + ũ l+1,m ] b l,m [ũ l,m g l] = 0, (3.10) with the termial coditio ũ l,m = g l ad Dirichlet boudary coditios. Similar to (3.7), b is defied by bl,m =(1 θ)p ũ l 1,m+1 +λ t +(1 (1 θ)p 0 )ũ l,m+1 [ (1 θ)(ĩũ 1) l,m+1 +θ(ĩũ 1) l,m]. +(1 θ)p + ũ l+1,m+1 (3.11) For each, we will solve the sparse liear system of equatios (3.10) usig the projected PSOR method (see eg. [20]). Remark 3.1. We will iterate (3.10) to approximate the solutio of (3.5), which ca be see as a global fixed poit iteratio algorithm. This global fixed poit algorithm is differet from the local fixed poit algorithm i [11], where d Hallui et al. implemeted the Crak Nicolso time steppig of a o-liear itegro-partial differetial equatio comig from a alterative represetatio (due to the pealty method) of the America optio price fuctio. Also see [12] for the case of Europea optios. Note that discretizig the o-liear PDE that arises from the pealized formulatio itroduces a extra error. We work with the variatioal formulatio directly. Each ũ approximates u, which itself is the value fuctio of a optimal stoppig problem, ad as we have discussed i Remark 2.1 provides a good imitatio of the America optio price fuctio. Each of these iteratios provide strictly decreasig free boudary curves with the same regularity ad jump properties as the free boudary curve for the America optio price fuctio, see Remarks 2.1 ad 4.2. The approximatig sequece i [11] does ot carry the same meaig, it is a techical step to carry out the Crak Nicolso time steppig of their o-liear itegro-pde Covergece of the Numerical Algorithm. I the followig, we will show the covergece of the umerical algorithm for the completely implicit Euler scheme (θ = 1). We first show that {ũ } 0 is a mootoe icreasig sequece. Extra care has to be give to make the approximatig sequece mootoe i the pealty formulatio of [11] (see Remark 4.3 o page 341), but the mootoicity comes out aturally i our formulatio. Next, we prove that the sequece {ũ } 0 is uiformly bouded above by the strike price K ad coverges to ũ at a expoetial rate. At last, we will argue that as the mesh sizes x ad t go to zero ũ coverges to the America optio value fuctio u. I the followig four propositios, we let t ad x to be sufficietly small so that costats p ad p + defied i (3.6) are positive. Propositio 3.1. The sequece {ũ } 0 is a mootoe icreasig sequece.

7 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 7 Proof. Whe θ = 1, subtractig the third equality for -th iteratio i (3.10) from the equality for ( + 1)-th iteratio, we obtai [ ][ ] p ũ l 1,m +(1+p 0 )ũ l,m p + ũ l+1,m b l,m ũ l,m +1 ũl,m ( ) ( ( ) + { p ũ l 1,m +1 ũ l 1,m +(1+p 0 ) ũ l,m +1 ũl,m ) p + ũ l+1,m +1 ũ l+1,m (3.12) ( ) ) } l,m [ ũ l,m+1 +1 ũ l,m+1 λ t (Ĩ(ũ ũ 1 ) ũ l,m +1 gl] = 0. i which we used the liearity of the operator Ĩ. Let us defie the vectors ( ) T e m +1 = ũ 0,m +1 ũ0,m,,ũl,m +1 ũl,m, ([ ) ] 0,m [ f+1 m = (ũ 0,m+1 +1 ũ 0,m+1 )+λ t (Ĩ(ũ ũ 1 ) ũ 0,m +1 g0],, Equatio (3.12) ca be represeted as i which the matrix A s etries are a l,j = [ ) ] L,m [ T (ũ L,m+1 +1 ũ L,m+1 )+λ t (Ĩ(ũ ũ 1 ) ũ L,m +1 gl]). Ae m +1 = fm +1, (3.13) p ( ũ l,m +1 gl ) ( ( (1+p 0 ) ũ l,m +1 )+ gl p ũ l 1,m +(1+p 0 )ũ l,m p + ũ l+1,m p + ( ũ l,m +1 gl ) ) b l,m j = l 1 j = l j = l+1 0 others. O the other had, usig the first ad secod iequalities i (3.10) ad the fact that p ad p + are positive, we see that A is a M-matrix, i.e. A has positive diagoals, o-positive off-diagoals ad the row sums are positive. As a result all etries of A 1 are oegative. Now we ca prove the propositio by iductio. Note that ũ 1 ũ 0 = g, as a result of the secod iequality i (3.10) ad the defiitio of ũ 0. Assumig ũ ũ 1, we will show that ũ +1 ũ, i.e. ũ l,m +1 ũl,m 0 for all l ad m, i the followig. First, the termial coditio of ũ gives us ũ l,m +1 ũl,m = 0. Secod, ) l,m (Ĩ(ũ ũ 1 ) is oegative from the assumptio (3.9). Assumig ũ l,m+1 +1 ũ l,m+1 oegative, we have f+1 m i (3.13) as a oegative vector. Combiig with the fact that all etries of A 1 are oegative, the oegativity of ũ l,m +1 ũl,m follows from multiplyig A 1 o both sides of (3.13). The the result follows from a iductio m. Propositio 3.2. {ũ } 0 are uiformly bouded above by the strike price K. Proof. Whe θ = 1, i the third equality of (3.10), there are some (l,m) such that ũ l,m = g l. Otherwise we have ) l,m. (1+p 0 )ũ l,m = p ũ l 1,m +p + ũ l+1,m +ũ l,m+1 +λ t(ĩũ 1 However, i both cases, we obtai the followig iequality (1+p 0 ) ũ l,m p B m +p +B m +Bm+1 +λ tb 1 +r tk, 0 l L,0 m M 1, (3.14) i which we defie B (max m = l ũ l,m K, ) B = max m Bm.

8 8 ERHAN BAYRAKTAR AND HAO XING Note that the right had side of (3.14) is idepedet of l. Moreover, (1+p 0 )K is also less tha or equal to the right had side of (3.14). Therefore, (3.14) gives us (1+(r +λ) t)b m Bm+1 +λ tb 1 +r tk. (3.15) Give B m+1 K ad B 1 K, it clear from (3.15) that B m K. Now the propositio follows from double iductio o m ad with iitial steps ũ M = g K ad ũ 0 = g K. As a result of Propositios 3.1, we ca defie ũ l,m = lim + ũl,m, 0 l L,0 m M. (3.16) It follows from Propositio 3.2 that ũ l,m K. Lettig go to +, we ca see from (3.10) that ũ satisfies the differece equatio (3.5). Therefore, I the followig, we will study the covergece rate of {ũ } 0. Propositio 3.3. ũ coverges to ũ uiformly ad max l,m where η = 1 1+(λ+r) t (0,1), K is a positive costat. ũ = ũ. (3.17) (ũl,m ũ l,m ) ( 1 η M ) ( λ λ+r Proof. Let us defie e l,m = ũ l,m ũ l,m, Em = max e l,m, E = max l m Em. Propositio 3.1 ad (3.17) esure that e l,m is oegative. Moreover e satisfies [ ] p ũ l 1,m +(1+p 0 )ũ l,m p +ũ l+1,m b l,m e l,m + { p e l 1,m +(1+p 0 )e l,m p +e l+1,m e l,m+1 λ t ) K, (3.18) (Ĩe l,m 1) l,m } [ũl,m g l] = 0. We ca drop the first term o the left-had-side of (3.19) because of the first iequality i (3.10) ad e l,m oegative. It gives us the iequality (3.19) beig (1+p 0 )e l,m[ũl,m g l] [ p e l 1,m +p + e l+1,m +e l,m+1 +λ te ][ũl,m 1 g l], (3.20) i which we also used the assumptio (3.9) to derive the upper boud for the itegral term. If there are some (l,m) such that ũ l,m = g l, sice {ũ } 0 is a icreasig sequece from Propositio 3.1, we have ũ l,m = ũ l,m for all. Therefore, e l,m = 0 for these (l,m). O the other had, if ũ l,m > g l for some (l,m), we ca divide ũ l,m g l o both sides of (3.20) to get (1+p 0 )e l,m p e l 1,m +p + e l+1,m +e l,m+1 +λ te 1 p E m +p+ E m +Em+1 +λ te 1. Sice the right-had-side of (3.21) does ot deped o l, we ca write E m ηe m+1 (3.21) λ +(1 η) λ+r E 1, (3.22) 1 i which η = 1+(λ+r) t (0,1). Note that (3.22) is also satisfied for all m, because eve if ũl,m = g l for some (l,m), e l,m = 0 as we proved above. If follows from (3.22) that E m ηm m Em M +(1 η)(1+η λ + +ηm m 1 ) λ+r E 1. (3.23)

9 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 9 Sice the termial coditio of ũ, we have E M = 0. Now maximizig the right-had-side of (3.23) over m, we obtai that E (1 η M λ ) λ+r E 1. As a result, E ( 1 η M) ( ) λ E 0 0, as +. (3.24) λ+r Remark 3.2. As M + ( ) M 1 η M 1 = 1 1 e (r+λ)t, 1+(λ+r)T/M which agree with the coverget rate (2.10) i the cotiuous case. Propositio 3.4. u(x k,m t) ũ(x k.m t) 0, (3.25) as x, t, z 0. Proof. Usig the triagle iequality, let us write u(x k,m t) ũ(x k,m t) u(x k,m t) u (x k,m t) + u (x k,m t) ũ (x k,m t) + ũ (x k,m t) ũ(x k,m t) ( K 1 e (r+λ)(t m t)) ( ) λ + O ( ( t)+( x) 2 +( z) 2) λ+r + K ( 1 η M) ( ) λ, λ+r (3.26) for some positive costats K ad K. The first ad third terms o the right-had-side of the secod iequality are due to (2.10) ad (3.18). The secod term arises sice the order of error from discretizig a PDE usig implicit Euler scheme is O(( t) + ( x) 2 ), the iterpolatio ad discretizatio error from umerical itegral are of order ( x) 2 ad ( z) 2 ad the total error made at each step propagates at most liearly i whe we sequetially discretize (3.4). Lettig t, x, z 0 i (3.26), we obtai that ( lim u(x k,m t) ũ(x k,m t) K + K ) ( ) ( λ 1 e (r+λ)t), t, x, z 0 λ+r i which we used (3.25). Sice is arbitrary the result follows. Remark 3.3. I Propositios , we have show the covergece of the algorithm for completely implicit Euler scheme (θ = 1). I order to have the time discretizatio error as O(( t) 2 ), we will choose Crak-Nicolso scheme with θ = 1/2 i the umerical experimets i the ext sectio. From umerical results i Table 4, we shall see that Crak-Nicolso Scheme is also stable ad the covergece is fast.

10 10 ERHAN BAYRAKTAR AND HAO XING 4. The Numerical Performace of the Proposed Numerical Algorithm I this sectio, we preset the umerical performace of the algorithm proposed i the previous sectio. First, we compare the prices we obtai to the prices obtaied i the literature. To demostrate our competitiveess we also list the time it takes to obtai the prices for certai accuracy. We will use either the PSOR or the Brea- Schwartz algorithm to solve the sparse liear system i (3.10); see Remark 4.1. All our computatios are performed with C++ o a Petium IV, 3.0 GHz machie. I Table 1, we take the jump distributio F to be the double expoetial distributio F(dz) = ( pη 1 e η1z 1 {z 0} +(1 p)η 2 e η2z 1 {z<0} ) dz. (4.1) We compare our performace with that of [17] ad [16]. [17] obtai a approximate America optio price formula, for by reducig the itegro-pde equatio V satisfied to a itegro-ode followig [5]. This approximatio is accurate for small ad large maturities. Also, they do ot provide error bouds, the magitude of which might deped o the parameters of the problem, therefore oe might ot be able to use this price approximatio without the guidace of aother umerical scheme. A more accurate umerical scheme usig a approximatio to the exercise boudary ad Laplace trasform was later developed by [16]. Our performace has the same order of magitude as theirs. Our method s advatage is that it works for a more geeral jump distributio ad we do ot have to assume a double expoetial distributio for jumps as [17] ad [16] do. I Table 2 we compute the prices of America ad Europea optios i a Merto jump diffusio model, i which the jump distributio F is specified to be the Gaussia distributio F(dz) = ( ) 1 (z µ) 2 exp 2π σ 2 σ 2 dz. (4.2) We list the accuracy ad time characteristics of the proposed umerical algorithm algorithm. We compare our prices to the oes obtaied by [11, 12]. [11] used a pealty method to approximate the America optio price, while we aalyze the variatioal iequalities directly (see (3.5) ad (3.10)). Moreover, our approximatig sequece is mootoe (see Propositio 3.1). I Table 3, We also list the approximated prices of Barrier optios. We compare the prices we obtai with [19] where a Mote Carlo method is used. We do ot list the time it takes for the alterative algorithms i Tables 2 ad 3 either because they are ot listed i the origial papers or they take ureasoably log time. I Table 4, we list the umerical covergece of the proposed algorithm algorithm with respect to grid sizes. We choose Crak-Nicolso scheme with θ = 1/2 i (3.10) ad solve the sparse liear system by either the Bera- Schwartz algorithm or the PSOR. Remark 4.1. Here we will aalyze the complexity of our algorithm. Let us fix x/ t as a costat ad choose the umber of grid poit i x to be N. For each time step, usig the FFT to calculate the itegral term i (3.10) costs O(N log N) computatios. O the other had, the Brea-Schwartz algorithm, which uses the LU decompositio to solve the sparse liear system i (3.10) (see [15] pp. 283), eeds 2N computatios for each time step. However, PSOR eeds C N computatios for each time step to solve (3.10) at each time step. Here, C is the umber of iteratios PSOR requires to coverge to a fixed small error tolerace ǫ. We will see i the followig that PSOR is umerically more expesive tha the Brea-Schwartz algorithm.

11 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 11 For PSOR, the umber of iteratios C icreases with respect to N. To see this, we start from the tri-diagoal matrix o the left-had-side of (3.10) 1+θp 0 θp +. θp 1+θp.. 0 A = θp+ θp 1+θp 0 For the SOR (without projectio), the optimal relaxatio parameter ω is give by (see [22]) ω = ρ 2 J, where ρ J is the spectral radius of the Jacobi iteratio matrix J = D 1 (A D) with D as the diagoal matrix of A. Sice ρ J J = θ(p + +p )/(1+θp 0 ), we have 2 ω ω 0 =. (4.3) 1+ 1 J 2 + We will use ω 0 as the optimal relaxatio parameter i our umerical experimets. O the other had, sice the largest eigevalue λ max of the SOR iteratio matrix is bouded above by ω 1, usig (3.6) ad (4.3) we obtai that C = mi{c 0 (λ max ) c ǫ} = O( N). (4.4) Sice O(N 3/2 ) domiates O(N logn), the complexity of the Iterated PSOR algorithm at each time step will be O(N 3/2 ). Therefore, with O(N) time steps, the complexity for Iterated PSOR algorithm is O(N 5/2 ). O the other had, for the Iterated Brea-Schwartz algorithm, sice O(N log N) domiates O(N), the complexity at each time step will be O(N logn). Therefore, the complexity of the Iterated Brea-Schwartz algorithm is O(N 2 logn) sice we have O(N) time steps. Please refer to Tables 1, 2, 3 ad 4 for umerical performace of both algorithms. Next, we illustrate the behavior of the sequece of fuctios {v (S,t)} 0 ad its limit V i Figures 1, 2 ad 3. All the figures are obtaied for a America put optio i the case of the double expoetial jump with K = 100, S 0 = 100, T = 0.25, r = 0.05, σ = 0.2, λ = 3, p = 0.6, η 1 = 25 ad η 2 = 25 (the same parameters are used i the 8th row of Table 1) at a sigle ru. Remark 4.2. (i) I Figure 1, we show, how V(S,0) depeds o the time to maturity, ad that it fits smoothly to the put-pay-off fuctio at s(0) (the exercise boudary). The y-axis is the differece betwee the optio price ad the pay-off fuctio. As the time to maturity icreases, the optio price V(S,0) icreases while the exercise boudary s(0) decreases. Eve though the stock price process has jumps, the optio price smoothly fits the pay-off fuctio at s(0), as i the classical Black-Scholes case without the jumps. (ii) I Figure 2, we illustrate the covergece of the exercise boudaries t s (t), 1. We ca see from the figure that all s (t) are covex fuctios. Also, the sequece {s } 1 is a mootoe decreasig sequece, which implies that the cotiuatio regio is gettig larger, ad that the covergece of the free boudary sequece is fast. Moreover, we otice that, whe the parameters are chose such that (2.12) is satisfied, the free boudaries are discotiuous at the maturity time. I additio, we have s(t ) = s (T ) = S < K, where S

12 12 ERHAN BAYRAKTAR AND HAO XING is the uique solutio of (2.14). Furthermore, if F is the double expoetial distributio as i (4.1), the itegral equatio (2.14) ca be solved aalytically to obtai ( ) 1/η1 S (η1 1)r = K. (4.5) λp With the parameters we choose, we get from (4.5) that S = It is close to our umerical result as oe ca see from Figure 2. (iii) I Figure 3, we illustrate the covergece of the sequece of prices {v (S,0)} 0. Observe that this is a mootoically icreasig sequece ad it coverges to its limit V(S, 0) very fast. Ackowledgmet We are grateful to the two aoymous referees for their detailed commets that helped us improve our paper. Refereces [1] F. Aitsahlia ad A. Ruemo. A caoical optimal stoppig problem for America optios uder a doubleexpoetial jump-diffusio model. Joural of Risk, 10:85 100, [2] A. Almedral ad C. Oosterlee. O America optios uder the variace gamma process. Applied Mathematical Fiace, 14(2): , [3] K. I. Ami. Jump diffusio optio valuatio i discrete time. Joural of Fiace, 48: , [4] L. Aderse ad J. Adrease. Jump-diffusio processes: Volatility smile fittig ad umerical methods for optio pricig. Review of Derivatives Research, 4(3): , [5] G. Baroe-Adesi ad R. E. Whaley. Efficiet aalytic approximatio of America optio values. Joural of Fiace, 42: , [6] E. Bayraktar. A proof of the smoothess of the fiite time horizo America put optio for jump diffusios. To appear i the SIAM Joural o Cotrol ad Optimizatio, Available at [7] E. Bayraktar ad H. Xig. Aalysis of the optimal exercise boudary of America optios for jump diffusios. Techical report, Uiversity of Michiga, Available at [8] M. J. Brea ad E. S. Schwartz. The valuatio of America put optios. Joural of Fiace, 32(2): , [9] R. Cot ad E. Voltchkova. A fiite differece scheme for optio pricig i jump diffusio ad expoetial Lévy models. SIAM Joural o Numerical Aalysis, 43(4): , [10] Rama Cot ad Peter Takov. Fiacial modellig with jump processes. Chapma & Hall/CRC Fiacial Mathematics Series. Chapma & Hall/CRC, Boca Rato, FL, [11] Y. d Hallui, P. A. Forsyth, ad G. Labah. A pealty method for America optios with jump diffusio processes. Numerische Mathematik, 97(2): , [12] Y. d Hallui, P. A. Forsyth, ad K. R. Vetzal. Robust umerical methods for cotiget claims uder jump diffusio processes. IMA Joural of Numerical Aalysis, 25(1):87 112, [13] A. Hirsa ad D. Mada. Pricig America optios uder variace gamma. Joural of Computatioal Fiace, 7(2):63 80, 2004.

13 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 13 [14] Keeth R. Jackso, Sebastia Jaimugal, ad Vladimir Surkov. Fourier space time steppig for optio pricig with Lévy models. To appear i the Joural of Computatioal Fiace, [15] P. Jaillet, D. Lamberto, ad B. Lapeyre. Variatioal iequalities ad the pricig of America optios. Acta Applicadae Mathematicae, 21(3): , [16] S. G. Kou, G. Petrella, ad H. Wag. Pricig path-depedet optios with jump risk via laplace trasforms. Kyoto Ecoomic Review, 74:1 23, [17] S. G. Kou ad H. Wag. Optio pricig uder a double expoetial jump diffusio model. Maagemet Sciece, 50: , [18] R. C. Merto. Optio pricig whe the uderlyig stock returs are discotiuous. Joural of Fiacial Ecoomics, 3: , [19] S. A. K. Metwally ad A. F. Atiya. Fast mote carlo valuatio of barrier optios for jump diffusio processes. Proceesigs of the Computatioal Itelligece for Fiacial Egieerig, pages , [20] Paul Wilmott, Sam Howiso, ad Jeff Dewye. The mathematics of fiacial derivatives. Cambridge Uiversity Press, Cambridge, A studet itroductio. [21] C. Yag, L. Jiag, ad B. Bia. Free boudary ad America optios i a jump-diffusio model. Europea Joural of Applied Mathematics, 17(1):95 127, [22] D. M. Youg. Iterative solutio of large liear system. Academic Press, New York, [23] X. L. Zhag. Valuatio of America optios i a jump-diffusio model. I Numerical methods i fiace, Publ. Newto Ist., pages Cambridge Uiv. Press, Cambridge, 1997.

14 14 ERHAN BAYRAKTAR AND HAO XING Table 1. Compariso betwee the proposed iterated jump algorithm with the method i [17] ad [16], where the parameters were chose as r = 0.05, S(0) = 100 ad p = 0.6. Ami s price is calculated i [17] usig the ehaced biomial tree method as i [3]. The accuracy of Ami s price is up to about a pey. The KPW 5EXP price from [16] is calculated o a Petium IV, 1.8 GHz, while the iterated price is calculated o Petium IV, 3.0GHz, both usig C++ implemetatio. Ru times are i secods. For umerical algorithm we propose, the umber of grid poits i x is chose as 2 6 ad t = x. The optio prices from both Iterated Brea-Schwartz ad Iterated PSOR are the same. Below B-S stads for the Brea-Schwartz. America Put Double Expoetial Jump Diffusio Model Parameter Values Ami s KW KPW 5EXP Proposed Algorithm K T σ λ η 1 η 2 Price Value Error Value Error Time Value Error B-S Time PSOR Time

15 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 15 Table 2. Optio price i Merto jump-diffusio model K=100, T=0.25, r=0.05, σ = 0.15, λ = 0.1. Stock price has logormal jump distributio with µ = 0.9 ad σ = For the iterated jump schemes, the umber of grid poits i x is chose as 2 7 ad t = x. Below B-S stads for the Brea-Schwartz. Optio Type a S(0) dflv b Proposed Algorithm Value Error LU(B-S) Time PSOR Time America Put c Europea Put Europea Call a The optio prices (for the same kid of optio) for differet S(0) are obtaied from a sigle ru. b The dflv price comes from [11, 12]. c the optio price is usig the iterated Brea-Schwartz scheme. Table 3. Europea dow-ad-out barrier call optio with Merto jump-diffusio model K=110, S(0)=100, T=1, r=0.05, σ = 0.25, λ = 2, rebate R=1, the Stock price has logormal jump distributio with µ = 0 ad σ = 0.1. For the algorithm we propose the umber of grid poits i x is chose as 2 6 ad t = x. Below we use the acroyms LU or SOR to tell wheher we use the LU factorizatio or the SOR to solve for the sparse liear systems at each time step. a The MA price comes from [19] Barrier H MA Price a Proposed Algorithm Value Error LU Time SOR Time

16 16 ERHAN BAYRAKTAR AND HAO XING Table 4. Covergece of the umerical algorithm with respect to grid sizes K=100, T=0.25, r=0.05, σ = 0.15, λ = 0.1, stock price has logormal jump distributio with µ = 0.9 ad σ = 0.45 (the same parameters that are used i [11]). The differetial equatio is discretized by the Crak-Nicolso scheme as (3.10) with θ = 1/2. The logarithmic variable x = logs is equally spaced discretized o a iterval [x mi,x max] with x = t. The umerical itegral is trucated o the smallest iterval [z mi,z max], such that [x + µ 4 σ,x + ũ + 4 σ] will be iside [z mi,z max] for ay x [x mi,x max]. The step legth for the umerical itegral is chose the same as the step legth i x, i.e. z = x. The umber of grid poits for to implemet the FFT is chose as a itegral power of 2. The error tolerace for PSOR method is 10 8 ad for the global iteratio is Ru times are i secods. Each row i the Differece colum of the followig table is v PSOR(L,M) v PSOR(L/2,M/2). B-S stads for the Brea-Schwartz algorithm. The umber of global iteratio is 3 for all the followig umerical experimets. S(0) No. of grid No. of time B-S Value B-S PSOR Value Differece PSOR Max. No. of poits i x ( L ) steps ( M ) v B S Time v PSOR Time PSOR iteratios a a a Usig (4.4), the umber of SOR iteratios ca be calculated. The calculatio gives 11,16,22 ad 31. Comparig with the last colum of above table, the maximum umbers of PSOR iteratio are slightly larger tha these theoretical predicted SOR iteratio times. Moreover, whe L = 512 the the ratio betwee the maximum umber of PSOR iteratio ad L is This cofirms the aalysis i Remark 4.1 that the maximal PSOR iteratio time grows as the order of L.

17 ON THE PRICING OF AMERICAN OPTIONS FOR JUMP DIFFUSIONS 17 The parameters for the followig three figures are K = 100, S 0 = 100, T = 0.25, r = 0.05, σ = 0.2, λ = 3, the stock price has double expoetial jump with p = 0.6, η 1 = 25 ad η 2 = 25 (the same parameters used i the 8th row of Table 1). Figure 1. The optio price fuctio S V(S,0) smoothly fits the pay-off fuctio (K S) + at s(0). V(S,0) icreases ad s(0) (V(S,0) (K S) + = 0 at s(0)) decreases as time to maturity T icreases. 4 T: time to maturity T = T = T = 0.15 Optio Price ( K S ) T = 0.10 T = T = Stock Price ( S )

18 18 ERHAN BAYRAKTAR AND HAO XING Figure 2. Iteratio of the Exercise Boudary: s (t) s(t), t [0,T). Both s (t) ad s(t) will coverge to S < K as t T. 100 S * approx = Cotiuatio Regio 95 Stock Price ( S ) s 1 ( t ) 90 s 2 ( t ) s i ( t ), i = 4, 5, 6 s 3 ( t ) t ( i years ) Stoppig Regio 10 Figure 3. Iteratio of the price fuctios: v (S,0) V(S,0), S v 3 ( S, 0 ) 7 Optio Price v i ( S, 0 ), i = 4, 5, 6 3 ( K S ) + 2 v 1 ( S, 0 ) 1 v 2 ( S, 0 ) Stock Price ( S )