ANA Advances in Numerical Analysis 1687-9570 1687-9562 Hindawi Publishing Corporation 868279 10.1155/2012/868279 868279 Research Article Solution of Nonlinear Volterra-Fredholm Integrodifferential Equations via Hybrid of Block-Pulse Functions and Lagrange Interpolating Polynomials Marzban Hamid Reza Hoseini Sayyed Mohammad Bermudez De Castro Alfredo 1 Department of Mathematical Sciences Isfahan University of Technology P.O. Box 8415683111, Isfahan Iran iut.ac.ir 2012 10 12 2012 2012 22 05 2012 08 11 2012 2012 Copyright © 2012 Hamid Reza Marzban and Sayyed Mohammad Hoseini. 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.

An efficient hybrid method is developed to approximate the solution of the high-order nonlinear Volterra-Fredholm integro-differential equations. The properties of hybrid functions consisting of block-pulse functions and Lagrange interpolating polynomials are first presented. These properties are then used to reduce the solution of the nonlinear Volterra-Fredholm integro-differential equations to the solution of algebraic equations whose solution is much more easier than the original one. The validity and applicability of the proposed method are demonstrated through illustrative examples. The method is simple, easy to implement and yields very accurate results.

1. Introduction

Integral and integrodifferential equations have many applications in various fields of science and engineering such as biological models, industrial mathematics, control theory of financial mathematics, economics, electrostatics, fluid dynamics, heat and mass transfer, oscillation theory, queuing theory, and so forth .

It is well known that it is extremely difficult to analytically solve nonlinear integrodifferential equations. Indeed, few of these equations can be solved explicitly. So it is required to devise an efficient approximation scheme for solving these equations. So far, several numerical methods are developed. The solution of the first order integrodifferential equations has been obtained by the numerical integration methods such as Euler-Chebyshev  and Runge-Kutta methods .

Moreover, a differential transform method for solving integrodifferential equations was introduced in . Shidfar et al.  applied the homotopy analysis method for solving the nonlinear Volterra and Fredholm integrodifferential equations. As a concrete example, we can express the mathematical model of cell-to-cell spread of HIV-1 in tissue cultures considered by Mittler et al. . Yalcinbas and Sezer  proposed an approximation scheme based on Taylor polynomials for solving the high-order linear Volterra-Fredholm integrodifferential equations of the following form: (1.1)j=0qfj(t)y(j)(t)=f(t)+λ10tk1(t,s)y(s)ds+λ20tfk2(t,s)y(s)ds,0t,stf. Maleknejad and Mahmoudi  developed a numerical method by using Taylor polynomials to solve the following type of nonlinear Volterra-Fredholm integrodifferential equations: (1.2)j=0qfj(t)y(j)(t)=f(t)+λ10tk1(t,s)[y(s)]pds+λ20tfk2(t,s)y(s)ds,0t,stf. Darania and Ivaz  suggested an efficient analytical and numerical procedure for solving the most general form of nonlinear Volterra-Fredholm integrodifferential equations (1.3)j=0qfj(t)y(j)(t)=f(t)+λ10tk1(t,s)g1(s,y(s))ds+λ20tfk2(t,s)g2(s,y(s))ds, under the mixed conditions (1.4)j=0q-1(αijy(j)(0)+βij  y(j)(tf))=μi,i=0,1,,q-1, where λ1 and λ2 are constants and 0t,stf. Moreover, f(t), k1(t,s), k2(t,s), g1(s,y(s)), g2(s,y(s)), and fj(t), j=0,1,,q, are functions that have suitable derivatives on the interval 0t,stf.

These kinds of equations can be found in numerous applications such as electrodynamics, electromagnetic, biomechanics, and elasticity . Akyüz and Sezer  presented a Taylor polynomial approach for solving high-order linear Fredholm integrodifferential equations in the most general form. Streltsov  developed an effective numerical method based on the use of Chebyshev and Legendre polynomials for solving Fredholm integral equations. Babolian et al.  suggested an effective direct method to determine the numerical solution of the specific nonlinear Volterra-Fredholm integrodifferential equations. Their approach was based on triangular functions. In [16, 17], the variational iteration method (VIM) was employed for solving integral and integrodifferential equations. In addition, iterative and noniterative methods for the solution of nonlinear Volterra integrodifferential equations were presented and their local convergence was proved. The iterative methods provide a sequence solution and make use of fixed-point theory whereas the noniterative ones result in series solutions and also make use of the fixed-point principles . In recent years, the meshless methods have gained more attention not only by mathematicians but also in the engineering community. In , the meshless moving least square method was employed for solving nonlinear Fredholm integrodifferential equations.

In the present paper, we introduce a new numerical method for solving (1.3) with the mixed conditions (1.4). The method is based on a hybrid of block-pulse functions and Lagrange interpolating polynomials. The associated operational matrices of integration and product together with the Kronecker property of hybrid functions are then used to reduce the solution of (1.3) and (1.4) to the solution of nonlinear algebraic equations whose solution is much more easier than the original one.

The paper is organized as follows. In Section 2, we describe the basic properties of hybrid functions. Section 3 is devoted to the solution of the nonlinear Volterra-Fredholm integrodifferential equations. Finally, in Section 4, various types of integrodifferential equations are given to demonstrate the efficiency and the accuracy of the proposed method.

2. Hybrid Functions 2.1. Properties of Hybrid Functions

In order to define hybrid functions, we first divide the interval [0,tf) into N equidistant subintervals [((n-1)/N)  tf,(n/N)  tf), for n=1,2,,N. The hybrid of block-pulse functions and Lagrange interpolating polynomials, bnm(t), n=1,2,,N, m=0,1,,M-1, are defined on the interval [0,tf) as, , (2.1)bnm(t)={Lm(2Ntf  t-2n+1),t[(n-1N)  tf,nN  tf),0,  otherwise, where n and m are the orders of block-pulse functions and Lagrange interpolating polynomials, respectively. Here, Lm(t) are defined in  as (2.2)Lm(t)=i=0imM-1(t-τiτm-τi),m=0,1,,M-1, in which τi, i=0,1,,M-1 are the zeros of Legendre polynomial of order M, with the Kronecker property (2.3)Lm(τi)=δmi={1,if  i=m,0,if  im. where δmi is the Kronecker delta function. No explicit formulas are known for the points τi. However, they can be computed numerically using existing subroutines . Since bnm(t) consists of block-pulse functions and Lagrange interpolating polynomials, which are both complete and orthogonal, the set of hybrid of block-pulse functions and Lagrange interpolating polynomials is a complete orthogonal set in the Hilbert space 2[0,tf). A more detailed error analysis concerning this type of hybrid functions was discussed in .

2.2. Function Approximation

A function f(t)2[0,tf) may be expanded as (2.4)f(t)=n=1m=0fnmbnm(t).

If the infinite series in (2.4) is truncated, then (2.4) can be written as (2.5)f(t)n=1Nm=0M-1fnmbnm(t)=FTB(t), where (2.6)F=[f10,,f1(M-1),,fN0,,fN(M-1)]T,B(t)=[b10(t),,b1(M-1)(t),,bN0(t),,bN(M-1)(t)  ]T. In (2.6), fnm, n=1,2,,N, m=0,1,,M-1 are the expansion coefficients of the function f(t) in the nth subinterval [((n-1)/N)tf,(n/N)tf), and bnm(t), n=1,2,,N, m=0,1,,M-1 are defined in (2.1). In order to evaluate the coefficients fnm, we use Gaussian nodes which are the zeros of Legendre polynomial of order M. Let τm, m=0,1,,M-1 be the Gaussian nodes which are defined in [-1,1] and let tnm, n=1,2,,N, m=0,1,,M-1 be the corresponding Gaussian nodes in the nth subinterval [((n-1)/N)  tf,(n/N)  tf]. The relation between τm and tnm is given by (2.7)tnm=tf2N(τm+2  n-1). With the aid of (2.3) we have (2.8)bnm(tij)={1,if  i=n,j=m,0,otherwise. Using (2.8), the coefficients fnm, can be obtained as (2.9)fnm=f(tnm).

Now, let g(t,s) be a function of two independent variables defined for t[0,tf] and s[0,sf]. Then g(t,s) can be expanded in terms of hybrid functions as follows: (2.10)g(t,s)BT(t)GB(s), where G is a matrix of order NM×NM whose elements can be calculated for i,n=1,2,,N, and j,m=0,1,,M-1 as follows: (2.11)G(n-1)M+m,(i-1)M+j=g(tnm,sij), in which sij can be determined by the following relation: (2.12)sij=sf2N(τj+2i-1). The integration of the vector B(t) can also be expanded in terms of hybrid functions as (2.13)0tB(s)  dsPB(t), where P is the NM×NM operational matrix of integration and given by, , (2.14)P=[EHHH0EHH00EH000E], in which E and H can be obtained in the following manner.

Let (2.15)E=(eij),H=(hij), then, for i,j=0,1,,M-1, we have (2.16)eij=tf2N-1τjLi(t)dt,hij=tf2N-11Li(t)dt, where τj, j=0,1,,M-1 are the zeros of Legendre polynomial of order M. It is noted that E is the associated operational matrix of integration for Lagrange interpolating polynomials on interval [((n-1)/N)tf,(n/N)  tf].

Integrating (2.5) from 0 to t and using (2.13) imply (2.17)0tf(s)dsFTPB(t). From (2.13), by using the interpolating property of hybrid functions we get (2.18)0tnmf(s)dsFTPB(tnm)=i=1n-1j=0M-1fijhjm+j=0M-1fnjejm.

2.3. The Integration of the Cross-Product of Two Hybrid Function Vectors

The following property of the cross-product of two hybrid function vectors will also be used. Integrating the square matrix B(t)BT(t) from 0 to tf, and using the orthogonality of hybrid functions, we get (2.19)W=0tfB(t)BT(t)dt, where W is a block diagonal matrix of order NM×NM given by (2.20)W=[A000A000A], moreover A is also a diagonal matrix of order M×M given as (2.21)A=diag[a0,a1,,aM-1]. According to the orthogonality of hybrid functions, the values of am, for m=0,1,,M-1, can be computed as follows: (2.22)am=tf2N-11Lm2(t)dt.

2.4. The Operational Matrix of Derivative

Let d/dt denote the derivative operator. By expanding (d/dt)B(t) in terms of B(t) we get (2.23)ddtB(t)=D*B(t). In (2.23), D* is called the operational matrix of derivative for hybrid functions which is a block diagonal matrix of order NM×NM and is given by (2.24)D*=[D000D000D], where D is also a matrix of order M×M and the (i,j) th element of this matrix denoted by [D]ij for i,j=0,1,,M-1 can be calculated as follows: (2.25)[D]ij=dLidt    (tij)=2Ntf  Li(2Ntf  tij-2j+1)=2Ntf  Li(τj). Let f(1)(t) denote the first derivative of f(t) with respect to t. By expanding f(1)(t) in terms of hybrid functions and using (2.5) and (2.23) we obtain (2.26)f(1)(t)n=1Nm=0M-1f(1)(tnm)bnm(t)=FTD*B(t). According to (2.8) we get (2.27)f(1)(tnm)=FTD*B(tnm)=l=0M-1fnl[D]lm. Let f(j)(t) denote the jth derivative of f(t) with respect to t. In order to obtain the relation between the expansion coefficients of f(j)(t) and f(t), we use the following fact: (2.28)djdtjB(t)=(D*)jB(t), where (D*)j is the jth power of the matrix D* given in (2.23). Equation (2.28) can easily be obtained from (2.23) by differentiating f(t), j times with respect to t, j=2,3,,q. It should be noted that the matrix (D*)j has the following structure: (2.29)(D*)j=[Dj000Dj000Dj]. By expanding f(j)(t) in terms of hybrid functions and using (2.5) and (2.28) we get (2.30)f(j)(t)n=1Nm=0M-1f(j)(tnm)  bnm(t)=FT(D*)jB(t). According to (2.8) it immediately follows that (2.31)f(j)(tnm)=FT(D*)jB(tnm)=l=0M-1fnl[Dj]lm. where [Dj]lm denotes the (l,m) th element of the matrix Dj.

3. Solution of the Nonlinear Volterra-Fredholm Integrodifferential Equations

Consider the nonlinear Volterra-Fredholm integrodifferential equation given in (1.3). Since the set of hybrid functions is a complete orthogonal set in the Hilbert space 2[0,tf), we can expand any function in this space in terms of hybrid functions.

Let (3.1)y(t)=n=1Nm=0M-1ynmbnm(t)=YTB(t), where (3.2)ynm=y(tnm). Using (2.30), for j=1,2,,q, we have (3.3)y(j)(t)=n=1Nm=0M-1y(j)(tnm)bnm(t). where (3.4)y(j)(tnm)=l=0M-1ynl[Dj]lm. Moreover, by expanding k2(t,s) in terms of hybrid functions we get (3.5)k2(t,s)=BT(t)K2B(s), where K2 is a matrix of order NM×NM and can be calculated similar to matrix G in (2.10). From (2.5) and (2.9) we obtain (3.6)g2(t,y(t))=n=1Nm=0M-1g2(tnm,y(tnm))bnm(t)=n=1Nm=0M-1g2(tnm,ynm)bnm(t)=BT(t)G2, where G2 is a vector of order NM×1. From (2.19), (3.5) and (3.6), we have (3.7)01k2(t,s)g2(s,y(s))ds=01BT(t)K2B(s)BT(s)G2ds=BT(t)V, where (3.8)V=K2WG2. It should be noted that V is a vector of order NM×1. Substituting (3.1)–(3.7) in (1.3) yields (3.9)j=0qfj(t)y(j)(t)=f(t)+λ10  tk1(t,s)g1(s,y(s))ds+λ2BT(t)V. We now collocate (1.3) at NM points. For a suitable collocation points, we use Gaussian nodes tnm,n=1,2,,N, m=0,1,,M-1, which are introduced in (2.7). Therefore we have (3.10)j=0qfj(tnm)y(j)(tnm)=f(tnm)+λ10  tnmk1(tnm,s)g1(s,y(s))ds+λ2  BT(tnm)V. Also, using (2.8) we get (3.11)BT(tnm)V=V(n-1)M+m, where V(n-1)M+m denotes the ((n-1)M+m)th element of vector V. Consequently, with the aid of (2.18) and (3.10) we obtain (3.12)j=0qfj(tnm)y(j)(tnm)=f(tnm)+λ2V(n-1)M+m+λ1(i=1n-1j=0M-1k1(tnm,sij)g1(sij,yij)hjmasdasdaasdasda+j=0M-1k1(tnm,snj)g1(snj,ynj)ejm), where y(j)(tnm) for n=1,2,,N, m=0,1,,M-1 can be calculated from (3.4). Equation (3.12) gives NM nonlinear equations. Finally, we approximate the mixed conditions given in (1.4). From (1.4), (3.1), (3.3), and the definition of hybrid functions we obtain (3.13)j=0q-1(αijm=0M-1y(j)(t1m)b1m(0)+βijm=0M-1y(j)(tNm)bNm(tf))=μi,i=0,1,,q-1.

In order that the approximate solution obtained by the present method be continuous, we impose the continuity condition at the points (n/N)tf, n=1,2,,N-1. Equations. (3.12)-(3.13) together with the continuity conditions give N(M+1)+q-1, nonlinear equations that can be solved by using the well-known Tau method  for the unknown coefficients ynm, n=1,2,,N, m=0,1,M-1.

4. Illustrative Examples

In this section, four examples are given to demonstrate the efficiency, the accuracy, and the applicability of the proposed method.

Example 4.1.

Consider the nonlinear Volterra-Fredholm integrodifferential equation (4.1)t2y′′(t)+2y(t)=2-56  t+12te-t2+0  ttse-y2(s)ds+01ty2(s)ds,y(0)=0,y(0)=1, whose exact solution is given in  as y(t)=t. We applied the present method to solve this problem.

Let us define (4.2)E=(0tf(ye(t)-y(t))2dt)1/2, where y(t) and ye(t) denote the approximate solution obtained by the present method and the exact solution, respectively. In Table 1, the computational results of the 2-norm error denoted by E between the approximate solution and the exact solution for different values of N and M are given. The simulation results reported in Table 1, indicate that there is an excellent agreement between the approximate and exact solution.

The computational results of 2-norm error for Example 4.1.

Present method E
N = 2 , M=8 9.28 e - 12
N = 2 , M=9 6.32 e - 12
N = 2 , M=10 1.85 e - 13
N = 4 , M=8 1.11 e - 13
N = 4 , M=9 6.26 e - 15
N = 4 , M=10 9.29 e - 17
Example 4.2.

As the second example, consider the nonlinear Volterra integrodifferential equation (4.3)y(t)=2t-12sin(t4)+0tt2scos(t2y(s))ds,y(0)=0,whose exact solution is given in  as y(t) = t2. We applied the present method to solve this problem. In Table 2, the computational results of the 2-norm error between the approximate solution and the exact solution for different values of N and M are summarized. The simulation results obtained by the present method are superior to those reported in .

The computational results of 2-norm error for Example 4.2.

Present method E
N = 2 , M=8 8.83 e - 12
N = 2 , M=9 3.37 e - 12
N = 2 , M=10 2.15 e - 13
N = 4 , M=8 7.00 e - 14
N = 4 , M=9 4.63 e - 15
N = 4 , M=10 9.70 e - 17
Example 4.3.

Let us consider the nonlinear Volterra integrodifferential equation (4.4)y(4)(t)=1+0te-sy2(s)ds,t[0,1],y(0)=1,y(0)=1,y(1)=e,y(1)=e,whose exact solution is given in  as y(t)=et. We solved this problem by using hybrid functions. In Table 3, the computational results of y(t) obtained by the present method with N=2 and M=10 together with the exact solution and the solution determined by the Tau method  are summarized. The computational results obtained by the present method are superior to that reported in . Moreover in Table 4, the 2-norm error between the approximate solution obtained by the present method and the exact solution for different values of N and M is given. It should be pointed out that only a very small numbers of hybrid basis functions are needed to obtain a quite satisfactory approximation to the exact solution.

The computational results of y(t) for Example 4.3.

t Tau method Present method Exact
n = 10 N = 2 , M=10
0.0 1.0000000000 1.0000000000 1.0000000000
0.2 1.2214027513 1.2214027582 1.2214027582
0.4 1.4918246782 1.4918246976 1.4918246976
0.6 1.8221187743 1.8221188004 1.8221188004
0.8 2.2255409119 2.2255409285 2.2255409285
1.0 2.7182818285 2.7182818284 2.7182818284

The computational results of 2-norm error for Example 4.3.

Present method E
N = 2 , M=8 4.64 e - 08
N = 2 , M=9 2.94 e - 10
N = 2 , M=10 3.68 e - 11
N = 4 , M=8 3.01 e - 09
N = 4 , M=9 4.52 e - 12
N = 4 , M=10 5.96 e - 13
Example 4.4.

As the last example, consider the nonlinear Volterra-Fredholm integrodifferential equation (4.5)y′′′(t)-ty′′(t)=13t3-2t+1-0ts  y(s)ds+012tesy(s)-sds,y(0)=1,y(0)=1,y(1)=1+e, which has the exact solution y(t)=t+et. To measure the accuracy of proposed approach, the 2-norm error between the approximate solution obtained by the present method for different values of N and M and the exact solution is given in Table 5. The computational results reveal that the method is convergent.

The computational results of 2-norm error for Example 4.4.

Present method E
N = 2 , M=8 2.18 e - 08
N = 2 , M=9 6.15 e - 10
N = 2 , M=10 1.31 e - 11
N = 4 , M=8 6.42 e - 10
N = 4 , M=9 9.28 e - 12
N = 4 , M=10 9.82 e - 14
5. Conclusion

In the present work, a hybrid approximation method was developed for solving the high-order nonlinear Volterra-Fredholm integrodifferential equations. The nice properties of hybrid functions together with the associated operational matrices of integration, derivative, and cross-product are used to reduce the solution of problem to the solution of nonlinear algebraic equations whose solution is much more easier than the original one. The matrices P, C, and D* are sparse, hence making the method computationally attractive without sacrificing the accuracy of the solution. Moreover, It was shown that small values for N and M are needed to achieve high accuracy and a satisfactory convergence. The numerical results support this claim. The method is general, easy to implement and yields the desired accuracy in a few terms of hybrid basis functions. The simulation results indicate the convergence and the effectiveness of the proposed approach.

Acknowledgments

The authors are very grateful to the referees for their comments and valuable suggestions which improved the paper.

Polyanin A. D. Manzhirov A. V. Handbook of Integral Equations 2008 2nd Boca Raton, Fla, USA Chapman & Hall/CRC 10.1201/9781420010558 2404728 ZBL1217.35166 Van der Houwen P. J. Sommeijer B. P. Euler-Chebyshev methods for integro-differential equations Applied Numerical Mathematics 1997 24 2-3 203 218 10.1016/S0168-9274(97)00021-4 1464724 ZBL0881.65141 Enright W. H. Hu M. Continuous Runge-Kutta methods for neutral Volterra integro-differential equations with delay Applied Numerical Mathematics 1997 24 2-3 175 190 10.1016/S0168-9274(97)00019-6 1464722 ZBL0876.65089 Darania P. Ebadian A. A method for the numerical solution of the integro-differential equations Applied Mathematics and Computation 2007 188 1 657 668 10.1016/j.amc.2006.10.046 2327392 ZBL1121.65127 Shidfar A. Molabahrami A. Babaei A. Yazdanian A. A series solution of the nonlinear Volterra and Fredholm integro-differential equations Communications in Nonlinear Science and Numerical Simulation 2010 15 2 205 215 10.1016/j.cnsns.2009.03.015 2572561 ZBL1221.65343 Mittler J. E. Sulzer B. Neumann A. U. Perelson A. S. Influence of delayed viral production on viral dynamics in HIV-1 infected patients Mathematical Biosciences 1998 152 2 143 163 2-s2.0-0032168279 10.1016/S0025-5564(98)10027-5 Yalcinbas S. Sezer M. The approximate solution of high-order linear Volterra-Fredholm integro-differential equations in terms of Taylor polynomials Applied Mathematics and Computation 2000 112 2-3 291 308 10.1016/S0096-3003(99)00059-4 1758715 ZBL1023.65147 Maleknejad K. Mahmoudi Y. Taylor polynomial solution of high-order nonlinear Volterra-Fredholm integro-differential equations Applied Mathematics and Computation 2003 145 2-3 641 653 10.1016/S0096-3003(03)00152-8 2009290 ZBL1032.65144 Darania P. Ivaz K. Numerical solution of nonlinear Volterra-Fredholm integro-differential equations Computers & Mathematics with Applications 2008 56 9 2197 2209 10.1016/j.camwa.2008.03.045 2466741 ZBL1165.65404 Bloom F. Asymptotic bounds for solutions to a system of damped integro-differential equations of electromagnetic theory Journal of Mathematical Analysis and Applications 1980 73 2 524 542 10.1016/0022-247X(80)90297-8 564002 Holmåker K. Global asymptotic stability for a stationary solution of a system of integro-differential equations describing the formation of liver zones SIAM Journal on Mathematical Analysis 1993 24 1 116 128 10.1137/0524008 1199530 ZBL0767.45005 Forbes L. K. Crozier S. Doddrell D. M. Calculating current densities and fields produced by shielded magnetic resonance imaging probes SIAM Journal on Applied Mathematics 1997 57 2 401 425 10.1137/S0036139995283110 1438760 Akyüz A. Sezer M. A Taylor polynomial approach for solving high-order linear Fredholm integro-differential equations in the most general form International Journal of Computer Mathematics 2007 84 4 527 539 10.1080/00207160701227848 2326310 ZBL1118.65129 Streltsov I. P. Application of Chebyshev and Legendre polynomials on discrete point set to function interpolation and solving Fredholm integral equations Computer Physics Communications 2000 126 1-2 178 181 10.1016/S0010-4655(99)00520-2 1760663 ZBL0963.65143 Babolian E. Masouri Z. Hatamzadeh-Varmazyar S. Numerical solution of nonlinear Volterra-Fredholm integro-differential equations via direct method using triangular functions Computers & Mathematics with Applications 2009 58 2 239 247 10.1016/j.camwa.2009.03.087 2535790 ZBL1189.65306 Bildik N. Konuralp A. Yalçınbaş S. Comparison of Legendre polynomial approximation and variational iteration method for the solutions of general linear Fredholm integro-differential equations Computers & Mathematics with Applications 2010 59 6 1909 1917 10.1016/j.camwa.2009.06.022 2595965 ZBL1189.65307 Wang S.-Q. He J.-H. Variational iteration method for solving integro-differential equations Physics Letters A 2007 367 3 188 191 10.1016/j.physleta.2007.02.049 2570733 ZBL1209.65152 Ramos J. I. Iterative and non-iterative methods for non-linear Volterra integro-differential equations Applied Mathematics and Computation 2009 214 1 287 296 10.1016/j.amc.2009.03.067 2541068 ZBL1202.65179 Dehghan M. Salehi R. The numerical solution of the non-linear integro-differential equations based on the meshless method Journal of Computational and Applied Mathematics 2012 236 9 2367 2377 10.1016/j.cam.2011.11.022 2879705 ZBL1243.65154 Marzban H. R. Hoseini S. M. Razzaghi M. Solution of Volterra's population model via block-pulse functions and Lagrange-interpolating polynomials Mathematical Methods in the Applied Sciences 2009 32 2 127 134 10.1002/mma.1028 2478910 ZBL1156.65106 Canuto C. Hussaini M. Y. Quarteroni A. Zang T. A. Spectral Methods in Fluid Dynamics 1987 New York, NY, USA Springer Marzban H. R. Tabrizidooz H. R. Razzaghi M. A composite collocation method for the nonlinear mixed Volterra-Fredholm-Hammerstein integral equations Communications in Nonlinear Science and Numerical Simulation 2011 16 3 1186 1194 10.1016/j.cnsns.2010.06.013 2736626 ZBL1221.65340 Ortiz E. L. The tau method SIAM Journal on Numerical Analysis 1969 6 480 492 0258287 10.1137/0706044 ZBL0195.45701 Ebadi G. Rahimi-Ardabili M. Y. Shahmorad S. Numerical solution of the nonlinear Volterra integro-differential equations by the tau method Applied Mathematics and Computation 2007 188 2 1580 1586 10.1016/j.amc.2006.11.024 2335652 ZBL1119.65123