The Legendre Spectral-Collocation method for a class of ... - arXiv

0 downloads 0 Views 502KB Size Report
[1] A.A. Kilbas, H.M. Srivastava and J.J. Trujillo, Theory and Applications of ... [2] I. Podlubny, Fractional Differential Equations, Academic Press, San Diego,Β ...
The Legendre Spectral-Collocation method for a class of fractional integral equations A. Yousefi, S. Javadi, E. Babolian Department of Computer Science, Faculty of Mathematical Sciences and Computer, Kharazmi University, Tehran, Iran [email protected], [email protected], [email protected]

Abstract In this paper, we consider spectral-collocation method base on Legendre-Gauss-Lobatto point. We present a computational method for solving a class of fractional integral equation of the second kind. Then based on Legendre-Gauss-Lobatto point and using, we derive a system of algebraic equations. The method is illustrated by applications and the results obtained are compared with the exact solutions in open literature. The obtained numerical results show that our proposed method is efficient and accurate for fractional integral equations of second kind. In addition, we prove that the error of the approximate solution decay exponentially in 𝐿2 norm. Keyword. Legendre Spectral-Collocation method, Fractional integral equation, Fractional order.

1. 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 has been many studies of fractional integral and differential equations (for instance, see [3-13]). 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 ([14-19]). In this paper, we are concerned with the numerical solutions of the following equation: (1) 𝑦(π‘₯) = π‘Ž(π‘₯) 𝐼 𝛼 {𝑏(π‘₯)𝑦(π‘₯)} + 𝑓(π‘₯), where 𝐼 Ξ± is the fractional integral of order 0 < 𝛼 < 1, π‘Ž, 𝑏 and 𝑓 are the continuous real functions on [0, 𝑇], 𝑦 is the unknown function that must be determined. Existence and uniqueness of the solution of the Eq.(1), in additional to some analytical properties and important inequalities, are investigated [20]. The authors in [20] applied the product trapezoidal method based on Picard iteration for approximating solution at the mesh point. 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, 𝑇], and therefore global method, 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 [21], 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 𝐿2 and 𝐿∞ norms. The results obtained in this paper may be considered as a generalization of the result in [20]. The structure of this paper is as follows: In section 2, some necessary definitions and mathematical tools of the fractional calculus which are required for our subsequent developments are introduced. In section 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, some numerical experiments are presented in Section 5 until efficiency of Legendre spectral-collocation method. Some conclusions are given in section 6. 2. Basic Definitions and Notations For the concept of fractional integral, we will adopt Riemann-Liouville definition. In whole of this paper, integral operator will be Riemann-Liouville fractional integral operator. Definition 1.2 A real function 𝑓 on [0, 𝑇] is said to be in the space πΆπœ‡ , πœ‡ ∈ ℝ, if there exists a real number 𝑝(> πœ‡), such that 𝑓(π‘₯) = π‘₯ 𝑝 𝑓1 (π‘₯), where 𝑓1 ∈ 𝐢[π‘œ, 𝑇], and it is said to be in the space πΆπœ‡π‘š iff 𝑓 (π‘š) ∈ πΆπœ‡ , π‘š ∈ β„•. Definition 2.2 The Riemann-Liouville fractional integral operator of order 𝛼 β‰₯ 0, of a function 𝑓 ∈ πΆπœ‡ , πœ‡ β‰₯ βˆ’1, is defined as π‘₯ 1 (2) 𝐼 𝛼 𝑓(π‘₯) = ∫ (π‘₯ βˆ’ 𝑑)π›Όβˆ’1 𝑓(𝑑)𝑑𝑑, 𝛼 > 0, π‘₯ ∈ [0, 𝑇], Ξ“(𝛼) 0 𝐼 0 𝑓(π‘₯) = 𝑓(π‘₯), where Ξ“(. ) is the Gamma function where

∞

Ξ“(𝛼) = ∫ 𝑑 π›Όβˆ’1 𝑒 βˆ’π‘‘ 𝑑𝑑. 0

Some properties of the Riemann-Liouville fractional integral are (𝑖) 𝐼 𝛼 𝐼𝛽 = 𝐼 𝛼+𝛽 ,

(𝑖𝑖) 𝐼 𝛼 𝐼𝛽 = 𝐼𝛽 𝐼 𝛼 , 𝛀(𝜈+1)

(𝑖𝑖𝑖) 𝐼 𝛼 (π‘₯ βˆ’ π‘Ž)𝜈 = (π‘₯ βˆ’ π‘Ž)𝛼+𝜈 , 𝛀(𝛼+𝜈+1) where 𝛼, 𝛽 β‰₯ 0 and 𝜈 > βˆ’1. Let πœ”π‘ž1 ,π‘ž2 (π‘₯) = (1 βˆ’ π‘₯)π‘ž1 (1 + π‘₯)π‘ž2 be a weight function in the usual sense, for π‘ž1 , π‘ž2 > βˆ’1. π‘ž ,π‘ž2 ∞ }𝑛=0

The Jacobi polynomials {𝐽𝑛 1

are orthogonal with respect to πœ”π‘ž1,π‘ž2 over interval (βˆ’1,1).

The set of Jacobi polynomials is a complete 𝐿2πœ”π‘ž1,π‘ž2 (βˆ’1,1) βˆ’orthogonal system, namely [17] 1

π‘ž ,π‘ž2 π‘ž1 ,π‘ž2 π‘ž1 ,π‘ž2 (π‘₯) π½π‘š πœ”

∫ 𝐽𝑛 1 βˆ’1

where π›Ώπ‘›π‘š is the Kronecker function, and

π‘ž ,π‘ž2

𝑑π‘₯ = 𝛾𝑛 1

π›Ώπ‘›π‘š ,

(3)

π‘ž ,π‘ž 𝛾𝑛 1 2

2π‘ž1 +π‘ž2 +1 Ξ“(π‘ž1 + 1)Ξ“(π‘ž2 + 1) , Ξ“(π‘ž1 + π‘ž2 + 2) = 2π‘ž1 +π‘ž2 +1 Ξ“(𝑛 + π‘ž1 + 1)Ξ“(𝑛 + π‘ž2 + 1) , {(2𝑛 + π‘ž1 + π‘ž2 + 1)n! Ξ“(𝑛 + π‘ž1 + π‘ž2 + 1)

𝑛 = 0, 𝑛 β‰₯ 1.

In particular, we find π‘ž ,π‘ž2

𝐽0 1

1 1 (π‘₯) = 1, 𝐽1π‘ž1 ,π‘ž2 (π‘₯) = (π‘ž1 + π‘ž2 + 2)π‘₯ + (π‘ž1 βˆ’ π‘ž2 ). 2 2

Let 𝑆𝑁 [βˆ’1,1] the set of all polynomials of degree ≀ 𝑁 (𝑁 β‰₯ 0). Thus, if we denote by π‘ž ,π‘ž2

{π‘₯𝑗 1

π‘ž ,π‘ž2 𝑁

, πœ”π‘— 1

}

𝑗=0

the nodes and the corresponding Christoffel numbers of the standard Jacobi-

Gauss interpolation on Ξ› = (βˆ’1,1), then we have 1

π‘ž ,π‘ž2

βˆ«βˆ’1 𝜈(π‘₯)πœ”π‘— 1

π‘ž1 ,π‘ž2 π‘ž ,π‘ž (π‘₯)𝑑π‘₯ = βˆ‘π‘ )πœ”π‘— 1 2 , 𝑗=0 𝜈(π‘₯𝑗

βˆ€πœˆ ∈ 𝑆2π‘βˆ’1 [βˆ’1,1].

(4)

π‘ž ,π‘ž

1 2 For any 𝑒 ∈ 𝐢(Ξ›), we denote by 𝑝π‘₯,𝑁 ∢ 𝐢(Ξ›) β†’ 𝑆𝑁 the Jacobi-Gauss interpolation operator, such that

π‘ž ,π‘ž

π‘ž ,π‘ž2

1 2 𝑝π‘₯,𝑁 𝑒(π‘₯𝑗 1

π‘ž ,π‘ž2

) = 𝑒(π‘₯𝑗 1

),

(5)

0 ≀ 𝑗 ≀ 𝑁.

It is clear that 𝑁 π‘ž1 ,π‘ž2 𝑝π‘₯,𝑁

π‘ž ,π‘ž2

𝑒(π‘₯) = βˆ‘ 𝑒𝑗 1

π‘ž ,π‘ž2

𝐽𝑗 1

(π‘₯),

(6)

𝑗=0

where π‘ž ,π‘ž 𝑒𝑗 1 2

=

𝑁

1 π‘ž ,π‘ž2

𝛾𝑗 1

π‘ž ,π‘ž2

βˆ‘ 𝑒(π‘₯𝑗 1

π‘ž ,π‘ž2

) 𝐽𝑗 1

π‘ž ,π‘ž2

(π‘₯𝑗 1

π‘ž ,π‘ž2

)πœ”π‘— 1

.

𝑗=0

In special case, if π‘ž1 = π‘ž2 = 0, then the Jacobi polynomial is reduced to the Legendre polynomial. Thus, we recall that the Legendre polynomials [17-18] {𝐿𝑖 (π‘₯); 𝑖 = 0,1, … } are defined on the Ξ› with the following recurrence formula: 2𝑖 + 1 𝑖 𝐿𝑖+1 (π‘₯) = π‘₯ 𝐿𝑖 (π‘₯) βˆ’ 𝐿 (π‘₯), 𝑖 = 1,2, …, 𝑖+1 𝑖 + 1 π‘–βˆ’1 with 𝐿0 (π‘₯) = 1 and 𝐿1 (π‘₯) = π‘₯. By using the Rodrigues formula 𝐿𝑖 (π‘₯) =

(βˆ’1)𝑖 2𝑖 𝑖!

𝐷 π‘˜ ((1 βˆ’ π‘₯ 2 )𝑖 ),

the Legendre polynomial has the expansion as follow 𝑖 [ ] 2

𝐿𝑖 (π‘₯) =

(2𝑖 βˆ’ 2π‘˜)! 1 βˆ‘(βˆ’1)π‘˜ π‘˜ π‘₯ π‘–βˆ’2π‘˜ . 𝑖 2 2 π‘˜! (𝑖 βˆ’ π‘˜)! (𝑖 βˆ’ 2π‘˜)! π‘˜=0

The set of {𝐿𝑖 (π‘₯); 𝑖 = 0,1, … } is a complete orthogonal system in 𝐿2 (Ξ›), and we have

(7)

1

(𝐿𝑗 , πΏπ‘˜ ) = ∫ 𝐿𝑗 (π‘₯) πΏπ‘˜ (π‘₯) 𝑑π‘₯ = β„Žπ‘— π›Ώπ‘—π‘˜ ,

(8)

βˆ’1

2

where π›Ώπ‘—π‘˜ is the Kronecker delta symbol and β„Žπ‘— = 2𝑗+1. Thus, for any 𝑔 ∈ 𝐿2 (Ξ›), we obtain ∞

𝑔(π‘₯) = βˆ‘ π‘Žπ‘– 𝐿𝑖 (π‘₯) , 𝑖=0

1 1 π‘Žπ‘– = ∫ 𝑔(π‘₯) 𝐿𝑖 (π‘₯) 𝑑π‘₯ . β„Žπ‘– βˆ’1

(9)

0,0 For simplicity, we let π‘₯𝑗 = π‘₯𝑗0,0 , πœ”π‘— = πœ”π‘—0,0 and 𝑝π‘₯,𝑁 = 𝑝𝑁π‘₯, . Thus, according to Eq.(4), we have 𝑁

1

∫ 𝜈(π‘₯)𝑑π‘₯ = βˆ‘ 𝜈(π‘₯𝑗 )πœ”π‘— , βˆ’1

βˆ€πœˆ ∈ 𝑆2π‘βˆ’1 [βˆ’1,1],

(10)

𝑗=0

where π‘₯𝑗 (0 ≀ 𝑗 ≀ 𝑁), i.e. are the zero of (1 βˆ’ π‘₯ 2 )𝐿′ 𝑁 (π‘₯), and πœ”π‘— =

2 𝑁 (𝑁+1) (𝐿𝑁 (π‘₯𝑗 ))

2

πœ”π‘— (0 ≀ 𝑗 ≀ 𝑁), i.e.

, are denoted to the nodes and Christoffel numbers of Legendre-Gauss-

Lobatto interpolation on the classical interval (βˆ’1,1), repectively. The norm and discrete inner product in 𝐿2 (Ξ›) are defined as 𝑁

(𝑒, 𝜈)𝑁 = βˆ‘ 𝑒(π‘₯𝑗 ) 𝜈(π‘₯𝑗 ) πœ”π‘— ,

1

‖𝑒‖𝑁 = (𝑒, 𝑒)2 .

(11)

𝑗=0

In the next section, we use some special techniques for converting Eq. (1) to a problem can be solved by the standard Legendre spectral collocation scheme. These topics will be presented in the next section. 3. The Spectral Collocation Method In this section, we shall propose a Legendre spectral collocation method for solving the Eq. (1). The Legendre spectral-collocation method takes advantage of both the Legendre polynomials and Legendre-Gauss-Lobatto interpolation nodes. The main idea is to use the suggested method consists 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 definitions 2.2 as follow π‘Ž(π‘₯) π‘₯ 𝑦(π‘₯) = ∫ 𝑏(𝑠)(π‘₯ βˆ’ 𝑠)π›Όβˆ’1 𝑦(𝑠) 𝑑𝑠 + 𝑓(π‘₯), Ξ“(𝛼) 0

π‘₯ ∈ 𝐼 = [0, 𝑇].

Next, in enjambment we let Ξ› = [βˆ’1,1]. For ease of analysis, we transfer the problem (1) to an equivalent problem defined in Ξ›. More specifically, we apply the change of variable 𝑑 =

2π‘₯ 𝑇

βˆ’1

𝑇

or π‘₯ = 2 (𝑑 + 1), such that 𝑑 ∈ Ξ›. Then, we have

𝑇 𝑦 ( (𝑑 + 1)) = 2

𝑇 π‘Ž ( 2 (𝑑 + 1)) Ξ“(𝛼)

∫

𝑇 (𝑑+1) 2

0

𝑇 +𝑓 ( (𝑑 + 1)) , 2

π›Όβˆ’1 𝑇 𝑏(𝑠) ( (𝑑 + 1) βˆ’ 𝑠) 𝑦(𝑠) 𝑑𝑠 2

(12)

𝑑 ∈ Ξ›.

𝑇

Moreover, to transfer the integral interval (0, 2 (𝑑 + 1)) to (βˆ’1, 𝑑), we make the transformation 𝑇

𝑠 = 2 (πœ‡ + 1). Then Eq. (12) can be written

𝑇 𝑇 𝑦 ( (𝑑 + 1)) = ( ) 2 2

𝛼

𝑇 π‘Ž (2 (𝑑 + 1)) Ξ“(𝛼)

𝑇 +𝑓 ( (𝑑 + 1)) , 2

𝑑 𝑇 𝑇 ∫ 𝑏 ( (πœ‡ + 1)) (𝑑 βˆ’ πœ‡)π›Όβˆ’1 𝑦 ( (πœ‡ + 1)) π‘‘πœ‡ 2 2 βˆ’1

(13)

𝑑 ∈ Ξ›.

Further, if we let 𝑇 π‘Œ(𝑑) = 𝑦 ( (𝑑 + 1)) , 2 𝑇 𝐡(𝑑) = 𝑏 ( (𝑑 + 1)) , 2

𝑇 𝐴(𝑑) = π‘Ž ( (𝑑 + 1)), 2 𝑇 𝐹(𝑑) = 𝑓 ( (𝑑 + 1)), 2

(14)

then, Eq. (13) can be reduced to 𝑇 𝛼 𝐴(𝑑) 𝑑 π‘Œ(𝑑) = ( ) ∫ 𝐡(πœ‡) (𝑑 βˆ’ πœ‡)π›Όβˆ’1 π‘Œ(πœ‡) π‘‘πœ‡ + 𝐹(𝑑), 𝑑 ∈ Ξ› = [βˆ’1,1]. 2 Ξ“(𝛼) βˆ’1

(15)

At last, under the following linear transformation πœ‡ = πœ‡(𝑑, πœƒ) = the Eq. (15) becomes

𝑑+1 π‘‘βˆ’1 πœƒ+ , 2 2

πœƒ ∈ Ξ›,

(16)

𝑇 𝛼 𝑑 + 1 𝛼 𝐴(𝑑) 1 π‘Œ(𝑑) = ( ) ( ) ∫ 𝐡(πœ‡(𝑑, πœƒ)) (1 βˆ’ πœƒ)π›Όβˆ’1 π‘Œ(πœ‡(𝑑, πœƒ)) π‘‘πœƒ + 𝐹(𝑑). 2 2 Ξ“(𝛼) βˆ’1

(17)

The Legendre spectral collocation method for Eq. (17) is to seek π‘ˆ ∈ 𝑆𝑁 (Ξ›) (𝑁 β‰₯ 1), such that 1 𝑇𝛼 π›Όβˆ’1,0 𝛼 π‘ˆ(𝑑) = Ξ± 𝑝𝑑,𝑁 {(𝑑 + 1) 𝐴(𝑑) ∫ (1 βˆ’ πœƒ)π›Όβˆ’1 π‘πœƒ,𝑁 (𝐡(πœ‡(𝑑, πœƒ)) π‘ˆ(πœ‡(𝑑, πœƒ))) π‘‘πœƒ} 4 Ξ“(𝛼) βˆ’1

(18)

+𝑝𝑑,𝑁 (𝐹(𝑑)).

Now, we want to describe a numerical implementation for Eq. (18). To this end, we use the following relations 𝑁

π‘ˆ(𝑑) = βˆ‘ 𝑒𝑖 𝐿𝑖 (𝑑), 𝑖=0 𝑁 π›Όβˆ’1,0 𝑝𝑑,𝑁 (π‘πœƒ,𝑁 ((𝑑

𝑁

+ 1)𝛼 𝐴(𝑑)𝐡(πœ‡(𝑑, πœƒ)) π‘ˆ(πœ‡(𝑑, πœƒ)))) βˆ‘ βˆ‘ 𝑑𝑖𝑗 𝐿𝑖 (𝑑) π½π‘—π›Όβˆ’1,0 (πœƒ) .

(19)

𝑖=0 𝑗=0

Then by using (19) and then Eq. (3), a direct computation leads to the following equation as 1 𝑇𝛼 π›Όβˆ’1,0 ∫ (1 βˆ’ πœƒ)π›Όβˆ’1 𝑝π‘₯,𝑁 (π‘πœƒ,𝑁 ((𝑑 + 1)𝛼 𝐴(𝑑)𝐡(πœ‡(𝑑, πœƒ)) π‘ˆ(πœ‡(𝑑, πœƒ)))) π‘‘πœƒ Ξ± 4 Ξ“(𝛼) βˆ’1 𝑁

𝑁

1 𝑇𝛼 (𝑑) = Ξ± βˆ‘ βˆ‘ 𝑑𝑖𝑗 𝐿𝑖 ∫ (1 βˆ’ πœƒ)π›Όβˆ’1 π½π‘—π›Όβˆ’1,0 (πœƒ) π‘‘πœƒ 4 Ξ“(𝛼) βˆ’1 𝑖=0 𝑗=0 𝑁

𝑇𝛼 = Ξ± βˆ‘ 𝑑𝑖0 𝐿𝑖 (𝑑) . 2 Ξ“(𝛼 + 1)

(20)

𝑖=0

Applying Eqs. (3) to (6), di0 can be readily obtained as 𝑁

𝑁

𝑙1

𝑙2

𝛼 (2𝑖 + 1) 𝛼 𝑑𝑖0 = )) π‘ˆ (πœ‡(π‘₯𝑙1 , πœƒπ‘™π›Όβˆ’1,0 )) 𝐿𝑖 (π‘₯𝑙1 )πœ”π‘™1 πœ”π‘™π›Όβˆ’1,0 . βˆ‘ βˆ‘(π‘₯𝑙1 + 1) 𝐴(π‘₯𝑙1 ) 𝐡 (πœ‡(π‘₯𝑙1 , πœƒπ‘™π›Όβˆ’1,0 𝛼+1 2 2 2 2

Hence, by inserting Eqs. (19) and (20) in Eq. (18), we conclude 𝑁

𝑁

𝑇𝛼 βˆ‘ 𝑒𝑖 𝐿𝑖 (𝑑) = Ξ± βˆ‘ 𝑑𝑖0 𝐿𝑖 (𝑑) . 2 Ξ“(𝛼 + 1) 𝑖=0

𝑖=0

(21)

Finally, by comparing the expansion coefficients of eq. (21), we can get 𝑇𝛼 𝑒𝑖 = Ξ± 𝑑 , 2 Ξ“(𝛼 + 1) 𝑖0

(22)

0 ≀ 𝑖 ≀ 𝑁.

By replacing the above relation, we obtain the expression of π‘ˆ(𝑑) accordingly.

4. Convergence Analysis In this section, we want to carry out the error analysis for the numerical method (18) under L2 (Ξ›) and L∞ (Ξ›). Here, we present some preparations as follow [14,17]:

L2πœ”π‘ž1,π‘ž2 (Ξ›): the Jacobi-weighted L2 Hilbert space with the scalar product 1

(u, v) = ∫ 𝑒(t) 𝑣(t) πœ”(t)dt , βˆ€u, v ∈ L2Ο‰q1,q2 (Ξ›). βˆ’1

Therefore, the corresponding norm is 1

‖𝑒‖𝐿2 π‘ž

πœ” 1 ,π‘ž2

π»πœ”π‘šπ‘ž1,π‘ž2 (𝛬) = {𝑒 ∈ 𝐿2πœ”π‘ž1,π‘ž2 (𝛬) :

𝑑𝑖 𝑒 𝑑π‘₯ 𝑖

= (𝑒, 𝑒)2 .

∈ 𝐿2πœ”π‘ž1,π‘ž2 (𝛬), 𝑖 = 0,1, … π‘š}: the Jacobi-weighted

Sobolev spaces with the norm and semi norm respectively π‘š

‖𝑒‖2𝐻 π‘šπ‘ž ,π‘ž πœ” 1 2

2

= βˆ‘β€–πœ•π‘‘π‘˜ 𝑒‖𝐿2

πœ”π‘ž1 ,π‘ž2

π‘˜=0

,

and

|𝑒|𝐻 π‘šπ‘ž ,π‘ž = β€–πœ•π‘‘π‘˜ 𝑒‖ π‘ž1+π‘˜,π‘ž2+π‘˜ πœ” 1 2 πœ”

π»πœ”π‘šπ‘ž1,π‘ž2 by its inner product is Hilbert space. For the purpose of convenience, we write 0

𝐿2 (Ξ›) = π»πœ”0,0 , β€–. β€–2 = β€–. ‖𝐿2(Ξ›) and β€–. β€–βˆž = β€–. β€–πΏβˆž(Ξ›) . m

Lemma 4.1 [17] Assume that 𝑒 ∈ Hπœ”π‘ž1,π‘ž2 with π‘ž1 , π‘ž2 > βˆ’1, and π‘š β‰₯ 1. Then for any integer number π‘˜ such that 0 ≀ π‘˜ ≀ π‘š ≀ 𝑁 + 1, exists the constant 𝐢 > 0 that the following relation holds: π‘ž1 ,π‘ž2 (𝑒))β€– β€–πœ•π‘‘π‘˜ (𝑒 βˆ’ 𝑝𝑑,𝑁

πœ”π‘ž1 +π‘˜ ,π‘ž2 +π‘˜

≀ 𝐢𝑁 π‘˜βˆ’π‘š β€–πœ•π‘‘π‘š π‘’β€–πœ”π‘ž1+π‘˜ ,π‘ž2+π‘˜ .

(23)

In particular, for any u ∈ H m (Ξ›) with 1 ≀ m ≀ N + 1, the above relation reduced to the following result 3

‖𝑒 βˆ’ 𝑝𝑑,𝑁 (𝑒)‖𝐻 1 (Ξ›) ≀ 𝐢𝑁 2βˆ’π‘š β€–πœ•π‘‘π‘š 𝑒‖2 .

(24)

Theorem 4.2 For sufficiently large 𝑁, the Legendre spectral-collocation approximation converges

to exact solution in 𝐿2 -norm, i.e. ‖𝑒𝑁 β€–2 = β€–π‘Œ βˆ’ π‘ˆβ€–2 β†’ 0. Proof. Assume that π‘ˆ(x) is obtained by using the Legendre spectral-collocation method Eq. (1). Then, according to definition of 𝑝𝑑,𝑁 in last section, we will have ‖𝑒𝑁 β€–2 ≀ β€–π‘Œ βˆ’ 𝑝𝑑,𝑁 (π‘Œ)β€– + ‖𝑝𝑑,𝑁 (π‘Œ) βˆ’ π‘ˆβ€– . 2 2

(25)

By using Lemma 4.1, we can get for any integer 1 ≀ π‘š ≀ 𝑁 + 1,

β€–π‘Œ βˆ’ 𝑝𝑑,𝑁 (π‘Œ)β€–2 ≀ 𝐢𝑁 βˆ’π‘š β€–πœ•π‘‘π‘š π‘Œβ€–πœ”π‘š ,π‘š .

(26)

The rest of our proof is to show ‖𝑝𝑑,𝑁 (π‘Œ) βˆ’ π‘ˆβ€–2 β†’ 0. To this end, along with definition π‘Œ in Eq. (17), for 𝑁 β‰₯ 1 we have 𝑇𝛼

𝑑

𝑝𝑑,𝑁 π‘Œ(𝑑) = 2𝛼 Ξ“(𝛼) 𝑝𝑑,𝑁 (𝐴(𝑑) βˆ«βˆ’1(𝑑 βˆ’ πœ‡)𝛼 𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) π‘‘πœ‡) + 𝑝𝑑,𝑁 (𝐹),

(27)

where, 𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) = 𝐡(πœ‡(𝑑, πœƒ))π‘Œ(πœ‡(𝑑, πœƒ)). Similarly, we can write 𝑑 𝑇𝛼 π‘ˆ(𝑑) = Ξ± 𝑝𝑑,𝑁 { 𝐴(𝑑) ∫ (𝑑 βˆ’ πœ‡)π›Όβˆ’1 π‘π›Όβˆ’1,0 (𝐹(𝐡(πœ‡), π‘ˆ(πœ‡))) π‘‘πœ‡} + 𝑝𝑑,𝑁 (𝐹(𝑑)), πœ‡,𝑁 2 Ξ“(𝛼) βˆ’1

(28)

where π›Όβˆ’1,0 π›Όβˆ’1,0 π‘πœ‡,𝑁 𝑣(πœ‡) = π‘πœƒ,𝑁 𝑣(πœ‡(𝑑, πœƒ))| 2πœ‡ π‘‘βˆ’1 . πœƒ= βˆ’

(29)

𝑑+1 𝑑+1

Before we proceed to the rest of proof, we present some relations that will be used later. By virtue of above definition and then using standard Jacobi-Gauss quadrature formula (4), we obtain 𝑑 1+𝑑 𝛼 1 ∫ (𝑑 βˆ’ πœ‡)π›Όβˆ’1 π‘π›Όβˆ’1,0 𝑣(πœ‡)π‘‘πœ‡ = ( ) ∫ (1 βˆ’ πœƒ)π›Όβˆ’1 π‘π›Όβˆ’1,0 𝑣(πœ‡(𝑑, πœƒ))π‘‘πœƒ πœ‡,𝑁 πœƒ,𝑁 2 βˆ’1 βˆ’1 𝑁

1+𝑑 𝛼 =( ) βˆ‘ 𝑣 (πœ‡(𝑑, πœƒπ‘—π›Όβˆ’1,0 )) πœ”π‘—π›Όβˆ’1,0 2 𝑗=0 𝑁

1+𝑑 𝛼 =( ) βˆ‘ 𝑣(πœ‡π‘—π›Όβˆ’1,0 )πœ”π‘—π›Όβˆ’1,0 . 2 𝑗=0

(30)

Again, in the same scheme we have 𝑑

∫ (𝑑 βˆ’ πœ‡) βˆ’1

π›Όβˆ’1

2

(π‘π›Όβˆ’1,0 𝑣(πœ‡)) πœ‡,𝑁

𝑁

1+𝑑 𝛼 π‘‘πœ‡ = ( ) βˆ‘ 𝑣 2 (πœ‡π‘—π›Όβˆ’1,0 )πœ”π‘—π›Όβˆ’1,0 . 2

(31)

𝑗=0

At last, by applying Lemma 4.1 and Eq. (29) together, for any integer π‘š (1 ≀ π‘š ≀ 𝑁 + 1), we conclude that 2

𝑑

π›Όβˆ’1,0 βˆ«βˆ’1(𝑑 βˆ’ πœ‡)π›Όβˆ’1 (𝑣(πœ‡) βˆ’ π‘πœ‡,𝑁 𝑣(πœ‡)) π‘‘πœ‡

=(

2 1+𝑑 𝛼 1 ( ))) ) ∫ (1 βˆ’ πœƒ)π›Όβˆ’1 (𝑣(πœ‡(𝑑, πœƒ)) βˆ’ π‘π›Όβˆ’1,0 𝑣 πœ‡ 𝑑, πœƒ π‘‘πœƒ ( πœƒ,𝑁 2 βˆ’1

2 1+𝑑 𝛼 1 π‘š ≀ 𝐢𝑁 βˆ’π‘š ( ) ∫ (πœ•πœƒ 𝑣(πœ‡(𝑑, πœƒ))) (1 βˆ’ πœƒ)𝛼+π‘šβˆ’1 (1 + πœƒ)π‘š π‘‘πœƒ 2 βˆ’1 π‘₯

2

= 𝐢𝑁 βˆ’π‘š ∫ (πœ•πœ‡π‘š 𝑣(πœ‡)) (π‘₯ βˆ’ πœ‡)𝛼+π‘šβˆ’1 (1 + πœ‡)π‘š π‘‘πœ‡ .

(32)

βˆ’1

Now, let us explain the rest of proof. By subtracting Eq. (28) from (27), we derive that 𝑝𝑑,𝑁 π‘Œ(𝑑) βˆ’ π‘ˆ(𝑑) =

𝑑 𝑇𝛼 π›Όβˆ’1,0 (𝑑 βˆ’ πœ‡)π›Όβˆ’1 (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ π‘πœ‡,𝑁 Γ— 𝑝 ∫ (𝐹(𝐡(πœ‡), π‘ˆ(πœ‡)))) π‘‘πœ‡). (𝐴(𝑑) 𝑑,𝑁 2𝛼 𝛀(𝛼) βˆ’1

The above formula can be rewritten as 𝑝𝑑,𝑁 π‘Œ(𝑑) βˆ’ π‘ˆ(𝑑) 𝑑 𝑇𝛼 π›Όβˆ’1,0 = 𝛼 Γ— 𝑝𝑑,𝑁 (𝐴(𝑑) ∫ (𝑑 βˆ’ πœ‡)π›Όβˆ’1 (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ π‘πœ‡,𝑁 (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)))) π‘‘πœ‡) 2 𝛀(𝛼) βˆ’1

+

𝑇𝛼 2𝛼 𝛀(𝛼) 𝑑

π›Όβˆ’1,0 Γ— 𝑝𝑑,𝑁 (𝐴(𝑑) ∫ (𝑑 βˆ’ πœ‡)π›Όβˆ’1 π‘πœ‡,𝑁 (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ 𝐹(𝐡(πœ‡), π‘ˆ(πœ‡))) π‘‘πœ‡ ).

(33)

βˆ’1

We let 𝐸1 (𝑑) =

𝑑 𝑇𝛼 π›Όβˆ’1,0 (𝑑 βˆ’ πœ‡)π›Όβˆ’1 (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ π‘πœ‡,𝑁 Γ— 𝑝 ∫ (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)))) π‘‘πœ‡), (𝐴(𝑑) 𝑑,𝑁 2𝛼 𝛀(𝛼) βˆ’1

𝐸2 (𝑑) =

𝑑 𝑇𝛼 π›Όβˆ’1,0 Γ— 𝑝𝑑,𝑁 (𝐴(𝑑) ∫ (𝑑 βˆ’ πœ‡)π›Όβˆ’1 π‘πœ‡,𝑁 (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ 𝐹(𝐡(πœ‡), π‘ˆ(πœ‡))) π‘‘πœ‡). 2𝛼 𝛀(𝛼) βˆ’1

Next, by using standard Legendre-Gauss quadrature formula we obtain ‖𝐸1 β€–2 𝑁

𝛼

=

2𝛼

1 2 2

π‘₯𝑗

𝑇 π›Όβˆ’1 π›Όβˆ’1,0 {βˆ‘ πœ”π‘— 𝐴(π‘₯𝑗 ) (∫ (π‘₯𝑗 βˆ’ πœ‡) (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ π‘πœ‡,𝑁 | (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)))) π‘‘πœ‡) } . π‘₯𝑗 𝛀(𝛼) βˆ’1 𝑗=0

By using the Cauchy-Schwartz inequality and Eq. (32), we get ‖𝐸1 β€–2 𝑁

π‘₯𝑗 π‘₯𝑗 𝑇𝛼 π›Όβˆ’1 π›Όβˆ’1 ≀ max 𝐴(𝑑) 𝛼 π‘‘πœ‡ ∫ (π‘₯𝑗 βˆ’ πœ‡) {βˆ‘ πœ”π‘— ∫ (π‘₯𝑗 βˆ’ πœ‡) (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) π‘‘βˆˆΞ› 2 𝛀(𝛼) βˆ’1 βˆ’1 𝑗=0

1 2

2 π›Όβˆ’1,0 βˆ’ π‘πœ‡,𝑁 | (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)))) π‘‘πœ‡} π‘₯𝑗

= max 𝐴(𝑑) π‘‘βˆˆΞ›

𝑇𝛼 2𝛼 𝛀(𝛼 + 1)

𝑁

𝛼

π‘₯𝑗

Γ— {βˆ‘ πœ”π‘— (π‘₯𝑗 + 1) ∫ (π‘₯𝑗 βˆ’ πœ‡)

π›Όβˆ’1

βˆ’1

𝑗=0

2

1 2

π›Όβˆ’1,0 (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ π‘πœ‡,𝑁 | (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)))) π‘‘πœ‡} π‘₯𝑗

𝑇𝛼 ≀ 𝐢 max 𝐴(𝑑) 𝛼 𝑁 βˆ’π‘š π‘‘βˆˆΞ› 2 𝛀(𝛼 + 1) 𝑁

𝛼

π‘₯𝑗

2

Γ— {βˆ‘ πœ”π‘— (π‘₯𝑗 + 1) ∫ (πœ•πœ‡π‘š 𝐹(𝐡(πœ‡), π‘Œ(πœ‡))) (π‘₯𝑗 βˆ’ πœ‡)

1 2 𝛼+π‘šβˆ’1

(1 + πœ‡)π‘š π‘‘πœ‡}

βˆ’1

𝑗=0

= 𝐢1 𝑁 βˆ’π‘š β€–πœ•πœ‡π‘š 𝐹(𝐡(. ), π‘Œ(. ))β€–πœ”π›Ό+π‘šβˆ’1,π‘š .

(34)

Similarly, by applying Eqs. (4) and (31) and the Cauchy-Schwartz inequality, we can get ‖𝐸2 β€–2 𝛼

=

2𝛼

𝑁

1 2 2

π‘₯𝑗

𝑇 π›Όβˆ’1 π›Όβˆ’1,0 π‘πœ‡,𝑁 | (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ 𝐹(𝐡(πœ‡), π‘ˆ(πœ‡))) π‘‘πœ‡] } {βˆ‘ πœ”π‘— [ 𝐴(π‘₯𝑗 ) ∫ (π‘₯𝑗 βˆ’ πœ‡) π‘₯𝑗 𝛀(𝛼) βˆ’1

≀ max 𝐴(𝑑) π‘‘βˆˆΞ›

𝑗=0

𝑇𝛼 𝛼

2 𝛀(𝛼)

π‘₯𝑗

𝑁

{βˆ‘ πœ”π‘— ∫ (π‘₯𝑗 βˆ’ πœ‡) 𝑗=0

Γ— ∫ (π‘₯𝑗 βˆ’ πœ‡) βˆ’1

π‘₯𝑗

π›Όβˆ’1

π‘‘πœ‡

βˆ’1

2 π›Όβˆ’1

1 2

(π‘π›Όβˆ’1,0 | (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ 𝐹(𝐡(πœ‡), π‘ˆ(πœ‡)))) π‘‘πœ‡} πœ‡,𝑁 π‘₯𝑗

= max 𝐴(𝑑) π‘‘βˆˆΞ›

𝑁

𝑇𝛼 2𝛼 𝛀(𝛼)

{βˆ‘ πœ”π‘—

( π‘₯ 𝑗 + 1) 𝛼

𝑗=0

π›Όβˆ’1,0

βˆ’ 𝐹 (𝐡(πœ‡π‘–

𝛼

𝑁

βˆ‘ πœ”π‘–π›Όβˆ’1,0 (𝐹 (𝐡(πœ‡π‘–π›Όβˆ’1,0 ), π‘Œ(πœ‡π‘–π›Όβˆ’1,0 )) 𝑖=0

2

1 2

), π‘ˆ(πœ‡π‘–π›Όβˆ’1,0 ))) } .

(35)

On the other hand, Jenson's inequality for any convex function 𝑔 on interval (π‘Ž, 𝑏) is 𝑔(πœ†π‘‘1 + (1 βˆ’ πœ†)𝑑2 ) ≀ πœ†π‘”(𝑑1 ) + (1 βˆ’ πœ†)𝑔(𝑑2 ),

βˆ€πœ† ∈ [0,1].

βˆ€π‘‘1 , 𝑑2 ∈ (π‘Ž, 𝑏),

Thus, since 𝑑2 𝑑 𝑑 (π‘₯𝑗 + 1) = (π‘₯𝑗 + 1) ln2 (π‘₯𝑗 + 1) > 0, 2 𝑑𝑑

π‘₯𝑗 is constant,

𝑑

then, for any given π‘₯𝑗 ∈ (βˆ’1,1), 𝑔(𝑑) = (π‘₯𝑗 + 1) is convex function and according to Jenson's inequality on interval [0,1], we will have 𝑔(𝑑) = 𝑔((1 βˆ’ 𝑑) Γ— 0 + 𝑑 Γ— 1) ≀ (1 βˆ’ 𝑑)𝑔(0) + 𝑑𝑔(1).

(36)

Therefore, by applying the previous inequality, we yield the following inequality 𝑁

𝑁

𝛼

βˆ‘ πœ”π‘— (π‘₯𝑗 + 1) ≀ βˆ‘ πœ”π‘— {(1 βˆ’ 𝛼) + 𝛼(π‘₯𝑗 + 1)} 𝑗=0

𝑗=0 𝑁

𝑁

= (1 βˆ’ 𝛼) βˆ‘ πœ”π‘— + 𝛼 βˆ‘ πœ”π‘— (π‘₯𝑗 + 1) 𝑗=0

𝑗=0

= (1 βˆ’ 𝛼) ∫ 𝑑𝑑 + 𝛼 ∫ (π‘₯ + 1)𝑑𝑑 Ξ›

= 2,

Ξ›

βˆ€π›Ό ∈ [0,1].

(37)

Next, by using relation (37), Eq. (35) will be reduced to ‖𝐸2 β€–2 ≀ max 𝐴(𝑑) π‘‘βˆˆΞ›

𝑇𝛼 2π›Όβˆ’1 𝛀(𝛼 + 1)

π‘₯𝑗

Γ— max { ∫ (π‘₯𝑗 βˆ’ πœ‡) 0≀𝑗≀𝑁

π›Όβˆ’1

βˆ’1

≀ max 𝐴(𝑑) max 𝐡(𝑑) π‘‘βˆˆΞ›

π‘‘βˆˆΞ›

2

1 2

π›Όβˆ’1,0 (π‘πœ‡,𝑁 | (𝐹(𝐡(πœ‡), π‘Œ(πœ‡)) βˆ’ 𝐹(𝐡(πœ‡), π‘ˆ(πœ‡)))) π‘‘πœ‡} π‘₯𝑗

𝑇𝛼 2π›Όβˆ’1 𝛀(𝛼 + 1)

π‘₯𝑗

Γ— max { ∫ (π‘₯𝑗 βˆ’ πœ‡) 0≀𝑗≀𝑁

π›Όβˆ’1

π›Όβˆ’1,0 (π‘πœ‡,𝑁 | (π‘Œ(πœ‡) βˆ’ π‘ˆ(πœ‡))) π‘‘πœ‡} π‘₯𝑗

βˆ’1

π‘‘βˆˆΞ›

𝑇𝛼 2π›Όβˆ’1 𝛀(𝛼 + 1)

π‘₯𝑗

π›Όβˆ’1

≀ max 𝐴(𝑑) max 𝐡(𝑑) π‘‘βˆˆΞ›

Γ— max {(∫ (π‘₯𝑗 βˆ’ πœ‡) 0≀𝑗≀𝑁

1 2

2

2

1 2

π›Όβˆ’1,0 (π‘πœ‡,𝑁 | (π‘Œ(πœ‡)) βˆ’ π‘Œ(πœ‡)) π‘‘πœ‡) π‘₯𝑗

βˆ’1

π‘₯𝑗

π›Όβˆ’1

+ (∫ (π‘₯𝑗 βˆ’ πœ‡)

2

1 2

(38)

(π‘Œ(πœ‡) βˆ’ π‘ˆ(πœ‡)) π‘‘πœ‡) }.

βˆ’1

At last, along with Lemma 4.1 and relations (32), (37), the previous result yields ‖𝐸2 β€–2 π‘₯𝑗

2

≀ 𝐢2 {𝐢𝑁 βˆ’π‘š max (∫ (πœ•πœ‡π‘š π‘Œ(πœ‡)) (π‘₯𝑗 βˆ’ πœ‡) 0≀𝑗≀𝑁

𝛼+π‘šβˆ’1

(1 + πœ‡)π‘š π‘‘πœ‡)

1 2

βˆ’1

1

1 2 2 2 + (∫ (π‘Œ(πœ‡) βˆ’ π‘ˆ(πœ‡)) π‘‘πœ‡) } 𝛼 βˆ’1

≀ 𝐢2 𝐢𝑁 βˆ’π‘š β€–πœ•πœ‡π‘š π‘Œβ€–πœ”π›Ό+π‘šβˆ’1,π‘š +

2𝐢2 ‖𝑒𝑁 β€–. 𝛼

(39)

Finally, relations (33) βˆ’ (34), (39) prove that the approximation is convergent in 𝐿2 βˆ’norm. Hence the theorem is proved. 5. Numerical example To show efficiency of our numerical method, the following examples which have exact solution are considered. Example 5.1. Consider the following fractional integral equation [20] 0.01

5

𝑑

1

3

𝑑3

𝑦(𝑑) = Ξ“(0.5) 𝑑 2 ∫0 (𝑑 βˆ’ 𝑠)βˆ’2 𝑦(𝑠) 𝑑𝑠 + βˆšπœ‹(1 + 𝑑)βˆ’2 βˆ’ 0.02 1+𝑑 ,

𝑑 ∈ [0,1].

3

The corresponding exact solution is given by 𝑦(𝑑) = βˆšπœ‹ (1 + 𝑑)βˆ’2 . Figure 1 present the approximate and exact solutions which are found in very good agreement. In Figure 2, the numerical errors are plotted in 𝐿2 βˆ’norm for 2 ≀ 𝑁 ≀ 24. As expected, an exponential rate of convergence is observed which confirmed our theoretical predictions in comparison of results in [20].

Figure 1. Comparison between exact solution and approximate solution of Example 5.1 for N = 10

Figure 2. The error function of Example 5.1 versus the number of interpolation operator Example 5.2. Our second example is the following fractional integral equation [20] 𝑦(𝑑) =

1

𝑑

1 2 1 8 ∫ 𝑠(𝑑 βˆ’ 𝑠)βˆ’3 𝑦(𝑠) 𝑑𝑠 + Ξ“ ( ) 𝑑 βˆ’ 𝑑 3, 2 0 3 40 27 Ξ“ ( ) 3

𝑑 ∈ [0,1].

2

The exact solution is 𝑦(π‘₯) = Ξ“ (3) 𝑑. We get the numerical solution of the above fractional integral equation as 2 𝑦1 (𝑑) = Ξ“ ( ) 𝑑, 3 2 𝑦2 (𝑑) = 1.21 Γ— 10βˆ’2 + Ξ“ ( ) 𝑑, 3 2 𝑦4 (𝑑) = 2.15 Γ— 10βˆ’9 + Ξ“ ( ) 𝑑 + 3.12 Γ— 10βˆ’10 𝑑 2 + 2.39 Γ— 10βˆ’7 𝑑 3 + 4.81 Γ— 10βˆ’9 . 3 The numerical results of proposed method can be seen from Figure 3 and Figure 4. These results indicate that the spectral accuracy is obtained for this problem, although the given function 𝑓(𝑑) is not very smooth.

Figure 3. Comparison between exact solution and approximate solution of Example 5.2

for 𝑁 = 6

Figure 4. The error function of Example 5.1 versus the number of interpolation operator

These results indicate that the spectral accuracy is obtained by using Legendre-Collocation method for this problem. Example 5.3. Consider the following fractional integral equation [20] 𝑑 1 1 𝑑 𝑦(𝑑) = βˆ’ ∫ (𝑑 βˆ’ 𝑠)βˆ’2 𝑦(𝑠) 𝑑𝑠 + 2√ , Ξ“(0.5) 0 πœ‹

𝑑 ∈ [0,1].

The exact solution is 𝑦(𝑑) = 1 βˆ’ et π‘’π‘Ÿπ‘“π‘(βˆšπ‘‘), such that π‘’π‘Ÿπ‘“π‘ is the complementary error function defined by π‘’π‘Ÿπ‘“π‘(𝑑) = 1 βˆ’ π‘’π‘Ÿπ‘“(𝑑) =

2 βˆšπœ‹

∞

2

∫ 𝑒 βˆ’π‘  𝑑𝑠. 𝑑

We list the numerical results in Table 1 to be able to test the convergence rate of the suggested. In Figure 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.

t

Proposed method at 𝑁 = 8 .0000000000 .2764215215 .3562117262 .4079815809 .4463937357 .4768434165 .5019754248 .5232972359 .5417539715 .5579785802 .5724164231

Proposed method at 𝑁 = 10 .0000000000 .2764215613 .3562117272 .4079815886 .4463937452 .4768434157 .5019754315 .5232972687 .5417539771 .5579785884 .5724164243

Exact solution

0 .0000000000 0.1 .2764215616 0.2 .3562117278 0.3 .4079815886 0.4 .4463937458 0.5 .4768434159 0.6 .5019754315 0.7 .5232972687 0.8 .5417539772 0.9 .5579785885 1 .5724164241 Table 1. Comparison of approximate solutions and exact solutions for various numbers N

Figure 2. The error function of Example 5.1 versus the number of interpolation operator for 𝑁 = 16

6. 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 𝐿2 βˆ’ π‘›π‘œπ‘Ÿπ‘š. We prove that our proposed method

is effective and has high convergence rate and better than some new methods which be define in this paper. During three numerical applications, we ensure that the proposed method is efficient and powerful numerical scheme for solving many fractional problems.

Reference [1] A.A. Kilbas, H.M. Srivastava and J.J. Trujillo, Theory and Applications of Fractional Differential Equations, North-Holland Math. Stud., vol. 204, Elsevier Science B.V., Amsterdam, (2006). [2] I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, (1999). [3] A. H. Bhrawy and M. A. Alghamdi, A shifted Jacobi-Gauss-Lobatto collocation method for solving nonlinear fractional Langevin equation involving two fractional orders in different intervals, Boundary Value Problems, vol. 2012, article 62, 13 pages, 2012. [4] E. Rawashdeh, Numerical solution of fractional integro-differential equations by collocation method, Appl. Math. Comput. 176, 1-6 (2006). [5] Y. Nawaz, Variational iteration method and Homotopy perturbation method for fourth-order fractional integro-differential equations, Comput. Math. Appl. 61, 2330-2341 (2011). [6] S. Momani and M. A. Noor, Numerical method for fourth-order fractional integro-differential equations, Appl. Math. Comput. 182,754-760 (2006). [7] H. Saeedi and F. Samimi, He's homotopy perturbation method for nonlinear Fredholm integrodifferential equation of fractional order, International J. of Eng. Res. And App., Vol. 2, no. 5, pp. 52-56, 2012. [8] A. Ahmed and S. A. H. Salh, Generalized Taylor matrix method for solving linear integro-fractional differential equations of Volterra type, App. Math. Sci., Vol. 5, no. 33-36, pp. 1765-1780, 2011. [9] A. Arikoglo and L. Ozkol, Solution of fractional integro-differential equation by using fractional differential transform method, Chaos Solitons Fractals 34, 1473-1481 (2007). [10] J. Wei and T. Tian, Numerical solution of nonlinear Volterra integro-differential equations of fractional order by the reproducing kernel method, Appl. Math. Model. 39, 4871-4876 (2015). [11] Li Zhu and Qibin Fan, Solving fractional nonlinear Fredholm integro-differential equations by the second kind Chebyshev wavelet, Common. Nonlinear Sci. Numer. Simul. 17, 2333–2341(2012). [12] H. Saeedi, M. Mohseni Moghadam, N. Mollahasani and GN. Chuev, A CAS wavelet method for solving nonlinear Fredholm integro-differential equation of fractional order, Common. Nonlinear Sci. Numer. Simul. 16, 1154-1163 (2011). [13] Z. Meng, L. Wang, H. Li and W. Zhang, Legendre wavelets method for solving fractional integrodifferential equations, Int. J. Comput. Math. 92, 1275-1291 (2015). [14] C. Canuto, M. Y. Hussaini, A. Quarteroni and T. A. Zang, Spectral methods: Fundamentals in single domain, Springer, Berlin (2006). [15] P. E. Bjorstad and B. P. Tjostheim, Efficient algorithm for solving a fourth-order equation with the spectral Galerkin method, SIAM J. Sci. Comput., 18, pp. 621-632 (1997). [16] J. Shen, Efficient spectral-Galerkin method I. Direct solvers for second- and fourth-order equations by using Legendre polynomials, SIAM J. Sci. Comput., 15, pp. 1489-1505 (1994). [17] J. Shen, T. Tang and Li-Lian Wang, Spectral Methods: Algorithms, Analysis and Applications, Springer-Verlag, Berlin Heidelberg (2011). [18] A. Yousefi, T. Mahdavi-Rad and S. G. Shafiei, A quadrature Tau method for solving Fractional integro-differential equations in the Caputo sense, J. of Math. Computer Sci. 16, 97-107 (2015).

[19] W. G. El-Sayed and A. M.A. El-Sayed, On the fractional integral equations of mixed type integrodifferential equations of fractional orders, App. Math. Comput. 154, 461-467 (2004). [20] S. Micula, An iterative numerical method for fractional integral equations of the second kind, J. of. Comput. and Appl. Math. (2017) DOI: http://doi.org/10.1016/j.cam.2017.12.006. [21] Y. Yang, Y. Chen and Y. Huang, Convergence analysis of the Jacobi spectral-collocation method for fractional integro-differential equations, J. of Acta Math. Sci. 34B (3), 673-690 (2014). [22] G. Mastroianni and D. Occorsto, Optimal systems of nodes for Lagrange interpolation on bounded intervals: a survey, J. Comput. Appli. Math. 134, 325-341 (2010).