Abstract and Applied Analysis
Volume 2018 (2018), Article ID 3569139, 6 pages
https://doi.org/10.1155/2018/3569139

Research Article

A Novel Method for Solving Nonlinear Volterra Integro-Differential Equation Systems

Mohammad Hossein Daliri Birjandi, Jafar Saberi-Nadjafi, and Asghar Ghorbani

Department of Applied Mathematics, School of Mathematical Sciences, Ferdowsi University of Mashhad, Mashhad, Iran

Correspondence should be addressed to Jafar Saberi-Nadjafi; najafi141@gmail.com

Received 12 December 2017; Accepted 2 April 2018; Published 3 June 2018

Academic Editor: Patricia J. Y. Wong

Copyright © 2018 Mohammad Hossein Daliri Birjandi 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.

An efficient iteration method is introduced and used for solving a type of system of nonlinear Volterra integro-differential equations. The scheme is based on a combination of the spectral collocation technique and the parametric iteration method. This method is easy to implement and requires no tedious computational work. Some numerical examples are presented to show the validity and efficiency of the proposed method in comparison with the corresponding exact solutions.

1. Introduction

Systems of integro-differential equations and their solutions play a pivotal role in the fields of science, industrial mathematics, control theory of financial mathematics, and engineering [1–3]. Physical systems, such as biological applications in population dynamics, and genetics where impulses arise naturally or are caused by control are modeled by a system of integro-differential equations [4, 5]. The initial value problem for a nonlinear system of integro-differential equations were used to model the competition between tumor cells and the immune system [6]. In [7], two systems of specific inhomogeneous integro-differential equations are studied in order to examine the noise term phenomenon. Thus applications of numerical methods for solving these equations are attractive. This has led to a great deal of research in recent years with the use of numerical methods such as the variational iteration method [8], differential transform method [9], Bezier curves method [10], radial basis function networks [11], biorthogonal systems [12], the block pulse functions method [13], and a collocation method in combination with operational matrices of Bernstein polynomials [14].

The parametric iteration method (PIM) is an analytic approximate method that provides the solution of linear and nonlinear problems as a sequence of iterations. In fact, the PIM as a fixed- point iteration method is a reconstruction of the variational iteration method. Since the implementation of the PIM generally leads to the calculation of unneeded terms, where more time is consumed in repeated calculations for series solutions, so to overcome these shortcomings, a useful improvement of the PIM was proposed in [15].

The aim of this work is to present an effective algorithm, requiring no tedious computational work, based on the improved PIM and the spectral collocation technique to obtaining an accurate solution for the following system of Volterra integro-differential equations as follows:

for and . The functions are given real valued functions, ,  , are kernels of the integral equations, and are linear or nonlinear functions of .

To demonstrate the utility of the proposed method, some examples of system of Volterra integro-differential equations are given, which are solved using the established method. The obtained results are compared with the exact solutions. In all cases, the present algorithm performed excellently.

2. Parametric Iteration Method

The PIM gives rapidly convergence by using successive approximations of the exact solution if such a solution exists; otherwise the approximations can be used for numerical aims. To convey the basic idea of this method, we first consider (1) as below:

where denotes the auxiliary linear operator with respect to . In (2) is a nonlinear continuous operator with respect to and is the source term.

Next, we construct a family of explicit iterative processes for (2) [15, 16]

where

with the initial conditions .

Also, we can construct a family of the implicit PIM for (2) as follows:

with the above initial condition.

is the initial guess which can be freely found from solving its corresponding linear equation ( or ) and the subscript denotes the th iteration. The parameter denotes the so-called auxiliary parameter, which can be identified easily and efficiently by the technique proposed in [15]. Also we are free to choose the auxiliary linear operator , the auxiliary parameter , and the initial approximation , which is fundamental to the validity and flexibility of the PIM. Accordingly, the successive approximations of , for the PIM iterative relation will be obtained readily in the auxiliary parameter . Consequently, the exact solution can be obtained by using the following:

When the original PIM fails, then the presence of the parameter in (3) or (5) can play an important role in the frame of the PIM. However, we can always discover a valid region of for every physical problem by plotting the solution or its derivatives versus the parameter in some points. An approximate optimal value of the convergence accelerating parameter can be determined at the order of approximation by the residual error [17]:

Now, one can minimize (7) by imposing the requirement .

3. Description of the Method

The PIM procedure provides the solution of the system of Volterra integro-differential equations as a sequence of iterates; its successive iterations may be very complex so that the resulting integrals in its iterative relation may not be performed analytically. Here, we will overcome this shortcoming of the PIM for solving (1) by proposing a spectral collocation PIM. As will be shown in this paper later, the new method will be very simple to implement and save time and calculations.

Assume that ,  , are the th component of the solution and also the th equation of system (1). Consider basis functions , which are polynomials of degree satisfying for the shifted Chebyshev nodes (note that and ):

The unknown function   () is approximated as a truncated series of polynomials. The polynomial

interpolates the points ,  , that is, , where and are vectors. The values of the interpolating polynomial of first derivative at the nodes are ,  , where the th element of the differentiation matrix is . Also the values of the integrals at the nodes are defined by ,  , where is the Volterra integration matrix [18, 19] (for more details observe [19]).

Generally, in order to solve system (1) using a spectral collocation scheme, the interpolating polynomials   () are required to satisfy the equations of the system at the interior nodes. The values of the interpolating polynomials at the interior nodes are , where denotes the row of the identity matrix and the derivative values are . The initial condition that involves the interpolating polynomials can be handled by using the formulas , where this notation denotes the last row of the identity matrix.

For the interpolating polynomial to satisfy the th differential equation of the system of Volterra integro-differential equations (1) at each interior node, the collocation equation

should be satisfied. Substituting the differentiation and the integration matrix relations into (10), we get

Now, in view of (3) and the definitions of and , by substituting the differentiation and integration matrix relations, we will have the following explicit PIM for solving (1) which is called the spectral PIM (SPIM):

If we define ,  , and , then we will have the following explicit iterative relation for finding the solution vector

Here the vector is defined as

Now one can start with the initial guess for obtaining the approximations .

4. Illustrative Examples

In this section, we give several test examples to confirm our analysis. To examine the accuracy of the results, and are employed to assess the efficiency of the method SPIM. All the computations were performed using software Matlab and terminated when the current iterate satisfies , where is the solution vector of the th SPIM iteration.

Example 1. Consider the following system of Volterra integro-differential equations [20]:

with the initial conditions and . The exact solution of this system is .

Table 1 illustrates the and errors for different values of as well as the number of iterations to reach the above-mentioned stopping criteria. Also Figure 1 shows the absolute error of the proposed method for . As expected, the exponential rate of convergence is observed for the system of nonlinear Volterra integro-differential equations, which confirmed our theoretical predictions.

Table 1: The and errors for Example 1.

No. of Itre.
23
21
31
17
21
25

Figure 1: Errors of and for Example 1 when .

Example 2. Consider the following system of Volterra integro-differential equations with the exact solution [8]:

subjected to initial conditions and .

The results of and errors of and for the different values of as well as the number of iterations to reach the stopping criteria are given in Table 2. Figure 2 depicts the absolute error of the presented method for . Again, the exponential rate of convergence is observed for the system of nonlinear Volterra integro-differential equations.

Table 2: The and errors for Example 2.

No. of Itre.
28
22
22
29
24
21

Figure 2: Errors of and for Example 2 when .

Example 3. Finally, let us consider the following system of nonlinear Volterra integro-differential equations [20]:

with the initial conditions of and whose exact solution is given by and .

The errors and number of iterations for the different values of are presented in Table 3. The other results are presented in Figure 3 with similar to Examples 1 and 2.

Table 3: The and errors for Example 3 when .

No. of Itre.
23
18
17
16
15
15

Figure 3: Errors of and for Example 3 when .

5. Conclusion

In this paper, we presented a powerful numerical approach based on a combination of the Chebyshev spectral collocation technique and the parametric iteration method for solving the linear and nonlinear system of Volterra integro-differential equations. This method inherits the strengths of the PIM and is easy to implement and is accurate when applied to the linear and nonlinear system of Volterra integro-differential equations. The comparison of the approximate solution and the exact solution reveals that the present method is very accurate and convenient for solving the linear and nonlinear system of Volterra integro-differential equations.

Conflicts of Interest

Mohammad Hossein Daliri Birjandi, Jafar Saberi-Nadjafi, and Asghar Ghorbani declare that there are no conflicts of interest regarding the publication of this paper.

References

M. A. Jaswon and G. T. Symm, Integral Equation Methods in Potential Theory and Elastostatics, Academic Press, London, UK, 1977.

S. Hatamzadeh-Varmazyar, M. Naser-Moghadasi, E. Babolian, and Z. Masouri, “Numerical approach to survey the problem of electromagnetic scattering from resistive strips based on using a set of orthogonal basis functions,” Progress in Electromagnetics Research, vol. 81, pp. 393–412, 2008.

A. Kyselka, “Properties of systems of integro-differential equations in the statistics of polymer chains,” Polymer Science U.S.S.R., vol. 19, no. 11, pp. 2852–2858, 1977.

S. Kheybari, M. T. Darvishi, and A. M. Wazwaz, “A semi-analytical algorithm to solve systems of integro-differential equations under mixed boundary conditions,” Journal of Computational and Applied Mathematics, vol. 317, pp. 72–89, 2017.

S. Kheybari, M. T. Darvishi, and A. M. Wazwaz, “A semi-analytical approach to solve integro-differential equations,” Journal of Computational and Applied Mathematics, vol. 317, pp. 17–30, 2017.

N. Bellomo, B. Firmani, and L. Guerri, “Bifurcation analysis for a nonlinear system of integro-differential equations modelling tumor-immune cells competition,” Applied Mathematics Letters, vol. 12, no. 2, pp. 39–44, 1999.

A.-M. Wazwaz, “The existence of noise terms for systems of inhomogeneous differential and integral equations,” Applied Mathematics and Computation, vol. 146, no. 1, pp. 81–92, 2003.

J. Saberi-Nadjafi and M. Tamamgar, “The variational iteration method: a highly promising method for solving the system of integro-differential equations,” Computers & Mathematics with Applications, vol. 56, no. 2, pp. 346–351, 2008.

A. Arikoglu and I. Ozkol, “Solutions of integral and integro—differential equation systems by using differential transform method,” Computers & Mathematics with Applications, vol. 56, no. 9, pp. 2411–2417, 2008.

F. Ghomanjani, A. Kiliçman, and S. Effati, “Numerical solution for IVP in volterra type linear integrodifferential equations system,” Abstract and Applied Analysis, vol. 2013, Article ID 490689, 2013.

A. Golbabai, M. Mammadov, and S. Seifollahi, “Solving a system of nonlinear integral equations by an {RBF} network,” Computers & Mathematics with Applications. An International Journal, vol. 57, no. 10, pp. 1651–1658, 2009.

M. I. Berenguer, A. I. Garralda-Guillem, and M. Ruiz Guiz Galán, “An approximation method for solving systems of Volterra integro-differential equations,” Applied Numerical Mathematics, vol. 67, pp. 126–135, 2013.

K. Maleknejad, H. Safdari, and M. Nouri, “Numerical solution of an integral equations system of the first kind by using an operational matrix with block pulse functions,” International Journal of Systems Science, vol. 42, no. 1, pp. 195–199, 2011.

K. Maleknejad, B. Basirat, and E. Hashemizadeh, “A Bernstein operational matrix approach for solving a system of high order linear Volterra-Fredholm integro-differential equations,” Mathematical and Computer Modelling, vol. 55, no. 3-4, pp. 1363–1372, 2012.

A. Ghorbani and J. Saberi-Nadjafi, “A piecewise-spectral parametric iteration method for solving the nonlinear chaotic Genesio system,” Mathematical and Computer Modelling, vol. 54, no. 1-2, pp. 131–139, 2011.

S. A. Alavi, A. Heydari, and F. Khellat, “On convergence and error analysis of the parametric iteration method,” International Journal of Applied Mathematical Research, vol. 4, no. 1, p. 129, 2014.

A. Ghorbani and M. Gachpazan, “A Spectral Quasilinearization Parametric Method for Nonlinear Two-Point Boundary Value Problems,” Bulletin of the Malaysian Mathematical Sciences Society.

J. A. Weideman and S. C. Reddy, “A {MATLAB} differentiation matrix suite,” ACM Transactions on Mathematical Software, vol. 26, no. 4, pp. 465–519, 2000.

T. Damian, Matrix Based Operatorial Approach to Differential and Integral Problems, Babes Bolyai University of Cluj-Napoca, Romania, 2011.

M. I. Berenguer, A. I. Garralda-Guillem et al., “An approximation method for solving systems of Volterra integro-differential equations,” Applied Numerical Mathematics, vol. 67, pp. 126–135, 2013.