The Role of High Precision Arithmetic in Calculating Numerical Laplace and Inverse Laplace Transforms
- 1 Department of Applied Mathematics, Western University, London, Ontario, Canada
- 2 Department of Applied Mathematics, Western University, London, Ontario, Canada
- 3 Department of Applied Mathematics, Western University, London, Ontario, Canada
Abstract
In order to find stable, accurate, and computationally efficient methods for performing the inverse Laplace transform, a new double transformation approach is proposed. To validate and improve the inversion solution obtained using the Gaver-Stehfest algorithm, direct Laplace transforms are taken of the numerically inverted transforms to compare with the original function. The numerical direct Laplace transform is implemented with a composite Simpson’s rule. Challenging numerical examples involving periodic and oscillatory functions, are investigated. The numerical examples illustrate the computational accuracy and efficiency of the direct Laplace transform and its inverse due to increasing the precision level and the number of terms included in the expansion. It is found that the number of expansion terms and the precision level selected must be in a harmonious balance in order for correct and stable results to be obtained.
- Cengel, Y. and Palm, W. (2013) Differential Equations for Engineers and Scientists. McGraw-Hill, New York.
- Edwards, C. and Penney, C. (2015) Differential Equations and Boundary Value Problems. 5th Edition, Pearson, London.
- Polyanin, A.D. (2016) Handbook of Linear Partial Differential Equations for Engineers and Scientists. 2nd Edition, Chapman and Hall/CRC, UK. https://doi.org/10.1201/b19056
- Cohen, A.M. (2007) Numerical Methods for Laplace Transform Inversion. Springer, Berlin.
- Kreyszig, E. (2011) Advanced Engineering Mathematics. 10th Edition, Wiley, Hoboken.
- Abate, J., Chudhuru, G. and Whitt, W. (2000) An Introduction to Numerical Transform Inversion and Its Application to Probability Models. In: Grassmann, W., Ed., Computational Probability, Springer, Berlin, 257-323.
- Badescu, A., Breuer, L., Da Silva Soares, A., Latouche, G., Remiche, M.A. and Stanford, D.A. (2005) Risk Processes Analyzed as Fluid Queues. Scandinavian Actuarial Journal, 2005, 127-141. https://doi.org/10.1080/03461230410000565
- Kao, E. (1997) An Introduction to Stochastic Processes. Duxbury Press, Pacific Grove, CA.
- Nadarajah, S. and Kotz, S. (2006) On the Laplace transform of the Pareto Distribution. Queueing Systems, 54, 243-244. https://doi.org/10.1007/s11134-006-0299-1
- Davison, M. and Doeschl, A. (2004) A Hyperbolic PDE with Parabolic Behavior. SIAM Reviev, 46, 115-127. https://doi.org/10.1137/S0036144502409007
- Cruz-Baez, D.I. and Gonazalex-Rodriguez, J.M. (2008) A Different Approach for Pricing Asian Options. Applied Mathematics Letters, 21, 303-306.
- Deakin, A. and Davison, M. (2010) An Analytic Solution for a Vasicek Interest Rate Convertible Bond Model. Journal of Applied Mathematics, 2010, Article ID: 263451. https://doi.org/10.1155/2010/263451
- Fusai, G. and Roncoroni, A. (2008) Implementing Models in Quantitative Finance: Methods and Cases. In: Avellaneda, M., Barone-Adesi, G., Broadie, M., Davis, M., Derman, E., Klüppelberg, C. and Schachermayer, W., Eds., Springer Finance, Springer-Verlag, Berlin, Heidelberg.
- Kuhlman, K.L. (2013) Review of Inverse Laplace Transform Algorithms for Laplace-Space Numerical Approaches. Numerical Algorithms, 63, 339-355. https://doi.org/10.1007/s11075-012-9625-3
- Stehfest, H. (1970) Algorithm 368: Numerical Inversion of Laplace Transform. Communications of the ACM, 13, 47-49. https://doi.org/10.1145/361953.361969