Journal of Applied Mathematics

Volume 2017 (2017), Article ID 1510267, 5 pages

https://doi.org/10.1155/2017/1510267

## Solutions of First-Order Volterra Type Linear Integrodifferential Equations by Collocation Method

^{1}Department of Mathematics and Statistics, Federal Polytechnic, Ilaro, Ogun State, Nigeria^{2}Department of Mathematics, College of Science and Technology, Covenant University, Ota, Ogun State, Nigeria

Correspondence should be addressed to Timothy A. Anake

Received 28 November 2016; Revised 22 February 2017; Accepted 2 March 2017; Published 20 March 2017

Academic Editor: Mehmet Sezer

Copyright © 2017 Olumuyiwa A. Agbolade and Timothy A. Anake. 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

The numerical solutions of linear integrodifferential equations of Volterra type have been considered. Power series is used as the basis polynomial to approximate the solution of the problem. Furthermore, standard and Chebyshev-Gauss-Lobatto collocation points were, respectively, chosen to collocate the approximate solution. Numerical experiments are performed on some sample problems already solved by homotopy analysis method and finite difference methods. Comparison of the absolute error is obtained from the present method and those from aforementioned methods. It is also observed that the absolute errors obtained are very low establishing convergence and computational efficiency.

#### 1. Introduction

Integrodifferential equation is a hybrid of integral and differential equations which have found extensive applications in sciences and engineering since it was established by Volterra [1]. A special class of these equations are the Volterra type which have been used to model heat and mass diffusion processes, biological species coexisting together with increasing and decreasing rate of growth, electromagnetic theory, and ocean circulations, among others [2].

First-order integrodifferential equation (IDE) of the Volterra type is generally of the formwhereIn solving (1), we seek the unknown function given the kernel , a nonsingular function defined on with . This kernel determines the nature of the solutions of integral equation (2) depending on its type [3]. In this paper, only separable or degenerate kernels have been considered.

The theory and application of integrodifferential equations are important subjects in applied mathematics. The existence and uniqueness of the solutions of integrodifferential equations, usually discussed in terms of their kernel, had been established already in Linz [1]. Generally, methods for solving integrodifferential equations combine methods of solving both integral and differential equations. Also, since closed form solutions may not be tractable for most applications, numerical methods are employed to obtain approximations to the exact solutions.

Some numerical approaches in literature include iterative methods [4], successive approximation methods [5], and standard integral collocation approximation methods [6]. Other methods such as power series methods, where Chebyshev and Legendre’s polynomials are used as basis functions, have been applied to obtain solutions of some higher order IDE of linear type. Akyaz and Sezer [7], for instance, presented Chebyshev collocation method for solving linear integrodifferential equations by truncated Chebyshev series. Recently, Gegele et al. [8] used power and Chebyshev series approximation methods to find numerical solution to higher order linear Fredholm integrodifferential equations using collocation methods. The result presented showed that the methods can give accurate results when compared with the exact solution. These methods proved efficient in the respective applications from the results provided but they seem yet to be applied to integrodifferential equations of Volterra type.

It is our aim here to extend the approach in Gegele et al. [8] to obtain approximate solutions for integrodifferential equations of Volterra type.

In the next section, we shall discuss the derivation of our methods; then the implementation using some sample problems is presented in Section 3. Finally, in Section 4 we shall present the results and draw our conclusions.

#### 2. Methodology

In the sequel, the combination of the power series approximation and collocation method is employed for the solution of IDE of Volterra type.

To proceed, (1) is reduced to the formwhere and . The initial conditions (3b) are required in order to find particular solutions of (3a).

Now, let the solution of Volterra type IDE, (3a) and (3b), be analytic and therefore possess the power serieswhere s are monomial bases and s are real coefficients to be determined.

Substituting equation (4) into both sides of (3a) givesHence,where and are known functions.

For an arbitrary choice of , (6) is obtained as a linear algebraic equation in unknowns as follows:We note that is given by the initial condition (3b) while the remaining , , are to be determined by collocation method.

To generate the collocation points, we shall consider two methods, namely, the standard and Chebyshev-Gauss-Lobatto Collocation Methods, respectively.

##### 2.1. Standard Collocation Method (SCM)

This method is used to determine the desired collocation points within an interval, say, , and is given by

##### 2.2. Chebyshev-Gauss-Lobatto Collocation Method (CGLCM)

The collocation points are obtained as follows:Interestingly, Chebyshev-Gauss-Lobatto points have also been used as collocation and interpolation points in the solutions of optimal control problems governed by Volterra integrodifferential equations [9, 10].

Using either of the two collocation points to collocate (7) together with the initial conditions given in (3b) will result in a system of linear algebraic equations in unknowns. Hence, the resultant matrix problem is as follows:where , , are polynomials evaluated at each collocation point . The values of the unknowns can be obtained using any convenient method of solving matrix equations of the form , where is invertible.

Substituting the values of the , , obtained from (4) yields the approximate solution. We note that the accuracy level desired for the approximate solution is determined by the degree of the approximating polynomial.

#### 3. Results

In this section, standard and Chebyshev-Gauss-Lobatto collocation points have been employed, respectively, to solve sample problems as described in Section 2. The numerical solutions obtained using the present method had been compared with the exact solutions of the sample problems. Similarly, absolute errors of results from this present method have been compared with those obtained in Behrouz [11] by homotopy analysis method (HAM) and finite difference method (FDM) for the same problems.

The absolute error of computation is defined in all cases as follows:

*Problem 1. *Exact solution: .

Using SCM, we obtained the following approximate solutions:Similarly using CGLCM, we obtained the approximate solution as follows:The solutions obtained from the implementation of the method for Problem 1 using SCM and CGLCM are compared with the exact solution in Table 1. Also absolute errors obtained are compared with absolute errors obtained from HAM and FDM in Table 2.