# Approximate Solution of Perturbed Volterra-Fredholm Integrodifferential Equations by Chebyshev-Galerkin Method.

1. IntroductionIntegrodifferential equations (IDEs) arise in many branches of mathematical sciences such as financial mathematics, mathematical modeling, and control theory, which are usually difficult to solve analytically, so numerical approach is required. Different methods have been used to obtain the solution of both linear and nonlinear IDEs such as Galerkin method [1-4], homotopy perturbation [5,6], Tau method [79], spline collocation method [10-12], Taylor collocation [1317], finite element approach [18-23], Legendre polynomials [24], variational iteration method [25], Haar wavelets [26], Krein's method [27], and Bessel collocation method [28], to mention a few.

The present work is motivated by the desire to obtain numerical solutions to initial value problems for integrodifferential equations via perturbed Chebyshev-Galerkin method. This paper is organized as follows. In Section 2, Chebyshev polynomial is discussed. In Section 3, preliminary steps towards application of the perturbed Chebyshev-Galerkin method are introduced. In Section 4, some numerical results are provided to demonstrate the efficiency and accuracy of using perturbed Chebyshev-Galerkin method and compared with those of [1,29] and lastly, Section 5 is the conclusion.

2. Chebyshev Polynomial

Chebyshev polynomials are widely used in applications in mathematics, mathematical physics, engineering, and computer science. Chebyshev polynomial is an orthogonal polynomial which satisfy the recurrence relation

[mathematical expression not reproducible]. (1)

In recent years, a lot of attention has been devoted to the study of Chebyshev methods to investigate various scientific models. Using these methods made it possible to solve differential equations of different forms [1, 7, 30-32], irrespective of the order of the differential equations.

In certain applications we need expressions for the products like [T.sub.i](x)[T.sub.j](x), which comes easily from

[T.sub.i](x)[T.sub.j](x) = 1/2[[T.sub.i+j](x) + [T.sub.i+j](x)]. (2)

3. Perturbed Chebyshev-Galerkin Method

We will consider the numerical solution of a class of linear Volterra-Fredholm integrodifferential initial value problems of the form

[mathematical expression not reproducible]. (3)

[[phi].sup.j](a) = [[gamma].sub.j], j = 0, 1, ..., m - 1, (4)

where [phi](x) is the unknown function, [p.sub.l](x), f(x), and [[tau].sub.k](x, t), k = 1, 2, are known functions, m is the order of (3), and [[lambda].sub.k], k = 1,2, are real numbers. Unless otherwise stated, x will always be the independent variable of the functions which appear throughout this paper and will be defined in a finite interval [a,b]. Moreover suppose that [[phi].sub.n](x) be the approximate solution of degree n to [phi](x), so we write

[[phi].sub.n] (x) = [n.summation over (i=0)] [[beta].sub.i][T.sub.i] (2x - (a + b)/b - a), (5)

where ft are determined by adding perturbation terms

[H.sub.n](x) = [m.summation over (l=0)] [[mu].sub.l+1][T.sub.n-m+l+1] (2x - (a + b)/b - a) (6)

to the right hand side of (3); we obtain

[mathematical expression not reproducible]. (7)

Applying Chebyshev-Galerkin approach discussed in [1], that is, multiplying both sides of (7) by Ts((2x - (a+b))/(b - a)), s = m, m + 1, ..., n + m + 1, and then integrating the resulting equation over interval [a, b], we obtain

[mathematical expression not reproducible]. (8)

From (8), we have

[theta][beta] = F, (9)

where [theta] is a matrix of (n+2) x (n + m + 2), [beta] and F are column matrices of (n + 2) x 1, and the other equations are derived from the initial conditions (4); that is,

[mathematical expression not reproducible]. (10)

Substitute the values of [[beta].sub.0], [[beta].sub.1], ..., [[beta].sub.n] obtained from (9) and (10) in (5) to obtain the approximate solution of degree n.

4. Numerical Examples

Example 1 (see [1, 33]). Consider the Fredholm integrodifferential equation

[phi]'(x) = [phi](x) - 1/2x 1/1 + x - log (1 + x) + 1/[(log(2)).sup.2] [[integral].sup.1.sub.0] x 1 + t [phi](t) dt, 0 [less than or equal to] x [less than or equal to] 1 (11)

subject to initial condition

[phi](0) = 0 and whose exact solution is [phi](x) = log (1 + x). (12)

The absolute errors are tabulated in Table 1 at different n. Table 5 exhibits a comparison between the errors obtained by using the perturbed Chebyshev-Galerkin and using Galerkin method [1]. Figure 1 shows the absolute errors at different n

and Figure 1(b) presents the perturbed Chebyshev-Galerkin method and exact solutions.

Example 2 (see [1, 34]). Consider the Volterra integrodifferential equation

[phi]'(x) - [[integral].sup.x.sub.0] [phi](t) dt = 1 - 2x sin (x), 0 [less than or equal to] x [less than or equal to] 1, [phi](0) = 0, (13)

whose exact solution is [phi](x) = x cos(x).

The numerical results for the absolute errors are displayed in Table 2 for different values of n; comparison of maximum

absolute errors is tabulated in Table 5 while Figure 2 exhibits the perturbed Chebyshev-Galerkin and exact solutions and the maximum absolute errors at different n.

Example 3 (see [1, 34]). Consider the Volterra integrodifferential equation

[phi]' (x) + [[integral].sup.x.sub.0] t[phi](t)dt = -1 + 1/2 [x.sup.2] - x exp (x), 0 [less than or equal to] x [less than or equal to] 1, [phi](0) = 0, (14)

whose exact solution is [phi](x) = 1 - exp(x).

The computational results for the absolute errors are summarized in Table 3 for different values of n and comparison of maximum absolute errors are tabulated in Table 5 while Figure 3 exhibits the perturbed Chebyshev-Galerkin and exact solutions and the maximum absolute error at different n.

Example 4 (see [13, 29]). Consider Fredholm-Volterra integrodifferential equation

[mathematical expression not reproducible]. (15)

The exact solution is [phi](x) = [e.sup.x].

For numerical results see Tables 4 and 6 while Figure 4 displays the approximant and the maximum errors as compared with the results from the literature [29].

5. Conclusion

This paper has discussed how the perturbed Chebyshev-Galerkin method can be applied for obtaining solutions of integral and integrodifferential equations. The formulation and implementation of the scheme are illustrated. The proposed method was tested using some problems with results. This paper discussed how the integrodifferential equation with variable and constant coefficients can be solved using the perturbed Chebyshev-Galerkin method. Maple and Matlab had been used to obtain the approximate solution and for plotting the graphs, respectively. Numerical results demonstrate that our method is an accurate and reliable numerical technique for solving mth order integrodifferential and integral equations. Finally, because of the accuracy and simplicity of the elegant method presented in this study, we recommend its application in finding the approximate solution to integrodifferential and integral equations.

http://dx.doi.org/ 10.1155/2017/8213932

Competing Interests

The authors declare that they have no competing interests.

References

[1] J. Biazar and F. Salehi, "Chebyshev Galerkin method for integrodifferential equations of the second kind," Iranian Journal of Numerical Analysis and Optimization, vol. 6, no. 1, pp. 31-42, 2016.

[2] T. Jangveladze, Z. Kiguradze, and B. Neta, "Galerkin finite element method for one nonlinear integro-differential model," Applied Mathematics and Computation, vol. 217, no. 16, pp. 6883-6892, 2011.

[3] M. Fathy, M. El-Gamel, and M. S. El-Azab, "Legendre-Galerkin method for the linear Fredholm integro-differential equations," Applied Mathematics and Computation, vol. 243, pp. 789-800, 2014.

[4] J. Shen, "Efficient spectral-galerkin method I. direct solvers of second- and fourth-order equations using legendre polynomials," SIAM Journal on Scientific Computing, vol. 15, no. 6, pp. 1489-1505, 1994.

[5] M. Dehghan and F. Shakeri, "Solution of an integro-differential equation arising in oscillating magnetic fields using he's homotopy perturbation method," Progress in Electromagnetics Research, vol. 78, pp. 361-376, 2008.

[6] J. Saberi-Nadjafi and A. Ghorbani, "He's homotopy perturbation method: an effective tool for solving nonlinear integral and integro-differential equations," Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2379-2390, 2009.

[7] S. M. Hosseini and S. Shahmorad, "Tau numerical solution of Fredholm integro-differential equations with arbitrary polynomial

bases," Applied Mathematical Modelling, vol. 27, no. 2, pp. 145-154, 2003.

[8] S. M. Hosseini and S. Shahmorad, "Numerical solution of a class of integro-differential equations by the tau method with an error estimation," Applied Mathematics and Computation, vol. 136, no. 2-3, pp. 559-570, 2003.

[9] A. K. Jimoh and K. Issa, "Comparison of some numerical methods for the solution of fourth order Integro-Differential equations," Journal of the Nigerian Association of Mathematical Physics, vol. 28, pp. 115-122, 2014.

[10] N. M. Chuong and N. V. Tuan, "Spline collocation methods for Fredholm Volterra integro differential equations of high order," Vietnam Journal of Mathematics, vol. 25, no. 1, p. 1524,1997

[11] H. Brunner, "The approximate solution of initial-value problems for general Volterra integro-differential equations," Computing. Archives for Scientific Computing, vol. 40, no. 2, pp. 125-137, 1988.

[12] N. Minh and N. Van, "Spline collocation methods for Fredholm integro-differential equations ofsecond order," Acta Mathematica Vietnamica, vol. 20, no. 1, pp. 85-98, 1995.

[13] S. Yalcinbas and M. Sezer, "The approximate solution of high-order linear Volterra-Fredholm integro-differential equations in terms of Taylor polynomials," Applied Mathematics and Computation, vol. 112, no. 2-3, pp. 291-308, 2000.

[14] Z. Avazzadeh and M. Heydari, "Chebyshev polynomials for solving two dimensional linear and nonlinear integral equations of the second kind," Computational & Applied Mathematics, vol. 31, no. 1, pp. 127-142, 2012.

[15] Y. Huang and X.-F. Li, "Approximate solution of a class of linear integro-differential equations by Taylor expansion method," International Journal of Computer Mathematics, vol. 87, no. 6, pp. 1277-1288, 2010.

[16] S. Yalcinbas and M. Sezer, "A Taylor collocation method for the approximate solution of general linear Fredholm-Volterra integro-difference equations with mixed argument," Applied Mathematics and Computation, vol. 175, no. 1, pp. 675-690, 2006.

[17] K. Maleknejad, S. Sohrabi, and Y. Rostami, "Numerical solution of nonlinear Volterra integral equations of the second kind by using Chebyshev polynomials," Applied Mathematics and Computation, vol. 188, no. 1, pp. 123-128, 2007.

[18] C. Chen and T. Shih, Finite Element Methods for Integro Differential Equations, World Scientific Publishing, Singapore, 1997.

[19] M. R. Crisci, E. Russo, and A. Vecchio, "Stability results for onestep discretized collocation methods in the numerical treatment of Volterra integral equations," Mathematics of Computation, vol. 58, no. 197, pp. 119-134,1992.

[20] A. Makroglou, "Convergence of a block-by-block method for nonlinear Volterra integro-differential equations," Mathematics of Computation, vol. 35, no. 151, pp. 783-796,1980.

[21] P. Linz, "Linear multistep methods for Volterra integrodifferential equations," Journal of the Association for Computing Machinery, vol. 16, pp. 295-301,1969.

[22] S. H. Chang, "On certain extrapolation methods for the numerical solution of integro-differential equations," Mathematics of Computation, vol. 39, no. 159, pp. 165-171,1982.

[23] Z. Avazzadeh, M. Heydari, and G. B. Loghmani, "Numerical solution of Fredholm integral equations of the second kind by using integral mean value theorem," Applied Mathematical Modelling, vol. 35, no. 5, pp. 2374-2383, 2011.

[24] A. H. Khater, A. B. Shamardan, D. K. Callebaut, and M. R. Sakran, "Numerical solutions of integral and integrodifferential equations using Legendre polynomials," Numerical Algorithms, vol. 46, no. 3, pp. 195-218, 2007.

[25] N. Bildik, A. Konuralp, and S. Yalcinbas, "Comparison of Legendre polynomial approximation and variational iteration method for the solutions ofgeneral linear Fredholm integro-differential equations," Computers & Mathematics with Applications, vol. 59, no. 6, pp. 1909-1917, 2010.

[26] A. Shahsavara, "Numerical solution of linear Volterra and Fredholm integro differential equations using Haar wavelets," Mathematical Sciences Journal, vol. 6, pp. 85-96, 2010.

[27] M. H. Fahmy, M. A. Abdou, and E. I. Deebs, "Fredholm-Volterra integral equation of the first kind with potential kernel," Le Matematiche, vol. 55, no. 1, pp. 55-74, 2000.

[28] S. Yuzbasi, "Improved Bessel collocation method for linear Volterra integro-differential equations with piecewise intervals and application of a Volterra population model," Applied Mathematical Modelling, vol. 40, no. 9-10, pp. 5349-5363, 2016.

[29] S. Shahmorad, "Numerical solution of the general form linear Fredholm-Volterra integro-differential equations by the Tau method with an error estimation," Applied Mathematics and Computation, vol. 167, no. 2, pp. 1418-1429, 2005.

[30] K. Issa and R. B. Adeniyi, "Extension of generalized recursive tau method to non-linear ordinary differential equations," Journal of the Nigerian Mathematical Society, vol. 35, no. 1, pp. 18-24, 2016.

[31] P. Deulfhard, Newton Methods for Nonlinear Problems: Affine Invariant and Adaptive Algorithms, Springer, Berlin, Germany, 2004.

[32] M. H. Aliabadi and E. L. Ortiz, "A Tau method based on nonuniform space-time elements for the numerical simulation of solitons," Computers & Mathematics with Applications, vol. 22, no. 9, pp. 7-19,1991.

[33] A. Golbabai and S. Seifollahi, "Radial basis function networks in the numerical solution of linear integro-differential equations," Applied Mathematics and Computation, vol. 188, no. 1, pp. 427-432, 2007.

[34] A. Wazwaz, Linear and Nonlinear Integral Equations, Methods and Applications, Springer, 2011.

K. Issa (1) and F. Salehi (2)

(1) Department of Statistics and Mathematical Sciences, Kwara State University, PMB1530, Malete, Ilorin, Kwara State, Nigeria

(2) Department of Mathematics, Islamic Azad University, Darab Branch, Darab, Iran

Correspondence should be addressed to K. Issa; issakazeem@yahoo.com

Received 4 October 2016; Accepted 22 December 2016; Published 12 January 2017

Academic Editor: Liwei Zhang

Caption: Figure 1: The approximate solution for Example 1 and its absolute errors for different n.

Caption: Figure 2: The approximate solution for Example 2 and its absolute errors for different n.

Caption: Figure 3: The approximate solution for Example 3 and its absolute errors for different n.

Caption: Figure 4: The approximate solution for Example 4 and its absolute errors for different n.

TABLE 1: Absolute errors for Example 1 at different n. x n = 4 n = 8 n =12 0 0 0 0 0.1 7.8696e - 06 2.9221e - 08 3.9291e - 11 0.2 1.3565e - 04 7.2475e - 08 8.0337e - 12 0.3 1.9086e - 04 2.5687e - 08 7.0680e - 11 0.4 1.4192e - 04 2.7660e - 08 1.8044e - 11 0.5 3.6209e - 05 5.6989e - 08 6.6114e - 11 0.6 5.5319e - 05 1.2399e - 07 9.7743e - 11 0.7 7.6735e - 05 7.0373e - 08 6.7821e - 11 0.8 1.4665e - 05 6.0573e - 08 1.4743e - 10 0.9 8.1616e - 05 1.5553e - 07 1.2997e - 10 1.0 8.5267e - 05 1.5711e - 07 1.8751e - 10 TABLE 2: Absolute errors for Example 2 at different n. X n = 4 n = 8 n=12 0 0 0 0 0.1 2.2012e - 05 2.3599e - 10 5.9771e - 17 0.2 1.5422e - 04 2.2695e - 10 1.0337e - 16 0.3 2.2010e - 04 2.8800e - 11 9.2897e - 17 0.4 1.7484e - 04 4.1892e - 10 9.1015e - 17 0.5 5.7582e - 05 5.7543e - 11 1.2640e - 17 0.6 5.4134e - 05 2.8830e - 10 6.0745e - 17 0.7 8.7363e - 05 1.1093e - 10 1.2000e - 16 0.8 1.3365e - 05 3.5180e - 10 7.3259e - 17 0.9 1.1663e - 04 9.5773e - 11 8.8650e - 17 1.0 1.4003e - 04 1.4225e - 10 3.1275e - 17 TABLE 3: Absolute errors for Example 3 at different n. X n = 4 n = 8 n =12 0 0 0 0 0.1 5.1479e - 06 4.8686e - 11 8.9783e - 18 0.2 5.4286e - 05 5.0102e - 11 1.4736e - 17 0.3 8.1338e - 05 1.4487e - 13 1.4206e - 17 0.4 6.5339e - 05 8.6061e - 11 1.2362e - 17 0.5 1.7563e - 05 9.9855e - 12 1.4139e - 18 0.6 3.2690e - 05 7.0793e - 11 1.0652e - 17 0.7 5.3816e - 05 1.2859e - 11 1.6736e - 17 0.8 2.9273e - 05 6.8704e - 11 1.2782e - 17 0.9 2.3223e - 05 3.3909e - 11 1.1424e - 17 1.0 3.0485e - 05 1.5662e - 11 2.1410e - 18 TABLE 4: Absolute errors for Example 4 at different n. X n = 4 n = 8 n=12 0 0 0 0 0.1 2.4216e - 05 5.5536e - 09 1.1520e - 12 0.2 5.6674e - 05 4.4620e - 08 9.1630e - 12 0.3 7.5702e - 05 1.4986e - 07 3.0638e - 11 0.4 1.0165e - 04 3.5091e - 07 7.1708e - 11 0.5 1.8402e - 04 6.7437e - 07 1.3786e - 10 0.6 3.8724e - 04 1.1437e - 06 2.3382e - 10 0.7 7.7496e - 04 1.7782e - 06 3.6348e - 10 0.8 1.3927e - 03 2.5927e - 06 5.2990e - 10 0.9 2.2486e - 03 3.5978e - 06 7.3533e - 10 1.0 3.2922e - 03 4.8012e - 06 9.8128e - 10 TABLE 5: Comparison of absolute maximum errors for Examples 1-3 at different n. n Example 1 Present method Galerkin method [1] 4 1.9086e - 04 4.629e - 3 8 1.5711e - 07 5.611e - 6 12 1.8751e - 10 5.930e - 09 n Example 2 Present method Galerkin method [1] 4 2.2010e - 04 2.300e - 3 8 4.1892e - 10 1.382e - 8 12 1.2000e - 16 8.087e - 15 n Example 3 Present method Galerkin method [1] 4 8.1338e - 05 7.204e - 4 8 8.6061e - 11 2.311e - 9 12 1.6736e - 17 9.235e - 16 TABLE 6: Comparison of absolute maximum errors for Example 4 at different n. n Present method Tau method [29] 5 7.1247e - 04 3.19e - 03 10 1.0014e - 09 2.10e - 06 15 1.4365e - 17 7.50e - 11

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Issa, K.; Salehi, F. |

Publication: | Journal of Mathematics |

Article Type: | Report |

Date: | Jan 1, 2017 |

Words: | 2903 |

Previous Article: | Chromatic Numbers of Suborbital Graphs for the Modular Group and the Extended Modular Group. |

Next Article: | Differential Calculus on N-Graded Manifolds. |

Topics: |