A Spectral Method Based on Hahn Polynomials for Numerical Solution of Fractional Integro-Differential Equations with Weakly Singular Kernel

Authors
Farideh Salehi
Habibollah Saeedi
Mahmoud Mohseni Moghadam
10.52547/mmr.6.1.65
Abstract
Introduction

Despite wide applications of constant order fractional derivatives, some systems require the use of derivatives whose order changes with respect to other parameters. Samko and Ross produced an extension of the classical fractional calculus with a continuously varying order for differential and integral operators. Variable-order fractional (V-OF) calculus has applications in optimal control, processing of geographical data, diffusion processes, description of anomalous diffusion, heat-transfer problems, etc. Due to the V-OF operators which are non-local with singular kernels, finding the exact solutions of V-OF problems is difficult. Therefore, efficient numerical techniques are necessary to be developed. The numerical solution of V-OF differential equation has been considered in some papers.

Recently, discrete orthogonal polynomials have been considered as basis functions instead of continuous orthogonal polynomials. Discrete orthogonal polynomials are orthogonal with respect to a weighted discrete inner product. These polynomials have important applications in chemical engineering, theory of random matrices, queuing theory and image coding. In this paper, we focus on a special class of discrete polynomials, called Hahn polynomials.

In this work, first, a new operational matrix is obtained for V-OF integral of Hahn polynomials. Then, we use a spectral collocation technique combined with the associated operational matrices of V-OF integral for solving weakly singular fractional integro-differential equations.

Material and methods

In this scheme, the operational matrix of fractional integration of Hahn polynomials is calculated. This method converts the weakly singular fractional integro-differential equations into an algebraic system which can be solved by a technique of linear algebra.

Results and discussion

In this paper, some numerical examples are provided to show the accuracy and efficiency of the presented method. By using a small number of Hahn polynomials, significant results are achieved which are compared to other methods. A comparison to the numerical solutions by CAS and Haar wavelets and Adomain decomposition method, shows that this technique is accurate enough to be known as a powerful device.

Conclusion

The following results are obtained from this research.

The operational matrix of fractional integration of Hahn polynomials is presented for the first time.
The main advantage of approximating a continuous function by Hahn polynomials is that they have a spectral accuracy at interval [0,N], where N is the number of bases.
Furthermore, for estimating the coefficients of the expansion of approximate solution, we only have to compute a summation which is calculated exactly.
Using Hahn polynomials, the numerical results achieved only by a small number of bases, are accurate in a larger interval and significant results are achieved../files/site1/files/61/7.pdf
Keywords
Weakly Singular Fractional Integro-Differential Equations
Hahn Polynomials
Operational matrix
Spectral method
1. Magin R., "Fractional Calculus in Bioengineering", Begell House Publishers (2006.( 2. S. A., L. J. L., P. E. M. J. A., "Fractional order calculus: historical apologia, basic concepts and some applications", Rev. Bras. Ensino Fs. ( 2011) 4302-4302.## 3. Jiang Y., Ma J., "Higher order finite element methods for time fractional partial differential equations", J. Comput. Appl. Math, (2011) 3285-3290. ## 4. Bhrawy A., Zaky M., Gorder R. V., "A space-time Legendre spectral tau method for the two-sided space-time Caputo fractional diffusion-wave equation", Numer. Algor (2016)151-180. ## 5. Ray S., "On Haar wavelet operational matrix of general order and its application for the numerical solution of fractional Bagley-Torvik equation", Appl. Math. Comput (2012) 5239-5248. ## 6. Mohseni Moghadam S. H. M., Mollahasani N., Chuev G., "A CAS wavelet method for solving nonlinear Fredholm integro-differential equations of fractional order", Commun. Nonlinear Sci. Numer. Simulat., (2011) 1154-1163. ## 7. R. P., O. Y., B. E., "Fractional-order Bernoulli wavelets and their applications", Appl. Math. Model., vol. 40, no. 17 (2016) 8087-8107. ## 8. Rahimkhani P., Ordokhani Y., Babolian E., "A new operational matrix based on Bernoulli wavelets for solving fractional delay differential equations", Numer. Algor., Vol. 74, No. 1, (2017) 223-245. ## 9. Rahimkhani P., Ordokhani Y., Babolian E., "Fractional-order Bernoulli functions and their applications in solving fractional Fredholem–Volterra integro-differential equations", Appl. Numer. Math., vol. 122 (2017) 66-81. ## 10. T. B. Q., L. X. F, "Solution of a class of Volterra integral equations with singular and weakly singular kernels", Appl. Math. Comput (2008) 406-413. ## 11. Zozulya V., Gonzalez-chi P., "Weakly singular, singular and hypersingular integrals in 3-Delas-ticity and fracture mechanics", J. Chin. Inst. Eng., vol. 22 (1999) 763-775. ## 12. K. P., Puri P., "Computational Methods for Linear Integral Equations", Boston: Birkhauser, (2002 (## 13. Niu H., Hao L., Yu Z., Yin J., "Adomian decomposition method for solving nonlinear integro-differential equations of fractional order", J. Liaoning Tech. Univ. Nat. Sci. ( 2013) 132-135. ## 14. M. H. J. Yi, "CAS wavelet method for solving the fractional integro-differential equation with a weakly singular kernel", Int. J. Comput. Math. (2015) 1715-1728. ## 15. Y. a. Z. L. Wang, "SCW method for solving the fractional integro-differential equations with a weakly singular kernel",,Appl. Math. Comput. (2016) 72-80. ## 16. Xiu G., Karniadakis G., "The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations", SIAM J. Sci. Comput (2002) 619-644. ## 17. Glaubitz J., Öffner P., Sonar T., "Application of modal filtering to a spectral difference method", Math. Comput (2016) 1-32. ## 18. McGregor K. S. J., "The Hahn polynomials, formulas and an application", Scripta Math. (1961) 33-46. ## 19. Koekoek R., Lesky P. , Swarttouw R., "Hypergeometric orthogonal polynomials and their q-analogues", Springer Science & Business Media (2010). ## 20. Öffner G. R., P., "On Hahn polynomial expansion of a continuous function of bounded variation", arXiv:1610.06748v1 (2016)1-26. ## 21. Kreyszig E., "Introductory Functional Analysis with Applications", New York: John Wiley and Sons Press (1978). ## 22. Abramowitz M., Stegun I., "Handbook of Mathematical Functions", New York: Dover Publications (1972). ## 23. Mokhtary P., "Numerical analysis of an operational Jacobi Tau method for fractional weakly singular integro-differential equations", Appl. Numer. Math., vol. 121 (2017) 52-67. ##