Numerical Solutions to an integro-differential parabolic problem arising in the pricing of financial options in a Levy market
|
|
- Bertram Carr
- 6 years ago
- Views:
Transcription
1 September 2, 2 7:6 Quantitative Finance methodsqfstyle Quantitative Finance, Vol., No., Sep 2, 5 Numerical Solutions to an integro-differential parabolic problem arising in the pricing of financial options in a Levy market Ionuţ Florescu, Maria Cristina Mariani and Granville Sewell 4 Dept. of Mathematical Sciences, Stevens Institute of Technology Department of Mathematical Sciences, The University of Texas at El Paso (September 2, 2) We study the numerical solutions for an integro-differential parabolic problem modeling a process with jumps and stochastic volatility in Financial Mathematics. We present two general algorithms to calculate numerical solutions. The algorithms are implemented in PDE2D, a general purpose, partial differential equation solver. Keywords: integro-differential equation, numerical methods, option valuation. Introduction In Financial Mathematics the old problem of finding the price of derivatives (options, futures, etc.) leads to the study of Partial Differential Equations. The standard type of equations obtained are of parabolic type. In recent years, the complexity of the models used has increased and in turn this lead to more and more complicated equations for the derivative prices. Of particular interest in a type of differential equations containing an integral term. These equations aptly named Partial Integro-Differential Equations (PIDE) are difficult to solve and numerical methods specially constructed for them are not easy to find. In Florescu and Mariani (2) we study these type of problems and we prove the existence of the solution under general hypotheses about the integral term. In the present work we are extending the work on PIDE by providing a completely novel algorithm which is suggested in the proof of existence of the solution. Additionally, we present another algorithm which is a more classical finite element scheme coupled with a discretization of the integral term. In our numerical applications the two schemes are convergent to the same solution. The work is structured as follows. In Sections and 2 we Corresponding author. mcmariani@utep.edu Quantitative Finance ISSN print/issn online c 2 Taylor & Francis DOI:.8/469768YYxxxxxxxx
2 September 2, 2 7:6 Quantitative Finance methodsqfstyle 2 I. Florescu, M.C. Mariani and G. Sewell introduce the problem, as well as previous results and definitions; in Section 3 we describe the two algorithms that we will use to find numerical solutions; in Section 4 we explain how do we find numerical solutions by using PDE2D, a general purpose partial differential equation solver which has been used to solve many mathematical finance applications. In Sections 4 and 5 we present and discuss the results obtained. 2. Problem Motivation and general PIDE results In financial mathematics, usually the Black-Scholes model (Black and Scholes 973), or variants of the Black-Scholes model (Duffie 2, Hull 28, Ikeda 989, Jarrow 22, Merton 992) have been used for pricing derivatives on the equity. By applying the fundamental theorem of asset pricing (Harrison and Pliska 98, Delbaen and Schachermayer 994) one obtains different types of backward parabolic partial differential equations. In all these models, an important quantity is the volatility which is a measure of the fluctuation (risk) in the asset prices, and corresponds to the diffusion coefficient in the Black-Scholes equation. 2.. Stochastic Volatility models In the standard Black-Scholes model, a basic assumption is that the volatility is constant. It was soon discovered that this assumption does not allow matching an entire option chain (option values for different strike values). Several models were proposed, allowing the volatility to be modeled as a stochastic variable the so called stochastic volatility models (Hull and White 987, Wiggins 987, Scott 987, Chesney and Scott 989, Stein and Stein 99, Hagan et al. 22, Heston 993). To exemplify, in Heston (993) model the underlying security S follows, ds t = µs t dt + σ t S t dz t, where Z = {Z t } t is a standard Brownian motion. Unlike the classical model, the variance v(t) = σ 2 (t) also follows a stochastic process given by dv t = κ(θ v(t)) dt + γ v t dw t, where W = {W t } t is another standard Brownian motion. The correlation coefficient between W and Z is denoted by ρ: Cov(dZ t, dw t ) = ρ dt This leads to a generalized PDE in two state variables and one temporal variable: 2 vs2 2 U S 2 + ργvs 2 U v S + 2 vγ2 2 U U + rs v2 S + [κ(θ v) λv] U v ru + U t =. We should mention that any stochastic volatility model (no matter how simple) has the capability of matching an entire option chain for a fixed maturity (Fouque et al. 2, Florescu and Viens 28).
3 September 2, 2 7:6 Quantitative Finance methodsqfstyle Numerical Solutions for PIDE s Jump models, Levy processes At about the same time as the stochastic volatility models were developed researchers argued that the bad fitting to real data was caused by the path continuity of the price process. Thus, the resulting model may have difficulties fitting financial data exhibiting large fluctuations. The necessity of taking into account large market movements and a great amount of information arriving suddenly (i.e., a jump) led researchers to propose models with jumps. Merton (992), Andersen and Andreasen (2) model the stock as a jump-diffusion model (a geometric Brownian motion plus a compound Poisson process multiplying the stock process). Under this modeling assumption one may derive differential equations for option prices but these equations will contain an integral term coming from the compensator of the Poisson process. This has led to the study of partial integro-differential equations (PIDE). As an example, in the above cited work, the derivative value F (S, t) solves the following PIDE under appropriate boundary conditions: 2 σ2 S 2 F SS + (r λk)sf S + F t rf + λe{f (SY, t) F (S, t)} =. (2.) Here r denotes the riskless rate, λ the jump intensity, and k = E(Y ), where E is the expectation operator and the random variable Y measures the percentage change in the stock price if the jump occurs The jump-diffusion component is a particular case of a Lévy process and indeed stock evolution was soon modeled using these more general processes (e.g., Barndorff-Nielsen et al. (998), Madan et al. (998), Cont and Tankov (23), Geman (22)). Similar with the jump-diffusion case, when using any of these more general models, the resulting partial differential equation for derivative prices contains an integral term coming from the associated Lévy measure. Practical applications confirm that a Lévy like stochastic process appears to be the best fit when modeling high frequency data (see Mariani et al. (29) and its references) Combining the stochastic volatility and the Levy component Each of the two proposed models are doing something different. One is allowing for varying fluctuations and the other cope well with large amount of information driving the price suddenly up or down. It is only natural to attempt and combine the two modeling assumptions. As we shall see this complicates the resulting equations for the stock price but they do remain of the same integro differential type. As a simple example we consider a jump-diffusion process with volatility replaced by a geometric Brownian motion: ds t = S t (µdt + σ t dz t + Y (g)dn t ) dσ t = σ t (αdt + βdw t ), where Z and W are two standard Brownian motion with correlation coefficient ρ, N t is a Poisson process with intensity λ, and Y (g) is the jump amplitude random variable with density g. The jump part may be written in a perhaps more traditional way as a compound Poisson process: Nt i= Y i. To obtain the price of a derivative in the model above we follow Merton (992) and we obtain
4 September 2, 2 7:6 Quantitative Finance methodsqfstyle 4 I. Florescu, M.C. Mariani and G. Sewell the following PIDE, F t + 2 σ2 S 2 2 F S σ2 β 2 2 F σ 2 + ρσ2 βs 2 F S σ + λ R + (r λk)s F S 2 ρσ2 β F σ [F (SY, σ, t) F (S, σ, t)]g(y )dy rf =. (2.2) Here, r denotes the riskless rate, and k = E (Y ), where E is the expectation operator, g is the density of the Y random variable. We mention that the quantity Y measures the percentage change in the stock price if the jump occurs, for more discussion see (Merton 992, Section 9.2). We note that the derivative price still solves a PIDE and this is a general characteristic of these type of processes A general existence result The previous discussion motivates to consider more general integro-differential parabolic problems. In recent work Florescu and Mariani (2) the authors prove the existence of solutions to a general partial integro-differential equation in an unbounded smooth domain (2.3) using the method of upper and lower solutions. Since we shall use some of the parts of the proof in the algorithm we state the problem and the essential results here as well. Let Ω R d an unbounded smooth domain, L a second order elliptic operator in non divergence form, i.e., Lu := d a ij (x, t)u xix j + i,j= d b i (x, t)u xi + c(x, t)u, with coefficients of L in the Hölder Space C δ,δ/2 ( Ω [, T ] ) and satisfying the following conditions: λ v 2 i= d a ij (x, t)v i v j Λ v 2 ( λ Λ) i,j= b i (x, t) C; c(x, t). These are classical conditions to insure the strict ellipticity of the operator. The problem we are solving is: Lu u t = G(t, u) in Ω (, T ) u(x, ) = u (x) on Ω {} u(x, t) = h(x, t) on Ω (, T ) (2.3) The operator G is a completely continuous integral operator as the ones defined in (2.) and (2.2), coming from the jump distribution. More precisely, we assume that G(t, u) = Ω g(x, t, u)dx, where g is any continuous function. The proof however is not constrained to integral operators, indeed any operator with the same properties will work. In this general model, the case in which g is increasing with respect to u and all jumps are positive corresponds to the evolution of a call option near a crash.
5 September 2, 2 7:6 Quantitative Finance methodsqfstyle Numerical Solutions for PIDE s 5 The existence result in Florescu and Mariani (2) reads as follows: Theorem 2. : Let L and G be the operators defined above. Assume that either: G is nonincreasing with respect to u, or there exist some continuous, and increasing one dimensional function f such that G(t, u) f(u) is nonincreasing with respect to u Furthermore, assume there exist α and β a lower and an upper solution of the problem with α β in Ω (, T ). Then, the problem (2.3) admits a solution u such that α u β in Ω (, T ). Remark : The assumptions above say that the result is applicable in particular for any G that is dominated by a polynomial. Remark 2 : We remind the interested reader that a smooth function u is called an upper (lower) solution of problem (2.3) if Lu u t ( ) G(t, u) in Ω (, T ) u(x, ) ( ) u (x) on Ω {} u(x, t) ( ) h(x, t) on Ω (, T ) The new algorithm we propose is inspired by the proof of the theorem and we outline here the main steps in the proof. We construct a series of regular PDE discretizations (for n {, 2,..., K}) for the problem (2.3), Lu n+ u n+ t = G(t, u n ) in U (, T ) u n+ (x, ) = u (x) in U {} u n+ (x, t) = φ U (x, t) in U (, T ) Note that the boundary conditions are identical with the original problem 2.3, but by using the known solution from the previous step in the G term the main equation loses its integral term and becomes a regular PDE. (2.4) Under the hypotheses given each of the systems (2.4) have a unique solution u n+ W 2, p (V ) where W 2, p (V ) = {v L p : v xi, v xi x j L p, v t L p } (see e.g., Krylov (996), Lieberman (996)) We prove that α u n β, n using the special properties of L and G We show that {u n } is a Cauchy sequence in Wp 2, (V ) using that G is a completely continuous operator, the previous step and the Lebesgue s dominated convergence theorem We conclude that u n u in the Wp 2, -norm, and then u is a strong solution of the problem. 3. Algorithms for numerical solutions We propose two algorithms to find the solution of the problem (2.3). The first algorithm is new and is inspired by the constructive proof of the theorem. It can be considered that this algorithm handles the integral part of the PIDE using an implicit method
6 September 2, 2 7:6 Quantitative Finance methodsqfstyle 6 I. Florescu, M.C. Mariani and G. Sewell because, after convergence, the u appearing in the integral term is evaluated at the current value of t. The second algorithm is a more standard explicit approximation, where the u appearing in the integral is evaluated at the previous time step. We mention Cont and Voltchkova (25) where the authors provide a similar approach for a more specific PIDE in only one spatial dimension and using a finite difference scheme (we are using finite elements). We have no proof of convergence for this algorithm and indeed we observed that for stability the time step had to be less than a threshhold value, which decreases as λ increases. This is not surprising, since while both algorithms use implicit methods for the partial differential equations, the second algorithm handles the integral part of the PIDE using an explicit scheme. The purpose of presenting this algorithm is twofold. One, we want to compare the solution provided by this algorithm with the solution provided by the other algorithm in hopes of showing that under reasonable parameter values and discretizations the two solutions are close. Two, this algorithm is much faster than the other algorithm (in our tests 5 times faster) and thus there is interest in applying it in the financial applications where computer time is very important. 3.. Algorithm The new algorithm we propose comes from the proof of theorem 2.. We first find α() lower and β() upper solutions of the problem and we ensure that the operators L and G verify the hypotheses of the theorem. We start with u = α the lower solution. For every n we solve the system (2.4). The PDE in this system cannot be solved analytically and we implement a finite element scheme for it. We take a three dimensional grid in t [, T ], S [, S max ], σ [, σ max ]. The upper bounds are suitably chosen large numbers. This grid is kept the same for all iterations. We approximate the integrals in G(t, u n ) using a midpoint rule, with the required values of u n (S, σ, t) interpolated from the values saved at all the points on the grid, using quadratic interpolation in S and σ and linear interpolation in t. Since the boundary condition is usually at T we solve the resulting partial differential equation backwards in time using a PDE solver (we will use PDE2D in our experiments). The result is u n+ calculated at all grid points and to be used in the next iteration. The u n sequence converges to the solution of the main system (2.3). Thus the algorithm stops when the maximum difference between two consecutive iterations at all points on the grid is small Algorithm 2 Algorithm 2 may be viewed as an explicit numerical scheme. It does not use the discretization (2.4), instead it works directly with the given system (2.3). It works backward in time from t = T. Specifically: Starts with the final condition at T. When a PDE solver (we will use PDE2D) is used to solve (2.3) backwards from t i+ to t i, G(t, u(t)) is replaced by G(t i+, u(t i+ )), the integral is approximated using a midpoint rule, and the required values of u at t i+ are approximated using quadratic interpolation to the values saved at the S, σ grid points. When t = we stop. The solution approximation is obtained for all grid points.
7 September 2, 2 7:6 Quantitative Finance methodsqfstyle Numerical Solutions for PIDE s 7 4. Numerical solutions using PDE2D In order to test the two algorithms we consider the PIDE resulting from the model combining stochastic volatility and jumps (system (2.2)). In this system we use the following parameter values. Parameters characterizing the stock: initial value S =, riskfree rate r =.5, volatility of volatility β =.4, correlation ρ =.3. The drift parameters µ and α do not enter in the PIDE. Parameters characterizing the jump component: jump frequency λ = 5 (expected five jumps per year) and the expected jump percentage k =.5, which can be calculated from the jump distribution. We are using a Gaussian mixture model for g(y ) =.75 g (Y ) +.25 g 2 (Y ), where g (Y ) and g 2 (Y ) are Normal densities with means.7 and.3 respectively, and variances equal to.. The values are chosen to insure that a jump is noticeable and separable from the stochastic variability of the model. Fan and Wang (27) show that in the presence of both jumps and stochastic volatility it is hard to separate the small sized jumps from the large stock variability. Both algorithms are solved for the value of an European Call in this model. Specifically, the option maturity was T =, and the strike price K =. The boundary conditions characterizing this option type were: F (S, σ, T ) = max(s K, ), F (, σ, t) =, and F S (S max, σ, t) =. We had no boundary conditions at σ = or σ = σ max. The algorithms were implemented using PDE2D, a general purpose partial differential equation solver (Sewell 25, 2) available from Visual Numerics, Inc. ( which has been used to solve many mathematical finance applications (Topper 25). PDE2D solves linear or nonlinear, steady-state, time-dependent and eigenvalue problems, in D intervals, general 2D regions (with curved boundaries), and a wide range of simple 3D regions. It has a sophisticated GUI interface which makes it extremely easy to use. It can solve D, 2D or 3D problems similar to the problems presented (2.),(2.2), using a collocation finite element method with cubic Hermite basis functions to discretize the spatial derivatives and adaptive finite difference methods to discretize the time derivatives. Without the integral term, a problem such as (2.2) is a very straightforward PDE2D application, and can be solved with little user effort. The integral term in (2.2), however, required special efforts. The two algorithms were implemented using PDE2D, and we present below the specific details of each implementation. For both algorithms we use S max = 4, σ max =. A equally spaced grid was constructed on [, T ] [, S max ] [, σ max ]. A simple modification of the example 2.3 in Florescu and Mariani (2) shows that the operators L and G would satisfy the hypotheses of the theorem if σ min >, and that upper and lower solutions exist with the lower solution α. Algorithm implementation. The iteration (2.4) is repeatedly solved backward in time from a final condition using PDE2D. The first iteration uses u = α =, so that no integral term appears, and the solution was saved for each point on the grid. Each iteration thereafter, (2.4) is solved with the F in the integral term approximated by interpolating the solution saved on the previous iteration, using quadratic interpolation in S and σ, and linear interpolation in time. For our examples, g(y ) was negligible outside of < Y < 2, so the integral limits in (2.4) were taken to be [, 2], and a midpoint rule was used to numerically approximate the integral. Since Y can be as large as 2, however, values of F are needed at points SY beyond S max, where they have not been calculated. For these values, a value of F is extrapolated using the boundary condition F S =, that is: F (SY, σ, t) F (S max, σ, t) + SY S max
8 September 2, 2 7:6 Quantitative Finance methodsqfstyle 8 I. Florescu, M.C. Mariani and G. Sewell σ 5 S 5 Figure. F(S,σ,) for the integral term equal to (Algorithm first iteration) Algorithm 2 implementation. This time (2.2) is only solved once using PDE2D, and at each time point t i+ = (i + ) t ( t = T/), PDE2D saves the solution on the S σ grid. Then when it integrates backward from t i+ to t i, the terms F (SY, σ, t) and F (S, σ, t) appearing in the integral are approximated by F (SY, σ, t i+ ) and F (S, σ, t i+ ), and these values are obtained using quadratic interpolation to the solution saved at the grid points at t = t i+. The numerical approximation to the integral is the same as in Algorithm. 4.. Results The solution is a function of three variables, so it is impossible to view on a regular 3D plot. For this reason we only plot the solutions when t = since usually this is the most relevant value. Figure shows the solution from Algorithm after the first iteration. This is the same as the homogeneous problem with λ =, that is, without the integral term, corresponding to the stochastic volatility model where the jump diffusion term does not exist. This figure is provided for comparison to view the effect of the integral term in the PIDE.
9 September 2, 2 7:6 Quantitative Finance methodsqfstyle Numerical Solutions for PIDE s σ 2 S 3 4 Figure 2. The final solution F(S,σ,) after convergence of Algorithm Algorithm showed significant oscillations near σ = for the first few iterations, but the oscillations damped after further iteration, and (2.4) converged after 3 iterations, using our stopping criterion of maximum difference < 3 5. The oscillations are believed to be due to the lack of a boundary condition at σ =. Figure 2 shows F (S, σ, ) after convergence, and Figure 3 shows the same solution in the vicinity of the origin. Algorithm 2 produced a solution which agreed with that produced by Algorithm closely (which is why we do not plot it), the maximum difference between the two solutions on the three dimensional grid was about.2. Since Algorithm 2 does only one iteration (only solves (2.2) once), it is expected to be the more efficient algorithm, and indeed it required about 5 times less computer time for this example. To give more details about the order of magnitude of the difference between the two solutions we calculate and plot in Figure 4 the Root Mean Squared Error defined as RMSE(t) = n (F (S, σ, t) F 2 (S, σ, t)) 2, and the Mean Absolute Deviation defined as MAD(t) = n F (S, σ, t) F 2 (S, σ, t). F and F 2 are the solutions obtained using Algorithm respectively Algorithm 2 and the sum runs over all the grid points. At t = we do not plot since the two solutions coincide with the boundary value. Finally, since at t = we know the S variable but the σ is unobservable we present in Figure 5 the option price evolution depending on the realized volatility value σ.
10 September 2, 2 7:6 Quantitative Finance methodsqfstyle I. Florescu, M.C. Mariani and G. Sewell σ 5 S 5 Figure 3. The final solution F(S,σ,) a closeup at origin 4.2. A discussion of the influence of the integral term The approach we present in the current article deals specifically with partial differential equations containing an integral term. As we may observe in any of the equations presented, λ, the parameter which governs the frequency of the jumps, ends multiplying the integral term in the equation. Thus we chose to present the results when λ increases in value and thus the integral term becomes dominant. In Figure 6 we present the graphs of the solution at t =, obtained when λ changes in value. We see in Figure 6(a) the result of the regular PDE with no integral term, and as λ increases the integral term becomes dominant. Note that when λ = (Figure 6(d)) the solution is almost F = S, which is to be expected since F = S is the steady-state solution when λ is infinite (and satisfies the boundary conditions).
11 September 2, 2 7:6 Quantitative Finance methodsqfstyle Numerical Solutions for PIDE s Root mean squared error Mean absolute deviation times Figure 4. Two measures of difference between the two solutions 5. Conclusions In this article we present a completely novel algorithm proven to converge to solve partial differential equations with an integral term. We compare the new algorithm with another more explicit algorithm that is not proven to converge. In our numerical example we found the difference between the two solutions to be minimal. This opens the way for numerical schemes to approximate more general PIDE produced by Lévy models of the stock price. The algorithms presented should work for these problems although a more detailed study of convergence for Algorithm 2 should be performed. In our numerical experiments we discovered that increasing λ produced instability which was fixed by using a smaller time step. This is very similar with
12 September 2, 2 7:6 Quantitative Finance methodsqfstyle 2 REFERENCES option Figure 5. The call option value depending on the volatility value at time t = vol the traditional Explicit Finite Difference scheme and a relation between parameters that would guarantee convergence would be very beneficial. References Andersen, L. and J. Andreasen (2). Jump-diffusion processes: Volatility smile fitting and numerical methods for option pricing. Review of Derivatives Research 4, Barndorff-Nielsen, O. E., J. L. Jensen, and M. Sorensen (998, Dec.). Some stationary processes in discrete and continuous time. Advances in Applied Probability 3 (4), Black, F. and M. Scholes (973). The valuation of options and corporate liability. Journal of Political Economy 8,
13 September 2, 2 7:6 Quantitative Finance methodsqfstyle REFERENCES F 4 F σ 5 S 5.5 σ 5 S 5 (a) λ = (b) λ = 2 5 F F 5.5 σ 5 S 5.5 σ 5 S 5 (c) λ = (d) λ = Figure 6. Changes in the solution at t = when the importance of the integral term increases Chesney, M. and L. Scott (989). Pricing european currency options: a comparison of the modified black-scholes model and a random variance model. J. Finan. Quant. Anal. 24, Cont, R. and P. Tankov (23). Financial modeling with jumps processes. CRC Financial mathematics series. Chapman & Hall. Cont, R. and E. Voltchkova (25). A finite difference scheme for option pricing in jump diffusion and. SIAM Journal of Numerical Analysis 43 (4), Delbaen, F. and W. Schachermayer (994). A general version of the fundamental theorem of asset pricing. Mathematische Annalen (3), Duffie, D. (2). Dynamic Asset Pricing Theory (3rd ed.). Princeton University Press. Fan, J. and Y. Wang (27). Multi-scale jump and volatility analysis for high-frequency financial data. Journal of the American Statistical Association 2, Florescu, I. and M. Mariani (2, May). Solutions to integro-differential parabolic problems arising in the pricing of financial options in a levy market. Electronic Journal of Differential Equations 2 (2)(62),. Florescu, I. and F. Viens (28, April). Stochastic volatility: option pricing using a multinomial recombining tree. Applied Mathematical Finance 5 (2), 5 8.
14 September 2, 2 7:6 Quantitative Finance methodsqfstyle 4 REFERENCES Fouque, J.-P., G. Papanicolaou, and K. R. Sircar (2). Derivatives in Financial Markets with Stochastic Volatility. Cambridge University Press. Geman, H. (22). Pure jump Levy processes for asset price modeling. Journal of Banking and Finance 26, Hagan, P., D. Kumar, A. Lesniewski, and D. Woodward (22). Managing smile risk. Wilmott Magazine. Harrison, M. and S. Pliska (98). Martingales and stochastic integrals in the theory of continuous trading. Stochastic Processes and their Applications, Heston, S. L. (993). A closed-form solution for options with stochastic volatility with applications to bond and currency options. Review of Financial Studies 6 (2), Hull, J. C. (28). Options, Futures, and other Derivatives (7th ed.). Prentice Hall. Hull, J. C. and A. D. White (987, June). The pricing of options on assets with stochastic volatilities. Journal of Finance 42 (2), Ikeda, N. (989). Stochastic Differential Equations and Diffusion Processes (2nd Revised ed.). North-Holland. Jarrow, R. A. (22). Modelling Fixed Income Securities and Interest Rate Options (2nd ed.). Stanford Economics and Finance. Krylov, N. (996). Lectures on elliptic and parabolic equations in Hölder Spaces, Volume 2 of Graduate Studies in Mathematics. American Mathematical Society. Lieberman, G. M. (996). Second order parabolic differential equations. World Scientific Publ. Madan, D. B., P. P. Carr, and E. C. Chang (998). The variance gamma process and option pricing. European Finance Review 2 (), Mariani, M., I. Florescu, M. B. Varela, and E. Ncheuguim (29, April). Long correlations and levy models applied to the study of memory effects in high frequency (tick) data. Physica A 388 (8), Merton, R. (992). Continuous-Time Finance. Wiley-Blackwell. Scott, L. O. (987). Option pricing when the variance changes randomly: Theory, estimation, and an application. The Journal of Financial and Quantitative Analysis 22 (4), Sewell, G. (25). The Numerical Solution of Ordinary and Partial Differential Equations (2nd ed.). John Wiley & Sons. Sewell, G. (2, May). Solving pdes in non-rectangular 3d regions using a collocation finite element method. Advances in Engineering Software 4, Stein, E. M. and J. C. Stein (99). Stock price distributions with stochastic volatility: An analytic approach. Review of Financial Studies 4 (4), Topper, J. (25). Financial Engineering with Finite Elements. John Wiley & Sons. Wiggins, J. B. (987, December). Option values under stochastic volatility: Theory and empirical estimates. Journal of Financial Economics 9 (2), Dr. Ionuţ Florescu ( Ionut.Florescu@stevens.edu) Department of Mathematical Sciences, Stevens Institute of Technology Castle Point on Hudson, Hoboken, NJ 73 USA Dr. Maria Christina Mariani ( mcmariani@utep.edu) Department of Mathematical Sciences, The University of Texas at El Paso Bell Hall 24, El Paso, Texas , USA
15 September 2, 2 7:6 Quantitative Finance methodsqfstyle REFERENCES 5 Dr. Granville Sewell ( sewell@utep.edu) Department of Mathematical Sciences, The University of Texas at El Paso Bell Hall 24, El Paso, Texas , USA
SOLUTIONS TO INTEGRO-DIFFERENTIAL PARABOLIC PROBLEMS ARISING IN THE PRICING OF FINANCIAL OPTIONS IN A LEVY MARKET
Electronic Journal of Differential Equations, Vol. 2010(2010), No. 62, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SOLUTIONS TO INTEGRO-DIFFERENTIAL
More informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationGreek parameters of nonlinear Black-Scholes equation
International Journal of Mathematics and Soft Computing Vol.5, No.2 (2015), 69-74. ISSN Print : 2249-3328 ISSN Online: 2319-5215 Greek parameters of nonlinear Black-Scholes equation Purity J. Kiptum 1,
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More informationTHE USE OF NUMERAIRES IN MULTI-DIMENSIONAL BLACK- SCHOLES PARTIAL DIFFERENTIAL EQUATIONS. Hyong-chol O *, Yong-hwa Ro **, Ning Wan*** 1.
THE USE OF NUMERAIRES IN MULTI-DIMENSIONAL BLACK- SCHOLES PARTIAL DIFFERENTIAL EQUATIONS Hyong-chol O *, Yong-hwa Ro **, Ning Wan*** Abstract The change of numeraire gives very important computational
More informationYoungrok Lee and Jaesung Lee
orean J. Math. 3 015, No. 1, pp. 81 91 http://dx.doi.org/10.11568/kjm.015.3.1.81 LOCAL VOLATILITY FOR QUANTO OPTION PRICES WITH STOCHASTIC INTEREST RATES Youngrok Lee and Jaesung Lee Abstract. This paper
More informationFX Smile Modelling. 9 September September 9, 2008
FX Smile Modelling 9 September 008 September 9, 008 Contents 1 FX Implied Volatility 1 Interpolation.1 Parametrisation............................. Pure Interpolation.......................... Abstract
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
More informationExact Sampling of Jump-Diffusion Processes
1 Exact Sampling of Jump-Diffusion Processes and Dmitry Smelov Management Science & Engineering Stanford University Exact Sampling of Jump-Diffusion Processes 2 Jump-Diffusion Processes Ubiquitous in finance
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationOption Valuation with Sinusoidal Heteroskedasticity
Option Valuation with Sinusoidal Heteroskedasticity Caleb Magruder June 26, 2009 1 Black-Scholes-Merton Option Pricing Ito drift-diffusion process (1) can be used to derive the Black Scholes formula (2).
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationOptimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models
Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University University of Missouri-Kansas City Department of Mathematics
More informationThe Numerical Valuation of Options with Underlying Jumps
The Numerical Valuation of Options with Underlying Jumps Gunter H. Meyer School of Mathematics Georgia Institute of Technology Atlanta, GA, 30332-060, U.S.A. Abstract A Black-Scholes type model for American
More informationPricing of a European Call Option Under a Local Volatility Interbank Offered Rate Model
American Journal of Theoretical and Applied Statistics 2018; 7(2): 80-84 http://www.sciencepublishinggroup.com/j/ajtas doi: 10.11648/j.ajtas.20180702.14 ISSN: 2326-8999 (Print); ISSN: 2326-9006 (Online)
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets Liuren Wu ( c ) The Black-Merton-Scholes Model colorhmoptions Markets 1 / 18 The Black-Merton-Scholes-Merton (BMS) model Black and Scholes (1973) and Merton
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationTime-changed Brownian motion and option pricing
Time-changed Brownian motion and option pricing Peter Hieber Chair of Mathematical Finance, TU Munich 6th AMaMeF Warsaw, June 13th 2013 Partially joint with Marcos Escobar (RU Toronto), Matthias Scherer
More informationQuadratic hedging in affine stochastic volatility models
Quadratic hedging in affine stochastic volatility models Jan Kallsen TU München Pittsburgh, February 20, 2006 (based on joint work with F. Hubalek, L. Krawczyk, A. Pauwels) 1 Hedging problem S t = S 0
More informationOptimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University
Optimal Hedging of Variance Derivatives John Crosby Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation at Baruch College, in New York, 16th November 2010
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More informationChapter 5 Finite Difference Methods. Math6911 W07, HM Zhu
Chapter 5 Finite Difference Methods Math69 W07, HM Zhu References. Chapters 5 and 9, Brandimarte. Section 7.8, Hull 3. Chapter 7, Numerical analysis, Burden and Faires Outline Finite difference (FD) approximation
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets (Hull chapter: 12, 13, 14) Liuren Wu ( c ) The Black-Scholes Model colorhmoptions Markets 1 / 17 The Black-Scholes-Merton (BSM) model Black and Scholes
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationAn IMEX-method for pricing options under Bates model using adaptive finite differences Rapport i Teknisk-vetenskapliga datorberäkningar
PROJEKTRAPPORT An IMEX-method for pricing options under Bates model using adaptive finite differences Arvid Westlund Rapport i Teknisk-vetenskapliga datorberäkningar Jan 2014 INSTITUTIONEN FÖR INFORMATIONSTEKNOLOGI
More informationPricing Implied Volatility
Pricing Implied Volatility Expected future volatility plays a central role in finance theory. Consequently, accurate estimation of this parameter is crucial to meaningful financial decision-making. Researchers
More informationApplication of Moment Expansion Method to Option Square Root Model
Application of Moment Expansion Method to Option Square Root Model Yun Zhou Advisor: Professor Steve Heston University of Maryland May 5, 2009 1 / 19 Motivation Black-Scholes Model successfully explain
More informationFinite Difference Approximation of Hedging Quantities in the Heston model
Finite Difference Approximation of Hedging Quantities in the Heston model Karel in t Hout Department of Mathematics and Computer cience, University of Antwerp, Middelheimlaan, 22 Antwerp, Belgium Abstract.
More information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More informationRough volatility models: When population processes become a new tool for trading and risk management
Rough volatility models: When population processes become a new tool for trading and risk management Omar El Euch and Mathieu Rosenbaum École Polytechnique 4 October 2017 Omar El Euch and Mathieu Rosenbaum
More informationOption pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard
Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard Indifference pricing and the minimal entropy martingale measure Fred Espen Benth Centre of Mathematics for Applications
More informationMulti-period mean variance asset allocation: Is it bad to win the lottery?
Multi-period mean variance asset allocation: Is it bad to win the lottery? Peter Forsyth 1 D.M. Dang 1 1 Cheriton School of Computer Science University of Waterloo Guangzhou, July 28, 2014 1 / 29 The Basic
More informationMSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013
MSc Financial Engineering 2012-13 CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL To be handed in by monday January 28, 2013 Department EMS, Birkbeck Introduction The assignment consists of Reading
More informationRichardson Extrapolation Techniques for the Pricing of American-style Options
Richardson Extrapolation Techniques for the Pricing of American-style Options June 1, 2005 Abstract Richardson Extrapolation Techniques for the Pricing of American-style Options In this paper we re-examine
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationAn Efficient Numerical Scheme for Simulation of Mean-reverting Square-root Diffusions
Journal of Numerical Mathematics and Stochastics,1 (1) : 45-55, 2009 http://www.jnmas.org/jnmas1-5.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 An Efficient Numerical Scheme for Simulation of
More informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More informationNon-semimartingales in finance
Non-semimartingales in finance Pricing and Hedging Options with Quadratic Variation Tommi Sottinen University of Vaasa 1st Northern Triangular Seminar 9-11 March 2009, Helsinki University of Technology
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationTHE OPTIMAL ASSET ALLOCATION PROBLEMFOR AN INVESTOR THROUGH UTILITY MAXIMIZATION
THE OPTIMAL ASSET ALLOCATION PROBLEMFOR AN INVESTOR THROUGH UTILITY MAXIMIZATION SILAS A. IHEDIOHA 1, BRIGHT O. OSU 2 1 Department of Mathematics, Plateau State University, Bokkos, P. M. B. 2012, Jos,
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationWeak Reflection Principle and Static Hedging of Barrier Options
Weak Reflection Principle and Static Hedging of Barrier Options Sergey Nadtochiy Department of Mathematics University of Michigan Apr 2013 Fields Quantitative Finance Seminar Fields Institute, Toronto
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationPricing Dynamic Solvency Insurance and Investment Fund Protection
Pricing Dynamic Solvency Insurance and Investment Fund Protection Hans U. Gerber and Gérard Pafumi Switzerland Abstract In the first part of the paper the surplus of a company is modelled by a Wiener process.
More informationCS 774 Project: Fall 2009 Version: November 27, 2009
CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS
More informationOvernight Index Rate: Model, calibration and simulation
Research Article Overnight Index Rate: Model, calibration and simulation Olga Yashkir and Yuri Yashkir Cogent Economics & Finance (2014), 2: 936955 Page 1 of 11 Research Article Overnight Index Rate: Model,
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2017 14 Lecture 14 November 15, 2017 Derivation of the
More informationA Study on Numerical Solution of Black-Scholes Model
Journal of Mathematical Finance, 8, 8, 37-38 http://www.scirp.org/journal/jmf ISSN Online: 6-44 ISSN Print: 6-434 A Study on Numerical Solution of Black-Scholes Model Md. Nurul Anwar,*, Laek Sazzad Andallah
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationOption Pricing Formula for Fuzzy Financial Market
Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,
More informationOption Pricing under Delay Geometric Brownian Motion with Regime Switching
Science Journal of Applied Mathematics and Statistics 2016; 4(6): 263-268 http://www.sciencepublishinggroup.com/j/sjams doi: 10.11648/j.sjams.20160406.13 ISSN: 2376-9491 (Print); ISSN: 2376-9513 (Online)
More informationPricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case
Pricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case Guang-Hua Lian Collaboration with Robert Elliott University of Adelaide Feb. 2, 2011 Robert Elliott,
More informationPDE Methods for the Maximum Drawdown
PDE Methods for the Maximum Drawdown Libor Pospisil, Jan Vecer Columbia University, Department of Statistics, New York, NY 127, USA April 1, 28 Abstract Maximum drawdown is a risk measure that plays an
More informationOption Pricing Modeling Overview
Option Pricing Modeling Overview Liuren Wu Zicklin School of Business, Baruch College Options Markets Liuren Wu (Baruch) Stochastic time changes Options Markets 1 / 11 What is the purpose of building a
More informationBasic Concepts in Mathematical Finance
Chapter 1 Basic Concepts in Mathematical Finance In this chapter, we give an overview of basic concepts in mathematical finance theory, and then explain those concepts in very simple cases, namely in the
More informationPricing Volatility Derivatives with General Risk Functions. Alejandro Balbás University Carlos III of Madrid
Pricing Volatility Derivatives with General Risk Functions Alejandro Balbás University Carlos III of Madrid alejandro.balbas@uc3m.es Content Introduction. Describing volatility derivatives. Pricing and
More informationPricing with a Smile. Bruno Dupire. Bloomberg
CP-Bruno Dupire.qxd 10/08/04 6:38 PM Page 1 11 Pricing with a Smile Bruno Dupire Bloomberg The Black Scholes model (see Black and Scholes, 1973) gives options prices as a function of volatility. If an
More informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationDynamic Hedging in a Volatile Market
Dynamic in a Volatile Market Thomas F. Coleman, Yohan Kim, Yuying Li, and Arun Verma May 27, 1999 1. Introduction In financial markets, errors in option hedging can arise from two sources. First, the option
More informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
More informationProbability in Options Pricing
Probability in Options Pricing Mark Cohen and Luke Skon Kenyon College cohenmj@kenyon.edu December 14, 2012 Mark Cohen and Luke Skon (Kenyon college) Probability Presentation December 14, 2012 1 / 16 What
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
More informationMarket Risk Analysis Volume I
Market Risk Analysis Volume I Quantitative Methods in Finance Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume I xiii xvi xvii xix xxiii
More informationInterest Rate Volatility
Interest Rate Volatility III. Working with SABR Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline Arbitrage free SABR 1 Arbitrage free
More informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More informationHedging Credit Derivatives in Intensity Based Models
Hedging Credit Derivatives in Intensity Based Models PETER CARR Head of Quantitative Financial Research, Bloomberg LP, New York Director of the Masters Program in Math Finance, Courant Institute, NYU Stanford
More informationA note on the existence of unique equivalent martingale measures in a Markovian setting
Finance Stochast. 1, 251 257 1997 c Springer-Verlag 1997 A note on the existence of unique equivalent martingale measures in a Markovian setting Tina Hviid Rydberg University of Aarhus, Department of Theoretical
More informationEstimation of Value at Risk and ruin probability for diffusion processes with jumps
Estimation of Value at Risk and ruin probability for diffusion processes with jumps Begoña Fernández Universidad Nacional Autónoma de México joint work with Laurent Denis and Ana Meda PASI, May 21 Begoña
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationSparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models
Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models Norbert Hilber Seminar of Applied Mathematics ETH Zürich Workshop on Financial Modeling with Jump Processes p. 1/18 Outline
More informationValuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005
Valuation of Volatility Derivatives Jim Gatheral Global Derivatives & Risk Management 005 Paris May 4, 005 he opinions expressed in this presentation are those of the author alone, and do not necessarily
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
More informationAn Adjusted Trinomial Lattice for Pricing Arithmetic Average Based Asian Option
American Journal of Applied Mathematics 2018; 6(2): 28-33 http://www.sciencepublishinggroup.com/j/ajam doi: 10.11648/j.ajam.20180602.11 ISSN: 2330-0043 (Print); ISSN: 2330-006X (Online) An Adjusted Trinomial
More informationarxiv: v1 [q-fin.cp] 1 Nov 2016
Essentially high-order compact schemes with application to stochastic volatility models on non-uniform grids arxiv:1611.00316v1 [q-fin.cp] 1 Nov 016 Bertram Düring Christof Heuer November, 016 Abstract
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationThe Black-Scholes PDE from Scratch
The Black-Scholes PDE from Scratch chris bemis November 27, 2006 0-0 Goal: Derive the Black-Scholes PDE To do this, we will need to: Come up with some dynamics for the stock returns Discuss Brownian motion
More informationEuropean option pricing under parameter uncertainty
European option pricing under parameter uncertainty Martin Jönsson (joint work with Samuel Cohen) University of Oxford Workshop on BSDEs, SPDEs and their Applications July 4, 2017 Introduction 2/29 Introduction
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU IFS, Chengdu, China, July 30, 2018 Peter Carr (NYU) Volatility Smiles and Yield Frowns 7/30/2018 1 / 35 Interest Rates and Volatility Practitioners and
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 19 11/20/2013. Applications of Ito calculus to finance
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.7J Fall 213 Lecture 19 11/2/213 Applications of Ito calculus to finance Content. 1. Trading strategies 2. Black-Scholes option pricing formula 1 Security
More informationA High-order Front-tracking Finite Difference Method for Pricing American Options under Jump-Diffusion Models
A High-order Front-tracking Finite Difference Method for Pricing American Options under Jump-Diffusion Models Jari Toivanen Abstract A free boundary formulation is considered for the price of American
More informationHedging with Life and General Insurance Products
Hedging with Life and General Insurance Products June 2016 2 Hedging with Life and General Insurance Products Jungmin Choi Department of Mathematics East Carolina University Abstract In this study, a hybrid
More informationCash Accumulation Strategy based on Optimal Replication of Random Claims with Ordinary Integrals
arxiv:1711.1756v1 [q-fin.mf] 6 Nov 217 Cash Accumulation Strategy based on Optimal Replication of Random Claims with Ordinary Integrals Renko Siebols This paper presents a numerical model to solve the
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
More informationAmerican Option Pricing Formula for Uncertain Financial Market
American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn
More informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
More informationAsset Pricing Models with Underlying Time-varying Lévy Processes
Asset Pricing Models with Underlying Time-varying Lévy Processes Stochastics & Computational Finance 2015 Xuecan CUI Jang SCHILTZ University of Luxembourg July 9, 2015 Xuecan CUI, Jang SCHILTZ University
More informationNUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE
Trends in Mathematics - New Series Information Center for Mathematical Sciences Volume 13, Number 1, 011, pages 1 5 NUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE YONGHOON
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 13 Lecture 13 November 15, 217 Derivation of the Black-Scholes-Merton
More informationOption Pricing for a Stochastic-Volatility Jump-Diffusion Model
Option Pricing for a Stochastic-Volatility Jump-Diffusion Model Guoqing Yan and Floyd B. Hanson Department of Mathematics, Statistics, and Computer Science University of Illinois at Chicago Conference
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationInfinite Reload Options: Pricing and Analysis
Infinite Reload Options: Pricing and Analysis A. C. Bélanger P. A. Forsyth April 27, 2006 Abstract Infinite reload options allow the user to exercise his reload right as often as he chooses during the
More informationUtility Indifference Pricing and Dynamic Programming Algorithm
Chapter 8 Utility Indifference ricing and Dynamic rogramming Algorithm In the Black-Scholes framework, we can perfectly replicate an option s payoff. However, it may not be true beyond the Black-Scholes
More informationRecent Advances in Fractional Stochastic Volatility Models
Recent Advances in Fractional Stochastic Volatility Models Alexandra Chronopoulou Industrial & Enterprise Systems Engineering University of Illinois at Urbana-Champaign IPAM National Meeting of Women in
More information