## Services on Demand

## Article

## Indicators

## Related links

- Similars in SciELO

## Share

## Computational & Applied Mathematics

*On-line version* ISSN 1807-0302

### Comput. Appl. Math. vol.30 no.3 São Carlos 2011

#### http://dx.doi.org/10.1590/S1807-03022011000300010

**Operational Tau approximation for a general class of fractional integro-differential equations **

**S. Karimi Vanani; A. Aminataei**

Department of Mathematics, K.N. Toosi University of Technology, P.O. Box: 16315-1618, Tehran, Iran, E-mails: solatkarimi@dena.kntu.ac.ir / ataei@kntu.ac.ir

**ABSTRACT**

In this work, an extension of the algebraic formulation of the operational Tau method (OTM) for the numerical solution of the linear and nonlinear fractional integro-differential equations (FIDEs) is proposed. The main idea behind the OTM is to convert the fractional differential and integral parts of the desired FIDE to some operational matrices. Then the FIDE reduces to a set of algebraic equations. We demonstrate the Tau matrix representation for solving FIDEs based on arbitrary orthogonal polynomials. Some advantages of using the method, errorestimation and computer algorithm are also presented. Illustrative linear and nonlinear experiments are included to show the validity and applicability of the presented method.

**Mathematical subject classification:** 65M70, 34A25, 26A33, 47Gxx.

**Key words:** spectral methods, operational Tau method, fractional integro-differential equations, error estimation, computer algorithm of the method.

**1 Introduction**

The main object of this work is to solve the fractional integro-differential equation of the following form:

where F and G are given smooth functions, D ^{α} is a fractional differential operator of order α in the Caputo's sense, B_{j}; j = 0,1,...,*m*-1, are *m* supplementary conditions and *u*(*x*) is a solution to be determined.

Differential and integral equations involving derivatives of non-integer order have shown to be adequate models for various phenomena arising in damping laws, diffusion processes, models of earthquake [1], fluid-dynamics traffic model [2], mathematical physics and engineering [3], fluid and continuum mechanics [4], chemistry, acoustics and psychology [5].

Some numerical methods for the solution of FIDEs are presented in the literature. We can point out to the collocation method [6], Adomian decomposition method (ADM) [7]-9], Spline collocation method [10], fractional differential transform method [11] and the method of combination of forward and central differences [12]. Out of the aforesaid methods, we desire to consider OTM for solving FIDEs.

Spectral methods provide a computational approach which achieved substantial popularity in the last three decades. Tau method is one of the most important spectral methods which is extensively applied for numerical solution of many problems. This method was invented by Lanczos [13] for solving ordinary differential equations (ODEs) and then the expansion of the method were done for many different problems such as partial differential equations (PDEs) [14]-[16], integral equations (IEs) [17], integro-differential equations (IDEs) [18] and etc. [19]-[22].

In this work, we are interested in solving FIDEs with an operational approach of the Tau method. Because in the Tau method, we are dealing with a system of equations wherein the matrix of unknown coefficients is sparse and can be easily invertible. Also, the differential and integral parts appearing in the equation is replaced by its operational Tau representation. Then, we obtain a system of algebraic equations wherein its solution is easy.

This work has been organized as follows. In Section 2, we briefly state the basic definitions of the fractional calculus. Section 3 is devoted to introduce the OTM and its application on FIDEs. In Section 4, the Legendre and Laguerre polynomials as base functions are considered. An efficient error estimation of the proposed method is presented in Section 5. The accuracy and efficiency of the scheme is investigated with three illustrative numerical experiments inSection 6. Finally, Section 7 consists of some obtained conclusions.

**2 Basic definitions of the fractional calculus**

In this section, we give some basic definitions and properties of the fractional calculus theory, which are used in this work [3,23].

**Definition 1.** A real function *u*(*x*), *x* > 0 is said to be in the space C _{µ}, µ ∈ , if there exists a real number p > µ, such that *u*(*x*) = *x*^{p}*v*(*x*), where *v*(*x*) ∈ C[0, ∞) and it is said to be in the space iff u^{(m)}(*x*) ∈ C _{µ}, *m* ∈ .

**Definition 2.** The Riemann-Liouville fractional integral operator of order α __>__ 0, of a function *u*(x) ∈ C _{µ}, µ __>__ -1, is defined as:

where Γ is the Gamma function. Some of the most important properties of operator J ^{α} for u(x) ∈ C _{µ}, µ __>__ -1, α, β __>__ 0 and γ > -1, are as follows:

**Definition 3.** The fractional derivative of *u*(*x*) in the Caputo's sense is defined as:

when m-1 < α __<__ m, m ∈ , x > 0, u(x) ∈ .

**3 Operational Tau method**

In this section, we state the structure of OTM and its application on FIDEs. Also, we present some preliminaries and notations using in this work.

For any integrable functions ψ(x) and ϕ(x) on [a,b], we define the scalar product ‹ , ∇ by

where and ω(x) is a weight function. Let be the space of all functions f:[a,b] →, with < ∞.

The main idea of the method is to seek a polynomial to approximate u(x) ∈ . Let ϕ_{x} = = **Φ****X**_{x} be a set of arbitrary orthogonal polynomial bases defined by a lower triangular matrix **Φ** and **X**_{x} = [1,x,x^{2},...]^{T}.

**Lemma 1.** Suppose that *u*(*x*) is a polynomial as u(x) = = **uX**_{x} , then we have:

and

where **u** = [u_{0},u_{1},...,u_{n},...], **X**_{a} = [1,a,a^{2},...]^{T}, a ∈ and **M**,**N** and **P** are infinite matrices with only nonzero elements **M**_{i+1,i} = i+1, **N**_{i,i+1} = 1,** **

**Proof. ** See [24].

□

Let us consider

to be an orthogonal series expansion of the exact solution of equations (1) and (2), where **u** = is a vector of unknown coefficients, **Φ****X**_{x} is an orthogonal basis for polynomials in .

In the Tau method, the aim is to convert equations (1) and (2) to an algebraic system using some operational matrices. First of all, we intend to describe an operational form of D ^{α} u(x). In order to do this, we assert the process as follows.

Using equations (3), (4) and (7); we get:

Using property iii) of Definition 2, we have:

where **Γ** is an infinite diagonal matrix with elements

and

By approximating x^{m-} ^{α}^{+} ^{γ}, γ = 0,1,...; as follows:

Substituting equation (10) in equation (9) and using equation (8); we obtain:

In next step, the aim is to linearize analytic functions *F*(*x*,*u*(*x*)) and *G*(*x*, *t*, *u*(*t*)). These functions can be written as:

and

Now, we state the following lemma and corollary.

**Lemma 2.** Let **X**_{x} = [1,x,x^{2},...]^{T}, **u** = [u_{0},u_{1},u_{2},...] be infinite vectors and \boldsymbol Φ = [ ϕ_{0}| ϕ_{1}| ϕ_{2}|...], ϕ_{i} are infinite columns of matrix Φ. Then, we have:

where **U** is an upper triangular matrix as:

In addition, if we suppose that u(x) = **u** **Φ****X**_{x} represents a polynomial, then for any positive integer p, the relation

is valid.

**Proof. ** We have:

therefore

if we call the last upper triangular coefficient matrix as **U**, then we have:

Now, in order to prove equation (16), we apply induction. For p = 1, it is obvious that u(x) = **u** **Φ****X**_{x}. For p = 2, we rewrite u^{2}(x) = **u** **Φ****X**_{x}**u** **Φ****X**_{x} = **u** **Φ**(**X**_{x}**u** **Φ****X**_{x}) and using equation (14), we have:

therefore, equation (16) is hold for p = 2. Now, suppose that equation (16) is hold for p = k, then we must prove that the relation is valid for s = k+1.Thus, u^k+1(x)

So, equation (16) is proved.

We rewrite equation (12) as:

The first summation can be considered as:

For the second summation, using equations (5) and (16) yields:

Therefore, we have:

In the same manner, we rewrite:

Using equation (5), the second summation is as:

Therefore,

The first part easily can be written as follows:

where

Using equations (5) and (6), the second part of equation is written as follows:

Thus, we have:

In addition, suppose that the supplementary conditions are generally as:

Therefore, using equation (4) we have:

Using equations (11), (18) and (27) we replace equation by the following operational form:

So, the residual R(x) of equation , can be written as:

where

Now, we set the residual matrix **R** = 0 or we use the following inner products,

Therefore an infinite system of algebraic equations is obtained. Since, somewhere we require finite terms of approximation, then we must truncate the series to finite number of terms. Thus, we choose n+1-m of the first equations in **R** and the following system is obtained:

or

Solving the aforesaid system yields the unknown vector **u** = [u_{0},u_{1},...,u_{n}]. This is the so-called * operational Tau method* which is applicable for finite, infinite, regular and irregular domains.

We summarize OTM in the following algorithm.

**4 Some orthogonal polynomials**

Orthogonal functions have received considerable attention in dealing with several problems. Their most important characteristic is reducing the computations and converting the problem to a system of algebraic equations. In this work, we use Legendre polynomials for finite domain [0,h] and Laguerre polynomials for infinite domain [0, ∞].

*4.1 Legendre polynomials*

The well-known Legendre polynomials are defined as follows:

They are orthogonal on the interval [-1, 1] with respect to the weight function w(x) = 1 and satisfy:

Since the interval of orthogonality of these polynomials may differ with the domain of the problem, we must shift the polynomial to the desired interval. Thus, to construct the shifted Legendre polynomials on arbitrary interval [0,h], it is sufficient to do the change of variable: . So, the shifted Legendre polynomials are defined as follows:

and satisfy:

*4.2 Laguerre polynomials*

Laguerre polynomials are defined as follows:

They are orthogonal on the interval [0, ∞) with respect to the weight function w(x) = e^{-x} and satisfy:

**5 Error estimation**

In this section, an error estimation for the approximate solution of equation with supplementary conditions (2) is obtained. Let us call e_{n}(x) = u(x)- u_{n}(x) as the error function of the approximate solution u_{n}(x) to u(x) where u(x) is the exact solution of equation (1). Hence, u_{n}(x) satisfies the following equations:

The perturbation term H_{n}(x) can be obtained by substituting the computedsolution u_{n}(x) into the equation:

We proceed to find an approximation e_{n,N}(x) to the error function e_{n}(x) in the same way as we did before for the solution of equation (1). Note that N is the degree of approximation of u_{n}(x). By subtracting equations (31) and (32) from equations (1) and (2), respectively we have:

or

It should be noted that in order to construct the approximate e_{n,N}(x) to e_{n}(x), only the related equations like as equations (7) through (31) needs to be recomputed and the structure of the method remains the same.

**6 Illustrative numerical experiments with some comments**

In this section, four experiments of linear and nonlinear FIDEs are given to illustrate the results. In all experiments, we consider the shifted Legendre polynomials as basis functions for finite domains and Laguerre polynomials for infinite domains. The computations associated with the experiments discussed above were performed in Maple 13 on a PC with a CPU of 2.4 GHz.

**Experiment 6.1. ** Consider the following FIDE [6,7]:

with the initial condition: u(0) = 0

The exact solution is: u(x) = x^{3}.

We have solved this experiment using OTM with Laguerre polynomials and some approximations are obtained as follows:

and so on.

Thus, we have u(x) = x^{3} which is the exact solution of the problem.

**Experiment 6.2. ** Consider the following nonlinear FIDE [8]:

with the boundary conditions:

u(0) = u'(0) = 1,

u(1) = u'(1) = e.

The only case which we know the exact solution for α = 4 is: u(x) = e^{x}.

We have solved this experiment for n = 7 with different α and have compared it with the closed form series solutions of the exact solution obtained by ADM [8]. The comparison is shown in Table 1.

From the numerical results in Table 1, it is easy to conclude that obtained results by OTM are in good agreement with those obtained using the ADM. In this experiment, the exact solution is not known, so a main question arises that which method is more accurate. We conclude that OTM is more accurate by considering the following notation and discussion.

**Note 1. ** In the theory of fractional calculus, it is obvious that when the fractional derivative α (m-1 < α __<__ m) tends to positive integer number m, then the approximate solution *continuously* tends to the exact solution of the problem with derivation m.

A closer look at the values obtained by ADM in Table 1 do not have this characteristic. The values for α = 3.25 are less than the values for α = 3.5, so the values for α = 3.5 must be less than the values for α = 3.75. But this fact have not occurred for ADM solutions but for OTM solutions, we have the reduction in the results. Therefore OTM is more reliable than ADM.

Moreover, for α = 4, using Legendre polynomials, the following sequences of approximate solution is obtained:

and so on.

Thus, we obtain:

This has the closed form u(x) = e^{x}, which is the exact solution of the problem. Thus, for positive integer derivatives, if the exact solution exists, then OTM produces its series solution.

**Experiment 6.3. ** Consider the following nonlinear FIDE:

with the initial conditions:

The only case which we know the exact solution for α = 2 is u(x) = sinx.

We have solved this experiment for n = 5 with different α and have used Laguerre polynomials as basis functions. The results are given in Table 2.

For α = 2, using Laguerre polynomials, the following sequence of approximate solution is obtained:

and so on.

Thus, we obtain:

This has the closed form u(x) = sinx, which is the exact solution of the problem.

**Experiment 6.4. ** Consider the following nonlinear FIDE:

with the initial conditions:

The only case which we know the exact solution for α = 2 is: u(x) = e^{-x}.

We have solved this experiment for n = 4 with different α and used shifted Legendre polynomials as basis functions. Figure 1 shows the approximate solutions and illustrates the aforesaid fact given by note 1.

The following figure illustrates the convergency of the method and the factthat the method tends continuously to the exact solution if fractional derivations tend to an integer order.

For α = 2, using shifted Legendre polynomials, the following sequence of approximate solution is obtained:

and so on.

Thus, we obtain:

This has the closed form u(x) = e^{-x}, which is the exact solution of theproblem.

**7 Conclusion**

In this work, operational Tau method is employed successfully to solve theFIDEs. Arbitrary orthogonal polynomial bases were applied as basis functions.Reducing the FIDEs to algebraic equations is the first characteristics of theproposed method. The applications of OTM on some problems including linear and nonlinear terms are considered and some useful results are obtained.The most important ones are the simplicity of the method, reducing the computations using orthogonal polynomials and having low run time of its algorithm. Furthermore, this method yields the desired accuracy only in a few terms in a series form of the exact solution. All of these advantages of the OTM to solve nonlinear problems assert the method as a convenient, reliable andpowerful tool.

**REFERENCES**

[1] J.H. He, *Nonlinear oscillation with fractional derivative and its applications*. In: International Conference on Vibrating Engineering, Dalian, China (1998), 288-291. [ Links ]

[2] J.H. He, *Some applications of nonlinear fractional differential equations and their approximations.* Bull. Sci. Technol., **15** (1999), 86-90. [ Links ]

[3] I. Podlubny, *Fractional Differential Equations.* Academic Press, New York (1999). [ Links ]

[4] F. Mainardi, *Fractals and Fractional Calculus Continuum Mechanics.* Springer Verlag, (1997), 291-348. [ Links ]

[5] W.M. Ahmad and R. El-Khazali, *Fractional-order dynamical models of love.* Chaos, Solitons & Fractals, **33** (2007), 1367-1375. [ Links ]

[6] E.A. Rawashdeh, *Numerical solution of fractional integro-differential equations by collocation method.* Applied Mathematics and Computation, **176** (2005), 1-6. [ Links ]

[7] R.C. Mittal and R. Nigam, *Solution of fractional integro-differential equations by adomian decomposition method.* Int. J. of Appl. Math. and Mech., **4** (2008),87-94. [ Links ]

[8] S. Momani and M.A. Noor, *Numerical methods for fourth-order fractional integro-differential equations.* Applied Mathematics and Computation, **182** (2006), 754-760. [ Links ]

[9] W.G. El-Sayed and A.M.A. El-Sayed, *On the functional integral equations of mixed type and integro-differential equations of fractional orders.* Applied Mathematics and Computation, **154** (2004), 461-467. [ Links ]

[10] A. Pedas and E. Tamme, *Spline collocation method for integro-differential equations with weakly singular kernels.* Journal of Computational and Applied Mathematics, **197** (2006), 253-269. [ Links ]

[11] D. Nazari and S. Shahmorad, *Application of the fractional differential transform method to fractional-order integro-differential equations with nonlocal boundary conditions.* Journal of Computational and Applied Mathematics, **234** (2010), 883-891. [ Links ]

[12] M.F. Al-Jamal and E.A. Rawashde, *The Approximate Solution of Fractional Integro-Differential Equations.* Int. J. Contemp. Math. Sciences, **4** (2009), 1067-1078. [ Links ]

[13] C. Lanczos, *Trigonometric interpolation of empirical and analytical functions.*J. Math. Phys., **17** (1938), 123-199. [ Links ]

[14] K.M. Liu and E.L. Ortiz, *Numerical solution of eigenvalue problems for partial differential equations with the Tau-lines method.* Comp. Math. Appl. B, **12** (1986), 1153-1168. [ Links ]

[15] E.L. Ortiz and K.S. Pun, *Numerical solution of nonlinear partial differential equations with Tau method.* J. Comp. Appl. Math., **12** (1985), 511-516. [ Links ]

[16] E.L. Ortiz and H. Samara, *Numerical solution of partial differential equations with variable coefficients with an operational approach to the Tau method.* Comp.Math. Appl., ** 10** (1984), 5-13. [ Links ]

[17] M.K. EL-Daou and H.G. Khajah, *Iterated solutions of linear operator equations with the Tau method.* Math. Comput., ** 66**(217) (1997), 207-213. [ Links ]

[18] J. Pour-Mahmoud, M.Y. Rahimi-Ardabili and S. Shahmorad, *Numerical solution of the system of Fredholm integro-differential equations by the Tau method.* Applied Mathematics and Computation, **168** (2005), 465-478. [ Links ]

[19] K.M. Liu and E.L. Ortiz, *Approximation of eigenvalues defined by ordinary differential equations with the Tau method.* Matrix Pencils, Springer, Berlin, (1983), 90-102. [ Links ]

[20] K.M. Liu and E.L. Ortiz, *Tau method approximation of differential eigenvalue problems where the spectral parameter enters nonlinearly.* J. Comput. Phys.,**72** (1987), 299-310. [ Links ]

[21] K.M. Liu and E.L. Ortiz, *Numerical solution of ordinary and partial functional-differential eigenvalue problems with the Tau method.* Computing, **41** (1989),205-217. [ Links ]

[22] E.L. Ortiz and H. Samara, *Numerical solution of differential eigenvalue problems with an operational approach to the Tau method.* Computing, **31** (1983), 95-103. [ Links ]

[23] S. Samko, A. Kilbas and O. Marichev, *Fractional Integrals and Derivatives.*Gordon and Breach, Yverdon (1993). [ Links ]

[24] K.M. Liu and C.K. Pan, *The automatic solution to systems of ordinary differential equations by the Tau method.* Computers Math. Applicat., **38** (1999), 197-210. [ Links ]

Received: 13/XI/10.

Accepted: 25/IV/11.

#CAM-285/10.