Research Article | Open Access

Volume 2017 |Article ID 6917828 | https://doi.org/10.1155/2017/6917828

Zhiqiang Zhou, Xuemei Gao, "Laplace Transform Methods for a Free Boundary Problem of Time-Fractional Partial Differential Equation System", Discrete Dynamics in Nature and Society, vol. 2017, Article ID 6917828, 9 pages, 2017. https://doi.org/10.1155/2017/6917828

# Laplace Transform Methods for a Free Boundary Problem of Time-Fractional Partial Differential Equation System

Revised24 May 2017
Accepted08 Jun 2017
Published10 Jul 2017

#### Abstract

We study the pricing of the American options with fractal transmission system under two-state regime switching models. This pricing problem can be formulated as a free boundary problem of time-fractional partial differential equation (FPDE) system. Firstly, applying Laplace transform to the governing FPDEs with respect to the time variable results in second-order ordinary differential equations (ODEs) with two free boundaries. Then, the solutions of ODEs are expressed in an explicit form. Consequently the early exercise boundaries and the values for the American option are recovered using the Gaver-Stehfest formula. Numerical comparisons of the methods with the finite difference methods are carried out to verify the efficiency of the methods.

#### 1. Introduction

Markov regime switching models were first introduced by Hamilton [1] and recently have become popular in financial applications including equity options [217], bond prices and interest rate derivatives [1820], portfolio selection [21], and trading rules [2226]. The Markov regime switching models allow the model parameters (drift and volatility coefficients) to depend on a Markov chain which can reflect the information of the market environments and at the same time preserve the simplicity of the models. However when the model parameters are governed by the Markov chain, the valuation of the options becomes complex.

American option pricing is a kind of classical free boundary problems, where the free boundary refers to early exercise boundary or optimal exercise boundary. Laplace transform methods have been developed to solve the free boundary problems arising in American option pricing under geometric Brownian motion (GBM) (see Mallier and Alobaidi [27], Zhu [28], and Zhu and Zhang [29]) and constant elasticity of variance (CEV) (see Wong and Zhao [30]). The essential idea of the Laplace transform methods for solving the American option pricing problems is described as follows. Applying Laplace transform to the governing free boundary partial differential equations (PDEs) with respect to the time variable results in a boundary value problem of second-order ordinary differential equations (ODEs). Unlike the case for the European option pricing, the ODE involves an unknown boundary (in Laplace space) which needs to be solved. If the solution of the ODE with the unknown boundary can be expressed in an explicit form, then a nonlinear algebraic equation for the unknown boundary can be derived. Using a simple solver for the nonlinear algebraic equation, we can get the value of the unknown boundary in Laplace space. Consequently the early exercise boundary and the value for the American option can be achieved using the inverse Laplace transform.

Chen et al. [31] studied a predictor-corrector finite difference methods for pricing American options under the FMLS model which is a kind of space-fractional derivative model. Liang et al. [32] used Fourier transform to solve a bifractional Black-Scholes-Merton differential equation. In this paper, we study the pricing of American options with time-fractional model which has essential difference to the space-fractional model. Following the model in [33], we assume that the underlying asset price still follows the classical Brownian motion, but the change in the option price is considered as a fractal transmission system. The price of such American option follows time-fractional partial differential equations (PDEs) with free boundaries. The solution of the time-fractional PDEs with two free boundaries is more challenging than solving the fractional PDEs with fixed boundary for European option pricing in [33] and much more complex than solving single Black-Scholes-Merton differential equation in [32].

In this paper, we develop the Laplace transform methods for pricing time-fractional American options under regime switching models. The fundamental difference to the GBM and CEV American option pricing is that the governing free boundary PDEs for the regime switching American option pricing are a coupled system. Therefore the Laplace transform in time for the system of free boundary PDEs leads to a coupled system of ODEs whose analytic solutions take much trouble to the approach. A careful derivation in different cases makes the approach successful. Numerical examples are provided to verify the efficiency of the approach. The generated option value and early exercise boundary are compared with the finite difference method which are usually used as the benchmark methods in the area of option pricing.

#### 2. Laplace Transform Methods for American Option Pricing

Let the underlying asset prices follow a two-state regime switching model under risk-neutral measure:where is a standard Brownian motion and is a continuous-time Markov chain with two states . Assume that, at each state , the interest rates and volatility for are nonnegative constants. Letbe the generator matrix of the Markov chain process with . Let be the option expiration date, be the time to maturity, represent the American put option prices with stock price , and be the corresponding optimal exercise boundaries at time . Then the American put options satisfy the following system of PDEs with free boundaries, for :where are the differential operators defined byHowever, it is argued that the time derivative should be replaced by the fractional derivative under the assumption that the change in the option price follows a fractal transmission system (see, e.g., [3234]). So, the time-fractional American put options satisfy the following system of PDEs:with the free boundary conditions and the initial condition:where is defined as the Caputo fractional derivative for :Condition (10) is called “high-contact condition” or “smooth-pasting condition.”

For , we define the Laplace-Carson transform (LCT) of the American put option price aswhich is essentially the Laplace transform (LT) asUsing the Laplace transform formula for the Caputo fractional derivative (see (2.253) in [35]),and the relationship (13), the LCT for is found asThe use of LCT simplifies notations in the later analysis and makes inverse Laplace solutions have a unique singularity . We note that for the case formula (15) is still true. Denoting and taking LCTs to (3)–(10), we obtain the following ODEs, for :The existence of Laplace transform and is guaranteed by the continuity and boundedness of and . To solve (16)–(20), we use the following variable and function transformations:Let . Then the calculations yield an equivalent form to (16)–(20) as follows:Compared to (16), the coefficients of the first and second derivatives in (22) become constant which will simplify the calculation. However, (22) includes two unknown variables and a nonsmooth term , which will make the solutions somewhat complicated.

From the definition of and monotone decreasing property of (see, e.g., [10]), we haveand thereforeMoreover and its first-order derivative are continuous at ( means that ). Without loss of generality, we assume . So . The solutions of (22)–(26) can be solved separately on four intervals , and .

Proposition 1. The solutions of (22)–(26) on interval are given by

Proof. The proof is straightforward by noticing (23).

The solutions of ODEs on the other intervals can be written as a general solution to the corresponding homogeneous part plus a particular one satisfying the nonhomogeneous system. More precisely, we discuss the solutions on the different intervals as follows.

Proposition 2. On interval , the solutions of (22)–(26) are given byand constant is uniquely determined by the boundary conditions (25) and (26); that is,Here the unspecified notations are defined in (34)–(37) below.

Proof. For , the solution and then satisfies the following equation:By solving ODE (32), we get the general solutionwhere is a particular solution and and are the fundamental increasing and decreasing solutions with explicit forms:Note that and . Denotingand using (23), we complete the proof.

Now we derive the solutions of (22)–(26) on interval . Let withThen (22) is equivalent to the following first-order ODE system:where

Lemma 3. For any real number , the homogeneous equation has real valued fundamental solutionswhere and are different negative real numbers, and are different positive real numbers, and are real valued column vectors.

Proof. It is well-known that has fundamental solution matrix as (41), with and representing characteristic values and corresponding characteristic vectors. So we only need to prove all are real numbers and different from each other. The characteristic polynomial of (39) isTherefore,whereSince and , for any we know that all are real numbers and and . Moreover, we haveThen there must exist four real roots solving characteristic equation , and we can list the ranges of these values:Note that the roots can be searched by the secant method on their corresponding ranges as above.

Proposition 4. Let be a submatrix extracting the first and the third rows from . Then the solutions of (22)–(26) on interval are given bywhere is given by (41) and by (52) below.

Proof. It can be verified that (22) has a particular solutionHere, for . Denote . Then the solution to (39) can be written as Using the continuity and smoothing-pasting of solutions at , we know that satisfiesSubstituting into the above equation and denotingwe obtain thatThus the proof is complete.

Proposition 5. Let be the submatrix extracting the first and the second columns from , where is given in Proposition 4. Then the solutions of (22)–(26) on interval are given bywith given by (56).

Proof. Since for , (22) becomes homogeneous. Furthermore using condition (see (24)), we can write the solution to the homogeneous part of (24) (or homogeneous part of (39)) aswhere is the submatrix containing the first and the second columns of and . Notations , and are , , and matrices, respectively. Since the first-order derivative of is also continuous at , by evaluating the first-order derivatives of (53) and (47) at , we obtain thatwhere and . Therefore, we obtain

Since is continuous at , by evaluating (53) and (47) at , we haveSubstituting (52) and (56) into (57) and simplifying the calculation, we obtain a system of nonlinear algebra equations with respect to and where constant matrices and are defined byWe summarize the above discussions in the following theorem.

Theorem 6. For any fixed , the optimal exercise boundaries and can be numerically computed from nonlinear equations (58) and the solutions of (22)–(26) are given by where , and are given by formulas (31), (52), and (56), respectively; the fundamental solution matrices , , and are defined in (37) and Propositions 4 and 5, respectively; the particular solutions and are given by (34) and (48), respectively.

Finally, the original optimal exercise boundaries and American put option prices can be expressed in terms of the Laplace inversion:

The numerical approximation for Laplace inversion can be computed using the most powerful Gaver-Stehfest algorithm (see, e.g., [36])withhere must be taken as even number.

#### 3. Numerical Examples

In this section, the Laplace transform method (LTM) is compared with the finite difference method (FDM). We use symbolic manipulation and multiprecision computing provided by Mathematica 9.0.

Table 1 lists the computational values of American put options under regime switching. Columns entitled “LTM" and “FDM” represent the option values obtained by Laplace transform method and finite different method, respectively. For the LTM the Gaver-Stehfest formula is applied with the number and for the FDM the number of time nodes and the number of space nodes . The relative error (RE) is the absolute value of the difference between the computed LTM option values and the FDM values divided by the FDM option values. In the implementation of the LTM, the work precision is set by 30 significant digits, while the numerical values are presented with 5 significant digits.

 Regime 1 Regime 2 LTM FDM RE (%) LTM FDM RE (%) 3.00 6.0000 6.0000 0.0000 6.0000 6.0000 0.0000 4.50 4.4924 4.5431 1.1160 4.5278 4.5117 0.3568 6.00 3.3692 3.4138 1.3065 3.3071 3.3502 1.2865 7.50 2.5536 2.5835 1.1573 2.4741 2.5026 1.1388 9.00 1.9512 1.9712 1.0146 1.8630 1.8816 0.9885 10.50 1.5047 1.5177 0.8566 1.4146 1.4265 0.8342 12.00 1.1717 1.1795 0.6613 1.0847 1.0915 0.6230 3.00 6.0000 6.0000 0.0000 6.0000 6.0000 0.0000 4.50 4.5084 4.5436 0.7747 4.5181 4.5149 0.0709 6.00 3.3787 3.3988 0.5914 3.3229 3.3420 0.5715 7.50 2.5258 2.5381 0.4846 2.4442 2.4537 0.3872 9.00 1.9016 1.9089 0.3824 1.7972 1.8037 0.3604 10.50 1.4633 1.4672 0.2658 1.3659 1.3691 0.2337 12.00 1.1517 1.1531 0.1214 1.0668 1.0678 0.0937 3.00 6.0000 6.0000 0.0000 6.0000 6.0000 0.0000 4.50 4.5275 4.5370 0.2094 4.5067 4.5112 0.0998 6.00 3.3630 3.3677 0.1396 3.3102 3.3146 0.1327 7.50 2.4738 2.4749 0.0444 2.3871 2.3877 0.0251 9.00 1.8271 1.8257 0.0767 1.7102 1.7085 0.0995 10.50 1.3915 1.3883 0.2305 1.2882 1.2849 0.2568 12.00 1.0940 1.0891 0.4499 1.0094 1.0046 0.4778 CPU time 1.98 s for LTM; 532.47 s for FDM

The relative errors shown in Table 1 are less than for all the cases, which illustrates that the LTM is competitively accurate for solving American option pricing problems. Since there are no exact solutions of FPDEs, we cannot judge which one of the two results is more accurate. In general, the relative error less than five percent is acceptable in financial engineering practice.

From the last row of Table 1, we see that the CPU time of LTM is much less than that of FDM. The computational cost of LTM is uniform at every time and the CPU consumption mainly comes from solving nonlinear algebraic equations with two variables and (see expression (58)). FDM is a time-stepping scheme which is needed to solve nonlinear systems at time or . Moreover, for the global dependence of the fractional differential operator, FDM needs to do some additional computation. We observe that the numbers of the quadrature nodes are much smaller than the numbers of time steps. In other words, to achieve the similar accuracy, there are much fewer nonlinear systems needed to be solved by the LTM than that by the FDM.

Figure 1 plots the early exercise boundaries obtained by LTM and FDM with the same values of parameters as in Table 1. It can be seen that the boundaries computed by FDM and LTM are pretty close. Again it shows that the LTM is robust in generating the early exercise boundaries.

#### 4. Conclusions

In this paper, we have developed Laplace transform methods to solve the time-fractional American option pricing under regime switching models. The value of the American option with regime switching is formulated as the solution to a free boundary problem of time-fractional partial differential equation system. The Laplace transform is executed for the time variables and the resulting system PDEs are solved analytically. Consequently a system of nonlinear algebraic equations for the free boundaries is obtained and solved using secant methods. Finally numerical Laplace inversion is applied to recover the early exercise boundaries and the option values. Comparisons between the LTM and the benchmark FDM are made via numerical examples, which shows that the LTM is efficient for pricing American options with regime switching. However, the LTM is still challenging for more complex models like the regime switching models with state-dependent jump diffusions. This will be left for the future studies.

#### Nomenclature of Notations

 : Stock price : Interest rate : Volatility : Free boundary : Strike price : Generator matrix of the Markov chain process : Order of the Caputo derivative : Option price.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### References

1. J. D. Hamilton, “A new approach to the economic analysis of nonstationary time series and the business cycle,” Econometrica, vol. 57, no. 2, pp. 357–384, 1989. View at: Publisher Site | Google Scholar | MathSciNet
2. D. D. Aingworth, S. R. Das, and R. Motwani, “A simple approach for pricing equity options with Markov switching state variables,” Quantitative Finance, vol. 6, no. 2, pp. 95–105, 2006. View at: Publisher Site | Google Scholar | MathSciNet
3. N. P. B. Bollen, “Valuing options in regime-switching models,” The Journal of Derivatives, vol. 6, no. 1, pp. 38–49, 1998. View at: Publisher Site | Google Scholar
4. P. P. Boyle and T. Draviam, “Pricing exotic options under regime switching,” Insurance: Mathematics and Economics, vol. 40, no. 2, pp. 267–282, 2007. View at: Publisher Site | Google Scholar | MathSciNet
5. J. Buffington and R. J. Elliott, “American options with regime switching,” International Journal of Theoretical and Applied Finance, vol. 5, no. 5, pp. 497–514, 2002. View at: Publisher Site | Google Scholar | MathSciNet
6. P. Eloe, R. H. Liu, and J. Y. Sun, “Double barrier option under regime-switching exponential mean-reverting process,” International Journal of Computer Mathematics, vol. 86, no. 6, pp. 964–981, 2009. View at: Publisher Site | Google Scholar | MathSciNet
7. X. Guo, “Information and option pricings,” Quantitative Finance, vol. 1, no. 1, pp. 38–44, 2001. View at: Publisher Site | Google Scholar | MathSciNet
8. X. Guo and Q. Zhang, “Closed-form solutions for perpetual American put options with regime switching,” SIAM Journal on Applied Mathematics, vol. 64, no. 6, pp. 2034–2049, 2004. View at: Publisher Site | Google Scholar | MathSciNet
9. M. R. Hardy, “A regime-switching model of long-term stock returns,” North American Actuarial Journal, vol. 5, no. 2, pp. 41–53, 2001. View at: Publisher Site | Google Scholar | MathSciNet
10. A. Q. Khaliq and R. H. Liu, “New numerical scheme for pricing American option with regime-switching,” International Journal of Theoretical and Applied Finance, vol. 12, no. 3, pp. 319–340, 2009. View at: Publisher Site | Google Scholar | MathSciNet
11. A. Q. Khaliq, B. Kleefeld, and R. H. Liu, “Solving complex PDE systems for pricing American options with regime-switching by efficient exponential time differencing schemes,” Numerical Methods for Partial Differential Equations, vol. 29, no. 1, pp. 320–336, 2013. View at: Publisher Site | Google Scholar | MathSciNet
12. R. H. Liu, “Regime-switching recombining tree for option pricing,” International Journal of Theoretical and Applied Finance, vol. 13, no. 3, pp. 479–499, 2010. View at: Publisher Site | Google Scholar | MathSciNet
13. R. H. Liu and J. L. Zhao, “A lattice method for option pricing with two underlying assets in the regime-switching model,” Journal of Computational and Applied Mathematics, vol. 250, pp. 96–106, 2013. View at: Publisher Site | Google Scholar | MathSciNet
14. D. D. Yao, Q. Zhang, and X. Y. Zhou, “A regime-switching model for European options,” in Stochastic processes, optimization, and control theory: applications in financial engineering, queueing networks, and manufacturing systems, H. M. Yan, G. Yin, and Q. Zhang, Eds., pp. 281–300, Springer, New York, NY, USA, 2006. View at: Publisher Site | Google Scholar | MathSciNet
15. F. L. Yuen and H. Yang, “Option pricing with regime switching by trinomial tree method,” Journal of Computational and Applied Mathematics, vol. 233, no. 8, pp. 1821–1833, 2010. View at: Publisher Site | Google Scholar | MathSciNet
16. F. L. Yuen and H. Yang, “Pricing Asian options and equity-indexed annuities with regime switching by the trinomial tree method,” North American Actuarial Journal, vol. 14, no. 2, pp. 256–277, 2010. View at: Publisher Site | Google Scholar | MathSciNet
17. J. Ma and T. Zhu, “Convergence rates of trinomial tree methods for option pricing under regime-switching models,” Applied Mathematics Letters. An International Journal of Rapid Publication, vol. 39, pp. 13–18, 2015. View at: Publisher Site | Google Scholar | MathSciNet
18. R. Bansal and H. Zhou, “Term structure of interest rates with regime shifts,” Journal of Finance, vol. 57, no. 5, pp. 1997–2043, 2002. View at: Publisher Site | Google Scholar
19. C. Landén, “Bond pricing in a hidden Markov model of the short rate,” Finance and Stochastics, vol. 4, no. 4, pp. 371–389, 2000. View at: Publisher Site | Google Scholar | MathSciNet
20. R. H. Liu, “A new tree method for pricing financial derivatives in a regime-switching mean-reverting model,” Nonlinear Analysis. Real World Applications. An International Multidisciplinary Journal, vol. 13, no. 6, pp. 2609–2621, 2012. View at: Publisher Site | Google Scholar | MathSciNet
21. X. Y. Zhou and G. Yin, “Markowitz's mean-variance portfolio selection with regime switching: a continuous-time model,” SIAM Journal on Control and Optimization, vol. 42, no. 4, pp. 1466–1482, 2003. View at: Publisher Site | Google Scholar | MathSciNet
22. P. Eloe, R. H. Liu, M. Yatsuki, G. Yin, and Q. Zhang, “Optimal selling rules in a regime-switching exponential Gaussian diffusion model,” SIAM Journal on Applied Mathematics, vol. 69, no. 3, pp. 810–829, 2008. View at: Publisher Site | Google Scholar | MathSciNet
23. G. Yin, R. H. Liu, and Q. Zhang, “Recursive algorithms for stock liquidation: a stochastic optimization approach,” SIAM Journal on Optimization, vol. 13, no. 1, pp. 240–263, 2002. View at: Publisher Site | Google Scholar | MathSciNet
24. G. Yin, Q. Zhang, F. Liu, R. H. Liu, and Y. Cheng, “Stock liquidation via stochastic approximation using NASDAQ daily and intra-day data,” Mathematical Finance, vol. 16, no. 1, pp. 217–236, 2006. View at: Publisher Site | Google Scholar | MathSciNet
25. Q. Zhang, “Stock trading: an optimal selling rule,” SIAM Journal on Control and Optimization, vol. 40, no. 1, pp. 64–87, 2001. View at: Publisher Site | Google Scholar | MathSciNet
26. Q. Zhang, G. Yin, and R. H. Liu, “A near-optimal selling rule for a two-time-scale market model,” Multiscale Modeling & Simulation. A SIAM Interdisciplinary Journal, vol. 4, no. 1, pp. 172–193, 2005. View at: Publisher Site | Google Scholar | MathSciNet
27. R. Mallier and G. Alobaidi, “Laplace transforms and American options,” Applied Mathematical Finance, vol. 7, no. 4, pp. 241–256, 2000. View at: Publisher Site | Google Scholar
28. S.-P. Zhu, “A new analytical approximation formula for the optimal exercise boundary of American put options,” International Journal of Theoretical and Applied Finance, vol. 9, no. 7, pp. 1141–1177, 2006. View at: Publisher Site | Google Scholar | MathSciNet
29. S.-P. Zhu and J. Zhang, “Using Laplace transform to price American puts,” Dynamics of Continuous, Discrete & Impulsive Systems. Series B. Applications & Algorithms, vol. 19, no. 4-5, pp. 447–469, 2012. View at: Google Scholar | MathSciNet
30. H. Y. Wong and J. Zhao, “Valuing American options under the {CEV} model by Laplace-Carson transforms,” Operations Research Letters, vol. 38, no. 5, pp. 474–481, 2010. View at: Publisher Site | Google Scholar | MathSciNet
31. W. Chen, X. Xu, and S.-P. Zhu, “A predictor-corrector approach for pricing American options under the finite moment log-stable model,” Applied Numerical Mathematics. An IMACS Journal, vol. 97, pp. 15–29, 2015. View at: Publisher Site | Google Scholar | MathSciNet
32. J.-R. Liang, J. Wang, W.-J. Zhang, W.-Y. Qiu, and F.-Y. Ren, “The solutions to a bi-fractional Black-Scholes-Merton differential equation,” International Journal of Pure and Applied Mathematics, vol. 58, no. 1, pp. 99–112, 2010. View at: Google Scholar | MathSciNet
33. W. Chen, X. Xu, and S.-P. Zhu, “Analytically pricing double barrier options based on a time-fractional Black-Scholes equation,” Computers & Mathematics with Applications. An International Journal, vol. 69, no. 12, pp. 1407–1419, 2015. View at: Publisher Site | Google Scholar | MathSciNet
34. Z. Zhou and X. Gao, “Numerical methods for pricing American options with time-fractional PDE models,” Mathematical Problems in Engineering, vol. 2016, Article ID 5614950, 8 pages, 2016. View at: Publisher Site | Google Scholar | MathSciNet
35. I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999. View at: MathSciNet
36. P. P. Valkó and J. Abate, “Comparison of sequence accelerators for the Gaver method of numerical Laplace transform inversion,” Computers & Mathematics with Applications, vol. 48, no. 3-4, pp. 629–636, 2004. View at: Publisher Site | Google Scholar | MathSciNet

Copyright © 2017 Zhiqiang Zhou and Xuemei Gao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.