Applying the Barycentric Jacobi Spectral Method to Price Options with Transaction Costs in a Fractional Black-Scholes Framework

B. F. Nteumagné, E. Pindza, E. Maré

Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria, Republic of South Africa.

DOI: 10.4236/jmf.2014.41004                

Abstract

The aim of this paper is to show how options with transaction costs under fractional, mixed Brownian-fractional, and subdiffusive fractional Black-Scholes models can be efficiently computed by using the barycentric Jacobi spectral method. The reliability of the barycentric Jacobi spectral method for space (asset) direction discretization is demonstrated by solving partial differential equations (PDEs) arising from pricing European options with transaction costs under these models. The discretization of these PDEs in time relies on the implicit Runge-Kutta Radau IIA method. We conducted various numerical experiments and compared our numerical results with existing analytical solutions. It was found that the proposed method is efficient, highly accurate and reliable, and is an alternative to some existing numerical methods for pricing financial options.

Keywords

Jacobi Spectral Method; European Options; Fractional Black-Scholes Model; Mixed Brownian-Fractional Brownian Black-Scholes Model; Transaction Cost; Subdiffusive Fractional Black-Scholes Model; Scaling

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	F. Black and M. Scholes, “Pricing of Options and Corporate Liabilities,” Journal of Political Economy, Vol. 81, No. 3, 1973, pp. 637-654. http://dx.doi.org/10.1086/260062
[2]	J. C. Hull, “Options, Futures and Other Derivatives,” 4th Edition, Prentice-Hall, Upper Saddle River, 2000.
[3]	H. E. Leland, “Option Pricing and Replication with Transaction Costs,” Journal of Finance, Vol. 40, No. 5, 1985, pp. 141-183. http://dx.doi.org/10.1111/j.1540-6261.1985.tb02383.x
[4]	P. Boyle and T. Vorst, “Option Replication in Discrete Time with Transaction Costs,” Journal of Finance, Vol. 47, No. 1, 1992, pp. 271-293. http://dx.doi.org/10.1111/j.1540-6261.1992.tb03986.x
[5]	M. Monoyios, “Option Pricing with Transaction Costs Using a Markov Chain Approximation,” Journal of Economic Dynamics and Control, Vol. 28, No. 5, 2004, pp. 889-913. http://dx.doi.org/10.1016/S0165-1889(03)00059-9
[6]	X. T. Wang, “Scaling and Long-Range Dependence in Option Pricing I: Pricing European Option with Transaction Costs under Fractional Black-Scholes Model,” Physics A, Vol. 389, No. 3, 2010, pp. 438-444. http://dx.doi.org/10.1016/j.physa.2009.09.041
[7]	M. Magdziarz, “Black-Scholes Formula in Subdiffusive Regime,” Journal of Statistical Physics, Vol. 136, No. 3, 2009, pp. 553-564. http://dx.doi.org/10.1007/s10955-009-9791-4
[8]	J. Wang, J. R. Liang, L. J. Lv, W. Y. Qui and F, Y. Ren, “Continuous Time Black-Scholes Equation with Transaction Costs in Subdifusive Fractional Brownian Motion Regime,” Physics A, Vol. 391, No. 3, 2012, pp. 750-759. http://dx.doi.org/10.1016/j.physa.2011.09.008
[9]	H. Gu and Y.-X. Z, “Pricing Option with Transaction Costs under the Subdiffusive Black-Scholes Model,” Journal of East China Normal University, Vol. 5, 2012, pp. 85-92.
[10]	X. T. Wang, E. H. Zhu, M. M. Tang and H. G. Yan, “Scaling and Long-Range Dependence in Option Pricing II: Pricing European Option with Transaction Costs under the Mixed Brownian-Fractional Brownian Model,” Physics A, Vol. 389, No. 3, 2010, pp. 445-451. http://dx.doi.org/10.1016/j.physa.2009.09.043
[11]	R. Baltensperger, J. P. Berrut and B. Noel, “Exponential Convergence of a Linear Rational Interpolant between Transformed Chebyshev Points,” Mathematics of Computation, Vol. 68, 1999, pp. 1109-1120. http://dx.doi.org/10.1090/S0025-5718-99-01070-4
[12]	E. Pindza and K. C. Patidar, “Spectral Method for Pricing Options in Illiquid Markets,” Proceedings of the International Conference on Numerical Analysis and Applied Mathematics, Vol. 1479, 2012, pp. 1403-1406.
[13]	E. Ngounda, K. C. Patidar and E. Pindza, “Contour Integral Method for European Options with Jumps,” Communications in Nonlinear Science and Numerical Simulation, Vol. 18, No. 3, 2013, pp. 478-492. http://dx.doi.org/10.1016/j.cnsns.2012.08.003
[14]	E. Ngounda, K. C. Patidar and E. Pindza, “A Robust Spectral Method for Solving Heston’s Volatility Model,” Journal of Optimization Theory and Applications, Springer US, 2013.
[15]	H. Wang and D. Huybrechs, “Explicit Barycentric Weights for Polynomial Interpolation in the Roots of Extrema of Classical Orthogonal Polynomials,” Report TW604, Department of Computer Science, Katholieke Universiteit Leuven, 2012.
[16]	E. Hairer and G. Wanner, “Solving Ordinary Differential Equations, II: Stiff and Differential-Algebraic Problems,” PringerVerlag, New York, 1991. http://dx.doi.org/10.1007/978-3-662-09947-6
[17]	J. L. Lagrange, “Lecons Elémentaires sur les Mathématiques, Données à l’Ecole Normale en 1795,” Journal de l’école Polytechnique, Vol. 7, 1795, pp. 183-287.
[18]	P. Henrici, “Essentials of Numerical Analysis with Pocket Calculator Demonstrations,” Wiley, New York, 1982.
[19]	W. Werner, “Polynomial Interpolation: Lagrange versus Newton,” Mathematics of Computation, Vol. 43, 1984, pp. 205-217. http://dx.doi.org/10.1090/S0025-5718-1984-0744931-0
[20]	J. P. Berrut and L. N. Trefethen, “Barycentric Lagrange Interpolation,” SIAM Review, Vol. 46, 2004, pp. 501-517. http://dx.doi.org/10.1137/S0036144502417715
[21]	T. W. Tee and L. N. Trefethen, “A Rational Spectral Collocation Method with Adaptively Transformed Chebyshev Grids Points,” Journal of Scientific Computing, Vol. 28, No. 5, 2006, pp. 1798-1811.