Minhyun Yoo, Darae Jeong, Seungsuk Seo, and Junseok Kim

Size: px
Start display at page:

Download "Minhyun Yoo, Darae Jeong, Seungsuk Seo, and Junseok Kim"

Transcription

1 Hoam Mathematical J. 37 (15), No. 4, pp A COMPARISON STUDY OF EXPLICIT AND IMPLICIT NUMERICAL METHODS FOR THE EQUITY-LINKED SECURITIES Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim Abstract. I this paper, we perform a compariso study of explicit ad implicit umerical methods for the equity-liked securities (ELS). The optio prices of the two-asset ELS are typically computed usig a implicit fiite differece method because a explicit fiite differece scheme has a restrictio for time steps. Nowadays, the three-asset ELS is gettig popularity i the real world fiacial market. I practical applicatios of the fiite differece methods i computatioal fiace, we typically use relatively large space steps ad small time steps. Therefore, we ca use a accurate ad efficiet explicit fiite differece method because the implemetatio is simple ad the computatio is fast. The computatioal results demostrate that if we use a large space step, the the explicit scheme is better tha the implicit oe. O the other had, if the space step size is small, the the implicit scheme is more efficiet tha the explicit oe. 1. Itroductio Equity-liked securities (ELS) are auto-callable optios whose retur o ivestmet is depedet upo the path of the uderlyig equities liked to the securities. ELS ca be made from a few umber of stocks or stock idexes such as the KOSPI i Korea. ELS is a derivative product i the market. ELS guaratees a debt ad is similar to a barrier optio. ELS comprises a large portio of exchage volume i Korea fiacial market. A distiguishig feature of ELS is the automatic earlyredemptio coditio before its maturity. Geerally, i order to get price Received April 15, 15. Accepted October 6, Mathematics Subject Classificatio. 91G6, 65M6. Key words ad phrases. Black Scholes partial differetial equatio, log trasformatio, explicit fiite differece method, equity-liked securities, o-uiform grid. *Correspodig author. Tel.: ; fax:

2 44 Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim of ELS, Mote Carlo simulatio (MCS) ad fiite differece method (FDM) are used. Typically, we use the implicit scheme with operator split method (OSM) or alteratig directio implicit (ADI) because they are stable. Although the implicit scheme with OSM has a advatage i stability, it is costly to solve the tridiagoal matrix implicitly. I particular, we eed to solve the system three times i each time step for the three-asset problems. I practice, we calculate oe time step with oe day, which is about 1/365 ad discretize the asset by oe uit. Therefore, the time step is small eough ad space step is large eough. Furthermore, if we wat to calculate the Greeks of the optio price, especially theta which is the rate of chage of the optio value with respect to chages i the time to maturity, the we have to use much smaller time step. For these cosideratios, it is better to use a explicit scheme with smaller time step, which gives much accurate solutios. The program implemetatio is simple ad fast. The mai purpose of this paper is to develop a explicit scheme o a o-uiform grid to solve value of ELS. The paper is orgaized as follows. We itroduce Black Scholes model i Sectio. I Sectio 3, we preset three-asset step-dow ELS. Numerical methods are preseted i Sectio 4. I Sectio 5, we perform umerical experimets. I Sectio 6, we take coclusio for this paper.. Black Scholes model To evaluate value of the ELS optio, we cosider the stadard Black Scholes model [1], which ca be writte as u(x, t) t = 1 d i,j=1 for (x, t) R d + [, T ]. ρ ij σ i σ j x i x j u(x, t) x i x j r d i=1 x i u(x, t) x i + ru(x, t), Here, u(x, t) is the value of the optio, where x = (x 1, x,..., x d ), d is the total umber of uderlyig assets, x i is the value of the i-th uderlyig assets, ad t is the time. Also, r represets the riskless iterest rate, σ i is the volatility of i-th the uderlyig assets, ρ ij is the correlatio coefficiet betwee x i ad x j, ad T is the maturity time of the optio. Switchig to the ew coordiate X = log x [4] ad usig the

3 A compariso study of umerical method for the equity-liked securities 443 trasformatio τ = T t, the stadard BS equatio ca be rewritte as (1) U(X, τ) τ = d i=1 (r σ i ru(x, τ), ) U(X, τ) X i + 1 d i,j=1 ρ ij σ i σ j U(X, τ) X i X j where U(X, τ) is the value of the optio ad X = (X 1, X,..., X d ). 3. Iitial coditio I this sectio, we briefly describe the cocept of step-dow ELS optio ad we itroduce a three-asset step-dow ELS optio as a example Step-dow ELS The payoff of ELS is determied by early redemptio or fial maturity redemptio. I oe-asset step-dow ELS, if a uderlyig asset is i the predetermied exercise price at certai maturity dates accordig to cotract, the ELS gives desigated retur ad is extermiated. However, if the uderlyig asset is ot i a certai price, the cotract is ot extermiated ad will be cotiued util ext maturity. If the cotract cotiues at fial maturity ad the uderlyig asset is ot i fial exercise barrier, the payoff is determied whether the cotract hit kock i barrier. If the uderlyig asset did ot hit the kock i barrier, the ELS gives predetermied retur, dummy. Otherwise, the ELS will make a loss i face value. The step-dow meas that the desigated strike price decreases. I this paper, we cosider a three-asset step-dow ELS which is similar to a oe asset step-dow ELS, as we explaied above. The differece with a oe asset ELS is that the base price is referred from the miimum of three uderlyig assets at certai maturity dates. The payoff structure of three-asset step-dow ELS is as follows [3]: Early redemptio occurs, ad the cotract is extermiated with predetermied retur if the value of the worst performer, which meas the miimum value of uderlyig assets is greater tha or equal to a give exercise price o a give date. If the early redemptio do ot occur util the fial maturity, the retur depeds upo whether Kock-I occurs or ot.

4 444 Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim 3.. Example for three-asset step-dow ELS To help readers uderstadig of the step-dow ELS optio, we iclude the example for three-asset step-dow ELS. We cosider the parameters as the referece price E = 1, the riskless iterest rate r =.3, the volatilities of the uderlyig assets σ 1 = σ = σ 3 =.3, the correlatios of uderlyig assets ρ 1 = ρ 13 = ρ 3 =.5, the face value F = 1, the Kock-I barrier level KIb =.65E, the dummy rate d =.3, ad the fial maturity time T = 1. At the early redemptio, observatio date (τ i ), exercise price (K i ), ad retur rate (c i ) are described i Table 1. Observatio date (τ i ) Exercise price (K i ) Retur Rate (c i ) τ 6 = T K 6 =.85E c 6 =.3 τ 5 = 5T/6 K 5 =.85E c 5 =.5 τ 4 = 4T/6 K 4 =.9E c 4 =. τ 3 = 3T/6 K 3 =.9E c 3 =.15 τ = T/6 K =.95E c =.1 τ 1 = T/6 K 1 =.95E c 1 =.5 Table 1. Observatio date (τ i ), exercise price (K i ), ad retur rate (c i ) used i example for three-asset step-dow ELS. Figure 1 illustrates the payoff of early obligatory redemptio before fial maturity. Payoff at early redemptio U(X 1,X,X 3,τ i ) V(X 1,X,X 3,τ i ) (1+c 5 )F (1+c 4 )F (1+c 3 )F (1+c )F (1+c 1 )F Observatio date at τ = τ 5 at τ = τ 4 at τ = τ 3 at τ = τ at τ = τ 1.85E.9E.95E S = mi(x 1,X,X 3 ) Figure 1. Payoff at early redemptio before fial maturity.

5 A compariso study of umerical method for the equity-liked securities 445 Alog with the characteristics of the step-dow ELS, we must take two cosideratios of the payoffs at maturity. Accordig to whether or ot the value of ELS hits the Kock-I barrier (KIb) durig the cotract, we defie two values V (X, τ) ad U(X, τ). The, the iitial coditios of U ad V are set to () ad U(X 1, X, X 3, ) = (1 + c 6 )F if S K 6 (1 + d)f if KIb < S < K 6 S F/E otherwise, (3) { (1 + c6 )F if S V (X 1, X, X 3, ) = K 6 S F/E otherwise, where S = mi(x 1, X, X 3 ) which is the value of the worst performer. I Fig., we ca see the correspodig payoffs of U ad V whe c 6 = d. U(X 1,X,X 3,) (1+d)F V(X 1,X,X 3,) (1+c 6 )F Payoff at maturity.65f Payoff at maturity.85f.65f KIb (a).85e S KIb (b).85e S Figure. Payoffs of (a) U ad (b) V at fial maturity. 4. Numerical solutio I this sectio, we describe the umerical discretizatio of Eq. (1) usig explicit scheme o the computatioal domai Ω = [1, S max ] 3, which is log-trasformed by the three-dimesioal fiite domai Ω = [e, e S max ]. Also, to prevet the spurious oscillatory solutio by explicit scheme, we derive the coditio for time step size.

6 446 Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim 4.1. Discretizatio of log-trasformed BS equatio Let X, Y, ad Z deote the log-trasform of x-, y-, ad z-variables, respectively. Now, we discretize the log-trasformed computatioal domai Ω = [1, S max ] 3 with o-uiform spatial step size h i 1 = X i X i 1 = Y i Y i 1 = Z i Z i 1 (see Fig. 3) ad temporal step size τ = T/N τ. Here, X = 1, X Nx = S max, where N x ad N τ are the umbers of grid poits i the X- ad τ-directios, respectively. 1 h i 1 h i X X 1 X i 1 X i X i+1 S max X Nx X Figure 3. Nouiform mesh o log-trasformed grid X. We deote the umerical solutio by U U(X i, Y j, Z k, τ ) for i = 1,,, N x, j = 1,,, N y, k = 1,,, N z, ad =, 1,, N τ. Applyig the explicit fiite differece scheme to Eq. (1) gives U +1 U ( = σ 1 ) ( U τ X + σ U (4) Y ( ) ( ) + r σ 1 U ( ) ( ) + r σ U X Y ( +ρ 1 σ1σ ) ( U + ρ 13 σ X Y 1σ3 U X Z ( +ρ 3 σσ 3 ) U ru Y Z. ) + σ 3 ) + ( r σ 3 ( U Z ) ) ( U Z Here, spatial differeces o the o-uiform grid are defied by ( ) U h i = X h i 1 (h i 1 + h i ) U i 1,jk + h i h i 1 U h i 1 h i ( ) U X = h i 1 + h i (h i 1 + h i ) U i+1,jk, h i 1 (h i 1 + h i ) U i 1,jk h i 1 h i U + h i (h i 1 + h i ) U i+1,jk, )

7 A compariso study of umerical method for the equity-liked securities 447 ( ) U = U i+1,j+1,k + U i 1,j 1,k U i+1,j 1,k U i 1,j 1,k. X Y h i h j + h i h j 1 + h i 1 h j + h i 1 h j 1 Other spatial differeces ca be similarly defied. Equatio (4) is rewritte as follows: [ U +1 = U + τ σ1 h ic 1 (5) h i 1 (h i 1 + h i ) U i 1,jk + σ 1 + (h i h i 1 )C 1 h i 1 h i U + σ 1 + h i 1C 1 h i (h i 1 + h i ) U i+1,jk + σ h jc h j 1 (h j 1 + h j ) U ij 1,k + σ + (h j h j 1 )C h j 1 h j + σ + h j 1C h j (h j 1 + h j ) U ij+1,k + σ 3 h kc 3 h k 1 (h k 1 + h k ) U 1 U + σ 3 + (h k h k 1 )C 3 U h k 1 h + σ 3 + h k 1 C 3 k h k (h k 1 + h k ) U +1 ru ( U +ρ 1 σ1σ i+1,j+1,k + Ui 1,j 1,k U i+1,j 1,k U i 1,j+1,k ) h i h j + h i h j 1 + h i 1 h j + h i 1 h j 1 ( U +ρ 13 σ1σ 3 i+1,jk+1 + Ui 1,jk 1 U i+1,jk 1 U i 1,jk+1 ) h i h k + h i h k 1 + h i 1 h k + h i 1 h k 1 ( U +ρ 3 σσ 3 ij+1,k+1 + Uij 1,k 1 U ij+1,k 1 U ij 1,k+1 ) ], h j h k + h j h k 1 + h j 1 h k + h j 1 h k 1 where C 1 = r.5σ 1, C = r.5σ, ad C 3 = r.5σ Coditio for the o-oscillatory solutio Now, we derive the coditios uder which the explicit scheme for Eq. (4) will ot make spurious oscillatios by usig the idea i referece []. The, we rewrite Eq. (4) as (6) U +1 ( σ = τ 1 + (h i h i 1 )C 1 + σ + (h j h j 1 )C h i h i 1 h j h j 1 ) U + σ 3 + (h k h k 1 )C 3 h k h k 1 r + 1 τ + τ(σ 1 h ic 1 ) h i 1 (h i + h i 1 ) U i 1,jk + τ(σ 1 + h i 1C 1 ) h i (h i + h i 1 ) + τ(σ h jc ) h j 1 (h j + h j 1 ) U ij 1,k + τ(σ + h j 1C ) h j (h j + h j 1 ) U i+1,jk U ij+1,k

8 448 Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim + τ(σ 3 h kc 3 ) h k 1 (h k + h k 1 ) U 1 + τ(σ 3 + h k 1C 3 ) U+1 h k (h k + h k 1 ) ( U i+1,j+1k + Ui 1,j 1,k + τρ 1 σ 1 σ U i+1,j 1,k U ) i 1,j+1,k h i h j + h i h j 1 + h i 1 h j + h i 1 h j 1 ( U i+1,jk+1 + Ui 1,jk 1 + τρ 13 σ 1 σ U i+1,jk 1 U ) i 1,jk+1 3 h i h k + h i h k 1 + h i 1 h k + h i 1 h k 1 + τρ 3 σ σ 3 ( U ij+1,k+1 + U ij 1,k 1 U ij+1,k 1 U ij 1,k+1 h j h k + h j h k 1 + h j 1 h k + h j 1 h k 1 Next, we substitute U +1 = β +1 (1 r τ) ito Eq. (6), where the superscript for (1 r τ) represets a expoet. The, we obtai (7) β +1 = τ [ σ 1 + (h i h i 1 )C 1 + σ + (h j h j 1 )C (1 r τ)h i h i 1 (1 r τ)h j h j 1 + σ 3 + (h k h k 1 )C 3 (1 r τ)h k h k 1 r (1 r τ) + 1 τ(1 r τ) ] β τ(σ1 + h ic 1 ) (1 r τ)h i 1 (h i + h i 1 ) β i 1,jk + τ(σ1 + h i 1C 1 ) (1 r τ)h i (h i + h i 1 ) β i+1,jk τ(σ + h jc ) (1 r τ)h j 1 (h j + h j 1 ) β ij 1,k + τ(σ + h j 1C ) (1 r τ)h j (h j + h j 1 ) β ij+1,k τ(σ3 + h kc 3 ) (1 r τ)h k 1 (h k + h k 1 ) β 1 + τ(σ3 + h k 1C 3 ) (1 r τ)h k (h k + h k 1 β +1 + τρ ( 1σ 1 σ β i+1,j+1k + βi 1,j 1,k β i+1,j 1,k ) β i 1,j+1,k (1 r τ) h i h j + h i h j 1 + h i 1 h j + h i 1 h j 1 + τρ ( 13σ 1 σ 3 β i+1,jk+1 + βi 1,jk 1 β i+1,jk 1 ) β i 1,jk+1 (1 r τ) h i h k + h i h k 1 + h i 1 h k + h i 1 h k 1 + τρ ( 3σ σ 3 β ij+1,k+1 + βij 1,k 1 β ij+1,k 1 ) β ij 1,k+1. (1 r τ) h j h k + h j h k 1 + h j 1 h k + h j 1 h k 1 Sice all coefficiets of β i Eq. (7) are positive, the followig coditios should be satisfied. σ1 + (h i h i 1 )C 1 + σ + (h j h j 1 )C (8) (1 r τ)h i h i 1 (1 r τ)h j h j 1 + σ 3 + (h k h k 1 )C 3 r (1 r τ)h k h k 1 (1 r τ) + 1 τ(1 r τ) >. Also, whe we solve the log-trasformed BS equatio o the adaptive grid, we will use o-decreasig spatial step size, that is, h i 1 h i for ).

9 A compariso study of umerical method for the equity-liked securities 449 all i. By settig the lower boud of spatial step size as h mi satisfyig h mi h 1, coditio (8) gives (9) r τh mi + h mi τσ 1 τσ τσ 3 >. Therefore, we have the followig restrictio coditio for time step size (1) τ < h mi rh mi + σ 1 + σ +. σ 3 5. Numerical experimets I this sectio, we implemet umerical tests with our umerical method. For umerical tests, we cosider three-asset step-dow ELS optio as the example described i sectio 3.. We compare o-oscillatory explicit FDM ad implicit FDM with respect to computatioal costs ad errors Numerical treatmet for three-asset step-dow ELS optio pricig Before evaluatig the optio value of the test problem which is stated i sectio. 3., we first have to cosider two cases accordig to the iitial payoff. Let U ad V be the umerical solutios with payoffs which kock-i evet does ot happe ad happe, respectively. With the iitial payoffs (3) ad () which are described i Fig., we solve Eq. (5). After solvig Eq. (5) oce, we replace the values of less tha KIb i U with the values of V [3]. Ad the, we update the value of U by usig the FDM scheme. These processes are repeated from to T every time step. Also, at the early redemptio before τ = T, U ad V follow the coditios which are described i Table 1 ad Fig Numerical test We perform umerical experimets for pricig three-asset step-dow ELS. The parameters we have used are listed i sectio 3.. Also, the computatio domai is used as Ω = [1, ] [1, ] [1, ]. I Figs. 4-5, we ca see the iitial payoff fuctio ad umerical result at τ = 1 of U ad V, respectively. Now, we compare the results from MCS ad the FDM, the implicit scheme with OSM ad o-oscillatory explicit scheme. We focus o the value of U(1, 1, 1). I this test, MCS is performed 1 6 samples with τ = 1/144 usig atithetic variates of variace reductio, ad results from MCS are used as a referece value [7]. I order to

10 45 Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim U(X1, X, 1) 13 5 V(X1,X,1) 13 5 X 1 1 X1 X 1 1 X1 Figure 4. Iitial payoff fuctios of U ad V. U(X1, X, 1) 13 5 U(X1, X, 1) 13 5 X 1 1 X1 (a) X 1 1 X1 (b) Figure 5. Fial solutio of U at X 3 = 1 ad τ = 1 usig the (a) explicit ad (b) implicit scheme. reduce errors of simulatio, we calculate the average of the 1 MCS cases. The optio value obtaied from MCS is as a referece value, ad the MCS takes 145 secods at a time. For FDM tests, a various of o-uiform mesh for each directio is used. We fix τ for 1/36 i the implicit scheme ad adjust τ to prevet from spurious oscillatio of explicit scheme satisfyig the umber of early redemptios. [ a : h : b ] i Tables 6 meas that computatioal mesh is [ a, a+h, a+h,, b h, b, ]. All tests were performed o Itel(R) Core(TM) Duo E84 CPU@3.GHZ with 3.46GB of RAM loaded MATLAB 14a [6]. The optio value of U(1, 1, 1) ad absolute relative percet error with the value of MCS ad each FDM are show i Tables 6, where h is spatial step size for o-uiform mesh. Table shows that o-oscillatory explicit scheme is more superior tha implicit scheme with OSM i terms of computatioal cost ad the error with MCS i every spatial step. Note that the more space steps are take, the faster computatioal time is take i our scheme. O the other had, Tables

11 A compariso study of umerical method for the equity-liked securities preset that explicit scheme has similar error with implicit scheme. Overall, the explicit scheme is better tha the implicit if we use a small umber of mesh poits. However, if the spatial step size is small, the the implicit scheme is more efficiet tha the explicit scheme. h N x N τ U(1, 1, 1) Error CPU time (s) Table. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h. h N x N τ U(1, 1, 1) Error CPU time (s) Table 3. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h. h N x N τ U(1, 1, 1) Error CPU time (s) Table 4. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h. Next, we similarly examie the optio value U(1, 1, 1) with respect to chagig volatility ad riskless iterest rate. We use the solutio of MCS with variace reductio. Tables 9 ad 1 show the optio price ad error with chagig the volatility ad riskfree iterest

12 45 Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim h N x N τ U(1, 1, 1) Error CPU time (s) Table 5. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h. h N x N τ U(1, 1, 1) Error CPU time (s) Table 6. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h. rate. Chages of volatility ad iterest rate require the differet umber of time steps o o-oscillatory explicit scheme. O the other had, chages of iterest rate have less of a effect o the umber of time step tha the chages of volatility. The referece values whe σ =. ad σ =.4 are ad , respectively. Also, ad are referece values for r =.1 ad r =.5. h N x N τ U(1, 1, 1) Error CPU time (s) Table 7. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h (σ 1 = σ = σ 3 =.).

13 A compariso study of umerical method for the equity-liked securities 453 h N x N τ U(1, 1, 1) Error CPU time (s) Table 8. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h (σ 1 = σ = σ 3 =.4). h N x N τ U(1, 1, 1) Error CPU time (s) Table 9. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h (r =.1). h N x N τ U(1, 1, 1) Error CPU time (s) Table 1. Error ad the value of U(1, 1, 1) with o-uiform mesh [1 6:h: ] ad various values of h (r =.5). 6. Coclusios I geeral, the implicit scheme is used with Thomas algorithm. However, there are several drawbacks for practical calculatio. First, the implicit scheme has more time complexity i terms of computatios for solvig optio values. Next, it is hard to use the methodology for multidimesioal problem sice it is ecessary to apply the OSM or ADI. Hece we compared the explicit scheme satisfyig o-oscillatory coditio ad the implicit method for multi-dimesioal fiacial optios.

14 454 Mihyu Yoo, Darae Jeog, Seugsuk Seo, ad Juseok Kim Throughout the paper, the o-oscillatory explicit scheme is simple to implemet ad superior if we use a small umber of mesh poits. O the other had, the implicit scheme is more efficiet tha the explicit scheme if the spatial step size is small. Refereces [1] F. Black, ad M. Scholes, The pricig of optios ad corporate liabilities, J. Polit. Eco., 81 (1973), [] R. Zva, P. A. Forsyth, ad K. R. Vetzal, Robust umerical methods for PDE models of Asia optios, J. Comput. Fiac., 1 (1998), [3] D. Jeog, I. Wee, ad J. Kim. A operator splittig method for pricig the ELS optio, J. Korea Soc. Id. Appl. Math., 14(3), (1), [4] L. Trigeorgis, A log-trasformed biomial umerical aalysis method for valuig complex multi-optio ivestmets, J. Fiac. Quat. Aal., 6, (1991), [5] C. Reisiger, ad G. Wittum, O multigrid for aisotropic equatios ad variatioal iequalities Pricig multi-dimesioal Europea ad America optios. Comput. Visual. Sci., 7, (4), [6] MathWorks, Ic., MATLAB: the laguage of techical computig, The MathWorks, Natick, MA., [7] R. Seydel, Tools for computatioal fiace, Spriger Sciece & Busiess Media, 1. Mihyu Yoo Departmet of Fiacial Egieerig, Korea Uiversity, Seoul , Korea. ymh1989@korea.ac.kr Darae Jeog Departmet of Mathematics, Korea Uiversity, Seoul , Korea. tiayoyo@korea.ac.kr Seugsuk Seo Garam Aalytics, Yosei Uiversity, Seoul 1-749, Korea. sseo@gaalytics.co.kr Juseok Kim Departmet of Mathematics, Korea Uiversity, Seoul , Korea.

15 A compariso study of umerical method for the equity-liked securities Homepage: cfdkim

Hopscotch and Explicit difference method for solving Black-Scholes PDE

Hopscotch and Explicit difference method for solving Black-Scholes PDE Mälardale iversity Fiacial Egieerig Program Aalytical Fiace Semiar Report Hopscotch ad Explicit differece method for solvig Blac-Scholes PDE Istructor: Ja Röma Team members: A Gog HaiLog Zhao Hog Cui 0

More information

Positivity Preserving Schemes for Black-Scholes Equation

Positivity Preserving Schemes for Black-Scholes Equation Research Joural of Fiace ad Accoutig IN -97 (Paper) IN -7 (Olie) Vol., No.7, 5 Positivity Preservig chemes for Black-choles Equatio Mohammad Mehdizadeh Khalsaraei (Correspodig author) Faculty of Mathematical

More information

Economic Computation and Economic Cybernetics Studies and Research, Issue 2/2016, Vol. 50

Economic Computation and Economic Cybernetics Studies and Research, Issue 2/2016, Vol. 50 Ecoomic Computatio ad Ecoomic Cyberetics Studies ad Research, Issue 2/216, Vol. 5 Kyoug-Sook Moo Departmet of Mathematical Fiace Gacho Uiversity, Gyeoggi-Do, Korea Yuu Jeog Departmet of Mathematics Korea

More information

Statistics for Economics & Business

Statistics for Economics & Business Statistics for Ecoomics & Busiess Cofidece Iterval Estimatio Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for the mea ad the proportio How to determie

More information

Binomial Model. Stock Price Dynamics. The Key Idea Riskless Hedge

Binomial Model. Stock Price Dynamics. The Key Idea Riskless Hedge Biomial Model Stock Price Dyamics The value of a optio at maturity depeds o the price of the uderlyig stock at maturity. The value of the optio today depeds o the expected value of the optio at maturity

More information

Subject CT1 Financial Mathematics Core Technical Syllabus

Subject CT1 Financial Mathematics Core Technical Syllabus Subject CT1 Fiacial Mathematics Core Techical Syllabus for the 2018 exams 1 Jue 2017 Subject CT1 Fiacial Mathematics Core Techical Aim The aim of the Fiacial Mathematics subject is to provide a groudig

More information

A Hybrid Finite Difference Method for Valuing American Puts

A Hybrid Finite Difference Method for Valuing American Puts Proceedigs of the World Cogress o Egieerig 29 Vol II A Hybrid Fiite Differece Method for Valuig America Puts Ji Zhag SogPig Zhu Abstract This paper presets a umerical scheme that avoids iteratios to solve

More information

A random variable is a variable whose value is a numerical outcome of a random phenomenon.

A random variable is a variable whose value is a numerical outcome of a random phenomenon. The Practice of Statistics, d ed ates, Moore, ad Stares Itroductio We are ofte more iterested i the umber of times a give outcome ca occur tha i the possible outcomes themselves For example, if we toss

More information

1 Estimating sensitivities

1 Estimating sensitivities Copyright c 27 by Karl Sigma 1 Estimatig sesitivities Whe estimatig the Greeks, such as the, the geeral problem ivolves a radom variable Y = Y (α) (such as a discouted payoff) that depeds o a parameter

More information

Mixed and Implicit Schemes Implicit Schemes. Exercise: Verify that ρ is unimodular: ρ = 1.

Mixed and Implicit Schemes Implicit Schemes. Exercise: Verify that ρ is unimodular: ρ = 1. Mixed ad Implicit Schemes 3..4 The leapfrog scheme is stable for the oscillatio equatio ad ustable for the frictio equatio. The Euler forward scheme is stable for the frictio equatio but ustable for the

More information

Pricing 50ETF in the Way of American Options Based on Least Squares Monte Carlo Simulation

Pricing 50ETF in the Way of American Options Based on Least Squares Monte Carlo Simulation Pricig 50ETF i the Way of America Optios Based o Least Squares Mote Carlo Simulatio Shuai Gao 1, Ju Zhao 1 Applied Fiace ad Accoutig Vol., No., August 016 ISSN 374-410 E-ISSN 374-49 Published by Redfame

More information

43. A 000 par value 5-year bod with 8.0% semiaual coupos was bought to yield 7.5% covertible semiaually. Determie the amout of premium amortized i the 6 th coupo paymet. (A).00 (B).08 (C).5 (D).5 (E).34

More information

CAPITAL ASSET PRICING MODEL

CAPITAL ASSET PRICING MODEL CAPITAL ASSET PRICING MODEL RETURN. Retur i respect of a observatio is give by the followig formula R = (P P 0 ) + D P 0 Where R = Retur from the ivestmet durig this period P 0 = Curret market price P

More information

The Time Value of Money in Financial Management

The Time Value of Money in Financial Management The Time Value of Moey i Fiacial Maagemet Muteau Irea Ovidius Uiversity of Costata irea.muteau@yahoo.com Bacula Mariaa Traia Theoretical High School, Costata baculamariaa@yahoo.com Abstract The Time Value

More information

Monetary Economics: Problem Set #5 Solutions

Monetary Economics: Problem Set #5 Solutions Moetary Ecoomics oblem Set #5 Moetary Ecoomics: oblem Set #5 Solutios This problem set is marked out of 1 poits. The weight give to each part is idicated below. Please cotact me asap if you have ay questios.

More information

Institute of Actuaries of India Subject CT5 General Insurance, Life and Health Contingencies

Institute of Actuaries of India Subject CT5 General Insurance, Life and Health Contingencies Istitute of Actuaries of Idia Subject CT5 Geeral Isurace, Life ad Health Cotigecies For 2017 Examiatios Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical techiques which

More information

Subject CT5 Contingencies Core Technical. Syllabus. for the 2011 Examinations. The Faculty of Actuaries and Institute of Actuaries.

Subject CT5 Contingencies Core Technical. Syllabus. for the 2011 Examinations. The Faculty of Actuaries and Institute of Actuaries. Subject CT5 Cotigecies Core Techical Syllabus for the 2011 Examiatios 1 Jue 2010 The Faculty of Actuaries ad Istitute of Actuaries Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical

More information

Proceedings of the 5th WSEAS Int. Conf. on SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 17-19, 2005 (pp )

Proceedings of the 5th WSEAS Int. Conf. on SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 17-19, 2005 (pp ) Proceedigs of the 5th WSEAS It. Cof. o SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 7-9, 005 (pp488-49 Realized volatility estimatio: ew simulatio approach ad empirical study results JULIA

More information

Anomaly Correction by Optimal Trading Frequency

Anomaly Correction by Optimal Trading Frequency Aomaly Correctio by Optimal Tradig Frequecy Yiqiao Yi Columbia Uiversity September 9, 206 Abstract Uder the assumptio that security prices follow radom walk, we look at price versus differet movig averages.

More information

A New Approach to Obtain an Optimal Solution for the Assignment Problem

A New Approach to Obtain an Optimal Solution for the Assignment Problem Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 231-7064 Idex Copericus Value (2013): 6.14 Impact Factor (2015): 6.31 A New Approach to Obtai a Optimal Solutio for the Assigmet Problem A. Seethalakshmy

More information

CHAPTER 2 PRICING OF BONDS

CHAPTER 2 PRICING OF BONDS CHAPTER 2 PRICING OF BONDS CHAPTER SUARY This chapter will focus o the time value of moey ad how to calculate the price of a bod. Whe pricig a bod it is ecessary to estimate the expected cash flows ad

More information

An Empirical Study of the Behaviour of the Sample Kurtosis in Samples from Symmetric Stable Distributions

An Empirical Study of the Behaviour of the Sample Kurtosis in Samples from Symmetric Stable Distributions A Empirical Study of the Behaviour of the Sample Kurtosis i Samples from Symmetric Stable Distributios J. Marti va Zyl Departmet of Actuarial Sciece ad Mathematical Statistics, Uiversity of the Free State,

More information

1 Random Variables and Key Statistics

1 Random Variables and Key Statistics Review of Statistics 1 Radom Variables ad Key Statistics Radom Variable: A radom variable is a variable that takes o differet umerical values from a sample space determied by chace (probability distributio,

More information

Introduction to Financial Derivatives

Introduction to Financial Derivatives 550.444 Itroductio to Fiacial Derivatives Determiig Prices for Forwards ad Futures Week of October 1, 01 Where we are Last week: Itroductio to Iterest Rates, Future Value, Preset Value ad FRAs (Chapter

More information

The Valuation of the Catastrophe Equity Puts with Jump Risks

The Valuation of the Catastrophe Equity Puts with Jump Risks The Valuatio of the Catastrophe Equity Puts with Jump Risks Shih-Kuei Li Natioal Uiversity of Kaohsiug Joit work with Chia-Chie Chag Outlie Catastrophe Isurace Products Literatures ad Motivatios Jump Risk

More information

AN OPERATOR SPLITTING METHOD FOR PRICING THE ELS OPTION

AN OPERATOR SPLITTING METHOD FOR PRICING THE ELS OPTION J. KSIAM Vol.14, No.3, 175 187, 21 AN OPERATOR SPLITTING METHOD FOR PRICING THE ELS OPTION DARAE JEONG, IN-SUK WEE, AND JUNSEOK KIM DEPARTMENT OF MATHEMATICS, KOREA UNIVERSITY, SEOUL 136-71, KOREA E-mail

More information

5. Best Unbiased Estimators

5. Best Unbiased Estimators Best Ubiased Estimators http://www.math.uah.edu/stat/poit/ubiased.xhtml 1 of 7 7/16/2009 6:13 AM Virtual Laboratories > 7. Poit Estimatio > 1 2 3 4 5 6 5. Best Ubiased Estimators Basic Theory Cosider agai

More information

Rafa l Kulik and Marc Raimondo. University of Ottawa and University of Sydney. Supplementary material

Rafa l Kulik and Marc Raimondo. University of Ottawa and University of Sydney. Supplementary material Statistica Siica 009: Supplemet 1 L p -WAVELET REGRESSION WITH CORRELATED ERRORS AND INVERSE PROBLEMS Rafa l Kulik ad Marc Raimodo Uiversity of Ottawa ad Uiversity of Sydey Supplemetary material This ote

More information

Non-Inferiority Logrank Tests

Non-Inferiority Logrank Tests Chapter 706 No-Iferiority Lograk Tests Itroductio This module computes the sample size ad power for o-iferiority tests uder the assumptio of proportioal hazards. Accrual time ad follow-up time are icluded

More information

Research Article The Probability That a Measurement Falls within a Range of n Standard Deviations from an Estimate of the Mean

Research Article The Probability That a Measurement Falls within a Range of n Standard Deviations from an Estimate of the Mean Iteratioal Scholarly Research Network ISRN Applied Mathematics Volume 0, Article ID 70806, 8 pages doi:0.540/0/70806 Research Article The Probability That a Measuremet Falls withi a Rage of Stadard Deviatios

More information

Productivity depending risk minimization of production activities

Productivity depending risk minimization of production activities Productivity depedig risk miimizatio of productio activities GEORGETTE KANARACHOU, VRASIDAS LEOPOULOS Productio Egieerig Sectio Natioal Techical Uiversity of Athes, Polytechioupolis Zografou, 15780 Athes

More information

Calculation of the Annual Equivalent Rate (AER)

Calculation of the Annual Equivalent Rate (AER) Appedix to Code of Coduct for the Advertisig of Iterest Bearig Accouts. (31/1/0) Calculatio of the Aual Equivalet Rate (AER) a) The most geeral case of the calculatio is the rate of iterest which, if applied

More information

Section 3.3 Exercises Part A Simplify the following. 1. (3m 2 ) 5 2. x 7 x 11

Section 3.3 Exercises Part A Simplify the following. 1. (3m 2 ) 5 2. x 7 x 11 123 Sectio 3.3 Exercises Part A Simplify the followig. 1. (3m 2 ) 5 2. x 7 x 11 3. f 12 4. t 8 t 5 f 5 5. 3-4 6. 3x 7 4x 7. 3z 5 12z 3 8. 17 0 9. (g 8 ) -2 10. 14d 3 21d 7 11. (2m 2 5 g 8 ) 7 12. 5x 2

More information

Chapter 8. Confidence Interval Estimation. Copyright 2015, 2012, 2009 Pearson Education, Inc. Chapter 8, Slide 1

Chapter 8. Confidence Interval Estimation. Copyright 2015, 2012, 2009 Pearson Education, Inc. Chapter 8, Slide 1 Chapter 8 Cofidece Iterval Estimatio Copyright 2015, 2012, 2009 Pearso Educatio, Ic. Chapter 8, Slide 1 Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for

More information

Lecture 4: Probability (continued)

Lecture 4: Probability (continued) Lecture 4: Probability (cotiued) Desity Curves We ve defied probabilities for discrete variables (such as coi tossig). Probabilities for cotiuous or measuremet variables also are evaluated usig relative

More information

EVEN NUMBERED EXERCISES IN CHAPTER 4

EVEN NUMBERED EXERCISES IN CHAPTER 4 Joh Riley 7 July EVEN NUMBERED EXERCISES IN CHAPTER 4 SECTION 4 Exercise 4-: Cost Fuctio of a Cobb-Douglas firm What is the cost fuctio of a firm with a Cobb-Douglas productio fuctio? Rather tha miimie

More information

Chapter Four 1/15/2018. Learning Objectives. The Meaning of Interest Rates Future Value, Present Value, and Interest Rates Chapter 4, Part 1.

Chapter Four 1/15/2018. Learning Objectives. The Meaning of Interest Rates Future Value, Present Value, and Interest Rates Chapter 4, Part 1. Chapter Four The Meaig of Iterest Rates Future Value, Preset Value, ad Iterest Rates Chapter 4, Part 1 Preview Develop uderstadig of exactly what the phrase iterest rates meas. I this chapter, we see that

More information

0.1 Valuation Formula:

0.1 Valuation Formula: 0. Valuatio Formula: 0.. Case of Geeral Trees: q = er S S S 3 S q = er S S 4 S 5 S 4 q 3 = er S 3 S 6 S 7 S 6 Therefore, f (3) = e r [q 3 f (7) + ( q 3 ) f (6)] f () = e r [q f (5) + ( q ) f (4)] = f ()

More information

A New Constructive Proof of Graham's Theorem and More New Classes of Functionally Complete Functions

A New Constructive Proof of Graham's Theorem and More New Classes of Functionally Complete Functions A New Costructive Proof of Graham's Theorem ad More New Classes of Fuctioally Complete Fuctios Azhou Yag, Ph.D. Zhu-qi Lu, Ph.D. Abstract A -valued two-variable truth fuctio is called fuctioally complete,

More information

AY Term 2 Mock Examination

AY Term 2 Mock Examination AY 206-7 Term 2 Mock Examiatio Date / Start Time Course Group Istructor 24 March 207 / 2 PM to 3:00 PM QF302 Ivestmet ad Fiacial Data Aalysis G Christopher Tig INSTRUCTIONS TO STUDENTS. This mock examiatio

More information

Chapter 13 Binomial Trees. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull

Chapter 13 Binomial Trees. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull Chapter 13 Biomial Trees 1 A Simple Biomial Model! A stock price is curretly $20! I 3 moths it will be either $22 or $18 Stock price $20 Stock Price $22 Stock Price $18 2 A Call Optio (Figure 13.1, page

More information

Combining imperfect data, and an introduction to data assimilation Ross Bannister, NCEO, September 2010

Combining imperfect data, and an introduction to data assimilation Ross Bannister, NCEO, September 2010 Combiig imperfect data, ad a itroductio to data assimilatio Ross Baister, NCEO, September 00 rbaister@readigacuk The probability desity fuctio (PDF prob that x lies betwee x ad x + dx p (x restrictio o

More information

ad covexity Defie Macaulay duratio D Mod = r 1 = ( CF i i k (1 + r k) i ) (1.) (1 + r k) C = ( r ) = 1 ( CF i i(i + 1) (1 + r k) i+ k ) ( ( i k ) CF i

ad covexity Defie Macaulay duratio D Mod = r 1 = ( CF i i k (1 + r k) i ) (1.) (1 + r k) C = ( r ) = 1 ( CF i i(i + 1) (1 + r k) i+ k ) ( ( i k ) CF i Fixed Icome Basics Cotets Duratio ad Covexity Bod Duratios ar Rate, Spot Rate, ad Forward Rate Flat Forward Iterpolatio Forward rice/yield, Carry, Roll-Dow Example Duratio ad Covexity For a series of cash

More information

arxiv: v1 [q-fin.pr] 31 Jul 2015

arxiv: v1 [q-fin.pr] 31 Jul 2015 Variable Auity with GMWB: surreder or ot, that is the questio X. Luo a ad P.V. Shevcheko a a CSIRO Risk Aalytics Group, Australia This versio: 31 July, 2015 Email: Xiaoli.Luo@csiro.au arxiv:1507.08738v1

More information

CAPITAL PROJECT SCREENING AND SELECTION

CAPITAL PROJECT SCREENING AND SELECTION CAPITAL PROJECT SCREEIG AD SELECTIO Before studyig the three measures of ivestmet attractiveess, we will review a simple method that is commoly used to scree capital ivestmets. Oe of the primary cocers

More information

Bayes Estimator for Coefficient of Variation and Inverse Coefficient of Variation for the Normal Distribution

Bayes Estimator for Coefficient of Variation and Inverse Coefficient of Variation for the Normal Distribution Iteratioal Joural of Statistics ad Systems ISSN 0973-675 Volume, Number 4 (07, pp. 7-73 Research Idia Publicatios http://www.ripublicatio.com Bayes Estimator for Coefficiet of Variatio ad Iverse Coefficiet

More information

NPTEL DEPARTMENT OF INDUSTRIAL AND MANAGEMENT ENGINEERING IIT KANPUR QUANTITATIVE FINANCE END-TERM EXAMINATION (2015 JULY-AUG ONLINE COURSE)

NPTEL DEPARTMENT OF INDUSTRIAL AND MANAGEMENT ENGINEERING IIT KANPUR QUANTITATIVE FINANCE END-TERM EXAMINATION (2015 JULY-AUG ONLINE COURSE) NPTEL DEPARTMENT OF INDUSTRIAL AND MANAGEMENT ENGINEERING IIT KANPUR QUANTITATIVE FINANCE END-TERM EXAMINATION (2015 JULY-AUG ONLINE COURSE) READ THE INSTRUCTIONS VERY CAREFULLY 1) Time duratio is 2 hours

More information

Optimizing of the Investment Structure of the Telecommunication Sector Company

Optimizing of the Investment Structure of the Telecommunication Sector Company Iteratioal Joural of Ecoomics ad Busiess Admiistratio Vol. 1, No. 2, 2015, pp. 59-70 http://www.aisciece.org/joural/ijeba Optimizig of the Ivestmet Structure of the Telecommuicatio Sector Compay P. N.

More information

Unbiased estimators Estimators

Unbiased estimators Estimators 19 Ubiased estimators I Chapter 17 we saw that a dataset ca be modeled as a realizatio of a radom sample from a probability distributio ad that quatities of iterest correspod to features of the model distributio.

More information

Simulation Efficiency and an Introduction to Variance Reduction Methods

Simulation Efficiency and an Introduction to Variance Reduction Methods Mote Carlo Simulatio: IEOR E4703 Columbia Uiversity c 2017 by Marti Haugh Simulatio Efficiecy ad a Itroductio to Variace Reductio Methods I these otes we discuss the efficiecy of a Mote-Carlo estimator.

More information

Department of Mathematics, S.R.K.R. Engineering College, Bhimavaram, A.P., India 2

Department of Mathematics, S.R.K.R. Engineering College, Bhimavaram, A.P., India 2 Skewess Corrected Cotrol charts for two Iverted Models R. Subba Rao* 1, Pushpa Latha Mamidi 2, M.S. Ravi Kumar 3 1 Departmet of Mathematics, S.R.K.R. Egieerig College, Bhimavaram, A.P., Idia 2 Departmet

More information

Valuation of Variable Annuities with Guaranteed Minimum Withdrawal and Death Benefits via Stochastic Control Optimization

Valuation of Variable Annuities with Guaranteed Minimum Withdrawal and Death Benefits via Stochastic Control Optimization Valuatio of Variable Auities with Guarateed Miimum Withdrawal ad Death Beefits via Stochastic Cotrol Optimizatio Xiaoli Luo 1, ad Pavel V. Shevcheko 2 arxiv:1411.5453v2 [q-fi.cp] 7 Apr 2015 Draft, 1st

More information

The pricing of discretely sampled Asian and lookback options: a change of numeraire approach

The pricing of discretely sampled Asian and lookback options: a change of numeraire approach The pricig of discretely sampled Asia ad lookback optios 5 The pricig of discretely sampled Asia ad lookback optios: a chage of umeraire approach Jesper Adrease This paper cosiders the pricig of discretely

More information

Dr. Maddah ENMG 624 Financial Eng g I 03/22/06. Chapter 6 Mean-Variance Portfolio Theory

Dr. Maddah ENMG 624 Financial Eng g I 03/22/06. Chapter 6 Mean-Variance Portfolio Theory Dr Maddah ENMG 64 Fiacial Eg g I 03//06 Chapter 6 Mea-Variace Portfolio Theory Sigle Period Ivestmets Typically, i a ivestmet the iitial outlay of capital is kow but the retur is ucertai A sigle-period

More information

5 Statistical Inference

5 Statistical Inference 5 Statistical Iferece 5.1 Trasitio from Probability Theory to Statistical Iferece 1. We have ow more or less fiished the probability sectio of the course - we ow tur attetio to statistical iferece. I statistical

More information

r i = a i + b i f b i = Cov[r i, f] The only parameters to be estimated for this model are a i 's, b i 's, σe 2 i

r i = a i + b i f b i = Cov[r i, f] The only parameters to be estimated for this model are a i 's, b i 's, σe 2 i The iformatio required by the mea-variace approach is substatial whe the umber of assets is large; there are mea values, variaces, ad )/2 covariaces - a total of 2 + )/2 parameters. Sigle-factor model:

More information

Overlapping Generations

Overlapping Generations Eco. 53a all 996 C. Sims. troductio Overlappig Geeratios We wat to study how asset markets allow idividuals, motivated by the eed to provide icome for their retiremet years, to fiace capital accumulatio

More information

AUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY

AUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY AUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY Dr. Farha I. D. Al Ai * ad Dr. Muhaed Alfarras ** * College of Egieerig ** College of Coputer Egieerig ad scieces Gulf Uiversity * Dr.farha@gulfuiversity.et;

More information

These characteristics are expressed in terms of statistical properties which are estimated from the sample data.

These characteristics are expressed in terms of statistical properties which are estimated from the sample data. 0. Key Statistical Measures of Data Four pricipal features which characterize a set of observatios o a radom variable are: (i) the cetral tedecy or the value aroud which all other values are buched, (ii)

More information

Monopoly vs. Competition in Light of Extraction Norms. Abstract

Monopoly vs. Competition in Light of Extraction Norms. Abstract Moopoly vs. Competitio i Light of Extractio Norms By Arkadi Koziashvili, Shmuel Nitza ad Yossef Tobol Abstract This ote demostrates that whether the market is competitive or moopolistic eed ot be the result

More information

The material in this chapter is motivated by Experiment 9.

The material in this chapter is motivated by Experiment 9. Chapter 5 Optimal Auctios The material i this chapter is motivated by Experimet 9. We wish to aalyze the decisio of a seller who sets a reserve price whe auctioig off a item to a group of bidders. We begi

More information

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

arxiv: v2 [q-fin.pr] 15 Jan 2017 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

More information

Building a Dynamic Two Dimensional Heat Transfer Model part #1

Building a Dynamic Two Dimensional Heat Transfer Model part #1 Buildig a Dyamic Two Dimesioal Heat Trasfer Model part #1 - Tis is te first alf of a tutorial wic sows ow to build a basic dyamic eat coductio model of a square plate. Te same priciple could be used to

More information

1 + r. k=1. (1 + r) k = A r 1

1 + r. k=1. (1 + r) k = A r 1 Perpetual auity pays a fixed sum periodically forever. Suppose a amout A is paid at the ed of each period, ad suppose the per-period iterest rate is r. The the preset value of the perpetual auity is A

More information

Analytical Approximate Solutions for Stochastic Volatility. American Options under Barrier Options Models

Analytical Approximate Solutions for Stochastic Volatility. American Options under Barrier Options Models Aalytical Approximate Solutios for Stochastic Volatility America Optios uder Barrier Optios Models Chug-Gee Li Chiao-Hsi Su Soochow Uiversity Abstract This paper exteds the work of Hesto (99) ad itegrates

More information

ON THE RATE OF CONVERGENCE

ON THE RATE OF CONVERGENCE ON THE RATE OF CONVERGENCE OF BINOMIAL GREEKS SAN-LIN CHUNG WEIFENG HUNG HAN-HSING LEE* PAI-TA SHIH This study ivestigates the covergece patters ad the rates of covergece of biomial Greeks for the CRR

More information

Online appendices from The xva Challenge by Jon Gregory. APPENDIX 10A: Exposure and swaption analogy.

Online appendices from The xva Challenge by Jon Gregory. APPENDIX 10A: Exposure and swaption analogy. APPENDIX 10A: Exposure ad swaptio aalogy. Sorese ad Bollier (1994), effectively calculate the CVA of a swap positio ad show this ca be writte as: CVA swap = LGD V swaptio (t; t i, T) PD(t i 1, t i ). i=1

More information

Faculdade de Economia da Universidade de Coimbra

Faculdade de Economia da Universidade de Coimbra Faculdade de Ecoomia da Uiversidade de Coimbra Grupo de Estudos Moetários e Fiaceiros (GEMF) Av. Dias da Silva, 65 300-5 COIMBRA, PORTUGAL gemf@fe.uc.pt http://www.uc.pt/feuc/gemf PEDRO GODINHO Estimatig

More information

First determine the payments under the payment system

First determine the payments under the payment system Corporate Fiace February 5, 2008 Problem Set # -- ANSWERS Klick. You wi a judgmet agaist a defedat worth $20,000,000. Uder state law, the defedat has the right to pay such a judgmet out over a 20 year

More information

14.30 Introduction to Statistical Methods in Economics Spring 2009

14.30 Introduction to Statistical Methods in Economics Spring 2009 MIT OpeCourseWare http://ocwmitedu 430 Itroductio to Statistical Methods i Ecoomics Sprig 009 For iformatio about citig these materials or our Terms of Use, visit: http://ocwmitedu/terms 430 Itroductio

More information

REVISIT OF STOCHASTIC MESH METHOD FOR PRICING AMERICAN OPTIONS. Guangwu Liu L. Jeff Hong

REVISIT OF STOCHASTIC MESH METHOD FOR PRICING AMERICAN OPTIONS. Guangwu Liu L. Jeff Hong Proceedigs of the 2008 Witer Simulatio Coferece S. J. Maso, R. R. Hill, L. Möch, O. Rose, T. Jefferso, J. W. Fowler eds. REVISIT OF STOCHASTIC MESH METHOD FOR PRICING AMERICAN OPTIONS Guagwu Liu L. Jeff

More information

Models of Asset Pricing

Models of Asset Pricing 4 Appedix 1 to Chapter Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

Estimating Forward Looking Distribution with the Ross Recovery Theorem

Estimating Forward Looking Distribution with the Ross Recovery Theorem roceedigs of the Asia acific Idustrial Egieerig & Maagemet Systems Coferece 5 Estimatig Forward Lookig Distributio with the Ross Recovery Theorem Takuya Kiriu Graduate School of Sciece ad Techology Keio

More information

Models of Asset Pricing

Models of Asset Pricing APPENDIX 1 TO CHAPTER 4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

Models of Asset Pricing

Models of Asset Pricing APPENDIX 1 TO CHAPTER4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

1 The Black-Scholes model

1 The Black-Scholes model The Blac-Scholes model. The model setup I the simplest versio of the Blac-Scholes model the are two assets: a ris-less asset ba accout or bod)withpriceprocessbt) at timet, adarisyasset stoc) withpriceprocess

More information

SUPPLEMENTAL MATERIAL

SUPPLEMENTAL MATERIAL A SULEMENTAL MATERIAL Theorem (Expert pseudo-regret upper boud. Let us cosider a istace of the I-SG problem ad apply the FL algorithm, where each possible profile A is a expert ad receives, at roud, a

More information

Sequences and Series

Sequences and Series Sequeces ad Series Matt Rosezweig Cotets Sequeces ad Series. Sequeces.................................................. Series....................................................3 Rudi Chapter 3 Exercises........................................

More information

Chapter Four Learning Objectives Valuing Monetary Payments Now and in the Future

Chapter Four Learning Objectives Valuing Monetary Payments Now and in the Future Chapter Four Future Value, Preset Value, ad Iterest Rates Chapter 4 Learig Objectives Develop a uderstadig of 1. Time ad the value of paymets 2. Preset value versus future value 3. Nomial versus real iterest

More information

Maximum Empirical Likelihood Estimation (MELE)

Maximum Empirical Likelihood Estimation (MELE) Maximum Empirical Likelihood Estimatio (MELE Natha Smooha Abstract Estimatio of Stadard Liear Model - Maximum Empirical Likelihood Estimator: Combiatio of the idea of imum likelihood method of momets,

More information

arxiv: v1 [q-fin.pr] 11 Mar 2019 Abstract

arxiv: v1 [q-fin.pr] 11 Mar 2019 Abstract Pricig Formulae of Power Biary ad Normal Distributio Stadard Optios ad Applicatios Hyog-Chol O, Dae-Sug Choe, Faculty of Mathematics, Kim Il Sug Uiversity, Pyogyag, D P R K arxiv:903.0406v q-fi.pr Mar

More information

Portfolio Optimization for Options

Portfolio Optimization for Options Portfolio Optimizatio for Optios Yaxiog Zeg 1, Diego Klabja 2 Abstract Optio portfolio optimizatio for Europea optios has already bee studied, but more challegig America optios have ot We propose approximate

More information

Inferential Statistics and Probability a Holistic Approach. Inference Process. Inference Process. Chapter 8 Slides. Maurice Geraghty,

Inferential Statistics and Probability a Holistic Approach. Inference Process. Inference Process. Chapter 8 Slides. Maurice Geraghty, Iferetial Statistics ad Probability a Holistic Approach Chapter 8 Poit Estimatio ad Cofidece Itervals This Course Material by Maurice Geraghty is licesed uder a Creative Commos Attributio-ShareAlike 4.0

More information

STRAND: FINANCE. Unit 3 Loans and Mortgages TEXT. Contents. Section. 3.1 Annual Percentage Rate (APR) 3.2 APR for Repayment of Loans

STRAND: FINANCE. Unit 3 Loans and Mortgages TEXT. Contents. Section. 3.1 Annual Percentage Rate (APR) 3.2 APR for Repayment of Loans CMM Subject Support Strad: FINANCE Uit 3 Loas ad Mortgages: Text m e p STRAND: FINANCE Uit 3 Loas ad Mortgages TEXT Cotets Sectio 3.1 Aual Percetage Rate (APR) 3.2 APR for Repaymet of Loas 3.3 Credit Purchases

More information

Fixed Income Securities

Fixed Income Securities Prof. Stefao Mazzotta Keesaw State Uiversity Fixed Icome Securities Sample First Midterm Exam Last Name: First Name: Studet ID Number: Exam time is: 80 miutes. Total poits for this exam is: 400 poits Prelimiaries

More information

Supersedes: 1.3 This procedure assumes that the minimal conditions for applying ISO 3301:1975 have been met, but additional criteria can be used.

Supersedes: 1.3 This procedure assumes that the minimal conditions for applying ISO 3301:1975 have been met, but additional criteria can be used. Procedures Category: STATISTICAL METHODS Procedure: P-S-01 Page: 1 of 9 Paired Differece Experiet Procedure 1.0 Purpose 1.1 The purpose of this procedure is to provide istructios that ay be used for perforig

More information

We learned: $100 cash today is preferred over $100 a year from now

We learned: $100 cash today is preferred over $100 a year from now Recap from Last Week Time Value of Moey We leared: $ cash today is preferred over $ a year from ow there is time value of moey i the form of willigess of baks, busiesses, ad people to pay iterest for its

More information

A DOUBLE INCREMENTAL AGGREGATED GRADIENT METHOD WITH LINEAR CONVERGENCE RATE FOR LARGE-SCALE OPTIMIZATION

A DOUBLE INCREMENTAL AGGREGATED GRADIENT METHOD WITH LINEAR CONVERGENCE RATE FOR LARGE-SCALE OPTIMIZATION A DOUBLE INCREMENTAL AGGREGATED GRADIENT METHOD WITH LINEAR CONVERGENCE RATE FOR LARGE-SCALE OPTIMIZATION Arya Mokhtari, Mert Gürbüzbalaba, ad Alejadro Ribeiro Departmet of Electrical ad Systems Egieerig,

More information

Basic formula for confidence intervals. Formulas for estimating population variance Normal Uniform Proportion

Basic formula for confidence intervals. Formulas for estimating population variance Normal Uniform Proportion Basic formula for the Chi-square test (Observed - Expected ) Expected Basic formula for cofidece itervals sˆ x ± Z ' Sample size adjustmet for fiite populatio (N * ) (N + - 1) Formulas for estimatig populatio

More information

FOUNDATION ACTED COURSE (FAC)

FOUNDATION ACTED COURSE (FAC) FOUNDATION ACTED COURSE (FAC) What is the Foudatio ActEd Course (FAC)? FAC is desiged to help studets improve their mathematical skills i preparatio for the Core Techical subjects. It is a referece documet

More information

x satisfying all regularity conditions. Then

x satisfying all regularity conditions. Then AMS570.01 Practice Midterm Exam Sprig, 018 Name: ID: Sigature: Istructio: This is a close book exam. You are allowed oe-page 8x11 formula sheet (-sided). No cellphoe or calculator or computer is allowed.

More information

Quantitative Analysis

Quantitative Analysis EduPristie www.edupristie.com Modellig Mea Variace Skewess Kurtosis Mea: X i = i Mode: Value that occurs most frequetly Media: Midpoit of data arraged i ascedig/ descedig order s Avg. of squared deviatios

More information

ON DIFFERENTIATION AND HARMONIC NUMBERS

ON DIFFERENTIATION AND HARMONIC NUMBERS ON DIFFERENTIATION AND HARMONIC NUMBERS ERIC MORTENSON Abstract. I a paper of Adrews ad Uchimura [AU, it is show how differetiatio applied to hypergeometric idetities produces formulas for harmoic ad q-harmoic

More information

CAUCHY'S FORMULA AND EIGENVAULES (PRINCIPAL STRESSES) IN 3-D

CAUCHY'S FORMULA AND EIGENVAULES (PRINCIPAL STRESSES) IN 3-D GG303 Lecture 19 11/5/0 1 CAUCHY'S FRMULA AN EIGENVAULES (PRINCIPAL STRESSES) IN 3- I II Mai Topics A Cauchy s formula Pricipal stresses (eigevectors ad eigevalues) Cauchy's formula A Relates tractio vector

More information

Math 312, Intro. to Real Analysis: Homework #4 Solutions

Math 312, Intro. to Real Analysis: Homework #4 Solutions Math 3, Itro. to Real Aalysis: Homework #4 Solutios Stephe G. Simpso Moday, March, 009 The assigmet cosists of Exercises 0.6, 0.8, 0.0,.,.3,.6,.0,.,. i the Ross textbook. Each problem couts 0 poits. 0.6.

More information

Estimating Volatilities and Correlations. Following Options, Futures, and Other Derivatives, 5th edition by John C. Hull. Chapter 17. m 2 2.

Estimating Volatilities and Correlations. Following Options, Futures, and Other Derivatives, 5th edition by John C. Hull. Chapter 17. m 2 2. Estiatig Volatilities ad Correlatios Followig Optios, Futures, ad Other Derivatives, 5th editio by Joh C. Hull Chapter 17 Stadard Approach to Estiatig Volatility Defie as the volatility per day betwee

More information

A universal difference method for time-space fractional Black-Scholes equation

A universal difference method for time-space fractional Black-Scholes equation Xiaozhog et al Advaces i Differece Equatios (016) 016:71 DOI 101186/s1366-016-079-8 R E S E A R C H Ope Access A uiversal differece method for time-space fractioal Black-Scholes equatio Yag Xiaozhog, Wu

More information

Annual compounding, revisited

Annual compounding, revisited Sectio 1.: No-aual compouded iterest MATH 105: Cotemporary Mathematics Uiversity of Louisville August 2, 2017 Compoudig geeralized 2 / 15 Aual compoudig, revisited The idea behid aual compoudig is that

More information

The Limit of a Sequence (Brief Summary) 1

The Limit of a Sequence (Brief Summary) 1 The Limit of a Sequece (Brief Summary). Defiitio. A real umber L is a it of a sequece of real umbers if every ope iterval cotaiig L cotais all but a fiite umber of terms of the sequece. 2. Claim. A sequece

More information

KEY INFORMATION DOCUMENT CFD s Generic

KEY INFORMATION DOCUMENT CFD s Generic KEY INFORMATION DOCUMENT CFD s Geeric KEY INFORMATION DOCUMENT - CFDs Geeric Purpose This documet provides you with key iformatio about this ivestmet product. It is ot marketig material ad it does ot costitute

More information