arxiv: v2 [q-fin.pr] 15 Jan 2017

Transcription

1 Valuatio of Variable Auities with Guarateed Miimum Withdrawal Beefit uder Stochastic Iterest Rate Pavel V. Shevcheko 1, ad Xiaoli Luo 2 Draft, 15 Jauary 217 1st versio 16 February 216) 1 Applied Fiace ad Actuarial Studies, Macquarie Uiversity, Australia; pavel.shevcheko@mq.edu.au 2 CSIRO Australia; Xiaoli.Luo@csiro.au arxiv: v2 [q-fi.pr] 15 Ja 217 Correspodig author Abstract This paper develops a efficiet direct itegratio method for pricig of the variable auity VA) with guaratees i the case of stochastic iterest rate. I particular, we focus o pricig VA with Guarateed Miimum Withdrawal Beefit GMWB) that promises to retur the etire iitial ivestmet through withdrawals ad the remaiig accout balace at maturity. Uder the optimal dyamic) withdrawal strategy of a policyholder, GMWB pricig becomes a optimal stochastic cotrol problem that ca be solved usig backward recursio Bellma equatio. Optimal decisio becomes a fuctio of ot oly the uderlyig asset but also iterest rate. Presetly our method is applied to the Vasicek iterest rate model, but it is applicable to ay model whe trasitio desity of the uderlyig asset ad iterest rate is kow i closed-form or ca be evaluated efficietly. Usig bod price as a uméraire the required expectatios i the backward recursio are reduced to two-dimesioal itegrals calculated through a high order Gauss-Hermite quadrature applied o a two-dimesioal cubic splie iterpolatio. The quadrature is applied after a rotatioal trasformatio to the variables correspodig to the pricipal axes of the bivariate trasitio desity, which empirically was observed to be more accurate tha the use of Cholesky trasformatio. Numerical compariso demostrates that the ew algorithm is sigificatly faster tha the partial differetial equatio or Mote Carlo methods. For pricig of GMWB with dyamic withdrawal strategy, we foud that for positive correlatio betwee the uderlyig asset ad iterest rate, the GMWB price uder the stochastic iterest rate is sigificatly higher compared to the case of determiistic iterest rate, while for egative correlatio the differece is less but still sigificat. I the case of GMWB with predefied static) withdrawal strategy, for egative correlatio, the differece i prices betwee stochastic ad determiistic iterest rate cases is ot material while for positive correlatio the differece is still sigificat. The algorithm ca be easily adapted to solve similar stochastic cotrol problems with two stochastic variables possibly affected by cotrol. Applicatio to umerical pricig of Asia, barrier ad other fiacial derivatives with a sigle risky asset uder stochastic iterest rate is also straightforward. Keywords: Variable auity, livig ad death beefits, stochastic iterest rate, optimal stochastic cotrol, Guarateed Miimum Withdrawal Beefit, Gauss-Hermite quadrature. 1

2 1 Itroductio The world populatio is gettig older fast with life expectacy raisig to above 9 years i some coutries. Logevity risk the risk of outlivig oe s savigs) became critical for retirees. Variable auity VA) with livig ad death beefit guaratees is oe the products that ca help to maage this risk. It takes advatage of market growth ad at the same time provides protectio of the savigs. VA guaratees are typically classified as guarateed miimum withdrawal beefit GMWB), guarateed miimum accumulatio beefit GMAB), guarateed miimum icome beefit GMIB), ad guarateed miimum death beefit GMDB). A good overview of VA products ad the developmet of their market ca be foud i Bauer et al. 28), Ledlie et al. 28) ad Kalberer ad Ravidra 29). Isurers started to sell these types of products from the 199s i Uited States. Later, these products became populari Europe, UK adjapa. The market of VAs is very large, for example, sales of these cotracts i Uited States betwee 211 ad 213 averaged about $16 billio per year accordig to the LIMRA Life Isurace ad Market Research Associatio) fact sheets. For clarity ad simplicity of presetatio, i this paper we cosider a VA cotract with a very basic GMWB guaratee that promises to retur the etire iitial ivestmet through cash withdrawals durig the policy life plus the remaiig accout balace at maturity, regardless of the portfolio performace. Thus eve whe the accout of the policyholder falls to zero before maturity, GMWB feature will cotiue to provide the guarateed cashflows. GMWB allows the policyholder to withdraw fuds below or at the cotractual rate without pealty ad above the cotractual rate with some pealty. If the policyholder behaves passively ad makes withdrawals at the cotractual rate defied at the begiig of the cotract, the the behavior of the policyholder is called static. I this case the paths of the wealth accout ca be simulated ad a stadard Mote Carlo MC) simulatio method ca be used for GMWB pricig. O the other had if the policyholder optimally decides the amout to withdraw at each withdrawal date, the the behavior of the policyholder is called dyamic. Uder the optimal withdrawal strategy, the pricig of variable auities with GMWB becomes a optimal stochastic cotrol problem. This problem caot be solved by a stadard simulatio-based method such as the well kow Least-Squares MC method itroduced i Logstaff ad Schwartz 21). This is because the paths of the uderlyig wealth process are altered by the optimal cash withdrawals that should be foud from the backward i time solutio ad the uderlyig wealth process caot be simulated forward i time. However, it should be possible to apply cotrol radomizatio methods extedig Least-Square MC to hadle optimal stochastic cotrol problems with cotrolled Markov processes recetly developed i Kharroubi et al. 214); though the accuracy ad robustess of this method for GMWB pricig has ot bee studied yet. It is importat to ote that the fair fee for the VA guaratee obtaied uder the assumptio that the policyholders behave optimally to maximise the value of the guaratee is a importat bechmark because it is a worst case sceario for the cotract writer. That is, uder the o-arbitrage 2

3 assumptio, if the guaratee is perfectly hedged the the issuer will receive a guarateed profit if the policyholder deviates from the optimal strategy. Pricig uder ay other strategy will lead to smaller fair fee. Of course, the strategy optimal i this sese may ot be optimal to the policyholder uder his circumstaces ad prefereces. O the other had, secodary markets for equity liked isurace products are growig ad fiacial third parties ca potetially geerate guarateed profit through hedgig strategies from VA guaratees which are ot priced accordig to the worst case assumptio about the optimal strategies. There are a umber of studies cosiderig these aspects ad we refer the reader to Shevcheko ad Luo 216) for discussio of this topic ad refereces therei. Pricig of VA with a GMWB feature assumig costat iterest rate has bee cosidered i may papers over the last decade. For example, Milevsky ad Salisbury 26) developed a variety of methods for pricig GMWB products. I their static withdrawal approach the GMWB product is decomposed ito a Quato Asia put optio plus a geeric term-certai auity. They also cosidered pricig whe the policyholder ca termiate surreder) the cotract at the optimal time, which leads to a optimal stoppig problem aki to pricig a America put optio. Bauer et al. 28) presets valuatio of variable auities with multiple guaratees via a multidimesioal discretizatio approach i which the Black-Scholes partial differetial equatio PDE) is trasformed to a oe-dimesioal heat equatio ad a quasi-aalytic solutio is obtaied through a simple piecewise summatio with a liear iterpolatio o a mesh. Dai et al. 28) developed a efficiet fiite differece algorithm usig the pealty approximatio to solve the sigular stochastic cotrol problem for a cotiuous time withdrawal model uder the optimal withdrawal strategy ad also fiite differece algorithm for discrete time withdrawal. Their results show that the GMWB values from the discrete time model coverge fast to those of the cotiuous time model. Huag ad Forsyth 212) did a rigorous covergece study of this pealty method for GMWB, ad Huag ad Kwok 214) deduce various asymptotes for the free boudaries that separate differet withdrawal regios i the domai of the GMWB pricig model. Che ad Forsyth 28) preset a impulse stochastic cotrol formulatio for pricig variable auities with GMWB uder the optimal policyholder behavior, ad develop a umerical scheme for solvig the Hamilto-Jacobi-Bellma variatioal iequality for the cotiuous withdrawal model as well as for pricig the discrete withdrawal cotracts. More recetly, Azimzadeh ad Forsyth 214) prove the existece of a optimal bag-bag cotrol for a Guarateed Lifelog Withdrawal Beefits GLWB) cotract. I particular, they fid that the holder of a GLWB ca maximize the cotract writer s losses by oly performig o-withdrawal, withdrawal at exactly the cotract rate or full surreder. This dramatically reduces the optimal strategy space. However, they also demostrate that the related GMWB cotract is ot covexity preservig, ad hece does ot satisfy the bag-bag priciple other tha i certai degeerate cases. GMWB pricig uder bag-bag strategy was studied i Luo ad Shevcheko 215c), ad Huag ad Kwok 215) have developed a regressio-based MC method for pricig GLWB. For GMWB uder the optimal withdrawal strategy, the umerical evaluatios have bee developed by Dai et al. 28) ad 3

4 Che ad Forsyth 28) usig fiite differece PDE methods ad by Luo ad Shevcheko 215a) usig direct itegratio method. Pricig of VAs with both GMWB ad death beefit both uder static ad dyamic regimes) has bee developed i Luo ad Shevcheko 215b). Some withdrawals from the VA type cotracts ca also attract coutry specific govermet additioal tax ad pealty. Recetly, Moeig ad Bauer 215) demostrated that icludig taxes sigificatly affects the value of withdrawal guaratees i variable auities producig results i lie with empirical market prices. These matters are ot cosidered i our paper but ca be hadled by the umerical methodology developed here. I the literature o pricig GWMB, iterest rate is typically assumed to be costat. Few papers cosidered the case of stochastic iterest rate. I particular, Peg et al. 212) cosidered pricig GMWB uder the Vasicek stochastic iterest rate i the case of static withdrawal strategy; they derived the lower ad upper bouds for the price because closed-form solutio is ot available due to withdrawals from the uderlyig wealth accout durig its stochastic evolutio. Baciello et al. 211) cosidered stochastic iterest rate ad stochastic volatility models uder the Cox-Igersoll-Ross CIR) models. They developed pricig i the case of static policyholder behavior via the ordiary MC method ad mixed valuatio where the policyholder is semiactive ad ca decide to surreder the cotract at ay time before the maturity) is performed by the Least-Squares MC. Forsyth ad Vetzal 214) cosidered modellig stochasticity i the iterest rate ad volatility via the Markov regime switchig models ad developed pricig uder the static ad dyamic withdrawal strategies. Uder this approach, the iterest rate ad volatility are assumed to have the fiite umber of possible values ad their evolutio i time is drive by the fiite state Markov chai variable represetig possible regimes of the ecoomy. I this paper, we develop direct itegratio method for pricig of VAs with guaratees uder the dyamic ad static withdrawal strategies whe the iterest rate follows the Vasicek stochastic iterest rate model. I the case of geeral stochastic processes for the uderlyig asset ad iterest rate, umerical pricig ca be accomplished by PDE methods that become slow ad difficult to implemet i the case of two ad more uderlyig stochastic variables. Our method is developed for the case whe the bivariate trasitio desity of the uderlyig asset ad iterest rate are kow i closed-form or ca be evaluated efficietly. That is, it should be possible to apply this method to the case of, for example, CIR stochastic iterest rate model. Usig chage of uméraire techique with bod price as a uméraire, the required expectatios i the backward recursio of the stochastic cotrol solutio are reduced to the two-dimesioal itegrals calculated through a high order Gauss-Hermite quadrature applied o a two-dimesioal cubic splie iterpolatio. The quadrature is applied after rotatioal trasformatio to the variables correspodig to the pricipal axes of the bivariate trasitio desity which appeared to be more efficiet tha the use of the stadard Cholesky trasformatio to the idepedet variables. For coveiece, hereafter we refer this ew algorithm as GHQC Gauss- Hermite quadrature o cubic splie). This allows us to get very fast ad accurate results for prices 4

5 of a typical GMWB cotract o the stadard desktop computer. Previously, i a similar spirit, we developed algorithm for the case of oe uderlyig stochastic risky asset ad o-stochastic iterest rate for pricig exotic optios i Luo ad Shevcheko 214) ad optimal stochastic cotrol problems for pricig GMWB i Luo ad Shevcheko 215a). For clarity of presetatio, we focus o pricig of a VA with a very basic GMWB structure. However, the developed methodology ca be easily applied to pricig other VA guaratees, see Shevcheko ad Luo 216) for geeral formulatio of these cotracts as the optimal stochastic cotrol problem. Fially we would like to metio that the preseted algorithm ca be easily adapted to solve similar stochastic cotrol problems with two state variables possibly affected by cotrol. Also, applicatios to pricig Asia, barrier ad other fiacial derivatives with a sigle uderlyig asset uder stochastic iterest rate are straightforward. I the ext sectio we preset the uderlyig stochastic model ad describe the GMWB cotract. Solutio as a optimal stochastic cotrol is preseted i Sectio 3. Sectio 4 gives a short descriptio of the well kow PDE approaches that ca be used for pricig. Sectio 5 presets our direct itegratio GHQC algorithm for pricig of GMWB cotracts uder both static ad dyamic policyholder behaviors. I Sectio 6, umerical results for the fair prices ad fair fees uder a series GMWB cotract coditios are preseted, i compariso with the results from the fiite differece method solvig correspodig two-dimesioal PDEs. The compariso demostrates that the ew algorithm produces results very close to those of the fiite differece PDE method, but at the same time it is sigificatly faster. Also, the results demostrate that stochastic iterest rate has sigificat impact o price aalysed i Sectio 6. Cocludig remarks are give i Sectio 7. Useful closed-form formulas for the required trasitio desities, bod ad vailla prices are derived i Appedix A. 2 Model Followig the existig literature, we assume o-arbitrage market with respect to the fiacial risk ad thus the price of the VA with GMWB ca be expressed as a expectatio with respect to the risk-eutral probability measure for the uderlyig risky asset. Also, there is o mortality risk i the evet of policyholder death, the cotract is maitaied by the beeficiary. Death beefit feature commoly offered to the policyholders i additio to GMWB ca be easily icluded ito pricig methodology as described i Luo ad Shevcheko 215b). Let Ω,F,Q) be a probability space with sample space Ω, filtratio F = {F t : t } sequece of σ-algebras F t icreasig with time t o Ω) ad risk-eutral probability measure Q such that all discouted asset price processes are Q-martigales, i.e. paymet streams ca be valuated as expected discouted values. Existece of measure Q implies that the fiacial market is arbitrage-free ad uiqueess of such measure implies that the market is complete. This meas that the cost of a portfolio replicatig VA cotract with guaratee is give by its expected discouted value uder Q. 5

6 This is a typical set up for pricig of fiacial derivatives, for a good textbook i this area we refer the reader to e.g. Björk 29). Cosider the joit dyamics for the referece portfolio of assets St), e.g. a mutual fud, uderlyig the cotract ad the stochastic iterest rate rt), uder the risk-eutral probability measure Q, govered by dst) St) = rt)dt+σ S ρdb 1 t)+ ) 1ρ 2 db 2 t), drt) = κθ rt))dt+σ r db 1 t). Here, B 1 t) ad B 2 t) are idepedet stadard Wieer processes, ρ is the correlatio coefficiet betwee St) ad rt) processes, ad σ S is the asset volatility parameter. The process for the iterest rate rt) is the well kow Vasicek model with costat parameters κ, θ ad σ r. For simplicity of otatio we assume that model parameters are costat i time though the results ca be geeralized to the case of time depedet parameters. We cosider time discretizatio = t < t 1 < < t N = T correspodig to the cotract withdrawal dates, where t = is today ad T is the cotract maturity. For this stochastic iterest rate model, the price of a zero coupo bod Pt,T) at time t with maturity T, ca be foud i closed-form B t,t = 1 κ Pt,T) := E Q t 1) [e T t ru)du ] = e A t,trt)b t,t, 2) 1e κtt)) ), A t,t = θ σ2 r 2κ 2 B t,t +tt) σ2 r 4κ B2 t,t, where E Q t [ ] deotes expectatio with respect to the probability measure Q coditioal o the iformatio available at time t. Correspodig stochastic dyamics is easily obtaied from 2) usig Itô s calculus to be dpt,t) Pt,T) = rt)dtσ rb t,t db 1 t). 3) Solutio for the process1), which is a bivariate Normal distributio forl St), rt)) gives), r)), ad the bod price formula are derived i Appedix A. Cosider the followig VA cotract with a basic GMWB ofte used i research studies which is coveiet for bechmarkig. The actual products may have extra features but these ca be easily icorporated i the model ad umerical algorithm developed i this paper. The premium paid by the policyholder upfrot at t is ivested ito the referece portfolio/risky asset St). The value of this portfolio hereafter referred to as wealth accout) at time t is deotedas Wt), sothat theupfrotpremiumpaidbythepolicyholderisw). GMWBguaratees the retur of the premium via the withdrawals γ allowed at times t, = 1,2,...,N. Let N w deote the umber of withdrawals per aum. The total of withdrawals caot exceed the guaratee W) ad withdrawals ca be differet from the cotractual guarateed) 6

7 withdrawal G = W)t t 1 )/T, with pealties imposed if γ > G. Deote the aual cotractual rate as g := 1/T. The the wealth accout Wt) evolves as Wt ) = Wt+ 1 ) St 1 ) St )e α, Wt + ) = max Wt )γ, ), = 1,2,...,N, 4) where = t t 1 ad α is the aual fee cotiuously charged by the cotract issuer. If the accout balace becomes zero or egative, the it will stay zero till maturity. The process for Wt) withi t 1,t ) is the same as the process for the uderlyig asset St) i 1) except that the drift term rt) is replaced by rt)α. Deote the value of the cotract guaratee at time t as At), hereafter referred to as guaratee accout, with A) = W). Hereafter, deote the time immediately before t i.e. before withdrawal) as t, ad immediately after t i.e. after withdrawal) as t + ad let all fuctios discotiuous at t be right-cotiuous with fiite left limit. The guaratee balace evolves as At + ) = At )γ = At + 1 )γ, = 1,2,...,N 5) with AT + ) =, i.e. W) = A) γ γ N ad At + 1 ) N k= γ k. The accout balace At) remais uchaged withi the iterval t 1, t ), = 1,2,...,N. The cashflow received by the policyholder at the withdrawal time t is give by C γ ) = { γ, if γ G, G +1β)γ G ), if γ > G, 6) where G is the cotractual withdrawal ad β [,1] is the pealty coefficiet applied to the portio of withdrawal above G. Let Q t W,r,A) be a price of the VA cotract with GMWB at time t, whe Wt) = W, rt) = r, At) = A. At maturity, the policyholder takes the maximum betwee the remaiig guaratee accout et of pealty charge ad the remaiig balace of the wealth accout, i.e. the fial payoff is Q t W,r,A) = maxw,c N A)). 7) N Durig the cotract, the policyholder receives cashflows C γ ), = 1,2,...,N 1 ad the fial payoff at maturity. Deote the Markov state vector at time t as V t = Wt),rt),At)) ad V = V t ) t T. Give the withdrawal strategy γ = γ 1,...,γ N1 ), the preset value of the total cotract payoff is H V,γ) = e T rτ)dτ max WT ),C N AT )) ) N1 + e rτ)dτ C γ ). 8) 7 =1

8 Uder the above assumptios/coditios, the fair o-arbitrage value of the cotract for the pre-defied static) withdrawal strategy γ ca be calculated as Q V ) = E Q t [H V,γ)]. 9) Uder the optimal dyamic) withdrawal strategy, where the decisio o withdrawal amout γ is based upo the iformatio F t available at time t, the fair cotract value is Q V ) = supe Q t [H V,γ)], 1) γ where γ 1,...,γ N1 are the cotrol variables withdrawals) chose to maximize the expected value of discouted cashflows ad supremum is take over all admissible strategies. Note that the withdrawal γ := γ V t ) at time t is a fuctio of the state variable V t, i.e. it ca be differet for differet realizatios of V t. Moreover the cotrol variable γ affects the trasitio law of the uderlyig wealth process from t to t +1. Ay strategy differet from optimal is sub-optimal ad leads to a smaller price. The today s value of the cotract Q V ) is a fuctio of a fee α charged by the issuer for GMWB guaratee. The fair fee value of α to becharged for providig GMWB feature correspodsto Q V ) = W). That is, oce a pricig of Q V ) for a give value of α is developed, the a umerical root search algorithm is required to fid the fair fee. It is importat to ote that the fair fee for the VA guaratee obtaied uder the assumptio that the policyholders behave optimally to maximise the value of the guaratee is a worst case sceario for the cotract writer. If the guaratee is perfectly hedged the the issuer will receive a guarateed profit if the policyholder deviates from the optimal strategy. Pricig uder ay other strategy will lead to smaller fair fee. Of course, the strategy optimal i this sese may ot be optimal to the policyholder uder his circumstaces ad prefereces but it is a importat bechmark. I practice, there will be a residual risk due to discrete i time hedgig ad icompleteesses of fiacial market that ca be hadled by addig extra loadig o the price uder the actuarial approach or adjustig risk premium uder the o-arbitrage fiacial mathematics approach so that the risk of hedgig error loss will ot exceed the required level. These adjustmets deped o the risk maagemet strategy for the product ad will ot be cosidered here; for discussio ad refereces, see e.g. Shevcheko ad Luo 216). Remark 2.1 Note that we started our modellig with assumptio of the stochastic model 1) uder the risk-eutral probability measure Q. For risk maagemet purposes, oe might be iterested to start with the process uder the real physical) probability measure P, dst) St) = µ t)dt+σ S ρdb1t)+ ) 1ρ 2 db2t), drt) = u r,t)dt+σ r db 1t), 8

9 with idepedet stadard Wieer processes B1 t) ad B 2 t), ad derive correspodig risk-eutral process 1) for the VA guaratee ad bod price valuatio. This ca be doe i the usual way by formig a portfolio Π t = U t W,r,A)+ S S+ P Pt,T), where U t W,r,A) := Q t W,r,A)W is the value of VA guaratee, S is the umber of uits of St) ad P is the umber of uits of bod Pt,T). The calculate the chage of portfolio dπ t usig Ito s lemma, ad set S = W S U t W,r,A) W ad P = U tw,r,a)/ r Pt,T)/ r to elimiate radom terms so that the portfolio ears risk free iterest rate dπ t = rπ t dt. This leads to a PDE 18) for Q t W,r,A) ad usig Feyma-Kac theorem oe ca establish that the process correspodig to this PDE is the risk-eutral process 1). For details, see e.g. Shevcheko ad Luo, 216, sectio 6.5) ad textbook Wilmott, 26, sectios 3.3 ad 33.6). It is importat to ote that this procedure will also itroduce the market price of iterest rate risk λr, t) such that the drift of the iterest rate risk-eutral process is u r,t)λr,t)σ r ; the uder the assumptio that λr,t) ad u r,t) are liear fuctios of r oe ca write the risk-eutral process for r as i 1). 3 Pricig GMWB as optimal stochastic cotrol Give that the discrete i time state vector V t = Wt ),rt ),At )), =,1,...,N is a Markov process, it is easy to recogize that the cotract valuatio uder the optimal withdrawal strategy 1) is the optimal stochastic cotrol problem for cotrolled Markov process that ca be solved recursively to fid the cotract value Q t ) at t, = N 1,..., via the well kow backward iductio Bellma equatio Q t Wt ),rt ),At )) = +E Q t + [ e +1 t rτ)dτ Q t +1 sup C γ ) γ At ) ]) Wt +1 ),rt +1 ),At +1 )) Wt + ),rt ),At + ) startig from the fial coditio Q t W,r,A) = maxw,c N A)). For a good textbook treatmet N of stochastic cotrol problem i fiace, see Bäuerle ad Rieder 211). Static pricig 9) uder the predefied strategy γ ca be also doe usig the above backward iductio with supremum removed. For each t, = 1,...,N 1, this backward recursio 11) ivolves calculatio of the expectatio Q t + W,r,A) = E [e Q +1 t + t rτ)dτ Q t Wt ),rt +1 ),At +1 W,r,A] )) 12) ad applicatio of the jump coditio across t Q t W,r,A) = max γ A [C γ )+Q t + maxw γ,),r,a γ )]. 13) Calculatig expectatio 12) is difficult as it would require three-dimesioal itegratio with respect to the joit distributio of three radom variables Wt +1 ), rt +1) ad Yt +1 ) = +1 t 9 11) ru)du

10 coditioal o Wt + ) ad rt ); ote that variable At) does ot chage withi t,t +1 ). Actually the required 3d distributio ca be foud i closed-form i the case of stochastic process 1) cosidered here, see Appedix A, which is useful for validatio of calculatios i the case of static withdrawals via direct simulatio of process 1). However if we chage uméraire from the moey market accout Mt) = e rτ)dτ to the bod Pt,t +1 ) with maturity t +1, i.e. chage probability measure with Rado-Nikodym derivative Z t = d Q = Mt ) Pt,t +1 ) dq Mt) Pt,t +1 ), t [t,t +1 ], 14) Ft the the expectatio 12) simplifies to the two-dimesioal itegratio [ e +1 t ru)du Q t Wt ),rt ), ) ] +1 E Q t + = Pt,t +1 )E Q t + [ Q t Wt ),rt ), ) ] +1, 15) where E Q t [ ] is expectatio uder the ew probability measure Q. The process for Z t is easily obtaied from the process 3) for the bod price Pt,t +1 ) as dz t = φt)z t db 1, φt) = σ r B t,t+1. The, usig Girsaov theorem the required trasformatio to the Wieer process is B 1 t) = φt)dt+ d B 1 t), ad the processes uder the ew measure Q for t t,t +1 ) are dst)/st) = rt)+σ S ρφt))dt+σ S ρd B 1 t)+ ) 1ρ 2 d B 2 t), ) drt) = κ θt)rt) dt+σ r d B 1 t); θt) = θ+ σ 16) r κ φt) with B 1 t) ad B 2 t) idepedet Wieer processes. Note that φt) is volatility of the bod Pt,t +1 ), see 3). Solutio for this process, which is a bivariate Normal distributio for l St), rt)) give S), r)), is derived i Appedix A. For a good textbook treatmet of chage of uméraire techique, see Björk, 29, chapter 26). It is importat to ote that for differet time steps, the chage of measure is based o bods of differet maturities. AssumigtheprobabilitydesityfuctioofWt +1 )adrt +1)atWt +1 ) = w adrt +1 ) = r coditioal o Wt + ) = w ad rt ) = r uder the ew probability measure Q is kow i closed-form p +1 w,r w,r), the required expectatio 15) ca be evaluated as Q t + w,r,a) = Pt,t +1 ) p +1 w,r w,r ) Q t w,r,a)dw dr. 17) +1 I the case of uderlyig stochastic process 1) the trasitio desity p +1 w,r w,r) is kow i closed-form ad we will use the Gauss-Hermite quadrature for evaluatio of the above itegratio over a ifiite domai. The required cotiuous fuctio Q t W,r,A) will be approximated by a 1

11 two-dimesioal cubic splie iterpolatio o a discretized grid i the W, r) space. Note, i geeral, a three-dimesioal iterpolatio i the W, r, A) space is required, but oe ca maage to avoid iterpolatio i A by smart umerical maipulatio settig the jump amouts i A spaced i such away thattheareducedafterjumpisalways oagridpoit. Below wediscussdetails ofthealgorithm of the umerical itegratio of 17) usig Gauss-Hermite quadrature o a cubic splie iterpolatio, followed by the applicatio of jump coditio 13). Note, for simple optios/cotracts where payoff depeds o the uderlyig asset oly ad is received at the cotract maturity, a chage of uméraire ca remove stochastic iterest rate dimesio from pricig, effectively reducig umerical problem to the determiistic iterest rate case. However, for pricig GMWB either static or dyamic cases, the additioal dimesio i the iterest rate r caot be avoided due to withdrawals jump coditios) durig the cotract life. 4 Numerical valuatio of GMWB via PDE I the case of cotiuous i time withdrawal, followig the procedure of derivig the Hamilto-Jacobi- Bellma HJB) equatios i stochastic cotrol problems, the value of the VA cotract with guaratee uder the optimal withdrawal is foud to be govered by a two-dimesioal PDE i the case of determiistic iterest rate; see Milevsky ad Salisbury26), Dai et al. 28) ad Che ad Forsyth 28), that will become three-dimesioal PDE i the case of stochastic iterest rate. For discrete withdrawals, the goverig PDE i the period betwee withdrawal dates is oe dimesio less tha the cotiuous case because the guaratee accout balace At) remais uchaged betwee withdrawals, similar to the Black-Scholes equatio, with jump coditios at each withdrawal date to lik the prices at the adjacet periods. I particular, the cotract value Q t W,r,A) at t t 1,t ) satisfies Q t t + σ2 S Q t 2 W2 2 W 2 +r α)w Q t W + σ2 r 2 Q t 2 r 2 +θ r) Q t r +ρσ Sσ r W 2 Q t W r rq t =, 18) that ca be solved umerically usig e.g. Crak-Nicholso fiite differece scheme for each A backward i time with jump coditio 13) applied at withdrawal dates t. A more efficiet ad very popular class of algorithms is the alteratig directio implicit ADI) method, amog which a stadout variatio is the so called hopscotch method, itroduced by Gourlay 197) by a reformulatio of a idea of Gordo 1965). It was show that hopscotch method was a ADI process with a ovel way of decomposig the problem ito simpler parts. The geeral idea is to solve alterative poits explicitly ad the employ a implicit scheme to solve for the remaiig poits explicitly. The origial hopscotch method caot be applied readily to equatios with mixed derivatives without itroducig a certai amout of implicitess. Gourlay ad McKee 1977) suggested two techiques for dealig with the mixed derivative ordered odd-eve hopscotch ad lie hopscotch. Numerical tests i Gourlay ad McKee 1977) showed that the lie hopscotch performed best for both costat ad variable coefficiet parabolic equatio cases, i compariso with the ordered odd-eve hopscotch ad 11

12 a locally oe dimesioal LOD) methods. I this work, for umerical validatio of our GHQC algorithm, we have implemeted the lie hopscotch method. I brief, assumig the parabolic equatio is discretized with the fiite differece grid poits x i,y j,t ), the lie hopscotch method first explicitly evaluates the solutio at those poits which have +j) eve, ad the solves implicitly for those poits with +j) odd. The alterative value of j for a give time step gives a tri-diagoal set of equatios, provided the fiite differece operators are chose i a certai maer. For details, see Gourlay ad McKee 1977). 5 GHQC direct itegratio method I this sectio we preset details of the algorithm for umerical itegratio 17) usig the Gauss- Hermite quadrature o a cubic splie iterpolatio, followed by the applicatio of jump coditio 13), referred to as GHQC. 5.1 Algorithm structure Our approach relies o computig expectatios 17) i a backward time-steppig betwee withdrawal dates through a high order Gauss-Hermite itegratio quadrature applied o a cubic splie iterpolatio. It is easier to implemet ad computatioally faster tha PDE method i the case of trasitio desity of uderlyig stochastic variables kow i closed-form. For a give guaratee accout variable A withi t,t +1 ), the price Q t + W,r,A) ca be umerically evaluated usig 17). For ow we leave details of computig 17) to the ext sectio ad assume it ca be doe with sufficiet accuracy ad efficiecy. Startig from a fial coditio at t = t N just immediately before the fial withdrawal), a backward time steppig usig 17) gives solutio at t = t + N1. Applyig jump coditio 13) to the solutio at t = t + N1 we obtai the solutio at t = t N1 from which further backward time steppig gives us solutio at t to fid Q W),r),W)). I order to apply the jump coditio at each withdrawal date, the solutio has to be foud for may differet levels of A. The umerical algorithm takes the followig key steps. Step 1. Geerate a auxiliary fiite grid = A 1 < A 2 < < A J = W) to track solutios for differet values of the guaratee accout A. Discretize the wealth accout W space as W,W 1,...,W M ad the iterest rate r space as r,r 1,...,r K. Step 2. At t = t N, iitialize Q t NW,r,A) with a give cotiuous payoff fuctio at maturity 7) required by the followig step of itegratio. Step 3. For t = t + N1, evaluate itegratio 17) for each ode poit W m,r k,a j ) ad usig a oedimesioal cubic splie iterpolatio i W to obtai the cotiuous fuctio Q t + W,r k,a j ) N1 required by the followig step of applyig the jump coditio. 12

13 Step 4. Apply the jump coditio 13) for all possible withdrawals γ N1 ad fid the withdrawal maximizig Q t W m,r k,a j ) for all grid poits j = 1,...,J, k =,...,K ad m =,...,M. N1 Use two-dimesioal cubic splie iterpolatio to obtai cotiuous fuctio Q t W,r,A j ) N1 required by the ext step of itegratio. Step 5. Repeat Step 3 ad Step 4 for t = t N2,t N3,...,t 1. Step 6. Evaluate itegratio 17) for the backward time step from t 1 to t for the sigle poit W),r),A)) to obtai solutio Q W),r),A)) for the cotract price at t = t. 5.2 Numerical evaluatio of the expectatio Similar to a fiite differece scheme, we discretize the wealth space domai [W mi,w max ] as W mi = W < W 1 < < W M = W max, wherew mi adw max arethelower adupperboudaryrespectively. Similarly, the iterest rate space is discretized as r mi = r < r 1 <... < r K = r max, where r mi ad r max are the bouds for the iterest rate. For pricig GMWB, due to the fiite reductio of W at each withdrawal date, we have to cosider the possibility of zero W, thus the lower boud W mi =. The upper boud is set sufficietly far from the iitial value at time zero W). I geeral for both W ad r dimesios, the proper choice of the lower ad upper bouds is guided by the joit distributio of lst) ad rt), derived i Appedix A, to esure that the probability for the radom process to go beyod the bouds is immaterial. The idea is to fid the cotract values at all grid poits at each time step t + 1,t ) through itegratio 17), startig at maturity t = t N. At each time step we evaluate the itegratio 17) for every grid poit by a high accuracy umerical quadrature. Uder the ew probability measure Q, the process for lwt) ad rt) betwee the withdrawal dates is a simple Gaussia process give by 4) ad 16), where the coditioal joit desity of lwt ),rt )) give lwt + 1 ) = x, rt 1 ) = r is a bivariate Normal desity fuctio, as show i Appedix A, with the mea, variace ad covariace give by ) µ r r ) : = meart )) = r e κ + θ σ2 r κ 2 b + σ2 r 2κ 2a ; 19a) τ 2 r : = varrt )) = σ2 r 2κ a ; µ x x,r ) : = mealwt )) = x + b κ r + b σr 2 2κ 2 ρσ Sσ r κ 2 κ b ) ) + θ σ2 r κ 2 α+ 1 2 σ2 S ) b ) κ ) ; 19b) 19c) 19d) τ 2 x : = varlwt )) = σ2 S + σ2 r 2κ 32κ 4b +a )+ 2ρσ Sσ r κ 2 κ b ); 19e) ρ xr : = covlwt ),rt )), covlwt τ x τ ),rt )) = ρσ Sσ r b + σ2 r r κ 2κ 22b a ), 19f) 13

14 where b = 1e κ, a = 1e 2κ ad = t t 1. For simplicity, here we omit time step idex i otatio for the meas ad covariaces. Thus the desity of Y 1 = lwt ) µ x )/τ x ad Y 2 = rt ) µ r )/τ r is the stadard bivariate Normal with zero meas, uit variaces, ad correlatio ρ xr. If we apply the chage of variables Z 1 = Y 1 21ρ 2 xr ), Z 2 = Y 2 21ρ 2 xr ) 2) the itegratio 17) becomes Q t + W,r,A) 1 = Pt 1,t ) 1ρ 2 xr π ) e z2 1 z2 2 2ρxrz 1z 2 Q Z) z t 1,z 2,A)dz 1 dz 2 which has the form suitable for itegratio usigthe Gauss-Hermite quadrature. Here, Q Z) t ) deotes Q t ) as a fuctio of Z 1 ad Z 2 after trasformatio from W ad r. For a arbitrary oe-dimesioal fuctio fx), the Gauss-Hermite quadrature is applied as + e x2 fx)dx q i=1 λ q) i fξ q) i ), 22) where q is the order of the Hermite polyomial, ξ q) i, i = 1,2,...,q are the roots of the Hermite polyomial H q x), ad the associated weights λ q) i are give by λ q) i = 2 q1 q! π q 2 H q1 ξ q) i ) This approximatio is exact whe fx) ca be represeted as polyomial of the order up to 2q 1. I geeral, the abscissas ξ q) i ad the weights λ q) i for the Gauss-Hermite quadrature for a give order q ca be readily computed, e.g. usig fuctios i Press et al. 1992), or available i precalculated tables. Decomposig the two-dimesioal itegratio i 21) ito ested oe-dimesioal itegratio ad applyig the oe-dimesioal Gauss-Hermite quadrature to each of the variable, we obtai ) e z2 1 z2 2 2ρxrz 1z 2 Q Z) z t 1,z 2,A)dz 1 dz 2 q 1,q 2 i=1,j=1 λ q 1) i λ q 2) j e 2ρxrξq 1) i ξ q 2 ) j Q Z) ξ q 1) t i,ξ q 2) j,a). 23) Note, i geeral, differet orders q 1 ad q 2 ca be used. For example we might let q 1 > q 2 to take ito cosideratio that the value fuctio Q Z) t z 1,z 2, ) chages more rapidly with z 1 tha with z 2, thus makig the quadrature poits more efficietly assiged to the variables. Ufortuately, theumerical itegratio 23) is ot efficiet dueto presece of the factor e 2ρxrz 1z 2, whe there is a o-zero correlatio betwee stock market ad the iterest rate. This ca be improved 14

15 by trasformatio to idepedet radom variables Z 1,Z 2 ) usig the stadard Cholesky trasformatio Y 1 = 2Z 1, Y 2 = 2ρ xr Z 1 + 1ρ 2 xrz 2 ). 24) However, we observed that more accurate results are obtaied by a trasformatio to idepedet variables Z 1, Z 2 ) correspodig to the pricipal axis of the joit desity, which ca be doe usig a matrix spectral decompositio Y 1 = 2aZ 1 +bz 2 ), Y 2 = 2bZ 1 +az 2 ), 25) where a = ρ xr + 1ρ xr ), b = ρ xr 1ρ xr ). Usig trasformatio 25), ot oly the cross term i the desity disappears, but also it is stadardized for applyig the Gauss-Hermite quadrature. Now, i terms of the ew variables Z 1,Z 2 ), the itegratio 23) chages to a simpler but more accurate approximatio + + e z2 1 z2 Z) 2 Q z t 1,z 2,A)dz 1 dz 2 q 1 q 2 i=1 j=1 ) λ q 1) i λ q 2) j Q Z) ξ q 1) t i,ξ q 2) j,a. 26) If we apply the chage of variable ad the Gauss-Hermite quadrature 26) as described above to every grid poit W m,r k,a j ), m =,1,...,M, k =,1,...,K ad j =,1,...,J, i.e. let Wt + 1 ) = W m, rt + 1 ) = r k ad A = A j, the the cotract values at time t = t + 1 for all the grid poits ca be evaluated. As is commoly practiced, we select the workig domai i the asset space to be i terms of X = lw/w)), i.e. we set X mi = lw mi /W)) ad X max = lw max /W)). The domai [X mi,x max ] is uiformly discretised with step δx = X max X mi )/M to yield the grid X m = X mi + mδx, m =,...,M. The grid poits W m, m =,1,...,M, are the give by W m = W)expX m ). The domai [r mi,r max ] is also uiformly discretised with step δr = r max r mi )/K to yield the grid r k = r mi +kδr, k =,...,K. For each grid poit X m,r k ), the cotract value at time t + 1 ca be expressed as the weighted sum of some cotract values at time t. Specifically, from 17), 25) ad 26) we have Q t + W m,r k,a) Pt 1,t ) 1 π w ijkm = exp 2τx aξ q 1) q 1,q 2 i=1,j=1 i +bξ q 2) j r ijk = 2τ r bξ q 1) i +aξ q 2) j λ q 1) i λ q 2) j Q t w ijkm,r ijk,a), 27) ) ) +µ x X m +lw),r k ), ) +µ r r k ). We foud that it is more efficiet to let q 1 > q 2. This is because a > b, i.e. Z 1 cotributes to W more tha Z 2, ad thus assigig a larger umber of quadrature poits to Z 1 is more efficiet. 28) 15

16 5.3 Cubic splie iterpolatios for itegratio ad jump coditio At time step t + 1,t ), the cotract value at t = t for ay give A is kow oly at the grid poits W m,r k ), m =,...,M, k =,...,K. I order to approximate the cotiuous fuctio Q t W,r, ) required for the itegratio, we propose to use the bi-cubic splie iterpolatio over the grid poits, which is smooth i the first derivative ad cotiuous i the secod derivative. The error of cubic splie is Oh 4 ), where h is the size for the spacig of the iterpolatig variable, assumig a uiform spacig. The cubic splie iterpolatio ivolves solvig a tri-diagoal system of liear equatios for the secod derivatives at all grid poits. For a fixed grid ad costat i time model parameters, the tri-diagoal matrix ca be iverted oce ad at each time step oly the back-substitutio i the cubic splie procedure is required. For uiform grids, the bi-cubic splie is about five times as expesive i terms of computig time as the oe-dimesioal cubic splie, as explaied below. Let Q X) t ) deote Q t ) as a fuctio of X = lw/w)). Suppose the itegratio requires the value Q X) X,r, ) at the poit X,r) located iside a grid: X t m X X m+1 ad r k r r k+1. Because the grid is uiform i both X ad r, the secod derivatives 2 Q X) / X 2 ad 2 Q X) / r 2 ca be accurately approximated by the three-poit cetral differece, ad cosequetly the oe-dimesioal cubic splie o a uiform grid ivolves oly four eighborig grid poits for ay sigle iterpolatio. I our bi-cubic splie case, we ca first obtai Q X) at four poits X,r k1 ), X,r k ), X,r k+1 ), X,r k+2 ) by applyig the oe-dimesioal cubic splie o the dimesio X for each poit ad the we use these four values to obtai Q X) X,r, ) through a oe-dimesioal cubic splie i r. Thus five oe-dimesioal cubic splie iterpolatios are required for a sigle poit, which ivolves sixtee grid poits eighborig the poit of iterest X,r). Needless to say, this aloe will make the pricig of GMWB uder stochastic iterest rate much more time cosumig tha GMWB uder determiistic iterest rate. Not oly the evaluatio per grid poit is more ivolved but also evaluatio has to be performed for larger umber of poits. To apply the jump coditios, oly oe-dimesioal cubic splie iterpolatio is ivolved, as the poitsbeforeadafter each jumpfall ogridpoitsir ada, olyiterpolatio ix isrequired. Let us itroducea auxiliary fiite grid = A 1 < A 2 < < A J = W) to track theremaiig guaratee balace A, where J is the total umber of odes i the guaratee balace amout coordiate. The upper limit W) is eeded because the remaiig guaratee balace caot exceed the target iitial accout value W). For each A j, we associate a cotiuous solutio Q t W,r,A) defied by the values at ode poits W m,r k ) ad a two-dimesioal cubic splie iterpolatio over these ode poits. At every jump we let A to be oe of the grid poits A j, 1 j J. Amog the ifiite umber of possible jumps, a most efficiet choice though ot ecessary) is to oly allow the guaratee balace to be equal to oe of the grid poits = A 1 < A 2 < < A J = W). This implies that, for a give balace A j at time t, the possible value after the withdraw at t + has to be oe of the grid poits equal to or less tha A j, i.e. A + j = A i, 1 i j. I other words, the withdrawal amout γ takes j 16

17 possible values: γ = A j A i, i = 1,...,j. Note the above restrictio that γ = A j A i, i = 1,...,j is ot ecessary. The oly real restrictio is γ A j. However, without the restrictio, the value of A + j after the jump falls betwee the grid poits ot exactly o a grid poit A k ) ad a costly two-dimesioal iterpolatio is required. The error due to this discretisatio restrictio ca be easily reduced to acceptable level by icreasig J. For ay ode W m,r k,a j ), m =,1,...,M, k =,1,...,K, j = 1,...,J, give that withdrawal amout ca oly take the pre-defied values γ = A j A i, i = 1,2,...,j, irrespective of time t ad accout value W m, the jump coditio 13) takes the followig discrete form [ ] Q t W m,r k,a j ) = max Q 1 i j t + maxw m A j +A i,),r k,a i )+C A j A i ). 29) For optimal strategy, we chose a value for 1 i j maximizig Q t W m,r k,a j ). The above jump has to be performed for every ode poit W m,r k,a j ), m M, k =,1,...,K, 1 j J at every withdrawal date. Obviously for every ode poit W m,r k,a j ) we have to attempt j jumps to fid the maximum value for Q t W m,r k,a j ). Figure 1 illustrates applicatio of the jump coditio. W W m m 1 t t A j Q t Wm k j, r, A ) A j W A j 1 A j 1 A Wm Wm 1 t t W A i Q t Wm k i, r, A) Ai A j A j 1 Aj 1 A W m W A A m j i Wm 1 Figure 1: Illustratio of the jump coditios applied o the fiite differece grids. Whe W m A j +A i >, the value Q t + W m A j +A i,r k,a i ) ca be obtaied by a oe-dimesioal cubic splie iterpolatio from the values at the M discrete grid poits. Iterpolatio scheme is importat, as show for example i a covergece study by Forsyth et al. 22), it is possible for a PDE 17

18 based umerical algorithm for discretely sampled path-depedet optio pricig to be o-coverget or coverget to a icorrect aswer) if the iterpolatio scheme is selected iappropriately. 6 Numerical Results I this sectio we first show results for a bechmark test where closed-form solutio exists. We the preset umerical results for pricig GMWB uder the static ad optimal policyholder strategies usig GHQC algorithm ad compare these with the MC ad PDE fiite differece results whe appropriate. 6.1 Vailla Europea optios I the case of stochastic dyamics 1) for the uderlyig asset ad iterest rate, there is a closedform solutio for Europea vailla optios, thus providig a valuable bechmark test for umerical algorithms. I particular, the formulas for prices of vailla call ad vailla put with strike K T ad maturity T, give S) ad r) at time t = are derived i Appedix A.5. For this test we set the iput as follows: asset volatility σ S = 2%, asset spot value S) = 1, iterest rate spot value r) = 5%, maturity T = 1, strike K T =.95, ad the Vasicek iterest rate model parameters κ =.349, θ = 5%, σ r = 1%, 3% ad ρ =.2,.,.2. We calculated the vailla prices usig the fiite differece ADI method solvig two-dimesioal PDE 18) ad our GHQC method developed i the previous sectios ad compared with the closedform solutio 37). Of course i the case of vailla optios, usig chage of uméraire to the bod P, T), the iterest rate dimesioality ca be removed from umerical pricig ad the required PDE ca be reduced to the oe-dimesioal PDE similar to determiistic iterest rate case. Here, we implemet ADI for the origial two-dimesioal PDE 18) for testig ad compariso purposes. For GHQC, we used q 1 = 12 ad q 2 = 3 quadrature, i.e. the total umber of quadrature poits for each itegratio is 36. The mesh for GHQC calculatios was fixed at M = 1 for X dimesio ad K = 2 for the iterest rate r dimesio, ad the total umber of time steps N = 5. Comparig with typical fiite differece calculatios, the above mesh ad time steps are quite coarse. Ideed, for ADI calculatios, i order to have a roughly compatible accuracy with GHQC, we had to set M = 2, K = 4 ad N = 3. Table 1 ad Table 2 show results for vailla call ad put prices for differet values of σ r ad ρ. The percetage umbers i the paretheses i both tables are the relative umerical errors of the price compared with the closed-form solutio. O average, for σ r = 1% the relative error is about.33% for ADI ad.42% for GHQC, while for σ r = 3% the relative error is about.77% for ADI ad.4% for GHQC. Sigificatly, the average relative error for ADI is more tha doubled whe the volatility of iterest rate is icreased from σ r = 1% to σ r = 3%, while the error for GHQC remais more or less the same. Both ADI ad GHQC took a fractio of a secod CPU to calculate oe of the optios i Table 1 ad Table 2. Averagig over 2 calculatios for call ad put optios with the same iputs as give 18

19 Table1: Vailla call optio price for differet values of σ r ad ρ. The other iput parameters are σ = 2%, S) = 1., r) = 5%, T = 1., K T =.95, κ =.349 ad θ = 5%. σ r ρ Closed-form ADI GHQC %) %) %) %) %) %) %) %) %) %) %) %) Table2: Vailla put optio price for differetvalues of σ r ad ρ. The other iput parameters are σ = 2%, S) = 1., r) = 5%, T = 1., K T =.95, κ =.349 ad θ = 5%. σ r ρ Closed-form ADI GHQC %) %) %) %) %) %) %) %) %) %) %) %) above, we foud the CPU time for each call or put calculatio is.55 secod for ADI ad.11 secod for GHQC, i.e. GHQC is about five times faster tha ADI i these tests. All the calculatios show i this study were performed o a desktop with a Itel Core i5-459 CPU@3.3GHz with a 4.GB RAM. It is worth commetig that for the vailla call ad put optios the fial payoff fuctio is oly piecewise liear, i.e. it is ot a polyomial fuctio ad it is ot smooth at the strike K T first derivative discotiuous at W = K T ). If we apply the Gauss-Hermite quadrature oly to the half domai W K T for call ad W K T for put), the the payoff fuctio at the maturity is a simple liear fuctio eedig o iterpolatio aywhere ad we foud the GHQC calculatios ca be made as accurate as desired, virtually limited oly by the machie accuracy. I the above tests we did ot take advatage of this specific feature because it is ot geerally applicable. 19

20 6.2 GMWB pricig results I the case of the static policyholder behavior, the withdrawal amouts are predetermied at the begiig of the cotract. I this case the paths of the wealth accout Wt) ca be simulated ad a stadardmc simulatio method ca beusedto calculate the priceof VA cotract with GMWB. Below we show results for both static ad optimal withdrawal cases. The static case allows a compariso betwee MC ad GHQC, further validatig the ew algorithm. We have also implemeted a efficiet fiite differece ADI algorithm for pricig of VA with GMWB uder the static policyholder behavior solvig correspodig two-dimesioal PDE. I what follows, results from the GHQC method will be compared with the MC ad ADI methods whe applicable GMWB with static policyholder behavior I a static case the withdrawal amout is pre-determied for each withdrawal date. I this case at each paymet date the jump coditio applies to the sigle solutio therefore o eed for a grid i the guaratee accout A dimesio). Sice the withdrawal amout is kow at every paymet date, the stochastic paths of the uderlyig W ca be simulated by MC method. Here we compare GHQC results with those of MC ad ADI methods. Table 3 shows the prices of the VA with GMWB uder the static withdrawal strategy as a fuctio of the correlatio ρ, comparig results betwee MC, ADI ad GHQC. The model iput parameters are σ S = 2%, σ r = 2%, S) = W) = 1., r) = 5%, g = 1%, T = 1/g = 1., N w = 4 quarterly withdraw frequecy), α =.6, κ =.349 ad θ = 5%. I Peg et al. 212), results for static GMWB pricig usig similar parameters were preseted for the lower ad upper boud estimates of the price comparig with the MC simulatio results. However, those results are for a cotiuous i time withdrawal case, while our results are for the discrete i time withdrawals at a quarterly frequecy. For the MC method, we have used oe millio sample paths simulated from the closed-form trasitio desity uder the ew measure Q, so that there is o time discretizatio error. For both ADI ad GHQC we have used two meshes, oe coarser ad oe fier, with the fier mesh doubles the umber of ode poits i both r ad X dimesios. Let M A deote the coarser mesh for ADI with M = 2 for X dimesio, K = 5 for r dimesio, ad MA deote fier mesh for ADI with M = 4 ad K = 1. I additio, the umber of time steps for each period betwee cosecutive withdrawal dates) was set N t = 1 for ADI whe usig the coarser mesh, ad it was doubled to N t = 2 whe usig the fier mesh MA. For GHQC, the coarser mesh, deoted as M G, has M = 5 ad K = 3, ad the fier mesh for GHQC, deoted as MG, has M = 1 ad K = 6. For both M G ad MG meshes, we have used a sigle time step betwee withdrawal dates, i.e. N t = 1. For the quadrature poits, we used q 1 = 5 ad q 2 = 3 for the coarser mesh, ad q 1 = 9 ad q 2 = 5 for the fier mesh. Note the fier mesh MG for GHQC is overall much coarser tha the coarser mesh M A for ADI. 2