A computational approach for solving fractional integral equations based on Legendre collocation method

In this paper, we propose a new method for solving second-order fractional integral equations. By using of spectral collocation method based on Gauss–Legendre–Lobatto points, a system of algebraic equations is obtained and then, we get a solution based on Legendre polynomials. If L2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^{2}$$\end{document} norm is used, this proposed method has a favorable convergence rate, which is shown by some examples.


Introduction
Fractional calculus has gained importance and popularity during the past three decades or so, due to mainly its demonstrated applications in numerous seemingly diverse fields of science and engineering. The advantage of fractional calculus becomes apparent in modeling mechanical and electrical properties of real material, as well as in the description of many fractional fields (see for details [1,2]). In recent years, there have been many studies of fractional integral and differential equations (for instance, see [3][4][5][6][7][8][9][10][11][12][13]).
The theory of integral equations is one of the most useful mathematical tools in both pure and applied mathematics. It has enormous applications in many physical problems. The Volterra integral equations were introduced in first by Vito Volterra and then studied by Traian Lalescu in his 1908 thesis. Throughout the last decade, some researchers have paid attention to the concept of fractional integral equations as Volterra kind [14][15][16]. In this paper, we are concerned with the numerical solution of the following fractional Volterra integral equation: where I a is the fractional integral of order 0\a\1; a; b and f are continuous real functions on 0; T ½ , y is the unknown function that must be determined. We directly observe that the fractional differential equation of a order [15] D a y x ð Þ ¼ H x; y x ð Þ ð Þ; a x b; can be written immediately as the fractional Volterra integral equation of second kind as y x ð Þ ¼ y a ð Þ þ I aþ1 H x; y x ð Þ ð Þ: Based on the above description, we intend to use numerical methods for solving the form of the Volterra fractional integral equation as from Eq. (1). Existence and uniqueness of the solution of Eq. (1), in additional to some analytical properties and important inequalities, are investigated in [17]. Some stability results for fractional integral equation are found in [18]. The authors in [17] applied the product trapezoidal method based on Picard iteration for approximating solution at the mesh points. Also, Atangana and et al. [19] use the Simpson's rules for solving special case of Volterra fractional integral equations. Spectral methods are an emerging area in the field of applied science and engineering. These methods provide a computational approach that has achieved substantial popularity over the last three decades. They have been applied successfully to numerical simulations of many problems in fractional calculus [20][21][22][23][24][25]. To the best of our knowledge, fractional derivative and integral are global, i.e., they are defined by an interval over the whole interval 0; T ½ , and therefore global methods, such as spectral methods, are better suited for this equation.
We introduce a collocation method based on the Gauss-Legendre interpolation for solving Eq. (1). Inspired by the work of [26], we extend the approach to Eq. (1) and provide a rigorous convergence analysis for the Legendre spectral collocation method. We show that an approximate solution is convergent in L 2 and L 1 norms. The results obtained in this paper may be considered as a generalization of the result in [17].
The structure of this paper is as follows: In Sect. 2, some necessary definitions and mathematical tools of the fractional calculus which are required for our subsequent developments are introduced. In Sect. 3, the Legendre spectral collocation method of the fractional integral equation of second kind is obtained. After this section, we discuss about convergence analysis and then, to show efficiency of Legendre spectral collocation method, some numerical experiments are presented in Sect. 5. Some conclusions are given in Sect. 6.

Basic definitions and notations
For the concept of fractional integral, we will adopt Riemann-Liouville definition. Throughout this paper, integral operator will be Riemann-Liouville fractional integral operator.
and it is said to be in the space C m l iff f m ð Þ 2 C l ; m 2 N: The Riemann-Liouville fractional integral operator of order a ! 0; of a function f 2 C l ; l ! À 1; is defined as ð2Þ where C : ð Þ is the Gamma function where C a ð Þ ¼ Some properties of the Riemann-Liouville fractional integral are where a; b ! 0 and m [ À 1: Some properties of the Riemann-Liouville fractional integral are where a; b ! 0 and m [ À 1: Þ q 2 be a weight function in the usual sense, for q 1 ; q 2 [ À 1: The Jacobi polynomials J q 1 ;q 2 n È É 1 n¼0 are orthogonal with respect to x q 1 ;q 2 over À1; 1 ½ : The set of Jacobi polynomials is a complete L 2 x q 1 ;q 2 À1; 1 ½ -orthogonal system, namely [23] where d nm is the Kronecker function, and In particular, we find Let S N À1; 1 ½ be the set of all polynomials of degree , the nodes and the corresponding Christoffel numbers of the standard Jacobi-Gauss interpolation on K ¼ À1; 1 ð Þ; then we have For any u 2 C K ð Þ; we denote by p q 1 ;q 2 x;N : C K ð Þ ! S N the Jacobi-Gauss interpolation operator, such that p q 1 ;q 2 x;N u x q 1 ;q 2 It is clear that In special case, if q 1 ¼ q 2 ¼ 0; then the Jacobi polynomial is reduced to the Legendre polynomial. Thus, we recall that the Legendre polynomials L i x ð Þ; i ¼ 0; 1; . . . f g are defined on K with the following recurrence formula [23,24]: , D is the differential operator, the Legendre polynomial has the expansion as following The set gis a complete orthogonal system in L 2 K ð Þ; and we have where d jk is the Kronecker delta symbol and h j ¼ 2 2jþ1 . Thus, for any g 2 L 2 K ð Þ; we have For simplicity, we let Þ are nodes and Christoffel numbers of Gauss-Legendre-Lobatto interpolation on the classical interval À1; 1 ½ ; respectively. The norm and discrete inner product in L 2 K ð Þ are defined as In the next section, we use some special techniques for converting Eq. (1) to a problem that can be solved by the standard Legendre spectral collocation scheme. These topics will be presented in the next section.

The spectral collocation method
In this section, we shall propose a Legendre spectral collocation method for solving Eq. (1). The Legendre spectral collocation method takes advantage of both the Legendre polynomials and the Gauss-Legendre-Lobatto interpolation nodes (for more details, see [27]). The main idea is to use the method of discretizing the fractional integral equation to obtain a system of algebraic equations with unknown coefficients. As a result, the solution of the obtained system can be easily solved.
Again, let us to consider Eq. (1) by applying Definition 2.2 as Next, in enjambment, we let K ¼ À1; 1 ½ : For ease of analysis, we transfer the problem (1) to an equivalent problem defined in K: More specifically, we apply the change of variable t ¼ 2x Moreover, to transfer the integral interval 0; (12) can be written as then, Eq. (13) can be reduced to At last, under the following linear transformation The Legendre spectral collocation method for Eq. (17) is Now, we want to describe a numerical implementation for Eq. (18). To this end, we use the following relations Then, by using (19) and then Eq. (3), a direct computation leads to the following equation Applying relations (3) to (6), d i0 can be readily obtained as Hence, by inserting Eqs. (19) and (20) in Eq. (18), we conclude Finally, by comparing the expansion coefficients of Eq. (21), we can get By replacing the above relation, we obtain the expression of U t ð Þ accordingly.

Convergence analysis
In this section, we want to carry out the error analysis for the numerical method (18) under L 2 K ð Þ and L 1 K ð Þ: Here, we present some preparations as follow [20,23]: Let L 2 x q 1 ;q 2 K ð Þ : the Jacobi-weighted L 2 Hilbert space with the scalar product Therefore, the corresponding norm is : the Jacobi-weighted Sobolev spaces with the norm and semi-norm, respectively, H m x q 1 ;q 2 by its inner product is a Hilbert space. For the purpose of convenience, we write Lemma 4.1 [23] Assume that u 2 H m x q 1 ;q 2 with q 1 ; q 2 [ À 1; and m ! 1: Then, for any integer number k such that 0 k m N þ 1; the constant C [ 0 exists such that the following relation holds: In particular, for any u 2 H m K ð Þ with 1 m N þ 1; the above relation is reduced to the following result u À p t;N u ð Þ Theorem 4.2 For sufficiently large N, the Legendre spectral collocation approximation converges to the exact solution in L 2 -norm, i.e., e N k k 2 ¼ Y À U k k 2 ! 0: Proof Assume that U x ð Þ is obtained by using the Legendre spectral collocation method Eq. (1). Then, according to definition of p t;N in last section, we will have By using Lemma 4.1, we can get for any integer 1 m N þ 1; The rest of our proof is to show p t;N Y ð Þ À U 2 ! 0 as N ! 1: To this end, along with definition Y in Eq. (17), for N ! 1, we have where, Similarly, we can write where Before we proceed to the rest of proof, we present some relations that will be used later. By virtue of (29) and then using the standard Jacobi-Gauss quadrature formula (4), we obtain Again, in the same scheme, we have At last, by applying Lemma 4.1 and Eq. (29) together, for any integer m 1 m N þ 1 ð Þ ; we conclude that Now, let us explain the rest of proof. By subtracting Eq. (28) from (27), we derive that The above formula can be rewritten as We let Next, by using the standard Legendre-Gauss quadrature formula, we obtain By using the Cauchy-Schwartz inequality and Eq. (32), we get Similarly, by applying Eqs. (4) and (31) and the Cauchy-Schwartz inequality, we can get On the other hand, Jenson's inequality for any convex function g on interval a; b ð Þ is Thus, since then, for any given x j 2 À1; 1 ð Þ, g t ð Þ ¼ x j þ 1 À Á t is convex function and according to Jenson's inequality on interval 0; 1 ½ , we will have Therefore, by applying the previous inequality, we yield the following inequality Next, by using relation (37), Eq. (35) will be reduced to At last, along with Lemma 4.1 and relations (32), (37), the previous result yields Finally, relations (33), (34), (39) prove that the approximation is convergent in L 2 -norm. Hence, the theorem is proved.

Numerical example
To show efficiency of our numerical method, the following examples which have exact solutions are considered.
Example 5.1 Consider the following fractional integral equation [17] y t ð Þ ¼ 0:01 The corresponding exact solution is given by y t ð Þ ¼ ffiffiffi p p 1 þ t ð Þ À 3 2 : Fig. 1 presents the approximate and exact solutions which are found in very good agreement. In Fig. 2, the numerical errors are plotted in L 2 -norm for 2 N 24: Example 5.2 Our second example is the following fractional integral equation [17] The exact solution is y x ð Þ ¼ C 2 3 À Á t: We get the numerical solution of the above fractional integral equation as t; t; Fig. 4 The error function of Example 5.1 versus the number of interpolation points t þ 3:12 Â 10 À10 t 2 þ 2:39 Â 10 À7 t 3 þ 4:81 Â 10 À9 t 4 : The numerical results of proposed method can be seen from Figs. 3 and 4. These results indicate that the spectral accuracy is obtained for this problem, although the given function f is not very smooth.
These results indicate that the spectral accuracy is obtained by using Legendre collocation method for this problem. The exact solution is y t ð Þ ¼ 1 À e t erfc ffi ffi t p À Á ; such that erfc is the complementary error function defined by erfc t ð Þ ¼ 1 À erf t ð Þ ¼ 2 ffiffiffi p p Z 1 t e Às 2 ds: We list the numerical results in Table 1 to be able to test the predicted convergence rate. The comparison of the results for N ¼ 8 and N ¼ 10 with the iterative numerical method [17] is investigated in Table 1. In Fig. 5, we can observe that our numerical solutions coincide closely with the exact ones. At a glance, we can find out the results of Legendre spectral collocation method are satisfactory in the few steps.

Conclusion
In this paper, we present a spectral Legendre collocation approximation of a class of fractional integral equations of the second kind. The most important contribution of this paper is that the errors of approximations decay exponentially in L 2norm. We prove that our proposed method is effective and has high convergence rate and is better than some new methods introduced in this paper. The numerical experiments show that the proposed method is efficient and powerful numerical scheme for solving many fractional problems.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creative commons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.