ON THE NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING TRANSFORMS AND QUADRATURE

Abstract

In this work, we extended the work of [12] to approximate the solution of fractional order dierential equations by an integral representation in the complex plane. The resultant integral is approximated to high order accuracy using quadrature. The accuracy of the method depends on the selection of optimal contour of integration. Several contour have been proposed in the literature for solving fractional dierential equations. In the present work, we will investigate the applicability of the recently developed optimal contour in [16] for solving fractional dierential equations. Various fractional order dierential equations are approximated and the results are compared with other methods to demonstrate the eciency and accuracy of the method for various optimal contour of integrations.

Keywords

• Fractional order ODE,Laplace Transform,Quadrature.

References

1. Daul, L., Klein, P., and Kempfle, S., (1991) Damping description involving fractional operators, Mech. systems Signal processing, 5, 81–88.
2. Diethelm, K., (1997), An algorithm for the numerical solution of differential equations of fractional order, Electronic Transactions on Numerical Analysis, 5, 1–6.
3. Diethelm, K. and Ford, N. J., (2002), Analysis of fractional differential equations, Journal of Mathe- matical Analysis and Applications, 265, 229–248.
4. Diethelm, K., (2015), Increasing the efficiency of shooting methods for terminal value problems of fractional order, Journal of Computational Physics, 293, 135–141.
5. Ford, N. J., Morgado, M. L., and Rebelo, M., High order numerical methods for fractional terminal value problems, Computational Methods in Applied Mathematics, 14, 55–70. Glockle, W. G., and Nonnenmacher, T. F., (1995), A fractional calculas approach to self-similar protein dynamics, Biophys J., 68, 46–53.
6. McLean, M., Thomee, V., (2007) Numerical solution via laplace transforms of a fractional order evolution equation, Journal of Integral Equations and Applications, 22, 57–94.
7. Metzler, R., Schick, W., Kilian, H. G., and Nonnenmacher, T. F., (1995) Relaxation in filled ploymers:A fractional calculus approach, J. Chem. Phys. 103, 7180–7186.
8. Oldham, K. B., and Spanier, J., (1974), The fractional calculas, Mathematics in Science and Engi- neering, Vol, 111, Academic Press, New York/London.

9. Rizzardi, M., (1995), A modification of Talbot’s method for the simultaneous approximation of several values of the inverse Laplace transform, ACM Trans. Math. Software, 21, 347–371.
10. Schwarz, H. R., (1989), Numerical Analysis, John Wiley and Sons Ltd., Chichester.
11. Sheen, D., Sloan, I. H., and Thomee, V., (2003), A parallel method for time discretization of parabolic equations based on Laplace transformations and quadrature, IMA J. Numer. Anal., 23, 269–299.
12. Talbot, A., (1079), The accurate numerical inversion of Laplace transform, J. Inst. Math. Appl., 23, –120.
13. Uddin, M., and Ahmad, Suleman., (2017), On the numerical solution of Bagley-Torvik equation via the Laplace transform, Tbilisi Mathematical Journal 10(1), 279-284.
14. Uddin, M., and Haq, S. (2011), RBFs approximation method for time fractional partial differential equations, Communications in Nonlinear Science and Numerical Simulation, 16(11), 4208–4214.
15. Weideman, J. A. C., and Trefethen, L. N., (2007) Parabolic and hyperbolic contours for computing the Bromwich integral, Mathematics of Computation, 76, 1341–1356.