Advances in Mathematical Physics
Volume 2017 (2017), Article ID 3821870, 12 pages
https://doi.org/10.1155/2017/3821870

Research Article

New Operational Matrix via Genocchi Polynomials for Solving Fredholm-Volterra Fractional Integro-Differential Equations

Jian Rong Loh, Chang Phang, and Abdulnasir Isah

Department of Mathematics and Statistics, Faculty of Science, Technology and Human Development, Universiti Tun Hussein Onn Malaysia, Johor, Malaysia

Correspondence should be addressed to Chang Phang; pchang@uthm.edu.my

Received 14 September 2016; Revised 10 November 2016; Accepted 8 December 2016; Published 16 January 2017

Academic Editor: Luigi C. Berselli

Copyright © 2017 Jian Rong Loh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

It is known that Genocchi polynomials have some advantages over classical orthogonal polynomials in approximating function, such as lesser terms and smaller coefficients of individual terms. In this paper, we apply a new operational matrix via Genocchi polynomials to solve fractional integro-differential equations (FIDEs). We also derive the expressions for computing Genocchi coefficients of the integral kernel and for the integral of product of two Genocchi polynomials. Using the matrix approach, we further derive the operational matrix of fractional differentiation for Genocchi polynomial as well as the kernel matrix. We are able to solve the aforementioned class of FIDE for the unknown function . This is achieved by approximating the FIDE using Genocchi polynomials in matrix representation and using the collocation method at equally spaced points within interval . This reduces the FIDE into a system of algebraic equations to be solved for the Genocchi coefficients of the solution . A few numerical examples of FIDE are solved using those expressions derived for Genocchi polynomial approximation. Numerical results show that the Genocchi polynomial approximation adopting the operational matrix of fractional derivative achieves good accuracy comparable to some existing methods. In certain cases, Genocchi polynomial provides better accuracy than the aforementioned methods.

1. Introduction

Fractional integro-differential equation (FIDE) is an equation that contains a fractional derivative term and an integral kernel operator term . An example is given as follows:

Here, is the unknown function to be solved and , are the Volterra integral kernel and Fredholm integral kernel, respectively. is Caputo’s fractional derivative and , are coefficient functions.

In the past few years, various wavelet operational matrices and polynomial operational matrices for both fractional differentiation and fractional integration have been derived one after another to solve various kinds of fractional differential equations (FDEs) and fractional integro-differential equations (FIDEs). Examples of polynomial operational matrices that have been proposed in the past are Shifted Legendre polynomials [1], shifted Chebyshev polynomials [2], and shifted Jacobi polynomials [3] to name a few. The operational matrix method is usually combined with methods such as collocation method, spectral tau method, and tau method to solve FDE and FIDE. In [4], Zhu and Fan have proposed Chebyshev wavelet operational matrix of fractional integration and applied it to solve a certain type of nonlinear fractional integro-differential equations (FIDEs). Apart from the operational matrix approach, there exist various approaches such as Adomian decomposition method [5], spectral tau method [6], Taylor expansion method [7], and Bernoulli matrix method [8]. In [9], Alshbool et al. have proposed a new Bernstein function to solve fractional order differential equations.

In this paper, different from the existing methods in numerically solving FIDE, we explore the application of a relatively new polynomial, namely, Genocchi polynomials, in the numerical solution of FIDE. Genocchi polynomial belongs to a larger family class of polynomials which is the Appell polynomial family. Throughout this paper, we denote Genocchi polynomials by . We apply these polynomials to solve the FIDEs given the advantage of Genocchi polynomials having smaller coefficients of each individual term and relatively lesser terms compared to classical orthogonal polynomials. Thus, this is expected to provide us with smaller computational errors.

The rest of the paper is organized as follows: Section 2 introduces preliminary definitions and properties of the Caputo fractional derivative. Section 3 gives basic definitions and properties of , function approximation by , and the analytical expression of which is the integral of the product of two Genocchi polynomials. Section 4.1 derives the Genocchi polynomial operational matrix of Caputo’s fractional derivative. Section 4.2 shows the way to approximate the integral kernel in terms of Genocchi polynomials and derives the analytical expression of the kernel matrix. Section 5 presents the general approach of approximating (1) into a system of algebraic equations and using collocation methods to solve (1). Section 6 describes the general procedure of using Genocchi polynomials approximation to solve (1). Section 7 shows results of the proposed method. Section 8 sums up the findings of this paper.

2. Preliminaries

2.1. Fractional Calculus: Definitions

Fractional differentiation comes in different versions [11, 12]. In this paper, we consider Caputo’s fractional differentiation which provides a more realistic physical interpretation in real-life applications. Caputo’s fractional derivative operator of a function is defined as follows:

Definition 1.

Below are some properties of Caputo’s fractional derivatives:

where denotes the smallest integer greater than or equal to and denotes the largest integer less than or equal to .

Caputo’s fractional differential operator is linear:

where and are arbitrary constants.

3. Genocchi Polynomials and Function Approximation

3.1. Genocchi Polynomials: Definitions and Basic Properties

The Genocchi polynomials are defined as follows [13–15].

Definition 2. The Genocchi polynomial is defined by the generating function :

where is the Genocchi number. are the Bernoulli number and Bernoulli polynomial, respectively.

Below are some of the important properties of Genocchi polynomials:

3.2. Function Approximation by Genocchi Polynomials

We may approximate a function in terms of Genocchi polynomials by the following infinite series:

where are the Genocchi polynomials and the Genocchi coefficients .

In practice, we truncate the infinite series up to number of Genocchi polynomials according to the desired accuracy of the problem into the following truncated Genocchi series:

where in matrix notation,

where is the Genocchi coefficient vector.

is the Genocchi vector.

In this case, the Genocchi coefficients may be computed as follows:

3.3. Computing Genocchi Coefficients by Matrix Approach

Equation (10) fails to work for functions that are not differentiable at the points and . An example is shown below.

Let ; ;

To avoid this problem that occurs in the case of differentiation, we may evaluate the Genocchi coefficients using Theorem 4. To prove Theorem 4, we first prove the following Theorem 3 which gives the analytical expression of the integral of the product of two Genocchi polynomials over an arbitrary interval , which will be used in the subsequent part of the paper.

Theorem 3. Given any two Genocchi polynomials , , for

where is the Genocchi number, , are the falling and rising factorial, respectively. In particular, we have the following for ,

Proof. Using the following expression:

We obtain

Continuing this relation recursively for times, we arrive at

Adopting the techniques used in [16] for function approximation, we prove the following theorem of function approximation using Genocchi polynomials.

Theorem 4. Let and be the set of Genocchi polynomials up to order . Let . Since is a finite dimensional closed subspace of , then is the unique best approximation in Genocchi polynomials such that can be approximated by unique coefficients :

The Genocchi coefficient matrix consisting of the unique coefficients is given by the following:

where and as derived in Theorem 3.

Proof.

Let . Then,

Thus, we have a system of equations which is written in matrix representation as follows:

Therefore, the Genocchi coefficients matrix is

where can be obtained from Theorem 3.

4. Main Result

4.1. Genocchi Polynomial Operational Matrix of Caputo’s Fractional Differentiation

In this section, we derive the analytical expression of the Genocchi polynomial operational matrix of Caputo’s fractional derivative, which is the matrix , where

To derive , we first prove the following Lemmas 5 and 6.

Lemma 5. Caputo’s fractional derivative of fractional order of a Genocchi polynomial of order is given by

Proof. For ,

where . Substituting

where is the Beta function for .

Lemma 6. The matrix defined as is given by the following:

where

Proof. From Lemma 5,

Now, we can prove the following theorem that computes the operational matrix .

Theorem 7. Given a set of , , of Genocchi polynomials, the Genocchi polynomial operational matrix of Caputo’s fractional derivative of order over the interval is the matrix and is given by

where is given in Lemma 6 and with elements given in Theorem 3.

Proof. From Lemma 5,

Here, we justify the better accuracy of the new operational matrix of Caputo’s fractional derivative due to Genocchi polynomials by comparing its errors and that obtained when using Shifted Legendre operational matrix. This may be due to the fact that Genocchi polynomials have some advantages over classical orthogonal polynomials in approximating function, such as lesser terms and smaller coefficients of individual terms. We define the absolute error and relative error respectively, as follows:

where is either the Genocchi polynomial or Shifted Legendre polynomial . In Figure 1, we compare the relative error obtained due to Genocchi polynomial operational matrix of fractional derivative (GPOMFD) and the Shifted Legendre operational matrix of fractional derivative (SLOMFD) derived in [17]. One clearly sees from the figure that the errors for GPOMFD become smaller than that of SLOMFD for and their difference gets larger as increases over the interval . Further, we also show the absolute error and relative error in Tables 1 and 2, respectively, for at different points.

Table 1: Comparison of absolute error between SLOMFD and GPOMFD with at points .

SLOMFD

GPOMFD

1
0.113121
0.001562
0.005656
2
0.113121
0.001306
0.005656

3
0.696638
0.010167
0.036679
4
0.003583
0.000429
0.000549

5
1.836044
0.030079
0.107732
6
0.008628
0.000272
0.001283

7
3.773422
0.077279
0.273317
8
0.048392
0.004996
0.007218

9
7.468363
0.257570
0.899849
10
0.440694
0.045400
0.065718

Table 2: Comparison of relative error between SLOMFD and GPOMFD with at points .

SLOMFD

GPOMFD

1
0.0061268
0.0009516
0.0025063
2
0.0061268
0.0009516
0.0025063

3
0.0095486
0.0256030
0.0081265
4
0.0016094
0.0000698
0.0003038

5
0.0235837
0.0171498
0.0159125
6
0.0014129
0.0000535
0.0002389

7
0.1314007
0.5765935
0.0302776
8
0.0014150
0.0000529
0.0002372

9
8.7021303
0.1871596
0.0797470
10
0.0014147
0.0000528
0.0002369

Figure 1: Comparison of relative error of operational matrix of fractional derivative (OMFD) with between Shifted Legendre polynomials and Genocchi polynomials.

4.2. Approximation of Integral Kernel by Genocchi Polynomials

The integral kernel function will be approximated by Genocchi polynomials for each variable as follows:

where we refer to as the kernel matrix by Genocchi polynomial. Theorem 8 provides the formula for obtaining the kernel matrix .

Theorem 8. Let be a two-variable continuous function in .
Then, can be approximated in terms of Genocchi polynomials up to order , that is, , where is the Genocchi polynomial basis matrix and is the integral kernel matrix in Genocchi basis given by

where

Proof. From ,

where .
Define . Thus,

5. Error Analysis

Following the analysis made in [18], we derive the error estimates of numerical approximation using Genocchi polynomials considering the two following scenarios.

(1) If the exact solution of (1) is known, and is the numerical approximate solution using Genocchi polynomials up to order , then we define the error estimation as follows:

where . In the following theorem, we derive a bound for the error estimation of our Genocchi polynomial method.

Theorem 9. Let be an arbitrary function. If is approximated by truncated Genocchi polynomial series and ( is finite), then the error in is norm; for order approximation is bounded above by

Proof.

Therefore,

Obviously, the estimation would become more precise if the order increases.

(2) In usual cases, the exact solution is unknown. So, first we obtain the approximate solution of order , . Then, we obtain the approximate solution of order , . Define the error estimation as follows:

For a given desired error , one is able to determine the smallest number which gives the desired error.

Theorem 10. Suppose and are numerical solutions for the unknown solution of (1) with order , respectively. The error estimation defined as is convergent.

Proof.

and are bounded by upper bound given in Theorem 9 and, in turn, are both convergent. Thus, is convergent too.

6. Solving Fractional Integro-Differential Equation by Genocchi Polynomial Approximation

As in any operational matrix method, we approximate each term in the following fractional integro-differential equation (FIDE) with Genocchi polynomials.

For the FIDE in (1), we approximate each term using a set of Genocchi polynomials up to order . The functions and , which are known a priori in the FIDE, will not be approximated since this adds additional computational error to the numerical solution.

After approximation by Genocchi polynomial, the FIDE written in the matrix representation becomes the following:

where and .

Then, we select equally spaced points as the collocation points in the interval ; that is,

After collocating (47), we obtain a system of algebraic equations consisting of Genocchi coefficients . Given initial conditions, , , we select equations out of the algebraic equations and solve numerically. Thus, we obtain the approximate solution of order .

7. Numerical Examples

We present some examples of FIDEs which are solved using the proposed method of Genocchi polynomial approximation and Genocchi polynomial operational matrix of fractional derivative (GPOMFD). As in [8], we compute the maximum error function of the approximate function of order for each numerical solution obtained as follows:

where is the exact solution and is the approximate solution of order . Also, we compute the absolute error of order approximation at a particular point as follows:

Example 1 (integer case). Given

with initial conditions , . . The exact solution is . Using , that is, , the Genocchi polynomial method gives us the exact solution .

Example 2 (integer case). As for Example 3 in [8], given

with initial condition , the exact solution is . Figure 2 shows both and over the interval . Tables 3 and 4 show the numerical results of the exact solutions and approximate solutions () by Genocchi polynomial for and approximate solutions () by Bernoulli matrix [8] together with absolute errors for both methods. For comparison, Genocchi polynomial approximation of (i.e., degree ) is compared with the Bernoulli matrix method of order . Therefore, the Genocchi method of order (i.e., degree ) is compared with the Bernoulli method of .

Table 3: Comparison of absolute errors for Example 2 using orders of Genocchi polynomials and Bernoulli method [8].

Exact
Genocchi
Genocchi
Bernoulli [8]

0.0
0.0000000000000000
0.0000000000000000

0.1
0.1105170918075647
0.1105170918090843

0.2
0.2442805516320338
0.2442805516333059

0.3
0.4049576422728009
0.4049576422741780

0.4
0.5967298790565080
0.5967298790578598

0.5
0.8243606353500640
0.8243606353514914

0.6
1.0932712802343040
1.0932712802357070

0.7
1.4096268952293330
1.4096268952309410

0.8
1.7804327427939730
1.7804327427938810

0.9
2.2136428000412540
2.2136427999488390

1.0
2.7182818284590450
2.7182818284628360

Table 4: Comparison of absolute errors for Example 2 using orders of Genocchi polynomials and Bernoulli method [8].

Exact
Genocchi
Genocchi
Bernoulli [8]

0.0
0.0000000000000000
0.0000000000000000

0.1
0.1105170918075647
0.1105170918075627

0.2
0.2442805516320338
0.2442805516320308

0.3
0.4049576422728009
0.4049576422727964

0.4
0.5967298790565080
0.5967298790565028

0.5
0.8243606353500640
0.8243606353500572

0.6
1.0932712802343040
1.0932712802342920

0.7
1.4096268952293330
1.4096268952293210

0.8
1.7804327427939730
1.7804327427939550

0.9
2.2136428000412540
2.2136428000412350

1.0
2.7182818284590450
2.7182818284590240

Figure 2: Comparison of exact solution and approximate solution for Example 2.

Example 3 (integer case). We consider the following problem as in [8]:

with the initial condition . The exact solution is . Table 5 shows the error of the Genocchi approximate solutions for order and that of Bernoulli method [8] for the corresponding order of .

Table 5: Comparison of error for Example 3 using orders of Genocchi polynomials and that of Bernoulli method [8].

Genocchi

Bernoulli [8]

7

—
—

8

7

9

8

Example 4 (fractional case). We consider the following problem in [10].

with the initial conditions and . The exact solution is . Table 6 shows the numerical results of the Genocchi approximate solutions, for , and Legendre wavelet method (LWM) approximate solution for [10]. Though the exact solution is a polynomial, one is unable to obtain the exact solution. This is because the operational matrix does not give exact approximation of fractional derivative , that is,

Table 6: Comparison of approximate solutions for Example 4 using order Genocchi polynomials and Legendre wavelet method (LWM) [10] for and .

Exact
Genocchi
Genocchi
LWM
LWM

[10]

0.000
0.000000
0.000003

0.000024

0.125
0.017578
0.015382

0.015822

0.250
0.078125
0.074814

0.075531

0.375
0.193359
0.189006

0.192880

0.500
0.375000
0.364240

0.361498

0.625
0.634765
0.608990

0.622950

0.750
0.984375
0.936805

0.930901

0.875
1.435540
1.365460

1.391340

Example 5 (fractional case). We consider the following problem in [10].

with the initial conditions and . The exact solution is . Table 7 shows the numerical results of the Genocchi approximate solutions for and for and Legendre wavelet method (LWM) approximate solution () in [10]. We manage to obtain comparable errors for Genocchi polynomials with that of Legendre wavelets which are known to be a more powerful tool with a higher accuracy than its counterpart, Legendre polynomials.

Table 7: Comparison of absolute errors for Example 5 using order of Genocchi polynomials with respect to Legendre wavelet method (LWM , ) of [10].

Genocchi
Genocchi
LWM [10]

0.000

0

0.125

0.250

0.375

0.500

0.625

0.750

0.875

Example 6 (fractional case and solution is unknown). We consider the following problem:

with the initial conditions and . The solution is unknown. Using Theorem 10 which follows the concept given in [18], we calculate for shown in Table 8. Table 9 shows , the so-called th level error estimate [18].
Based on the results of Table 8 which is calculated using Theorem 10, one should choose if is to be achieved, where is called the tolerance.

Table 8: of Example 6 for of Genocchi polynomials.

0.295460
0.000466
0.000461
0.000337
0.000318

Table 9: of Example 6 for of Genocchi polynomials.

0.00
0.000000
0.000000
0.000000
0.000000
0.000000

0.10
0.027424
0.000156
0.000099
0.000047
0.000004

0.20
0.092740
0.000364
0.000091
0.000073
0.000142

0.30
0.174125
0.000393
0.000120
0.000267
0.000237

0.40
0.254577
0.000185
0.000379
0.000340
0.000223

0.50
0.321910
0.000189
0.000507
0.000290
0.000243

0.60
0.368759
0.000577
0.000435
0.000276
0.000349

0.70
0.392579
0.000786
0.000268
0.000416
0.000406

0.80
0.395643
0.000630
0.000251
0.000605
0.000322

0.90
0.385042
0.000020
0.000653
0.000503
0.000363

1.00
0.372688
0.001202
0.001572
0.000118
0.000917

8. Conclusion

In this paper, FIDE is solved with a numerical approximation scheme using Genocchi polynomial and Genocchi polynomial operational matrix of fractional derivative (GPOMFD). To overcome the problem that occurs for the classical formula of Genocchi coefficients in the case of nondifferentiable function, we use the matrix approach to determine the Genocchi coefficients of an arbitrary integrable function. Also, we derive the analytical expression of which is the integral of the product of Genocchi polynomials, which is the integral kernel matrix, and which is the Genocchi polynomial operational matrix of Caputo’s fractional derivatives. The approximation of Genocchi polynomial operational matrix of fractional derivative (GPOMFD) is shown to achieve smaller error than that of Shifted Legendre polynomial operational matrix of fractional derivative (SLOMFD). This is due to the fact that Genocchi polynomials have smaller coefficients for each -term than the Shifted Legendre polynomial. Applying Genocchi polynomial approximation, the FIDE is transformed into a system of algebraic equations of Genocchi coefficients of the unknown function . Using collocation method of equally spaced points over the interval together with the initial conditions given, we are able to solve FIDE with good accuracy. Comparison with the Bernoulli matrix method [8] and the Legendre wavelet method [10] indicates that Genocchi polynomial can achieve comparable results as those of the aforementioned methods which are known to have better accuracy than Legendre polynomial. Estimation errors occur when one tries to solve the FIDE in the case when the derivative is of fractional order.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this manuscript.

Acknowledgments

This work was supported in part by Malaysia UTHM FRGS Grant (Vot 1433).

References

M. H. Akrami, M. H. Atabakzadeh, and G. H. Erjaee, “The operational matrix of fractional integration for shifted Legendre polynomials,” Iranian Journal of Science and Technology, vol. 37, no. 4, pp. 439–444, 2013.

A. H. Bhrawy and A. S. Alofi, “The operational matrix of fractional integration for shifted Chebyshev polynomials,” Applied Mathematics Letters. An International Journal of Rapid Publication, vol. 26, no. 1, pp. 25–31, 2013.

A. H. Bhrawy, M. M. Tharwat, and M. A. Alghamdi, “A new operational matrix of fractional integration for shifted Jacobi polynomials,” Bulletin of the Malaysian Mathematical Sciences Society, Series 2, vol. 37, no. 4, pp. 983–995, 2014.

L. Zhu and Q. Fan, “Solving fractional nonlinear Fredholm integro-differential equations by the second kind Chebyshev wavelet,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 6, pp. 2333–2341, 2012.

S. S. Ray and R. K. Bera, “Analytical solution of the Bagley Torvik equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 168, no. 1, pp. 398–410, 2005.

P. Mokhtary and F. Ghoreishi, “The L²-convergence of the Legendre spectral Tau matrix formulation for nonlinear fractional integro differential equations,” Numerical Algorithms, vol. 58, no. 4, pp. 475–496, 2011.

L. Huang, X.-F. Li, Y. Zhao, and X.-Y. Duan, “Approximate solution of fractional integro-differential equations by Taylor expansion method,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 1127–1134, 2011.

A. H. Bhrawy, E. Tohidi, and F. Soleymani, “A new Bernoulli matrix method for solving high-order linear and nonlinear Fredholm integro-differential equations with piecewise intervals,” Applied Mathematics and Computation, vol. 219, no. 2, pp. 482–497, 2012.

M. H. T. Alshbool, A. S. Bataineh, I. Hashim, and O. R. Isik, “Solution of fractional-order differential equations based on the operational matrices of new fractional Bernstein functions,” Journal of King Saud University Science, vol. 28, pp. 119–202, 2015.

Z. Meng, L. Wang, H. Li, and W. Zhang, “Legendre wavelets method for solving fractional integro-differential equations,” International Journal of Computer Mathematics, vol. 92, no. 6, pp. 1275–1291, 2015.

I. Podlubny, Fractional Differential Equations : An Introduction to Fractional Derivatives, Fractional Differential Equations, vol. 198 of Some Methods of Their Solution and Some of Their Applications, Academic Press, 1998.

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and applications of fractional differential equations, vol. 204 of North-Holland Mathematics Studies, Elsevier Science, Amsterdam, The Netherlands, 1st edition, 2006.

S. Araci, “Novel identities for q-genocchi numbers and polynomials,” Journal of Function Spaces and Applications, vol. 2012, Article ID 214961, 13 pages, 2012.

S. Araci, “Novel identities involving Genocchi numbers and polynomials arising from applications of umbral calculus,” Applied Mathematics and Computation, vol. 233, pp. 599–607, 2014.

A. F. Horadam, “Genocchi polynomials,” in Proceedings of the 4th International Conference on Fibonacci Numbers and Their Applications, pp. 145–166, Kluwer Academic, Winston-Salem, NC, USA, 1991.

E. Keshavarz, Y. Ordokhani, and M. Razzaghi, “A numerical solution for fractional optimal control problems via Bernoulli polynomials,” Journal of Vibration and Control, vol. 22, no. 18, 2015.

A. Saadatmandi and M. Dehghan, “A new operational matrix for solving fractional-order differential equations,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1326–1336, 2010.

Y. Chen, X. Ke, and Y. Wei, “Numerical algorithm to solve system of nonlinear fractional differential equations based on wavelets method and the error analysis,” Applied Mathematics and Computation, vol. 251, pp. 475–488, 2015.