On the goodness of t of Kirk s formula for spread option prices

Size: px
Start display at page:

Download "On the goodness of t of Kirk s formula for spread option prices"

Transcription

1 On he goodne of of Kirk formula for pread opion price Elia Alò Dp. d Economia i Emprea Univeria Pompeu Fabra c/ramon ria Farga, Barcelona, Spain Jorge A. León y Conrol Auomáico CINVESAV-IPN Aparado Poal México, D.F., Mexico Abrac In hi paper we inveigae he goodne of of he Kirk approximaion formula for pread opion price in he correlaed lognormal framework. oward hi end, we ue he Malliavin calculu echnique o nd an expreion for he hor-ime implied volailiy kew of opion wih random rike. In paricular, we obain ha hi kew i very pronounced in he cae of pread opion wih exremely high correlaion, which canno be reproduced by a conan volailiy approximaion a in he Kirk formula. hi fac agree wih he empirical evidence. Numerical example are given. Keyword: Spread opion, Kirk formula, Malliavin calculu, Skorohod inegral. JEL code: G, G3 Mahemaical Subjec Clai caion: 9B8, 9B70, 60H07. Inroducion A pread opion i a derivaive wrien a he di erence of wo underlying ae. Namely, he payo of a call pread opion wih rike K wih ime o mauriy i (S S K) ; where S and S denoe wo ae price. I i well-known ha if S and S are wo geomeric Brownian moion (ha may be correlaed) and K = 0; he correponding opion price i given by he Margrabe fomula (ee Margrabe (978)), which can be deduced from he fac ha S =S i a log-normal proce. hu, in hi cae, he pread opion value can be expreed a he claical Black-Schole call price wih iniial ae price S 0; where we ake he rike Suppored by he gran ECO and MEC FEDER MM y Parially uppored by a CONACy gran

2 equal o he expeced value of S q and volailiy equal o 0 ( 0 ). Here and 0 are he volailiy parameer of S and S, repecively, and denoe he correlaion. In he cae K 6= 0 he fracion S = S K i no log-normal and hen he argumen ued in he deducion of Margrabe formula canno be applied anymore. One propoed oluion i o aume he diribuion of he pread i able o be approximaed by he normal diribuion. hi lead u o he Bachelier mehod (ee Shimko (994)). Even hough he approximaion of he pread diribuion by he normal one i poor, hi mehod can be precie in ome range of parameer (ee, for example, Carmona and Durrleman (003)). A more ucceful mehod i uggeed by Kirk (995), who applied he Margrabe formula, aproaching S = S K by a log-normal random proce. Nowaday Kirk formula i he mo popular opion pricing approximaion expreion for pread opion due o i accuracy and i impliciy. Recenly, ome oher auhor have propoed new formula and mehod o eimae pread opion price. Among hem, we can menion Alexander and Venkaramanan (0), Bjerkund and Senland (0), Borovkova e al. (007), Carmona and Durrleman (003), and Deng e al. (008). Numerical evidence ha hown ha hee approache may improve Kirk eimae, pecially for high correlaed cae (ee, for inance, Bjerkund and Senland (0)). One inereing ool in he developmen of an accurae approximaion formula i he knowledge of he properie of he correponding implied volailiy. For example, in he cae of vainilla opion wih ochaic volailiy, we know ha implied volailiie exhibi mile and kew (ee Renaul and ouzi (996)). hen, we have o ake ino accoun ha accurae approximaion hould reproduce hem and, moreover, we can underand why ome formula fail (for example, we expec ha a conan volailiy expreion will fail for in or ou-of-he-money opion). Up o our knowledge, here are no imilar analyical udie in he cae of opion wih random rike (a for example pread opion). he main goal of hi paper i o conider analyically he relaion beween he implied volailiy and he ock price S for opion wih random rike. By mean of he Malliavin calculu hecnique, we develop an exenion of he Margrabe formula ha allow u o nd an expreion for he hor-ime kew lope for pread opion. hi analyi of he implied volailiy kew evidence ha he dependence beween he implied volailiy and he ock price S increae rongly when he correlaion parameer i cloe o. Hence we can expec ha Kirk approximaion reduce i accuracy for highly correlaed ae due o he fac ha in hi approximaion he volailiy i conan a a funcion of he ock price S. hi agree wih he numerical evidence (ee, for example, Baeva (0)). he organizaion of hi paper i a follow. In Secion we inroduce he framework of hi paper. Secion 3 i devoed o preen a exended Margrabe formula for he cae K 6= 0. hi formula i ued in Secion 4 o gure up an expreion for he derivaive of he implied volailiy wih repec o he logock price. he hor ime behaviour of hi derivaive i analyzed in Secion

3 5. Finally in Secion 6 we obain our reul on he goodne of he of Kirk formula a an applicaion. Saemen of he model and noaion In hi paper we conider he following model for he log-price of a ock under a rik-neural probabiliy meaure Q: dx = r d (dw p db ); [0; ]: () Here, r i he inananeou inere rae, W and B are independen andard Brownian moion and ( ; ). For he ake of impliciy, we aume ha he volailiy proce i a quare-inegrable deerminiic funcion which i righ-coninuou. In he following we denoe by F W and F B he lraion generaed by W and B; repecively. We de ne F := F W _ F B : In hi paper we conider European call opion wih payo h(x ) := e X K ; where we allow he rike K o be random. More preciely, we aume K i a quare-inegrable, poiive, coninuou, bounded and F W - meaurable proce. Noice ha hi choice include ome popular clae of opion a bake one. I i well-known ha he price of an European call wih rike K i given by he formula V = e r( ) E (e X K ) jf : () In he equel, we will make ue of he following noaion: M := E K j F W : Oberve ha, by he maringale repreenaion heorem, M = E (K ) 0 m(; )dw ; for ome F W meaurable and adaped proce m (; ) : v := Y ; wih Y := R a d; where a d := d dhm ;Xi M dhm ;M i (M ). Noe ha a = m(; ) M m (; ) (M ) = m(; ) M i a poiive quaniy. Alhough he righ-hand-ide of he la equaliy depend on, we denoe i by a in order o implify he noaion. I i eay o ee ha for all ( ; ), a C () : 3

4 BS(; x; K; ) denoe he price of an European call opion under he claical Black-Schole model wih conan volailiy, curren log ock price x, ime o mauriy ; rike price K and inere rae r: Remember ha in hi cae: BS(; x; K; ) = e x N(d ) Ke r( ) N(d ); where N denoe he cumulaive probabiliy funcion of he andard normal law and wih ~x := ln K r( ): d := x ~x p p ; L BS and for he Black-Schole di erenial operaor, in he log variable, wih volailiy : L xx (r )@ x r I i well known ha L BS BS(; ; ) = 0: Now we decribe ome baic noaion ha i ued in hi aricle. For hi, we aume ha he reader i familiar wih he elemenary reul of he Malliavin calculu, a given for inance in Nualar (006). Le u conider a andard Brownian moion = f ; [0; ]g de ned on a complee probabiliy pace (; F; P ): he e D ; i he domain of he derivaive operaor D in he Malliavin calculu ene. D ; i a dene ube of L () and D i a cloed and unbounded operaor from L () ino L ([0; ] ): We alo conider he ieraed derivaive D ;n ; for n > ; whoe domain i denoed by D n; W : he adjoin of he derivaive operaor D, denoed by ; i an exenion of he Iô inegral in he ene ha he e L a([0; ] ) of quare inegrable and adaped procee (wih repec o he lraion generaed by ) i included in Dom and he operaor rericed o L a([0; ] ) coincide wih he Iô inegral. We make ue of he noaion (u) = R 0 u d : We recall ha L n; := L ([0; ]; D n; W ) i included in he domain of for all n : 3 A decompoiion reul Before proving an exenion of he Hull and Whie formula, we ae he following reul, which i nedeed in he remaining of he paper. 4

5 Lemma Le K be bounded, 0 < and G := F _ F W. Aume ( ; ) :hen, for any n 0, here exi C = C(n; ) uch ha E (@ n )BS(; X ; M ; v )jg x n C d (n) : Proof: In order o how hi reul, we proceed a in he proof of Lemma 4. in Alò, León and Vive (007) and we ue he fac ha K i a bounded and F W meaurable and adaped proce o obain ha E (@ n )BS(; X ; M ; v )jg x n C ( ) d a d (n) : We know ha for all ( ; ), a C () for ome poiive conan C (). hen R a d C () R d; from where he reul follow. Now we are able o prove he main reul of hi ecion, he exended Hull and Whie formula. We will need he following hypohei: (H) he proce a L ; : heorem Conider he model () and aume ha hypohei (H) hold. hen i follow ha V = E BS(; X ; M ; v ) F E ( where W := e r( 3 xx BS(; X ; M ; v ) W d ) e r( x BS(; X ; M ; v ) W m(; )d F ; (3) h D W R i a (r)dr : Proof: hi proof i imilar o he one of he main heorem in Alò, León and Vive (007), o we only kech i. Noice ha BS(; X ; M ; v ) = V : hen, from (), we have e r V = E(e r BS(; X ; K ; v )jf ): Now, uing he Iô formula o he proce e r BS(; X ; M ; v ) 5

6 and proceeding a in Alò, León and Vive (007) (ee alo Alò and Nualar (998), Alò (006) or Nualar (006)), we can wrie e r BS(; X ; M ; v ) = e r BS(; X ; M ; v ) e r L BS (v )BS(; X ; M ; v )d e x BS(; X ; M ; v ) (dw p db ) e K BS(; X ; M ; v )dm e xkbs(; X ; M ; v )d M ; X e BS(; X ; M ; v ) v v ( a ) d e xbs(; X ; M W ; v ) v ( ) d e KBS(; X ; M ; v ) W m(; ) v ( ) d e x BS(; X ; M ; v ) v d e KKBS(; X ; M ; v )d M ; M : 6

7 Hence, he fac ha L BS (v)bs(; X ; M ; v ) = 0; muliplying by e r and aking condiional expecaion we can eablih E e r( ) BS(; X ; M ; v ) F ( = E BS(; X ; M ; v ) e r( e r( BS(; X ; M ; v ) v a v ( ) d e r( xbs(; X ; M W ; v ) v ( ) d e r( KBS(; X ; M ; v ) W m(; ) v ( ) d e r( xkbs(; X ; M ; v )d M ; X e r( x BS(; X ; M ; v ) v d KKBS(; X ; M ; v )d M ; M ) F : Conequenly, he claical relaionhip beween he x BS BS ( xkbs BS K( KKBS BS K ( ) give E e r( ) BS(; X ; M ; v ) F ( = E BS(; X ; M ; v ) e r( BS(; X ; M ; v ) M v ( e r( BS(; X ; M ; v ) v a v ( ) d e r( xbs(; X ; M W ; v ) v ( ) d e r( KBS(; X ; M ; v ) W m(; ) v ( ) d e r( BS(; X ; M ; v ) v v ( e r( BS(; X ; M ; v ) M v ( ) d ) d M ; M ) F : ) d M ; X 7

8 ha i, E e r( ) BS(; X ; M ; v ) F ( = E BS(; X ; M ; v ) " d M ; X M e r( BS(; X ; M ; v ) v ( ) v a d v d e r( xbs(; X ; M W ; v ) v ( ) d ) e r( KBS(; X ; M ; v ) W m(; ) v ( ) d F : d M ; M (M ) Since, a d := d dhm ;Xi dhm ;M i we obain M (M ) E e r( ) BS(; X ; M ; v ) F ( = E BS(; X ; M ; v ) e r( xbs(; X ; M W ; v ) v ( ) d ) e r( KBS(; X ; M ; v ) W m(; ) v ( ) d F ; a we waned o prove. Example 3 Aume he model () wih conan volailiy : We conider a call pread opion wih rike equal o K = S 0 K, where K i a non-negaive deerminiic conan and S 0 i anoher ock price of he form S 0 = exp (X) 0 ; where dx 0 ( 0 ) = r d 0 dw ; [0; ]; for ome poiive conan 0. hen we can eaily check ha m (; ) = exp(r( ))S 0 0 ; M = exp(r( ))S0 K; and hen a := 0 exp(r( ))S 0 exp(r( ))S 0 K (0 ) (exp(r( ))S 0 ) (exp(r( ))S 0 K) : Noice ha, if K = 0 and D W a := 0 ( 0 ) a () = 0: hen, he equaliy (3) reduce o V = BS ; X ; exp (r ( )) S; q 0 0 ( 0 ) ; from where we recover he well-known Margrabe formula (ee Margrabe (978)). # 8

9 Remark 4 Noice ha, in he conex of he previou example, when K i negaive, he call opion on he pread S S 0 i equivalen o he correponding pu opion on he pread S 0 S wih poiive rike K. hen, wihou lo of generaliy, we can aume ha he pread opion i wrien wih a poiive K. 4 Derivaive of he implied volailiy Le I (X ) denoe he implied volailiy proce, which ai e by de niion V = BS(; X ; M ; I (X )): In hi ecion we prove a formula for i a-hemoney derivaive ha we ue in Secion 5 o udy he hor-ime behavior of he implied volailiy. Propoiion 5 Aume ha he model () hold wih a L ; and ha, for R every xed [0; ) ; d <. hen i follow (x ) = E( R e r( ) (@ x F (; X ; M ; v ) F (; X ; M ; v BS(; x ; M ; I (x )) where X=x ; a.. F (; X ; M ; v ) 3 xx BS(; X ; M ; v ) x BS(; X ; M ; v ) W m(; ) and x = ln(m ) r( ): Proof: Uing heorem and he expreion V = BS(; X ; M ; I (X )) we x BS(; X ; M ; I (X BS(; X ; M ; I (X ): (4) V = E(BS(; X ; M ; v )jf ) E( e r( ) F (; X ; M ; v )djf ); which implie = E(@ x BS(; X ; M ; v )jf ) E( e r( x F (; X ; M ; v )djf (5) We can check ha he condiional expecaion E( R e r( x F (; X ; M ; v )djf )i R well de ned and nie a.. due o he fac ha d <. hu, (4) 9

10 and (5) imply (x ) BS(; x ; M ; I (x BS(; x ; M ; v )jf x BS(; x ; M ; I (x )) )) E( e r( x F (; X ; M ; v )djf )# : X=x Noice ha E(@ x BS(; x ; M ; v )jf ) x E(BS(; x; M ; v )jf ) x=x x BS(; x; M ; I 0 (x))j x=x ; (7) where, by he Hull and Whie formula, I 0 (X ) i he implied volailiy of call opion wih conan rike M, for a cerain ochaic volailiy model where = 0 and he volailiy proce i given by a. x (BS(; x; M ; I 0 (x)) x=x x BS(; x ; M ; I 0 (x BS(; x ; M ; I 0 (x ) : (8) From Renaul and ouzi (996) we (x ) = 0: hen, (6), (7) and (8) imply ha (x @ BS(; x ; M ; I (x BS(; x ; M ; I 0 (x x BS(; x ; M ; I (x )) )) E( e r( x F (; X ; M ; v )djf )# : X=x On he oher hand, raighforward calculaion lead u o and x BS(; x ; M ; ) = e x N( p ) BS(; x ; M ; ) = e x (N( p ) N( p x BS(; x ; M ; ) = (ex BS(; x ; M ; )); which yield 0

11 @ x BS(; x ; M ; I 0 (x x BS(; x ; M ; I (x )) = (BS(; x ; M ; I 0 (x )) BS(; x ; M ; I (x ))) = E(BS(; x ; M ; v ) V jf ) = E( e r( ) F (; X ; M ; v )djf ) hi, ogeher wih (9), implie ha he reul hold. 5 Shor-ime behaviour he pourpoe of hi ecion i o udy he limi (x ) a # : he following reul i par of he ool needed for our reul. : Lemma 6 Aume he model () i ai ed. hen I (x ) p : 0 a.. a Proof: Noice ha he fac ha K i a quare-inegrable and coninuou random proce and he dominaed convergence heorem lead o ge V j X=x = E(e r( ) (e X K ) jf ) X=x = E (e r( ) (e X X e r( ) M K ) jf ) X=x E ((e X X M K e r( ) ) jf ) X=x = E ((e X X e r( ) )M e r( ) (M K )) jf ) X=x E (je X X e r( ) jm jf ) X=x E (jm K je r( ) jf ) X=x M E (je X X e r( ) jjf ) X=x E (jm K je r( ) jf ) 0 a::; X=x a. Hence, aking ino accoun ha, in he a-he-money cae, V j X=x BS(; x ; M ; I (x )); we deduce ha = I(x BS(; x ; M ; I (x )) = M e N r( ) ) p 0 a::;

12 and hi allow u o complee he proof. Henceforh we conider he following hypohee: (H ) a L ; and, moreover, here exi a poiive conan C uch ha, for all 0 < < < r < ; D W a r D W D W a r C: Noice ha hi hypohee implie ha (H) hold. (H) here exi wo poiive conan c ; c uch ha for all r [0; ] c r c : Noice ha, a for all ( ; ), a C () for ome poiive conan C () ; hi hypohei implie ha a i lower bounded. (H3) he proce m(; ) L ; and moreover, here exi a poiive F proce C uch ha for all > > r > ; E jm(; r)j D W F E m(; r) F C : adaped Propoiion 7 Aume ha he model () and Hypohee (H )-(H3) hold. Alo aume ha here i a conan c > 0 uch ha c < K ; for all [0; ]. BS(; x ; M ; I (x (x = xxx xx BS(; x ; M ; v x BS(; x ; M ; v ) W a. O( ): Proof: Propoiion 5 give u BS(; x ; M ; I (x ) W d m(; )d F = E e r( ) (@ x xx BS(; X ; M ; v ) W e r( ) (@ x )@ x BS(; X ; M ; v ) W m(; )d F =: : X =x Now he proof i decompoed ino wo ep. d

13 Sep. Here we ee ha = E L(; x ; M ; v ) W djf O ( ) ; (0) where L(; X ; M ; v ) 3 xx BS(; X ; M ; v ): In fac, applying Iô formula o e r L(; X ; M ; v )( r W r a in he proof of heorem and aking condiional expecaion wih repec o F ; we obain ha dr) E( e r( ) L(; X ; M ; v ) W djf ) = E L(; X ; M ; v )( W d)jf 4 E( e r( ) (@ xxx xx)l(; X ; M ; v ) W r W r dr djf ) 4 E( e r( K (@ x )L(; X ; M ; v ) W m(; ) E( W r e r( r dr djf ) (D W W r ) r dr x L(; X ; M ; v ) djf )j E( e r( K L(; X ; M ; v ) m(; ) (D W W r ) r dr = E L(; X ; M ; v )( S S S 3 S 4 : djf ) W d)jf 3

14 Uing he proof of Lemma and Hypohee (H ) and (H), we can wrie js j = 4 E( e r( ) E (@ xxx xx)l(; X ; M ; v ) G ( W r r dr) W djf ) 6X k 3 C E 4 a d j( W r r dr) W jdjf 5 k=4 " 6X # C E ( ) k j( W r r dr) W jdjf k=4 Hence, uing Hypohee (H ), (H), and (H3), we can wrie Similarly, we have js j = 4 E ( js j C herefore, he BS(; x; K; K 6X ( ) k 4 = O( ): k=4 e r( ) K (@ x )L(; X ; M ; v ) G W r r dr) W m(; )djf : = x; K; BS(; x; K; ) x ogheer wih he hypohee of he Propoiion, implie 6X js j C E ( ) 3 jm(; )j F = O( ): In a imilar way, k=3 js 3 j = E( C C 4X E k=3 4X k=3 e r( x L(; X ; M ; v ) (D W W r ) r dr djf )j X=x ( ) (D W W r ) r dr d F ( ) = O( ): 4

15 Finally, he ame argumen give u ha js 4 j = O( ): Sep. In order o nih he proof we only need o proceed a in Sep. Here we ee ha = E P (; x ; M ; v ) W m(; )djf O ( ) ; () where P (; X ; M ; v ) = (@ x )@ x BS(; X ; M ; v ): In fac, applying Iô formula o e r P (; X ; M ; v )( m(; r) W r a in he proof of heorem and aking condiional expecaion wih repec o F ; we obain ha E e r( ) P (; X ; M ; v ) W m(; )d F = E P (; X ; M ; v )( m(; ) W d)jf dr) 4 E( e r( ) (@ xxx xx)p (; X ; M ; v ) W m(; r) W r dr djf ) 4 E( e r( k (@ x )P (; X ; M ; v ) W m(; ) W r E( e r( m(; r)dr djf )j x P (; X ; M ; v ) D W W r m(; r) dr djf ) E( e r( K P (; X ; M ; v ) m(; ) D W W r m(; r) dr djf ): Now, following he ame argumen a in Sep he proof i complee. 5

16 Remark 8 hi proof only need ome inegrabiliy and regulariy condiion. So, depending on he coe cien of he model () and he proce K, Hypohee (H )-(H3) can be ubiued by appropiae inegrabiliy condiion. Now we can ae he main reul of hi paper. oward hi end, we need o ae he following aumpion: (H4) Aume ha m(; ) ha coninou pah and ha, for each [0; ] xed, up E a r a <^^r< ~a F 0 a ; a.. and m(; ) up E m(; )a r a <^^r< ~a F 0 a ; a.. where, by convenion, ~a := m(; ) K m (; ) K (H5) here exi a F -meaurable random variable D a uch ha, for every xed > 0; up E D W a r D a F 0; a.. <<r< a : heorem 9 Conider he model () and uppoe ha Hypohee (H )-(H5) hold and here exi a poiive conan c uch ha c < K: (x ) = m(; ) K D a : () ~a Proof: We can BS(; x ; M ; I (x )) = M e r( p p ; ) e I(x ) ( ) x 3 xx BS(; x ; M ; v ) = M e r( ) = M e r( ) p e x BS(; x ; M ; v ) ( ) 8 v ( ) 6 p e v ( ) M v : ( ) 8 v 3 ( ) 3

17 hen we can wrie, due o Lemma 6 and Propoiion (x ) e I(x ) ( ) 8 ( ) E(e v ( ) 8 v 3 e I (x ) ( ) 8 ( ) E e v ( ) 8 v W djf ) M v ( ) W m(; )djf O( ) = : S S O( ) : By Lemma 6, we know ha I (x ) ( ) 0 a.. a : hen, " lim S = # lim ( ) E(e v ( ) 8 v 3 W djf ) and lim S = lim Now, le u ee ha In fac, we can eablih where and lim ( ) E(e v W m(; )djf ) ( ) 8 v # M v ( ) : (3) lim S ~a D a = 0 a..: (4) S ~a D a = lim E v A := exp ( ) 8 A B ~a D a F B := v ( ) a r D W a r drd: Conequenly lim E A B ~a D a F = lim E A B F lim E B D a ~a ~a ~a F = lim U ~a lim U : 7 v

18 Applying Schwarz inequaliy for condiional expecaion, i follow ha U " E A # ~a F E B F From he dominaed convergence heorem and (H), i i eay o ee ha E A ~a F end o zero a.. a ; and a imple calculaion give u ha (H ) and (H) imply ha E B F i bounded, from where we deduce ha lim U = 0: Oberve ha we alo have, ju j = ( ) E a r v D W a r D a drd ~a F C ( ) E a r v D W a r drd ~a F C ( ) E D W a r D a F drd = : ju ; j ju ; j : Uing Hypohee (H ) and (H) we obain ha C ju ; j ( ) E a r v ~a drd F C ( ) E a r v ~a drd F C = ( ) E a r a ~a ( ) d drd F C ( ) 3 E a r a ~a F ddrd; which end o zero, a.. a, becaue of Hypohei (H4). Similarly, C ju ; j ( ) E D W a r D a F drd ; which end o zero by Hypohei (H5). hu we have proved (4) i rue. On he oher hand, by (3) we can wrie m(; ) lim S K ~a D m(; ) a = lim E A B K ~a D a F : 8

19 bu now and v A := exp ( ) 8 M v B := v ( ) a r m(; )D W a r drd: Finally, proceeding imilarly a before, we have (3) yield ha S converge o m(;) K ~a D a, which, ogeher wih (4), implie ha () i ai ed. 6 Applicaion o he udy of pread opion hi ecion i devoed o apply he previou reul o udy he implied volailiy behaviour for pread opion. hi udy allow u o predic when he Kirk approximaion formula for pread opion may fail. Moreover, we will ee how hi analyi give u a ool o improve Kirk formula. 6. Shor-ime behaviour of he implied volailiy for pread opion Conider an pread opion wih K = S 0 K a in Example 3. For he ake of impliciy we will aume he inere rae r = 0. hen i i eay o ee ha We can eaily check ha, for < a := 0 S 0 S 0 K (0 ) (S 0 ) (S 0 K) : D W a = 0 K (S 0 K) (0 ) S 0 S 0 K = ( 0 ) S 0 0 Hence, we deduce ha S 0 K S 0 K (S 0 K) : K (S 0 K) 0 S 0 q D W a = D W a = DW r a p a S = p 0 0 a S 0 ( 0 ) SK 0 K (S 0 K) : 9

20 hen, from heorem 9, we ge = lim (x m(; ) K 0 S 0 S 0 K D a ~a p 3 ( 0 ) SK 0 a (S 0 K) : (5) Remark 0 Noice ha he above quaniy i allway poiive. In he following example we will udy i behaviour a a funcion of K and Example In Figure we plo (x ) a a funcion of K for = 0:9 (olid) and = (dah), and wih S = 00; = 0:5 and 0 = 0:4: We oberve he limi kew (x ) i zero in he cae K = 0. hi wa expeced from Example 3, where we found ha in hi cae he implied volailiy i conan, = 0: Noice alo ha, even hi kew increae wih K, hi incremen eem o be clearly bigger in he compleely correlaed cae = : implied volailiy kew Figure (x ) a a funcion of K for = 0:9 (olid) and = (dah). Here = 0:5; 0 = 0:4: Example In Figure we plo (x ) a a funcion of for K = 5 (olid) and K = 0 (dah) and for he ame parameer value a in Fig.. can oberve he limi kew (x ) ha i maximum in he compleely correlaed cae = : Noice ha hi mean ha he conan volailiy approximaion given by Kirk formula i expeced o be le accurae in hi cae. hi fac i conien numerical empirical evidence (ee for example Baeva (0) and Borovkova (007)). 0

21 iimplied volailiy kew Figure (x ) a a funcion of for K = 5 (olid) and K = 0 (dah). Here = 0:5; 0 = 0:4: Example 3 In Figure 3 we (x ) a a funcion of and K for he ame parameer value a in Fig. and Fig.. Noice ha hi limi kew i ubanially bigger near he cae = : ro z K correlaion Figure (x ) a a funcion of and K: 6. Applicaion o he udy of he accuracy of Kirk formula Wih he above noaion, he Kirk approximaion for an pread opion can be wrien a: q BS(; X ; M ; a ): I i well-known ha Kirk formula i a very accurae approximaion for pread opion given i impliciy (ee for example Baeva (0), Bjerkund and Senland (0) or Carmona and Durrleman (0)). Neverhele, i i well-known i may fail for highly correlaed ae (ee for example Baeva (0)). he reul in he above ecion give an analyical reaon for hi phenomenon.

22 In fac, noice ha p a (he volailiy parameer in he Kirk formula) i a proce ha doe no depend on X nor on he ime o mauriry : hen, Kirk formula may no reproduce he hor-ime volailiy kew ha we have een appear in he highly correlaed cae ( cloe o ) and we can expec i can fail when i near one. In he following example, we compare Kirk approximaion price wih he one obained wih a Mone-Carlo imulaion procedure wih 00,000 rial.we ue hee imulaion reul a he benchmark for he rue pread opion value Example 4 In he following able we can compare he price given by Kirk formula and by he Mone-Carlo imulaion, for di eren value for K and : Here X = ln(00); S 0 = 00; r = 0; = 0:5; = 0:5 and 0 = 0:4. he di erence i given in % of he Mone-Carlo price. A expeced from our analyical udy, he error increae rongly when he correlaion i cloe o. K= 0:60 0:98 0:99 0:999 5 Mone-Carlo Kirk Error 9; ; ; 40 % ; 890 ; 59 ; 30 % ; 8386 ; 8775 ; 7 % 0 Mone-Carlo Kirk Error 7; ; ; 545 % ; 74 ; ; 84 % ; 007 ; 05 7 ; 93 % ; 50 ; 540 ; 75 % 0; ; 8848 ; 56 % Noice ha our udy of he volailiy kew give no only a heoreical underanding of i goodne-of- of Kirk approximaion, bu alo give u ome hin o improve i. In fac, from (5) and by uing aylor expanion we can expec ha for mall ime o mauriie and for near-he-money opion, he expreion q ^I (X ) := a S 0 0 S 0 K p 3 ( 0 ) SK 0 a (S 0 K) (X x ) can be a reaonable approximaion for he implied volailiy I (X ). In fac, le u conider he modi ed Kirk approximaion given by BS(; X ; M ; ^I (X )): In he following example we will check numerically he goodne-of- of hi approximaion. Example 5 In he following able we can compare he price given by he modi ed Kirk formula and by he Mone-Carlo imulaion, for di eren value for K and and for he ame parameer a in Example 4. Noice ha we have reduced igni caively he error of approximaion wih repec o he reul in hi la example, pecially in he cae of highly correlaed ae.

23 K= 0:60 0:98 0:99 0:999 5 Mone-Carlo Modi ed Kirk Error 9; ; ; 37 % ; 890 ; ; 809 % ; 8386 ; ; 804 % 0 Mone-Carlo Modi ed Kirk Error 7 Concluion 7; ; ; 45 % ; 74 ; 888 ; 368 % ; 007 ; 0367 ; 660 % ; 50 ; ; 44 % 0; ; 80 ; 400 % We have ued he Malliavin calculu echnique o nd an expreion for he hor-ime implied volailiy kew of opion wih random rike. In paricular, we have een ha hi analyical udy give u a key o underand why ome approximaion formula may fail for ome e of parameer. A an applicaion, we have een ha hi kew i very pronounced for pread opion in he high correlaion cae, wich explain why a conan volailiy approximaion a Kirk formula canno be accurae in hi cae. Finally, our approach give u ome hin o improve hi eimae, by inroducing a kew in he approximaion of he correponding implied volailiy. Our preliminar numerical analyi how ha hi can be a naural and e cien way o improve Kirk formula. A precie developemen of hi improvemen i lef for fuure reearch. Reference [] C. Alexander and A. Venkaramanan: Cloed form approximaion for pread opion. Applied Mahemaical Finance 8 (5), (0) [] E. Alò: A generalizaion of he Hull and Whie formula wih applicaion o opion pricing approximaion. Finance Soch. 0, (006). [3] E. Alò, J.A. León and J. Vive: On he hor-ime behavior of he implied volailiy for jump-di uion model wih ochaic volailiy. Finance Soch., (007). [4] E. Alò and D. Nualar: An exenion of Iô formula for anicipaing procee. J. heor. Probab., (998) [5]. Baeva: On he pricing and eniiviy of pread opion on wo correlaed ae. Preprin (0). [6] P. Bjerkund and G. Senland: Cloed form pread opion valuaion. Quaniaive Finance ifir, -0 (0). [7] S. Borovkova, F. J. Permana and H. v.d.weide: A cloed form approach o he valuaion and hedging of bake and pread opion. he Journal of Derivaive 4 (4), 8-4 (007). 3

24 [8] R. Carmona andv. Durrleman: Pricing and hedging pread opion. SIAM Rev. 45 (4), (003). [9] S. Deng, M. Li and J. hou: Cloed-form approximaion for pread opion price and greek. Preprin (008) [0] J.-P. Fouque, G. Papanicolaou, R. Sircar and K. Solna: Mauriy cycle in implied volailiy. Finance Soch. 8 (4), (004). [] J. Jacod and P. Proer: Rik-neural compaibiliy wih opion price. Finance Soch. 4 (), (00). [] E. Kirk: Correlaion in he energy marke. In Managing Energy Price Rik (Fir Ediion). London: Rik Publicaion and Enron, pp (995) [3] W. Margrabe: he value of an opion o exchange one ae for anoher. he Journal of Finance 33 (), (978). [4] A. Medvedev and O. Scaille: Approximaion and calibraion of hor-erm implied volailiie under jump-di uion ochaic volailiy. Rev. Finance Sud. 0 (), (007). [5] D. Nualar: he Malliavin Calculu and Relaed opic. Second Ediion. Springer-Verlag, Berlin (006). [6] E. Renaul and N. ouzi: Opion hedging and implied volailiie in a ochaic volailiy model. Mah. Finance 6 (3) (996). [7] D. Shimko: Opion on fuure pread: hedging, peculaion and valuaion. Journal of Fuure Marke 4 (), 83 3 (994). 4

The Implied Volatility of Forward Starting Options: ATM Short-Time Level, Skew and Curvature

The Implied Volatility of Forward Starting Options: ATM Short-Time Level, Skew and Curvature The Implied Volailiy of Forward Saring Opion: ATM Shor-Time Level, Skew and Crvare Elia Alò Anoine Jacqier Jorge A León May 017 Barcelona GSE Working Paper Serie Working Paper nº 988 The implied volailiy

More information

ESTIMATES FOR THE DERIVATIVE OF DIFFUSION SEMIGROUPS

ESTIMATES FOR THE DERIVATIVE OF DIFFUSION SEMIGROUPS Elec. Comm. in Probab. 3 (998) 65 74 ELECTRONIC COMMUNICATIONS in PROBABILITY ESTIMATES FOR THE DERIVATIVE OF DIFFUSION SEMIGROUPS L.A. RINCON Deparmen of Mahemaic Univeriy of Wale Swanea Singleon Par

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 31 Signal & Syem Prof. Mark Fowler Noe Se #27 C-T Syem: Laplace Tranform Power Tool for yem analyi Reading Aignmen: Secion 6.1 6.3 of Kamen and Heck 1/18 Coure Flow Diagram The arrow here how concepual

More information

Fractional Ornstein-Uhlenbeck Bridge

Fractional Ornstein-Uhlenbeck Bridge WDS'1 Proceeding of Conribued Paper, Par I, 21 26, 21. ISBN 978-8-7378-139-2 MATFYZPRESS Fracional Ornein-Uhlenbeck Bridge J. Janák Charle Univeriy, Faculy of Mahemaic and Phyic, Prague, Czech Republic.

More information

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship Laplace Tranform (Lin & DeCarlo: Ch 3) ENSC30 Elecric Circui II The Laplace ranform i an inegral ranformaion. I ranform: f ( ) F( ) ime variable complex variable From Euler > Lagrange > Laplace. Hence,

More information

6.8 Laplace Transform: General Formulas

6.8 Laplace Transform: General Formulas 48 HAP. 6 Laplace Tranform 6.8 Laplace Tranform: General Formula Formula Name, ommen Sec. F() l{ f ()} e f () d f () l {F()} Definiion of Tranform Invere Tranform 6. l{af () bg()} al{f ()} bl{g()} Lineariy

More information

Introduction to SLE Lecture Notes

Introduction to SLE Lecture Notes Inroducion o SLE Lecure Noe May 13, 16 - The goal of hi ecion i o find a ufficien condiion of λ for he hull K o be generaed by a imple cure. I urn ou if λ 1 < 4 hen K i generaed by a imple curve. We will

More information

On the Exponential Operator Functions on Time Scales

On the Exponential Operator Functions on Time Scales dvance in Dynamical Syem pplicaion ISSN 973-5321, Volume 7, Number 1, pp. 57 8 (212) hp://campu.m.edu/ada On he Exponenial Operaor Funcion on Time Scale laa E. Hamza Cairo Univeriy Deparmen of Mahemaic

More information

Randomized Perfect Bipartite Matching

Randomized Perfect Bipartite Matching Inenive Algorihm Lecure 24 Randomized Perfec Biparie Maching Lecurer: Daniel A. Spielman April 9, 208 24. Inroducion We explain a randomized algorihm by Ahih Goel, Michael Kapralov and Sanjeev Khanna for

More information

Explicit form of global solution to stochastic logistic differential equation and related topics

Explicit form of global solution to stochastic logistic differential equation and related topics SAISICS, OPIMIZAION AND INFOMAION COMPUING Sa., Opim. Inf. Compu., Vol. 5, March 17, pp 58 64. Publihed online in Inernaional Academic Pre (www.iapre.org) Explici form of global oluion o ochaic logiic

More information

Generalized Orlicz Spaces and Wasserstein Distances for Convex-Concave Scale Functions

Generalized Orlicz Spaces and Wasserstein Distances for Convex-Concave Scale Functions Generalized Orlicz Space and Waerein Diance for Convex-Concave Scale Funcion Karl-Theodor Surm Abrac Given a ricly increaing, coninuou funcion ϑ : R + R +, baed on he co funcional ϑ (d(x, y dq(x, y, we

More information

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER #A30 INTEGERS 10 (010), 357-363 FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER Nahan Kaplan Deparmen of Mahemaic, Harvard Univeriy, Cambridge, MA nkaplan@mah.harvard.edu Received: 7/15/09, Revied:

More information

BSDE's, Clark-Ocone formula, and Feynman-Kac formula for Lévy processes Nualart, D.; Schoutens, W.

BSDE's, Clark-Ocone formula, and Feynman-Kac formula for Lévy processes Nualart, D.; Schoutens, W. BSD', Clark-Ocone formula, and Feynman-Kac formula for Lévy procee Nualar, D.; Schouen, W. Publihed: 1/1/ Documen Verion Publiher PDF, alo known a Verion of ecord (include final page, iue and volume number)

More information

Stability in Distribution for Backward Uncertain Differential Equation

Stability in Distribution for Backward Uncertain Differential Equation Sabiliy in Diribuion for Backward Uncerain Differenial Equaion Yuhong Sheng 1, Dan A. Ralecu 2 1. College of Mahemaical and Syem Science, Xinjiang Univeriy, Urumqi 8346, China heng-yh12@mail.inghua.edu.cn

More information

CHAPTER 7: SECOND-ORDER CIRCUITS

CHAPTER 7: SECOND-ORDER CIRCUITS EEE5: CI RCUI T THEORY CHAPTER 7: SECOND-ORDER CIRCUITS 7. Inroducion Thi chaper conider circui wih wo orage elemen. Known a econd-order circui becaue heir repone are decribed by differenial equaion ha

More information

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max ecure 8 7. Sabiliy Analyi For an n dimenional vecor R n, he and he vecor norm are defined a: = T = i n i (7.) I i eay o how ha hee wo norm aify he following relaion: n (7.) If a vecor i ime-dependen, hen

More information

Mathematische Annalen

Mathematische Annalen Mah. Ann. 39, 33 339 (997) Mahemaiche Annalen c Springer-Verlag 997 Inegraion by par in loop pace Elon P. Hu Deparmen of Mahemaic, Norhweern Univeriy, Evanon, IL 628, USA (e-mail: elon@@mah.nwu.edu) Received:

More information

EE202 Circuit Theory II

EE202 Circuit Theory II EE202 Circui Theory II 2017-2018, Spring Dr. Yılmaz KALKAN I. Inroducion & eview of Fir Order Circui (Chaper 7 of Nilon - 3 Hr. Inroducion, C and L Circui, Naural and Sep epone of Serie and Parallel L/C

More information

Rough Paths and its Applications in Machine Learning

Rough Paths and its Applications in Machine Learning Pah ignaure Machine learning applicaion Rough Pah and i Applicaion in Machine Learning July 20, 2017 Rough Pah and i Applicaion in Machine Learning Pah ignaure Machine learning applicaion Hiory and moivaion

More information

FIXED POINTS AND STABILITY IN NEUTRAL DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS

FIXED POINTS AND STABILITY IN NEUTRAL DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 136, Number 3, March 28, Page 99 918 S 2-9939(7)989-2 Aricle elecronically publihed on November 3, 27 FIXED POINTS AND STABILITY IN NEUTRAL DIFFERENTIAL

More information

Macroeconomics 1. Ali Shourideh. Final Exam

Macroeconomics 1. Ali Shourideh. Final Exam 4780 - Macroeconomic 1 Ali Shourideh Final Exam Problem 1. A Model of On-he-Job Search Conider he following verion of he McCall earch model ha allow for on-he-job-earch. In paricular, uppoe ha ime i coninuou

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

ANSWERS TO ODD NUMBERED EXERCISES IN CHAPTER

ANSWERS TO ODD NUMBERED EXERCISES IN CHAPTER John Riley 6 December 200 NWER TO ODD NUMBERED EXERCIE IN CHPTER 7 ecion 7 Exercie 7-: m m uppoe ˆ, m=,, M (a For M = 2, i i eay o how ha I implie I From I, for any probabiliy vecor ( p, p 2, 2 2 ˆ ( p,

More information

Fractional Brownian motion and applications Part I: fractional Brownian motion in Finance

Fractional Brownian motion and applications Part I: fractional Brownian motion in Finance Fracional Brownian moion and applicaion Par I: fracional Brownian moion in Finance INTRODUCTION The fbm i an exenion of he claical Brownian moion ha allow i dijoin incremen o be correlaed. Moivaed by empirical

More information

U T,0. t = X t t T X T. (1)

U T,0. t = X t t T X T. (1) Gauian bridge Dario Gabarra 1, ommi Soinen 2, and Eko Valkeila 3 1 Deparmen of Mahemaic and Saiic, PO Box 68, 14 Univeriy of Helinki,Finland dariogabarra@rnihelinkifi 2 Deparmen of Mahemaic and Saiic,

More information

On Delayed Logistic Equation Driven by Fractional Brownian Motion

On Delayed Logistic Equation Driven by Fractional Brownian Motion On Delayed Logiic Equaion Driven by Fracional Brownian Moion Nguyen Tien Dung Deparmen of Mahemaic, FPT Univeriy No 8 Ton Tha Thuye, Cau Giay, Hanoi, 84 Vienam Email: dungn@fp.edu.vn ABSTRACT In hi paper

More information

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility Saisics 441 (Fall 214) November 19, 21, 214 Prof Michael Kozdron Lecure #31, 32: The Ornsein-Uhlenbeck Process as a Model of Volailiy The Ornsein-Uhlenbeck process is a di usion process ha was inroduced

More information

Parameter Estimation for Fractional Ornstein-Uhlenbeck Processes: Non-Ergodic Case

Parameter Estimation for Fractional Ornstein-Uhlenbeck Processes: Non-Ergodic Case Parameer Eimaion for Fracional Ornein-Uhlenbeck Procee: Non-Ergodic Cae R. Belfadli 1, K. E-Sebaiy and Y. Ouknine 3 1 Polydiciplinary Faculy of Taroudan, Univeriy Ibn Zohr, Taroudan, Morocco. Iniu de Mahémaique

More information

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2) Laplace Tranform Maoud Malek The Laplace ranform i an inegral ranform named in honor of mahemaician and aronomer Pierre-Simon Laplace, who ued he ranform in hi work on probabiliy heory. I i a powerful

More information

Option pricing and implied volatilities in a 2-hypergeometric stochastic volatility model

Option pricing and implied volatilities in a 2-hypergeometric stochastic volatility model Opion pricing and implied volailiies in a 2-hypergeomeric sochasic volailiy model Nicolas Privaul Qihao She Division of Mahemaical Sciences School of Physical and Mahemaical Sciences Nanyang Technological

More information

Algorithmic Discrete Mathematics 6. Exercise Sheet

Algorithmic Discrete Mathematics 6. Exercise Sheet Algorihmic Dicree Mahemaic. Exercie Shee Deparmen of Mahemaic SS 0 PD Dr. Ulf Lorenz 7. and 8. Juni 0 Dipl.-Mah. David Meffer Verion of June, 0 Groupwork Exercie G (Heap-Sor) Ue Heap-Sor wih a min-heap

More information

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Discussion Session 2 Constant Acceleration/Relative Motion Week 03 PHYS 100 Dicuion Seion Conan Acceleraion/Relaive Moion Week 03 The Plan Today you will work wih your group explore he idea of reference frame (i.e. relaive moion) and moion wih conan acceleraion. You ll

More information

T-Rough Fuzzy Subgroups of Groups

T-Rough Fuzzy Subgroups of Groups Journal of mahemaic and compuer cience 12 (2014), 186-195 T-Rough Fuzzy Subgroup of Group Ehagh Hoeinpour Deparmen of Mahemaic, Sari Branch, Ilamic Azad Univeriy, Sari, Iran. hoienpor_a51@yahoo.com Aricle

More information

Chapter 7: Inverse-Response Systems

Chapter 7: Inverse-Response Systems Chaper 7: Invere-Repone Syem Normal Syem Invere-Repone Syem Baic Sar ou in he wrong direcion End up in he original eady-ae gain value Two or more yem wih differen magniude and cale in parallel Main yem

More information

Approximate Controllability of Fractional Stochastic Perturbed Control Systems Driven by Mixed Fractional Brownian Motion

Approximate Controllability of Fractional Stochastic Perturbed Control Systems Driven by Mixed Fractional Brownian Motion American Journal of Applied Mahemaic and Saiic, 15, Vol. 3, o. 4, 168-176 Available online a hp://pub.ciepub.com/ajam/3/4/7 Science and Educaion Publihing DOI:1.1691/ajam-3-4-7 Approximae Conrollabiliy

More information

On the Benney Lin and Kawahara Equations

On the Benney Lin and Kawahara Equations JOURNAL OF MATEMATICAL ANALYSIS AND APPLICATIONS 11, 13115 1997 ARTICLE NO AY975438 On he BenneyLin and Kawahara Equaion A Biagioni* Deparmen of Mahemaic, UNICAMP, 1381-97, Campina, Brazil and F Linare

More information

Research Article Existence and Uniqueness of Solutions for a Class of Nonlinear Stochastic Differential Equations

Research Article Existence and Uniqueness of Solutions for a Class of Nonlinear Stochastic Differential Equations Hindawi Publihing Corporaion Abrac and Applied Analyi Volume 03, Aricle ID 56809, 7 page hp://dx.doi.org/0.55/03/56809 Reearch Aricle Exience and Uniquene of Soluion for a Cla of Nonlinear Sochaic Differenial

More information

Analytical Pricing of An Insurance Embedded Option: Alternative Formulas and Gaussian Approximation

Analytical Pricing of An Insurance Embedded Option: Alternative Formulas and Gaussian Approximation Journal of Informaic and Mahemaical Science Volume 3 (0), Number, pp. 87 05 RGN Publicaion hp://www.rgnpublicaion.com Analyical Pricing of An Inurance Embedded Opion: Alernaive Formula and Gauian Approximaion

More information

Time Varying Multiserver Queues. W. A. Massey. Murray Hill, NJ Abstract

Time Varying Multiserver Queues. W. A. Massey. Murray Hill, NJ Abstract Waiing Time Aympoic for Time Varying Mulierver ueue wih Abonmen Rerial A. Melbaum Technion Iniue Haifa, 3 ISRAEL avim@x.echnion.ac.il M. I. Reiman Bell Lab, Lucen Technologie Murray Hill, NJ 7974 U.S.A.

More information

ON FRACTIONAL ORNSTEIN-UHLENBECK PROCESSES

ON FRACTIONAL ORNSTEIN-UHLENBECK PROCESSES Communicaion on Sochaic Analyi Vol. 5, No. 1 211 121-133 Serial Publicaion www.erialpublicaion.com ON FRACTIONAL ORNSTEIN-UHLENBECK PROCESSES TERHI KAARAKKA AND PAAVO SALMINEN Abrac. In hi paper we udy

More information

CHAPTER 7. Definition and Properties. of Laplace Transforms

CHAPTER 7. Definition and Properties. of Laplace Transforms SERIES OF CLSS NOTES FOR 5-6 TO INTRODUCE LINER ND NONLINER PROBLEMS TO ENGINEERS, SCIENTISTS, ND PPLIED MTHEMTICINS DE CLSS NOTES COLLECTION OF HNDOUTS ON SCLR LINER ORDINRY DIFFERENTIL EQUTIONS (ODE")

More information

Chapter 6. Laplace Transforms

Chapter 6. Laplace Transforms 6- Chaper 6. Laplace Tranform 6.4 Shor Impule. Dirac Dela Funcion. Parial Fracion 6.5 Convoluion. Inegral Equaion 6.6 Differeniaion and Inegraion of Tranform 6.7 Syem of ODE 6.4 Shor Impule. Dirac Dela

More information

Backward Stochastic Differential Equations and Applications in Finance

Backward Stochastic Differential Equations and Applications in Finance Backward Sochaic Differenial Equaion and Applicaion in Finance Ying Hu Augu 1, 213 1 Inroducion The aim of hi hor cae i o preen he baic heory of BSDE and o give ome applicaion in 2 differen domain: mahemaical

More information

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001 CS 545 Flow Nework lon Efra Slide courey of Charle Leieron wih mall change by Carola Wenk Flow nework Definiion. flow nework i a direced graph G = (V, E) wih wo diinguihed verice: a ource and a ink. Each

More information

CHAPTER. Forced Equations and Systems { } ( ) ( ) 8.1 The Laplace Transform and Its Inverse. Transforms from the Definition.

CHAPTER. Forced Equations and Systems { } ( ) ( ) 8.1 The Laplace Transform and Its Inverse. Transforms from the Definition. CHAPTER 8 Forced Equaion and Syem 8 The aplace Tranform and I Invere Tranform from he Definiion 5 5 = b b {} 5 = 5e d = lim5 e = ( ) b {} = e d = lim e + e d b = (inegraion by par) = = = = b b ( ) ( )

More information

Notes on cointegration of real interest rates and real exchange rates. ρ (2)

Notes on cointegration of real interest rates and real exchange rates. ρ (2) Noe on coinegraion of real inere rae and real exchange rae Charle ngel, Univeriy of Wiconin Le me ar wih he obervaion ha while he lieraure (mo prominenly Meee and Rogoff (988) and dion and Paul (993))

More information

Chapter 6. Laplace Transforms

Chapter 6. Laplace Transforms Chaper 6. Laplace Tranform Kreyzig by YHLee;45; 6- An ODE i reduced o an algebraic problem by operaional calculu. The equaion i olved by algebraic manipulaion. The reul i ranformed back for he oluion of

More information

Measure-valued Diffusions and Stochastic Equations with Poisson Process 1

Measure-valued Diffusions and Stochastic Equations with Poisson Process 1 Publihed in: Oaka Journal of Mahemaic 41 (24), 3: 727 744 Meaure-valued Diffuion and Sochaic quaion wih Poion Proce 1 Zongfei FU and Zenghu LI 2 Running head: Meaure-valued Diffuion and Sochaic quaion

More information

Research Article Stochastic Analysis of Gaussian Processes via Fredholm Representation

Research Article Stochastic Analysis of Gaussian Processes via Fredholm Representation Inernaional Journal of Sochaic Analyi Volume 216, Aricle ID 8694365, 15 page hp://dx.doi.org/1.1155/216/8694365 Reearch Aricle Sochaic Analyi of Gauian Procee via Fredholm Repreenaion ommi Soinen 1 and

More information

Lecture 26. Lucas and Stokey: Optimal Monetary and Fiscal Policy in an Economy without Capital (JME 1983) t t

Lecture 26. Lucas and Stokey: Optimal Monetary and Fiscal Policy in an Economy without Capital (JME 1983) t t Lecure 6. Luca and Sokey: Opimal Moneary and Fical Policy in an Economy wihou Capial (JME 983. A argued in Kydland and Preco (JPE 977, Opimal governmen policy i likely o be ime inconien. Fiher (JEDC 98

More information

Approximation for Option Prices under Uncertain Volatility

Approximation for Option Prices under Uncertain Volatility Approximaion for Opion Price under Uncerain Volailiy Jean-Pierre Fouque Bin Ren February, 3 Abrac In hi paper, we udy he aympoic behavior of he wor cae cenario opion price a he volailiy inerval in an uncerain

More information

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY Ian Crawford THE INSTITUTE FOR FISCAL STUDIES DEPARTMENT OF ECONOMICS, UCL cemmap working paper CWP02/04 Neceary and Sufficien Condiion for Laen

More information

GLOBAL ANALYTIC REGULARITY FOR NON-LINEAR SECOND ORDER OPERATORS ON THE TORUS

GLOBAL ANALYTIC REGULARITY FOR NON-LINEAR SECOND ORDER OPERATORS ON THE TORUS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 131, Number 12, Page 3783 3793 S 0002-9939(03)06940-5 Aricle elecronically publihed on February 28, 2003 GLOBAL ANALYTIC REGULARITY FOR NON-LINEAR

More information

18 Extensions of Maximum Flow

18 Extensions of Maximum Flow Who are you?" aid Lunkwill, riing angrily from hi ea. Wha do you wan?" I am Majikhie!" announced he older one. And I demand ha I am Vroomfondel!" houed he younger one. Majikhie urned on Vroomfondel. I

More information

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

More information

and Van Ne [3], where hey inroduced fbm a a cenered Gauian proce and gave he r repreenaion of i a an inegral wih repec o andard BM. he ur parameer i n

and Van Ne [3], where hey inroduced fbm a a cenered Gauian proce and gave he r repreenaion of i a an inegral wih repec o andard BM. he ur parameer i n Porfolio Opimizaion wih conumpion in a fracional Black-Schole marke Frederi G. Vien Λ ao Zhang y Yalοc n Sarol z November 4, 6 Abrac In hi paper we conider he claical Meron problem of nding he opimal conumpion

More information

MALLIAVIN CALCULUS FOR BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS AND APPLICATION TO NUMERICAL SOLUTIONS

MALLIAVIN CALCULUS FOR BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS AND APPLICATION TO NUMERICAL SOLUTIONS The Annal of Applied Probabiliy 211, Vol. 21, No. 6, 2379 2423 DOI: 1.1214/11-AAP762 Iniue of Mahemaical Saiic, 211 MALLIAVIN CALCULUS FOR BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS AND APPLICATION TO

More information

News-generated dependence and optimal portfolios for n stocks in a market of Barndor -Nielsen and Shephard type.

News-generated dependence and optimal portfolios for n stocks in a market of Barndor -Nielsen and Shephard type. New-generaed dependence and opimal porfolio for n ock in a marke of Barndor -Nielen and Shephard ype. Carl Lindberg Deparmen of Mahemaical Saiic Chalmer Univeriy of Technology and Göeborg Univeriy Göeborg,

More information

18.03SC Unit 3 Practice Exam and Solutions

18.03SC Unit 3 Practice Exam and Solutions Sudy Guide on Sep, Dela, Convoluion, Laplace You can hink of he ep funcion u() a any nice mooh funcion which i for < a and for > a, where a i a poiive number which i much maller han any ime cale we care

More information

CSC 364S Notes University of Toronto, Spring, The networks we will consider are directed graphs, where each edge has associated with it

CSC 364S Notes University of Toronto, Spring, The networks we will consider are directed graphs, where each edge has associated with it CSC 36S Noe Univeriy of Torono, Spring, 2003 Flow Algorihm The nework we will conider are direced graph, where each edge ha aociaed wih i a nonnegaive capaciy. The inuiion i ha if edge (u; v) ha capaciy

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull Chaper 14 Wiener Processes and Iô s Lemma Copyrigh John C. Hull 014 1 Sochasic Processes! Describes he way in which a variable such as a sock price, exchange rae or ineres rae changes hrough ime! Incorporaes

More information

Note on Matuzsewska-Orlich indices and Zygmund inequalities

Note on Matuzsewska-Orlich indices and Zygmund inequalities ARMENIAN JOURNAL OF MATHEMATICS Volume 3, Number 1, 21, 22 31 Noe on Mauzewka-Orlic indice and Zygmund inequaliie N. G. Samko Univeridade do Algarve, Campu de Gambela, Faro,85 139, Porugal namko@gmail.com

More information

SINGULAR FORWARD BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS AND EMISSIONS DERIVATIVES

SINGULAR FORWARD BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS AND EMISSIONS DERIVATIVES he Annal of Applied Probabiliy 213, Vol. 23, No. 3, 186 1128 DOI: 1.1214/12-AAP865 Iniue of Mahemaical Saiic, 213 SINGULAR FORWARD BACKWARD SOCHASIC DIFFERENIAL EQUAIONS AND EMISSIONS DERIVAIVES BY RENÉ

More information

Research Article An Upper Bound on the Critical Value β Involved in the Blasius Problem

Research Article An Upper Bound on the Critical Value β Involved in the Blasius Problem Hindawi Publihing Corporaion Journal of Inequaliie and Applicaion Volume 2010, Aricle ID 960365, 6 page doi:10.1155/2010/960365 Reearch Aricle An Upper Bound on he Criical Value Involved in he Blaiu Problem

More information

Research Article On Double Summability of Double Conjugate Fourier Series

Research Article On Double Summability of Double Conjugate Fourier Series Inernaional Journal of Mahemaic and Mahemaical Science Volume 22, Aricle ID 4592, 5 page doi:.55/22/4592 Reearch Aricle On Double Summabiliy of Double Conjugae Fourier Serie H. K. Nigam and Kuum Sharma

More information

EE Control Systems LECTURE 2

EE Control Systems LECTURE 2 Copyrigh F.L. Lewi 999 All righ reerved EE 434 - Conrol Syem LECTURE REVIEW OF LAPLACE TRANSFORM LAPLACE TRANSFORM The Laplace ranform i very ueful in analyi and deign for yem ha are linear and ime-invarian

More information

1 Motivation and Basic Definitions

1 Motivation and Basic Definitions CSCE : Deign and Analyi of Algorihm Noe on Max Flow Fall 20 (Baed on he preenaion in Chaper 26 of Inroducion o Algorihm, 3rd Ed. by Cormen, Leieron, Rive and Sein.) Moivaion and Baic Definiion Conider

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

A proof of Ito's formula using a di Title formula. Author(s) Fujita, Takahiko; Kawanishi, Yasuhi. Studia scientiarum mathematicarum H Citation

A proof of Ito's formula using a di Title formula. Author(s) Fujita, Takahiko; Kawanishi, Yasuhi. Studia scientiarum mathematicarum H Citation A proof of Io's formula using a di Tile formula Auhor(s) Fujia, Takahiko; Kawanishi, Yasuhi Sudia scieniarum mahemaicarum H Ciaion 15-134 Issue 8-3 Dae Type Journal Aricle Tex Version auhor URL hp://hdl.handle.ne/186/15878

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

Pricing the American Option Using Itô s Formula and Optimal Stopping Theory

Pricing the American Option Using Itô s Formula and Optimal Stopping Theory U.U.D.M. Projec Repor 2014:3 Pricing he American Opion Uing Iô Formula and Opimal Sopping Theory Jona Bergröm Examenarbee i maemaik, 15 hp Handledare och examinaor: Erik Ekröm Januari 2014 Deparmen of

More information

6.302 Feedback Systems Recitation : Phase-locked Loops Prof. Joel L. Dawson

6.302 Feedback Systems Recitation : Phase-locked Loops Prof. Joel L. Dawson 6.32 Feedback Syem Phae-locked loop are a foundaional building block for analog circui deign, paricularly for communicaion circui. They provide a good example yem for hi cla becaue hey are an excellen

More information

On a Fractional Stochastic Landau-Ginzburg Equation

On a Fractional Stochastic Landau-Ginzburg Equation Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

Energy Equality and Uniqueness of Weak Solutions to MHD Equations in L (0,T; L n (Ω))

Energy Equality and Uniqueness of Weak Solutions to MHD Equations in L (0,T; L n (Ω)) Aca Mahemaica Sinica, Englih Serie May, 29, Vol. 25, No. 5, pp. 83 814 Publihed online: April 25, 29 DOI: 1.17/1114-9-7214-8 Hp://www.AcaMah.com Aca Mahemaica Sinica, Englih Serie The Ediorial Office of

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

Local Risk-Minimization for Defaultable Claims with Recovery Process

Local Risk-Minimization for Defaultable Claims with Recovery Process Local Rik-Minimizaion for Defaulable Claim wih Recovery Proce Franceca Biagini Aleandra Crearola Abrac We udy he local rik-minimizaion approach for defaulable claim wih random recovery a defaul ime, een

More information

Pathwise description of dynamic pitchfork bifurcations with additive noise

Pathwise description of dynamic pitchfork bifurcations with additive noise Pahwie decripion of dynamic pichfork bifurcaion wih addiive noie Nil Berglund and Barbara Genz Abrac The low drif (wih peed ) of a parameer hrough a pichfork bifurcaion poin, known a he dynamic pichfork

More information

An introduction to the (local) martingale problem

An introduction to the (local) martingale problem An inroducion o he (local) maringale problem Chri Janjigian Ocober 14, 214 Abrac Thee are my preenaion noe for a alk in he Univeriy of Wiconin - Madion graduae probabiliy eminar. Thee noe are primarily

More information

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance. 1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard

More information

Clark s Representation of Wiener Functionals and Hedging of the Barrier Option

Clark s Representation of Wiener Functionals and Hedging of the Barrier Option aqarvelo mecnierebaa erovnuli akademii moambe, 8, #, 04 ULLEIN OF HE GEORGIAN NAIONAL ACADEMY OF SCIENCES, vol 8, no, 04 Mahemaic Clark Repreenaion of Wiener Funcional and Hedging of he arrier Opion Omar

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0.

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0. Advanced Financial Models Example shee 3 - Michaelmas 217 Michael Tehranchi Problem 1. Le f : [, R be a coninuous (non-random funcion and W a Brownian moion, and le σ 2 = f(s 2 ds and assume σ 2

More information

CONTROL SYSTEMS. Chapter 10 : State Space Response

CONTROL SYSTEMS. Chapter 10 : State Space Response CONTROL SYSTEMS Chaper : Sae Space Repone GATE Objecive & Numerical Type Soluion Queion 5 [GATE EE 99 IIT-Bombay : Mark] Conider a econd order yem whoe ae pace repreenaion i of he form A Bu. If () (),

More information

FUZZY n-inner PRODUCT SPACE

FUZZY n-inner PRODUCT SPACE Bull. Korean Mah. Soc. 43 (2007), No. 3, pp. 447 459 FUZZY n-inner PRODUCT SPACE Srinivaan Vijayabalaji and Naean Thillaigovindan Reprined from he Bullein of he Korean Mahemaical Sociey Vol. 43, No. 3,

More information

1 CHAPTER 14 LAPLACE TRANSFORMS

1 CHAPTER 14 LAPLACE TRANSFORMS CHAPTER 4 LAPLACE TRANSFORMS 4 nroducion f x) i a funcion of x, where x lie in he range o, hen he funcion p), defined by p) px e x) dx, 4 i called he Laplace ranform of x) However, in hi chaper, where

More information

Multidimensional Markovian FBSDEs with superquadratic

Multidimensional Markovian FBSDEs with superquadratic Mulidimenional Markovian FBSDE wih uperquadraic growh Michael Kupper a,1, Peng Luo b,, Ludovic angpi c,3 December 14, 17 ABSRAC We give local and global exience and uniquene reul for mulidimenional coupled

More information

Backward stochastic differential equations and Feynman±Kac formula for LeÂvy processes, with applications in nance

Backward stochastic differential equations and Feynman±Kac formula for LeÂvy processes, with applications in nance Bernoulli 7(5), 21, 761±776 Backward ochaic differenial equaion and Feynman±Kac formula for LeÂvy procee, wih applicaion in nance DAVID NUALAT 1 and WIM SCHOUTENS 2 1 Univeria de Barcelona, Gran Via de

More information

Heavy Tails of Discounted Aggregate Claims in the Continuous-time Renewal Model

Heavy Tails of Discounted Aggregate Claims in the Continuous-time Renewal Model Heavy Tails of Discouned Aggregae Claims in he Coninuous-ime Renewal Model Qihe Tang Deparmen of Saisics and Acuarial Science The Universiy of Iowa 24 Schae er Hall, Iowa Ciy, IA 52242, USA E-mail: qang@sa.uiowa.edu

More information

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS

A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS Theory of Sochasic Processes Vol. 14 3), no. 2, 28, pp. 139 144 UDC 519.21 JOSEP LLUÍS SOLÉ AND FREDERIC UTZET A FAMILY OF MARTINGALES GENERATED BY A PROCESS WITH INDEPENDENT INCREMENTS An explici procedure

More information

Lecture 4: Processes with independent increments

Lecture 4: Processes with independent increments Lecure 4: Processes wih independen incremens 1. A Wienner process 1.1 Definiion of a Wienner process 1.2 Reflecion principle 1.3 Exponenial Brownian moion 1.4 Exchange of measure (Girsanov heorem) 1.5

More information

Use of variance estimation in the multi-armed bandit problem

Use of variance estimation in the multi-armed bandit problem Ue of variance eimaion in he muli-armed bi problem Jean Yve Audiber CERTIS - Ecole de Pon 19, rue Alfred Nobel - Cié Decare 77455 Marne-la-Vallée - France audiber@cerienpcfr Rémi Muno INRIA Fuur, Grappa

More information

Heat semigroup and singular PDEs

Heat semigroup and singular PDEs Hea emigroup and ingular PDE I. BAILLEUL 1 and F. BERNICOT wih an Appendix by F. Bernico & D. Frey Abrac. We provide in hi work a emigroup approach o he udy of ingular PDE, in he line of he paraconrolled

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

Equivalence of K-andJ -methods for limiting real interpolation spaces

Equivalence of K-andJ -methods for limiting real interpolation spaces Available online a wwwciencedireccom Journal of Funcional Analyi 261 2011) 3696 3722 wwweleviercom/locae/jfa Equivalence of K-andJ -meho for limiing real inerpolaion pace Fernando Cobo a,,1, Thoma Kühn

More information

Fractional Brownian Bridge Measures and Their Integration by Parts Formula

Fractional Brownian Bridge Measures and Their Integration by Parts Formula Journal of Mahemaical Reearch wih Applicaion Jul., 218, Vol. 38, No. 4, pp. 418 426 DOI:1.377/j.in:295-2651.218.4.9 Hp://jmre.lu.eu.cn Fracional Brownian Brige Meaure an Their Inegraion by Par Formula

More information

Loss of martingality in asset price models with lognormal stochastic volatility

Loss of martingality in asset price models with lognormal stochastic volatility Loss of maringaliy in asse price models wih lognormal sochasic volailiy BJourdain July 7, 4 Absrac In his noe, we prove ha in asse price models wih lognormal sochasic volailiy, when he correlaion coefficien

More information

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5 Mah 225-4 Week 2 April 2-6 coninue.-.3; alo cover par of.4-.5, EP 7.6 Mon Apr 2:.-.3 Laplace ranform and iniial value problem like we udied in Chaper 5 Announcemen: Warm-up Exercie: Recall, The Laplace

More information