Next Article in Journal / Special Issue
Optimal B-Spline Bases for the Numerical Solution of Fractional Differential Problems
Previous Article in Journal / Special Issue
On Partial Cholesky Factorization and a Variant of Quasi-Newton Preconditioners for Symmetric Positive Definite Matrices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Collocation Methods for Volterra Integral and Integro-Differential Equations: A Review

by
Angelamaria Cardone
1,*,
Dajana Conte
1,
Raffaele D’Ambrosio
2 and
Beatrice Paternoster
1
1
Department of Mathematics, University of Salerno, 84084 Fisciano, Italy
2
Department of Engineering and Computer Science and Mathematics, University of L’Aquila, Via Vetoio, Loc. Coppito, 67100 L’Aquila, Italy
*
Author to whom correspondence should be addressed.
Axioms 2018, 7(3), 45; https://doi.org/10.3390/axioms7030045
Submission received: 27 April 2018 / Revised: 18 June 2018 / Accepted: 19 June 2018 / Published: 1 July 2018
(This article belongs to the Special Issue Advanced Numerical Methods in Applied Sciences)

Abstract

:
We present a collection of recent results on the numerical approximation of Volterra integral equations and integro-differential equations by means of collocation type methods, which are able to provide better balances between accuracy and stability demanding. We consider both exact and discretized one-step and multistep collocation methods, and illustrate main convergence results, making some comparisons in terms of accuracy and efficiency. Some numerical experiments complete the paper.

1. Introduction

It is the purpose of this paper to illustrate recent results on collocation methods for Volterra integral equations (VIEs) and Volterra integro-differential equations (VIDEs), mainly due to the authors. Such equations model evolutionary problems with memory in many applications, such as dynamics of viscoelastic materials with memory, electrodynamics with memory, heat conduction in materials with memory [1,2,3,4,5,6]. The numerical solution of these equations has a high computational cost due both to the nonlinearity of the advancing term and to the evaluation of the lag term, which contains the past history of the solution. Therefore, a crucial point is finding accurate and efficient numerical methods.
Collocation methods have several desirable properties. They provide an approximation over the entire integration interval to the solution of the equation, which reveals to be quite useful in a variable-stepsize implementation: indeed, it is easy to recover the missing past values when the stepsize is changed, by evaluating the collocation polynomial. Other good properties of collocation methods are their high order of convergence, strong stability properties and flexibility. As a matter of fact, if some information is known on the behavior of the exact solution, then it is possible to choose the collocation functions in order to better follow such behavior, so giving rise to mixed collocation methods, see for example [7] in the case of ordinary differential equations (ODEs), and [8] in the case of VIEs. It is also worthwhile mentioning that collocation also has an important theoretical relevance: in fact, many numerical methods are difficult to be analyzed as discrete schemes while, re-casted as collocation-based methods, their analysis is reasonably simplified and can be carried out in a very elegant way. There is, however, a remarkable drawback of one-step collocation methods: they suffer from order reduction phenomenon when applied to stiff problems [9,10,11], since the order of convergence is not uniform (for instance, in the case of s-stage collocation based Runge-Kutta methods on Gauss-Legendre collocation points, the order is p = 2 s in the grid points, but it degenerates to p = s for stiff problems, since the order is s in the internal stages). Such a drawback is successfully solved by two-step collocation methods [12], having high uniform order on the overall integration interval. On the side of computational cost, collocation methods are usually more expensive than other classes of methods. In fact, a collocation method with m collocation parameters requires at each time-step the solution of a nonlinear system of dimension m. To face this drawback, multistep collocation methods can be adopted which increase the order of convergence at the same computational cost of one-step ones. When a collocation method is applied to an integral equation, several integrals must be computed, thus suitable quadrature rules are needed to complete the discretization, with the introduction of an additional error. Lastly, a reliable error estimation for collocation methods for integral equations is still missing: there have been some advances (compare [13] and references therein contained), however considerable work needs to be done.
One-step collocation methods first appeared in the literature and main results are collected in the monographs [2,3]. Recently, we have proposed multistep collocation methods [13,14,15,16] and two step almost collocation methods [13,17,18], where the collocation polynomial depends on the approximate solution in a fixed number of previous time steps, with the aim of increasing the order of convergence of classical one–step collocation methods, without additional computational cost at each time step, and at the same time obtaining highly stable methods. This idea has been already proposed for the numerical solution of ODEs [19,20,21] (see also [11], Section V.3), and afterward modified in [12], by also using the inherent quadratic technique [22,23,24]. We also underline that they have high uniform order, thus they do not suffer from the order reduction phenomenon, well-known in the ODEs context [9]. Other approaches, based on multistep collocation, have been proposed in [25,26,27,28,29,30,31,32].
Here we briefly introduce one-step collocation methods and illustrate with more detail the construction and analysis of multistep collocation methods for VIEs and VIDEs, with the aim of giving a complete idea on the recent developments in this context. We give practical indications on how to choose the quadrature formulas in the discretized methods for an efficient implementation. In this review, we consider VIEs and VIDEs with smooth kernel and solution. We illustrate methods with a uniform mesh, however they could easily be applied to a non-uniform mesh (compare [2] for one-step collocation methods).
The paper is organized as follows. Section 2 and Section 3 deal with one-step and multistep collocation methods for VIEs, respectively. Section 4 illustrates two-step almost collocation methods for VIEs. Section 5 and Section 6 focus on one-step and multistep collocation methods for VIDEs, respectively.

2. One Step Collocation Methods for VIES

We consider VIEs of the form
y ( t ) = g ( t ) + 0 t k ( t , τ , y ( τ ) ) d τ , t I = [ 0 , T ] ,
where k C ( D × R ) , with D : = ( t , τ ) : 0 τ t T , and g C ( I ) . In the literature, many authors (see [2,3] and references therein contained) have analyzed one step collocation methods for VIEs. As it is well known, a collocation method is based on the idea of approximating the exact solution of a given integral equation with a suitable function belonging to a chosen finite dimensional space, usually a piecewise algebraic polynomial which satisfies the integral equation exactly on a certain subset of the integration interval (called the set of collocation points).
Let us discretize the interval I by introducing a uniform mesh
I h = t n : = n h , n = 0 , , N , h 0 , N h = T .
The Equation (1) can be rewritten, by relating it to this mesh, as
y ( t ) = F n ( t ) + Φ n ( t ) t [ t n , t n + 1 ] ,
where
F n ( t ) : = g ( t ) + 0 t n k ( t , τ , y ( τ ) ) d τ
and
Φ n ( t ) : = t n t k ( t , τ , y ( τ ) ) d τ
represent respectively the lag term and the increment function.
Collocation methods provide an approximation P ( t ) to the solution y ( t ) of (1) on [ 0 , T ] , such that its restriction to each interval ( t n , t n + 1 ] is a polynomial:
P ( t ) | ( t n , t n + 1 ] = P n ( t ) .

2.1. Exact One-Step Collocation Methods

Let us fix m collocation parameters 0 c 1 < < c m 1 and denote by t n j = t n + c j h the collocation points. The collocation polynomial, restricted to the interval [ t n , t n + 1 ] , is of the form:
P n ( t n + s h ) = j = 1 m L j ( s ) Y n j s [ 0 , 1 ] n = 0 , , N 1
where L j ( s ) is the j-th Lagrange fundamental polynomial with respect to the collocation parameters and Y n j : = P n ( t n j ) . Exact collocation methods are obtained by imposing that the collocation polynomial (2) exactly satisfies the VIE (1) in the collocation points t n , i and by computing y n + 1 = P n ( t n + 1 ) :
Y n i = F n i + Φ n i y n + 1 = j = 1 m L j ( 1 ) Y n j ,
where
F n i = g ( t n i ) + h ν = 0 n 1 0 1 k ( t n i , t ν + s h , P ν ( t ν + s h ) ) d s
Φ n i = h 0 c i k ( t n i , t n + s h , P n ( t n + s h ) ) d s ,
i = 1 , , m . Note that the first equation in (3) represents a system of m nonlinear equations in the m unknowns Y n i . We recall that generally P ( t ) is not continuous in the mesh points, as
P ( t ) S m 1 ( 1 ) ( I h ) ,
where
S μ ( d ) ( I h ) = v C d ( I ) : v ( t n , t n + 1 ] Π μ , n = 0 , 1 , , N 1 .
Here, Π μ denotes the space of (real) polynomials of degree not exceeding μ .
The classical collocation methods have uniform order O ( h m ) for any choice of the collocation parameters, and can achieve local superconvergence in the mesh points by opportunely choosing the collocation parameters, i.e., order 2 m 2 with m Lobatto points or m 1 Gauss points with c m = 1 and order 2 m 1 with m Radau II points. The optimal superconvergence order O ( h 2 m ) in the mesh points can be achieved with Gauss nodes in the iterated collocation methods [2,3].

2.2. Discretized One-Step Collocation Methods

The collocation Equation (3) is not yet in a form amenable to numerical computation: another discretization step, based on quadrature formulas F ˜ n i F n i and Φ ˜ n i Φ n i for the approximation of (4) and (5) are needed in order to obtain the fully discretised collocation schemes, thus leading to Discretized collocation methods.
The discretized collocation polynomial is of the form
P ˜ n ( t n + s h ) = j = 1 m L j ( s ) Y ˜ n j s [ 0 , 1 ] n = 0 , , N 1
where Y ˜ n j : = P ˜ n ( t n j ) . The m unknowns Y ˜ n j are determined by imposing that the collocation polynomial (7) satisfies exactly the integral equation at the collocation points and by using quadrature formulas of the form
Φ ˜ n ( t n i ) = h l = 0 μ 0 w i l k ( t n i , t n + d i l h , P ˜ n ( t n + d i l h ) )
F ˜ n ( t n i ) = g ( t n i ) + h ν = 0 n 1 l = 0 μ 1 b l k ( t n i , t ν + ξ l h , P ˜ ν ( t ν + ξ l h ) ) ,
i = 1 , , m , for approximating the lag term (4) and the increment function (5). The Formulas (8) and (9) are obtained by using quadrature formulas of the form
( ξ l , b l ) l = 1 μ 1 , ( d i l , w i l ) l = 1 μ 0 , i = 1 , , m ,
where the quadrature nodes ξ l and d i l satisfy 0 ξ 1 < < ξ μ 1 1 and 0 d i 1 < < d i μ 0 1 , μ 0 and μ 1 are positive integers and w i l , b l are suitable weights.
The numerical method is then of the form:
Y ˜ n i = F ˜ n ( t n i ) + Φ ˜ n ( t n i ) y ˜ n + 1 = j = 1 m L j ( 1 ) Y ˜ n j ,
where Φ ˜ n ( t n i ) and F ˜ n ( t n i ) are given by (8) and (9).
Note that the first equation in (9) represents a system of m nonlinear equations in the m unknowns Y ˜ n i .
Such methods preserve, under suitable hypothesis on the quadrature Formulas (8) and (9), the same order of the exact collocation methods [3].
A collocation method for VIEs is equivalent to an implicit Runge-Kutta method for VIEs (VRK method) if and only if c m = 1 (see Theorem 5.2.2 of [3]). As the lag–term computation is the most expensive part in the numerical solution of VIEs, fast collocation and Runge–Kutta methods have been constructed for convolution VIEs of Hammerstein type [33,34] in order to reduce the computational effort in the lag–term computation. The stability analysis of collocation methods for VIEs can be found in [3,35] and the related bibliography.

3. Multistep Collocation Methods for VIEs

Multistep collocation methods for VIEs have been introduced in [16] by adding interpolation conditions in r previous step points, with the aim of increasing the uniform order of convergence of one step collocation methods without increasing the computational cost. The multistep collocation polynomial, restricted to the interval [ t n , t n + 1 ] , is of the form
P n ( t n + s h ) = k = 0 r 1 φ k ( s ) y n k + j = 1 m ψ j ( s ) Y n j s [ 0 , 1 ] ,
n = r , , N 1 , where again
Y n j : = P n ( t n j )
and φ k ( s ) , ψ j ( s ) are the following polynomials of degree m + r 1
φ k ( s ) = i = 1 m s c i k c i · i = 0 i k r 1 s + i k + i , ψ j ( s ) = i = 0 r 1 s + i c j + i · i = 1 i j m s c i c j c i .
The collocation parameters are assumed to satisfy c i c j and c 1 0 .

3.1. Exact Multistep Collocation

The exact multistep collocation methods are obtained by imposing that the collocation polynomial (12) exactly satisfies the VIE (1) at the collocation points t n i , and by computing y n + 1 = P n ( t n + 1 ) :
Y n i = F n i + Φ n i , y n + 1 = k = 0 r 1 φ k ( 1 ) y n k + j = 1 m ψ j ( 1 ) Y n j ,
where the lag–term F n i and increment–term Φ n i are given by (4) and (5) respectively. The r-step m-point exact collocation method (12)–(15) has uniform convergence order of at least p = m + r , for any choice of distinct collocation abscissas 0 < c 1 < < c m 1 , as stated in the following theorem proved in [16].
Theorem 1.
Let ε ( t ) = y ( t ) P ( t ) be the error of the exact collocation method (12)–(15) and p = m + r . Suppose that
i. 
the given functions describing the VIE (1) satisfy k C ( p ) ( D × R ) , g C ( p ) ( I ) .
ii. 
the starting error is ε , [ 0 , t r ] = O ( h p ) .
iii. 
ρ ( A ) < 1 , where
A = [ c ] c | c 0 r 1 , 1 I r 1 φ r 1 ( 1 ) φ r 2 ( 1 ) , , φ 0 ( 1 )
and ρ denotes the spectral radius.
Then
ε = O ( h m + r ) .
Moreover, a suitable choice of collocation parameters can ensure superconvergence in the mesh points, as pointed out in the following theorem [16].
Theorem 2.
Let us suppose that
  • the hypothesis of the Theorem 1 hold with p = 2 m + r 1 .
  • the collocation parameters c 1 , , c m are the solution of the system
    c m = 1 1 i + 1 k = 0 r 1 β k ( k ) i j = 1 m γ j ( c j ) i = 0 , i = m + r , . . . , 2 m + r 2
    with
    β k = 0 1 φ k ( s ) d s , γ j = 0 1 ψ j ( s ) d s
then
max n = 0 , , N | ε ( t n ) | = O ( h 2 m + r 1 ) .

3.2. Discretized Multistep Collocation

The discretized multistep collocation methods are obtained by using quadrature formulas of the form (8) to (9) for approximating the lag term and the increment function. The discretized multistep collocation polynomial, denoted by P ˜ n ( t ) , is then of the form
P ˜ n ( t n + s h ) = k = 0 r 1 φ k ( s ) y ˜ n k + j = 1 m ψ j ( s ) Y ˜ n j , s [ 0 , 1 ]
n = 0 , , N 1 ,where the functions φ k ( s ) and ψ j ( s ) are given by (14), and Y ˜ n j : = P ˜ n ( t n j ) are determined by the solution of the following nonlinear system
Y ˜ n i = F ˜ n i + Φ ˜ n i , y ˜ n + 1 = k = 0 r 1 φ k ( 1 ) y ˜ n k + j = 1 m ψ j ( 1 ) Y ˜ n j .
The following theorem [16] shows that, as in the exact case, the r-step m-point discretized collocation method (19) and (20) has convergence order of at least p = m + r , for any choice of distinct collocation abscissas 0 < c 1 < < c m 1 .
Theorem 3.
Let ε ˜ ( t ) : = y ( t ) P ˜ ( t ) be the error of the discretized collocation method (19) and (20) and let p = m + r . Suppose that
i. 
the given functions describing the VIE (1) satisfy k C ( p ) ( D ) , g C ( p ) ( I ) ;
ii. 
the lag–term and increment–term quadrature Formulas (10) are of order respectively at least p + 1 and p;
iii. 
the starting error is ε ˜ , [ 0 , t r ] = O ( h p ) .
iv. 
ρ ( A ) < 1 , where A is given by (16).
Then
ε ˜ = O ( h m + r ) .
An analogous result holds concerning the local superconvergence:
Theorem 4.
Let us suppose that
  • the hypothesis of the Theorem 3 hold with p = 2 m + r 1 .
  • the collocation parameters c 1 , , c m are the solution of the system (17).
Then
max n = 0 , , N | ε ˜ ( t n ) | = O ( h 2 m + r 1 ) .

4. Two Step Almost Collocation Collocation Methods for VIEs

Within the class of multistep collocation methods, although methods with unbounded stability regions exist, no A-stable methods have been found [16]. In order to determine A-stable methods, two step almost collocation (TSAC) methods have been introduced in [18] and further analyzed in [13,17].
The collocation polynomial P n ( t ) for TSAC methods is computed by employing the information about the equation on two consecutive steps:
P n ( t n + s h ) = φ 0 ( s ) y n 1 + φ 1 ( s ) y n + j = 1 m χ j ( s ) Y j [ n ] + j = 1 m ψ j ( s ) ( F j [ n ] + Φ j [ n + 1 ] ) ,
where Y j [ n ] = P ( t n 1 , j ) . Then the method assumes the form:
Y i [ n + 1 ] = φ 0 ( c i ) y n 1 + φ 1 ( c i ) y n + j = 1 m χ j ( c i ) Y j [ n ] + j = 1 m ψ j ( c i ) F j [ n ] + Φ j [ n + 1 ] , y n + 1 = φ 0 ( 1 ) y n 1 + φ 1 ( 1 ) y n + j = 1 m χ j ( 1 ) Y j [ n ] + j = 1 m ψ j ( 1 ) F j [ n ] + Φ j [ n + 1 ] ,
where F j [ n ] and Φ j [ n + 1 ] are suitable sufficiently high order quadrature formulae for the discretization of F [ n ] ( t n j ) and Φ [ n + 1 ] ( t n j ) respectively, assuming the form
F j [ n ] = g ( t n j ) + h ν = 1 n l = 0 m + 1 b l k t n j , t ν 1 , l , Y l [ ν ] ,
and
Φ j [ n + 1 ] = h l = 0 m + 1 w j l k t n j , t n l , Y l [ n + 1 ] .
In the quadrature Formulas (23) and (24) we mean t ν 1 , 0 = t ν 1 , t ν 1 , m + 1 = t ν , Y 0 [ ν ] = P n ( t ν 1 ) , Y m + 1 [ ν ] = P n ( t ν ) and t n 0 = t n . We observe as the method (22) requires, at each step, the solution of a nonlinear system of ( m + 1 ) d equations in the stage values Y i [ n + 1 ] and y n + 1 .
The basis functions φ 0 ( s ) , φ 1 ( s ) , χ j ( s ) and ψ j ( s ) , j = 1 , 2 , , m , are polynomials of degree p, determined from the continuous order conditions, according to the following theorem [18]:
Theorem 5.
Assume that the kernel k ( t , η , y ) and the function g ( t ) in (1) are sufficiently smooth. Then the method (21) and (22) has uniform order p, i.e.,
η ( t n + s h ) = O ( h p + 1 ) , h 0 ,
for s [ 0 , 1 ] , if the polynomials φ 0 ( s ) , φ 1 ( s ) , χ j ( s ) and ψ j ( s ) , j = 1 , 2 , , m satisfy the system of equations
1 φ 0 ( s ) φ 1 ( s ) j = 1 m χ j ( s ) j = 1 m ψ j ( s ) = 0 , s k ( 1 ) k φ 0 ( s ) j = 1 m ( c j 1 ) k χ j ( s ) j = 1 m c j k ψ j ( s ) = 0 ,
s [ 0 , 1 ] , k = 1 , 2 , , p , where
η ( t n + s h ) = y ( t n + s h ) φ 0 ( s ) y ( t n h ) φ 1 ( s ) y ( t n ) j = 1 m χ j ( s ) y ( t n + ( c j 1 ) h ) + ψ j ( s ) y ( t n + c j h ) .
is the local truncation error.
As regards the global error, the method has uniform order of convergence p * = min { l + 1 , q , p + 1 } , where l and q are the order of the starting procedure (for the computation of the starting values y 1 and Y i [ 1 ] , i = 1 , 2 , , m ) and the order of the quadrature Formulas (23) and (24) respectively (see Theorem 2.5 in [18]). Then we use as starting procedure a one step collocation method having uniform order of convergence l = p .
Two-step collocation methods are obtained by solving the system of order conditions up to the maximum uniform attainable order p = 2 m + 1 , and, in this way, all the basis functions are determined as the unique solution of such system. However, as observed in [18], it is not convenient to impose all the order conditions because it is not possible to achieve high stability properties (e.g., A-stability) without getting rid of some of them. Therefore, almost collocation methods have been introduced by relaxing a specified number r of order conditions, i.e., by a priori opportunely fixing r basis functions, and determining the remaining ones as the unique solution of the system of order conditions up to p = 2 m + 1 r . Within the class of TSAC methods, A-stable methods have been constructed in [18] by fixing one (case r = 1 ) or both (case r = 2 ) of the polynomials φ 0 ( s ) and φ 1 ( s ) as
φ 0 ( s ) = k = 1 m ( s c k ) ( q 0 + q 1 s + + q p m s p m ) , φ 1 ( s ) = k = 1 m ( s c k ) ( p 0 + p 1 s + + p p m s p m ) ,
where α j and β j , j = 0 , 1 , , p m , are free parameters, which have to be determined in order to obtain desired stability properties.
A error estimation of the local discretization error for TSAC methods has been derived in [13].
Example 1.
Let us consider the methods with two stages m = 2 and order p = 2 m = 4 . Classes of A-stable methods were derived in [13,18] by considering
φ 0 ( s ) = s ( s c 1 ) ( s c 2 ) ( q 0 + q 1 s ) ,
where c 1 , c 2 , q 0 , q 1 are free parameters. The weights in (23) and (24) were computed in [18] as
b = 1 + 2 c 1 + 2 c 2 6 c 1 c 2 12 c 1 c 2 1 2 c 2 12 c 1 ( c 1 1 ) ( c 1 c 2 ) 2 c 1 1 12 c 2 ( c 2 1 ) ( c 2 c 1 ) 3 + 4 c 1 + 4 c 2 6 c 1 c 2 12 ( c 1 1 ) ( c 2 1 ) , W = c 1 2 3 c 1 c 2 6 c 2 c 1 ( 2 c 1 3 c 2 ) 6 ( c 1 c 2 ) c 1 3 6 c 2 ( c 1 c 2 ) 0 c 2 2 3 c 1 c 2 6 c 1 c 1 3 6 c 1 ( c 1 c 2 ) c 2 ( 2 c 2 3 c 1 ) 6 ( c 1 c 2 ) 0 .
An A-stable method is obtained by choosing for example q 0 = 15 / 10 , q 1 = 1 , c 1 = 0.9 , c 2 = 0.95 , see [13].

4.1. Diagonally Implicit TSAC Methods for VIEs

The computational cost associated to the solution of the nonlinear system (22) can be reduced by making the coefficient matrix have a structured shape, e.g., lower triangular or diagonal. This strategy, in the field of Runge–Kutta methods for ODEs, leads to the raise of the famous classes of Diagonally Implicit and Singly Diagonally Implicit Runge-Kutta methods (DIRK and SDIRK), see [10,11] and bibliography therein contained. Moreover, in the field of collocation-based methods for ODEs, an analogous strategy has been applied, obtaining TSAC methods having structured coefficient matrix [12].
In fact, a lower triangular matrix allows to solve the equations in m successive stages, with only a d-dimensional system to be solved at each stage. Moreover, if all the elements on the diagonal are equal, in solving the nonlinear systems by means of Newton-type iterations, one may hope to use repeatedly the stored L U factorization of the Jacobian. If the structure is diagonal, the problem reduces to the solution of m independent systems of dimension d, and can therefore be solved in a parallel environment.
Methods of this type have been derived in [17], where first of all it was assumed w j , m + 1 = 0 , j = 1 , , m , in such a way that (22) becomes a nonlinear system of dimension m d only depending on the stage values Y i [ n + 1 ] , i = 1 , , m , and assumes the following form
Y i [ n + 1 ] h j = 1 m l = 1 m ψ j ( c i ) w j l k ( t n j , t n l , Y l [ n + 1 ] ) = B i [ n ] , y n + 1 = P n ( t n + 1 ) ,
where
B i [ n ] = φ 0 ( c i ) y n 1 + φ 1 ( c i ) y n + j = 1 m χ j ( c i ) Y j [ n ] + j = 1 m ψ j ( c i ) F j [ n ] + h j = 1 m ψ j ( c i ) w j 0 k ( t n j , t n , y n ) .
By defining
Y [ n + 1 ] = Y 1 [ n + 1 ] , Y 2 [ n + 1 ] , , Y m [ n + 1 ] T , B [ n ] = B 1 [ n ] , B 2 [ n ] , , B m [ n ] T , Ψ = ψ j ( c i ) i , j = 1 m ,
W = w j l j , l = 1 m , K ( t n c , t n c , Y [ n + 1 ] ) = K ( t n i , t n j , Y j [ n + 1 ] ) i , j = 1 m ,
the nonlinear system in (28) takes the form
Y [ n + 1 ] h ( Ψ I ) ( W I ) · K ( t n c , t n c , Y [ n + 1 ] ) e = B [ n ] ,
where · denotes the usual Hadamard product, I is the identity matrix of dimension d and e is the unit vector of dimension m d . The tensor form (30) clearly shows as the matrices which determine the structure of the nonlinear system (28) are Ψ and W. In [17] a strategy was described to obtain lower triangular or diagonal structures for the matrices Ψ and W: in particular a quadrature formula of the form
0 c j f ( s ) d s w j 0 f ( 0 ) + l = 1 m w ˜ j l f ( c l 1 ) + l = 1 j w j l f ( c l ) ,
was proposed for the increment
Φ [ n + 1 ] ( t n j , P ( · ) ) = h 0 c j k ( t n j , t n + s h , P n ( t n + s h ) ) d s ,
in addition to the quadrature formula
0 1 f ( s ) d s b 0 f ( 0 ) + l = 1 m b l f ( c l ) + b m + 1 f ( 1 ) ,
for the approximation of the lag term
F [ n ] ( t n j , P ( · ) ) = g ( t n j ) + h ν = 1 n 0 1 k ( t n j , t ν 1 + s h , P ν 1 ( t ν 1 + s h ) ) d s .
We observe that in Formula (31), in case of triangular structure, w ˜ j l = 0 , l = 1 , , j while, in case of diagonal structure, w ˜ j 1 = 0 and w j l = 0 , l = 1 , , j 1 . The determination of the weights in Formulas (31) and (33) has been described in [17].
Assuming that Ψ and W are lower triangular, we obtain the diagonally implicit TSAC methods (DITSAC)
Y i [ n + 1 ] h ψ i ( c i ) w i i k ( t n i , t n i , Y i [ n + 1 ] ) = B i [ n ] + B ˜ i [ n ] + h l = 1 i 1 j = l i ψ j ( c i ) w j l k ( t n j , t n l , Y l [ n + 1 ] ) , y n + 1 = φ 0 ( 1 ) y n 1 + φ 1 ( 1 ) y n + j = 1 m χ j ( 1 ) Y j [ n ] + j = 1 m ψ j ( 1 ) F j [ n ] + Φ j [ n + 1 ] ,
where B i [ n ] is given by (29),
B ˜ i [ n ] = h j = 1 i l = 1 m ψ j ( c i ) w ˜ j l k ( t n j , t n 1 , l Y l [ n ] ) ,
and F j [ n ] , Φ j [ n + 1 ] are approximations of (34) by means of the quadrature Formulas (31) and (33).

4.2. Numerical Results

We present some numerical results which confirm that, differently from one step collocation methods, the TSAC methods do not suffer form the order reduction in the integration of stiff systems, as we expect from the uniform order of convergence stated in Theorem 5. In order to illustrate this phenomenon, we show the results obtained on both a non stiff and a stiff equation:
  • the non stiff VIE
    y ( t ) = 2 cos ( t ) 0 t sin ( t y ( τ ) τ ) d τ , t [ 0 , 3 ] ,
    with exact solution y ( t ) 1 ;
  • the stiff VIE
    y ( t ) = 0 t λ y ( τ ) sin ( τ ) + cos ( τ ) d τ , t [ 0 , 3 4 π ] ,
    with λ = 10 4 and exact solution y ( t ) = sin ( t ) . This is a stiff problem because it is equivalent to the Prothero-Robinson problem for ODEs.
We compare TSAC methods with superconvergent one step collocation methods of [2,3], where m denotes the number of collocation points and p denotes the order of the method:
  • G2: 1 point Gauss collocation, c 2 = 1 , m = 2 , p = 2 ;
  • R2: 2 points Radau collocation, m = 2 , p = 3 ;
  • TSAC2: 2 points TSAC method, m = 2 , p = 4 .
The method TSAC2 is the two-stage TSAC method described in Example 1. The accuracy is defined by the number of correct significant digits c d at the end point (the maximal absolute end point error is written as 10 c d ). For each test we plot in Figure 1 the number of c d versus the number of mesh points N. We observe as for non stiff Problem (37) the effective order of the all methods is coherent with the theoretical order, while for stiff Problem (38) the one step methods show order reduction as the effective order reduces to p = 2 .

5. One-Step Collocation Methods for VIDEs

We concentrate on VIDEs of type:
y ( t ) = g ( t , y ( t ) ) + 0 t k ( t , τ , y ( τ ) ) d τ , t I = [ 0 , T ] , y ( 0 ) = y 0 ,
where g ( t , y ) : I × I R d I R d , k ( t , s , y ) : S × I R d I R d , S = { ( t , s ) | 0 s t T } . For sake of completeness we report the theorem of existence and uniqueness of solution for (39) [3].
Theorem 6.
Let g ( t , y ) and k ( t , s , y ) be continuous functions and satisfy a uniform Lipschitz condition with respect to y. Then there exists a unique solution y C 1 ( [ 0 , T ] ) of the problem (39).
Let I h = { t n : 0 < t 0 < t 1 < < t N = T } be a partition of the time interval [ 0 , T ] with constant stepsize h = t n + 1 t n , n = 0 , , N 1 . The integro-differential Equation (39) can be written as follows:
y ( t ) = g ( t , y ( t ) ) + F n ( t , y ( · ) ) + Φ n ( t , y ( · ) ) , t [ t n , t n + 1 ] ,
where
F n ( t , y ( · ) ) = 0 t n k ( t , τ , y ( τ ) ) d τ , Φ n ( t , y ( · ) ) = t n t k ( t , τ , y ( τ ) ) d τ ,
represent respectively the lag term and the increment function.

5.1. Exact One-Step Collocation Methods

Here we briefly expose the classical one-step collocation methods for VIDEs and their main properties [2,3].
A one-step collocation method approximates y ( t ) by a piecewise polynomial P ( t ) , with P ( t ) = P n ( t ) , t [ t n , t n + 1 ] , n = 0 , , N 1 , where
P n ( t n + s h ) = y n + h j = 1 m β j ( s ) U n j , s [ 0 , 1 ] ,
with y n = P n ( t n ) , U n j = P n ( t n + c j h ) , β j ( s ) = 0 s L j ( τ ) d s , L j ( τ ) being the j-th Lagrange fundamental polynomial with respect to the collocation parameters.
The m unknowns U n j are found by imposing that P n ( t ) satisfies (39) at the collocation points t n j : = t n + c j h , j = 1 , , m , n = 0 , , N 1 , where 0 c 1 < < c m 1 are the collocation parameters.
The numerical approximation at the point t n + 1 is then given by y n + 1 = P n ( t n + 1 ) . The final form of an exact collocation method is
U n i = g ( t n i , P n ( t n i ) ) + F n ( t n i , P ( · ) ) + Φ n ( t n i , P ( · ) ) , i = 1 , , m , y n + 1 = y n + h i = 1 m β j ( 1 ) U n i ,
n = 0 , , N , where the lag term and the increment function can be written as
F n ( t n i , P ( · ) ) = h ν = 0 n 1 0 1 k t n i , t ν + τ h , P ν ( t n + τ h ) d τ ,
Φ n ( t n i , P ( · ) ) = h 0 c i k t n i , t n + τ h , P n ( t n + τ h ) d τ .
The first equation in (41) requires, at each time step, the solution of an m-dimensional nonlinear system in the unknowns { U n i } i = 1 m .
For every choice of the collocation parameters c 1 , , c m , the collocation polynomial P ( t ) is continuous on [ 0 , T ] and provides a uniform approximation of order O ( h m ) . Moreover, if c 1 , , c m are suitably chosen, the order of convergence at the mesh points increases (local superconvergence): is 2 m 2 for the Lobatto points, 2 m 1 for the Radau points and 2 m for the Gauss ones [2,3].

5.2. Discretized One-Step Collocation Methods

In the general case, the integrals appearing in (42) and (43) cannot be exactly evaluated, so a further approximation is needed in order to fully discretize the method. Let us suppose to approximate these integrals by quadrature formulae of the type:
F ˜ n ( t n i , P ( · ) ) = h ν = 0 n 1 l = 1 μ 1 w l k ( t n i , t ν + d l h , P ν ( t ν + d l h ) ) ,
Φ ˜ n ( t n i , P ( · ) ) = h l = 1 μ 0 w i l k ( t n i , t n + d i l h , P n ( t n + d i l h ) ) .
These formulae are then used to define the discretized collocation methods as
U ˜ n i = g ( t n i , P ˜ n ( t n i ) ) + F ˜ n ( t n i , P ˜ ( · ) ) + Φ ˜ n ( t n i , P ˜ ( · ) ) y ˜ n + 1 = y ˜ n + h i = 1 m β j ( 1 ) U ˜ n i ,
where the collocation polynomial is now of the form
P ˜ n ( t n + s h ) = y n + h j = 1 m β j ( s ) U ˜ n j , s [ 0 , 1 ] .
The discretized collocation methods are a special class of the Runge-Kutta extended methods and preserve the order of convergence and superconvergence of exact collocation methods, if the quadrature Formulae (44) and (45) are sufficiently accurate [3].
Some relevant stability results for one-step collocation methods are derived in [36,37].

6. Multistep Collocation for VIDEs

6.1. Exact Multistep Collocation

Recently, in order to obtain an higher order of convergence at the same computational effort, multistep collocation methods have been introduced: the solution y ( t ) is approximated by a piecewise algebraic polynomial P ( t ) :
P ( t n + s h ) = k = 0 r 1 φ k ( s ) y n k + h j = 1 m ψ j ( s ) U n j , s [ 0 , 1 ] ,
where again
U n j : = P ( t n j ) , j = 1 , , m ,
and the functions φ k ( s ) , ψ j ( s ) are polynomials of degree m + r 1 which are determined by imposing that the polynomial (48) satisfies (49) and the interpolation conditions:
P ( t n k ) = y n k , k = 0 , , r 1 .
For any fixed set of collocation parameters c 1 , , c m , conditions (49) and (50) lead to the following non linear system of ( r + m ) 2 equations, where the ( r + m ) 2 unknowns are the coefficients of the polynomials φ k ( s ) and ψ j ( s ) :
φ l ( k ) = δ l k , φ l ( c j ) = 0 , ψ i ( c j ) = δ i j , ψ i ( k ) = 0 ,
l , k = 0 , , r 1 , i , j = 1 , , m .
Exact multistep collocation methods are obtained by imposing that the collocation polynomial (48) satisfies the VIDE at the collocation points t n i , and by computing y n + 1 = P n ( t n + 1 ) :
U n i = g ( t n i , P ( t n i ) ) + F n ( t n i , P ( · ) ) + Φ n ( t n i , P ( · ) ) , i = 1 , , m y n + 1 = k = 0 r 1 φ k ( 1 ) y n k + h i = 1 m ψ j ( 1 ) U n i .
n = r 1 , , N , where now the lag term and the increment function can be written as
F n ( t n i , P ( · ) ) = h ν = 0 n 1 0 1 k t n i , t ν + τ h , k = 0 r 1 φ k ( τ ) y ν k + h j = 1 m ψ j ( τ ) U ν j d τ ,
Φ n ( t n i , P ( · ) ) = h 0 c i k t n i , t n + τ h , k = 0 r 1 φ k ( τ ) y n k + h j = 1 m ψ j ( τ ) U n j d τ .
We note that at each time step, the approximations y n k , k = 0 , , r 1 are already known, so only the unknowns { U n i } i = 1 m need to be computed, by solving the nonlinear system given by the first equation of (52).
Observe that we are able to give an approximate value P ( t ) of the solution y ( t ) at each point t of the integration interval, therefore we have a uniform approximation of the solution on [ 0 , T ] .
The classical one-step collocation methods described in the previous section can be seen as a particular case of multistep methods with r = 1 and
φ 0 ( s ) 1 , ψ j ( s ) = 0 s L j ( τ ) d τ ,
where L j ( τ ) is the j-th Lagrange fundamental polynomial with respect to the collocation parameters. We observe that, at each time step, both one-step and multistep collocation methods require the solution of a non linear system of dimension m for the stages U n i , i = 1 , , m . The multistep methods only need in addition the computation of the starting values y 1 , , y r 1 .

6.2. Discretized Multistep Collocation

As in the case of one-step collocation methods, it is evident that the exact multistep collocation methods (52) are not directly applicable for the implementation, since approximations of the integrals F n ( t n i , P ( · ) ) and Φ n ( t n i , P ( · ) ) are needed. With the aim of fully discretizing the multistep collocation methods we consider the following quadrature formulas
0 c i α ( x ) d x Q i ( α ( · ) ) : = l = 1 μ 0 w i l α ( d i l ) , 0 1 α ( x ) d x Q ( α ( · ) ) : = l = 1 μ 1 w l α ( d l ) ,
where the weights and nodes are suitably chosen, as it will be illustrated later.
The discretized multistep collocation method for the problem (39) approximates the solution y ( t ) with a piecewise polynomial P ˜ ( t ) , with
P ˜ ( t n + s h ) = k = 0 r 1 φ k ( s ) y ˜ n k + h j = 1 m ψ j ( s ) U ˜ n j , s [ 0 , 1 ] ,
where the polynomials { φ k ( s ) } k = 0 r 1 , { ψ j ( s ) } j = 1 m are the same as in the exact collocation, and can be computed by solving the system (51).
We impose that at the collocation points P ˜ ( t ) satisfies the VIDE (39), where the integrals appearing in both the lag term (53) and the increment function (54) are approximated by the quadrature formulae defined in (55), and we set y ˜ n + 1 = P ˜ ( t n + h ) . Thus the discretized multistep method is
U ˜ n i = g ( t n i , P ˜ ( t n i ) ) + F ˜ n ( t n i , P ˜ ( · ) ) + Φ ˜ n ( t n i , P ˜ ( · ) ) , i = 1 , , m y ˜ n + 1 = k = 0 r 1 φ k ( 1 ) y ˜ n k + h i = 1 m ψ i ( 1 ) U ˜ n i .
where F ˜ n ( t n i , P ˜ ( · ) ) and Φ ˜ n ( t n i , P ˜ ( · ) ) are of the form
F ˜ n ( t n i , P ˜ ( · ) ) = h ν = 0 n 1 l = 1 μ 1 w l k t n i , t ν + d l h , k = 0 r 1 φ k ( d l ) y ˜ ν k + h j = 1 m ψ j ( d l ) U ˜ ν j Φ ˜ n ( t n i , P ˜ ( · ) ) = h l = 1 μ 0 w i l k t n i , t n + d i l h , k = 0 r 1 φ k ( d i l ) y ˜ n k + h j = 1 m ψ j ( d i l ) U ˜ n j .
i.e., they are obtained by applying the quadrature Formula (55) to the integrals appearing in (53) and (54).

6.3. Convergence Analysis

The multivalue nature of the multistep methods imposes to analyze first the zero-stability of the methods. When h 0 , second equation of (52) reduces to
y n + 1 = k = 0 r 1 φ k ( 1 ) y n k .
Therefore, the method (48) and (52) is said to be zero-stable, if all of the roots of the polynomial
p ( λ ) = λ r k = 0 r 1 φ k ( 1 ) λ r k 1
have modulus less than or equal to unity, and those of modulus unity are simple.
On this basis, the following theorem studies the convergence of the method.
Theorem 7.
Consider the problem (39) with d = 1 . Let p = m + r 1 and assume that:
  • k C p ( S × I R ) and g C p ( [ 0 , T ] × I R ) and have bounded derivatives with respect to y;
  • the method (48) and (52) is zero-stable;
  • the starting error satisfies | e ( t ) | = O ( h p ) , for any t [ t 0 , t r 1 ] .
Then, the global error e ( t ) = y ( t ) P ( t ) of the exact MCM (48) and (52) satisfies
max [ 0 , T ] | e ( t ) | C h m + r 1 .
By a suitable choice of the collocation parameters, it is possible to increase the order of convergence at the mesh points (local superconvergence), following the lines of multistep methods for ODEs (compare [21], Section 3).
Theorem 8.
Assume that hypotheses of Theorem 7 hold with p = 2 m + r 1 and that the collocation parameters satisfy these conditions
k = 1 r 1 1 c i + k + 2 j = 1 j i m 1 c i c j = 0 , i = 1 , , m .
Then the order of the exact MCM (48) and (52) at the mesh points is p, i.e.,:
max 1 n N | e ( t n ) | = O ( h 2 m + r 1 ) .
Similar convergence and superconvergence results hold also for the discretized MCM (56) and (57). We can summarize them in the following theorem.
Theorem 9.
Assume that hypotheses of Theorem 7 hold. If quadrature formulae Q and Q i defined in (55) have order m + r and m + r 1 respectively, then the uniform order of the discretized method (56) and (57) is equal to m + r 1 .
Moreover, if hypotheses of Theorem 8 are fulfilled, Q and Q i defined in (55) have order 2 m + r and 2 m + r 1 respectively, then the order of the discretized method (56) and (57), at the mesh points, is 2 m + r 1 .
We observe that, at the same cost of one-step collocation methods with m collocation parameters, multistep collocation methods have an higher computational cost. A further improvement of the efficiency could be obtained by exploiting parallel techniques, as done for example in [38,39,40].
An extensive analysis of the stability properties on basic test equations is contained in [14]. A possible future development may regard new multistep methods with some relaxing order conditions, which leave some parameters free to perform a numerical search for the methods with optimal stability properties, as done in [22,23,41,42,43] in the context of ODEs.

6.4. Numerical Results

Now we give a short numerical illustration of discretized MCMs (56) and (57), on the linear test equation
y ( t ) = g ( t , y ) 0 t t 2 exp ( s t ) y ( s ) d s , t [ 0 , 1 ] , y ( 0 ) = 1 ,
with g ( t , y ) such that y ( t ) = exp ( t ) ; and on the nonlinear problem
y ( t ) = g ( t , y ) 0 t 2 t sin ( s ) exp ( y ( s ) ) d s , t [ 0 , 1 ] , y ( 0 ) = 1 ,
with g ( t , y ) such that y ( t ) = cos ( t ) .
We consider three methods
  • TS3: superconvergent discretized two-step collocation method, with r = 2 and m = 1 , with order p = 3 ;
  • TS3b: two-step discretized collocation method, with r = 2 and m = 2 , c 1 = 0.9 , c 2 = 1 , with uniform order 3;
  • TS5: superconvergent discretized two-step collocation method, with r = 2 and m = 2 , with order p = 5 .
Method TS3b has an unbounded stability region, while TS3 and TS5 have a bounded stability region. The exact expression of the methods and their stability region can be found in [14]. To confirm the theoretical order of convergence, in Figure 2 the error (in logarithmic scale) produced by methods TS3, TS3b and TS5 when applied to problems (61) and (62), and the slopes corresponding to order 3 and 5. We see that the effective order is equal to the theoretical one.

7. Conclusions

We have illustrated multistep collocation methods for VIEs and VIDEs and gave an overview of their convergence and superconvergence properties. This idea may be exploited to obtain high order methods for solving other types of equations as well. For example, recently two-step collocation methods have been proposed for fractional differential equations [44], and further developments may be achieved for other fractional models, as time fractional differential equations [45]. Further issues of this research will focus on oscillatory problems [46,47] and in particular on the application of multistep collocation methods to periodic integral equations [48,49]. Moreover, it seems reasonable to consider the possibility of employing collocation spaces based on functions other than polynomials, as in [50,51,52] and similarly as in the case of oscillatory problems [53], and merge into the numerical scheme as many known qualitative properties of the continuous problem as possible, in a structure-preserving perspective [54].
The literature on the numerical treatment of VIEs is quite rich and goes beyond the results considered in this review. Here we would like to mention some other results, at least. In [55] the modified Newton–Kantorovich method combined with collocation were applied non linear and nonlinear VIE with piecewise smooth kernels. Such VIE were introduced in [56] and asymptotic approximations to parametric families of solutions were constructed and the existence of continuous solutions was proved. The review of the numerical methods of optimal accuracy (spline-collocation technique) for multidimensional weakly singular VIEs is given in [57]. Some other interesting papers regard the distance between the approximate and exact solutions of various generalizations of the Volterra equations [58,59,60,61,62,63]. Lastly, we underline that in the practical applications of VIE based models it is extremely important to have the numerical method to be stable with respect to the measurement errors both in the source function and in the kernel. It is well known that the 1st kind of VIEs enjoy self-regularization property when the mesh step serves as the regularisation parameter. In addition, the Lavrentiev type regularisation is a good option [64,65]. These issues have been discussed in [66].

Funding

This research was supported by GNCS-INdAM.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Bonaccorsi, S.; Fantozzi, M. Volterra Integro-Differential Equations with Accretive Operators and Non-Autonomous Perturbations. J. Integral Equ. Appl. 2006, 18, 437–470. [Google Scholar] [CrossRef] [Green Version]
  2. Brunner, H. Collocation Methods for Volterra Integral and Related Functional Equations; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  3. Brunner, H.; van der Houwen, P.J. The Numerical Solution of Volterra Equations; CWI Monographs, 3; Elsevier Science Ltd.: Amsterdam, The Netherlands, 1986. [Google Scholar]
  4. Hoppensteadt, F.C.; Jackiewicz, Z.; Zubik-Kowal, B. Numerical solution of Volterra integral and integro-differential equations with rapidly vanishing convolution kernels. BIT Numer. Math. 2007, 47, 325–350. [Google Scholar] [CrossRef]
  5. Hrusa, W.J.; Renardy, M. A model problem in one-dimensional viscoelasticity with a singular kernel. In Volterra Integrodifferential Equations in Banach Spaces and Applications; Pitman Research Notes in Mathematics; da Prato, G., Iannelli, M., Eds.; Longman Scientific & Technical: Harlow, UK; Wiley: New York, NY, USA, 1989; Volume 190, pp. 221–230. [Google Scholar]
  6. Nohel, J.A.; Rogers, R.C.; Tzavaras, A.E. Hyperbolic conservation laws in viscoelasticity. In Volterra Integro Differential Equations in Banach Spaces and Applications; Da Prato, G., Iannelli, M., Eds.; Pitman Research Notes in Mathematics Series, 190; Longman Science Technology: Harlow, UK, 1989; pp. 320–338. [Google Scholar]
  7. Coleman, J.P.; Duxbury, S.C. Mixed collocation methods for y = f(x,y). J. Comput. Appl. Math. 2000, 126, 47–75. [Google Scholar] [CrossRef]
  8. Brunner, H.; Makroglou, A.; Miller, R.K. Mixed interpolation collocation methods for first and second order Volterra integro-differential equations with periodic solution. Appl. Numer. Math. 1997, 23, 381–402. [Google Scholar] [CrossRef]
  9. Braś, M.; Cardone, A.; Jackiewicz, Z.; Welfert, B. Order reduction phenomenon for general linear methods. J. Comput. Appl. Math. 2015, 290, 44–64. [Google Scholar] [CrossRef]
  10. Butcher, J.C. Numerical Methods for Ordinary Differential Equations, 2nd ed.; John Wiley & Sons: Chichester, UK, 2008. [Google Scholar]
  11. Hairer, E.; Wanner, G. Solving Ordinary Differential Equations. II. In Springer Series in Computational Mathematics; Springer: Berlin, Germany, 1991; Volume 14. [Google Scholar]
  12. D’Ambrosio, R.; Paternoster, B. Two-step modified collocation methods with structured coefficient matrices for ordinary differential equations. Appl. Numer. Math. 2012, 62, 1325–1334. [Google Scholar] [CrossRef]
  13. Capobianco, G.; Conte, D.; Paternoster, B. Construction and implementation of two-step continuous methods for Volterra Integral Equations. Appl. Numer. Math 2017, 119, 239–247. [Google Scholar] [CrossRef]
  14. Cardone, A.; Conte, D. Multistep collocation methods for Volterra integro-differential equations. Appl. Math. Comput. 2013, 221, 770–785. [Google Scholar] [CrossRef]
  15. Cardone, A.; Conte, D.; Paternoster, B. A family of multistep collocation methods for Volterra integro-differential equations. AIP Conf. Proc. 2009, 1168, 358–361. [Google Scholar] [CrossRef]
  16. Conte, D.; Paternoster, B. Multistep collocation methods for Volterra Integral Equations. Appl. Numer. Math. 2009, 59, 1721–1736. [Google Scholar] [CrossRef]
  17. Conte, D.; D’Ambrosio, R.; Paternoster, B. Two-step diagonally-implicit collocation-based methods for Volterra Integral Equations. Appl. Numer. Math. 2012, 62, 1312–1324. [Google Scholar] [CrossRef]
  18. Conte, D.; Jackiewicz, Z.; Paternoster, B. Two-step almost collocation methods for Volterra integral equations. Appl. Math. Comput. 2008, 204, 839–853. [Google Scholar] [CrossRef]
  19. Guillou, A.; Soulé, J.L. La résolution numérique des problèmes différentiels aux conditions initiales par des méthodes de collocation. Rev. Fr. Inform. Rech. Opér. 1969, 3, 17–44. [Google Scholar] [CrossRef] [Green Version]
  20. Lie, I. The stability function for multistep collocation methods. Numer. Math. 1990, 57, 779–787. [Google Scholar] [CrossRef]
  21. Lie, I.; Nørsett, S. Superconvergence for multistep collocation. Math. Comput. 1989, 52, 65–79. [Google Scholar] [CrossRef]
  22. Braś, M.; Cardone, A. Construction of Efficient General Linear Methods for Non-Stiff Differential Systems. Math. Model. Anal. 2012, 17, 171–189. [Google Scholar] [CrossRef]
  23. Braś, M.; Cardone, A.; D’Ambrosio, R. Implementation of explicit Nordsieck methods with inherent quadratic stability. Math. Model. Anal. 2013, 18, 289–307. [Google Scholar] [CrossRef]
  24. Jackiewicz, Z. General Linear Methods for Ordinary Differential Equations; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  25. Darania, P. Superconvergence analysis of multistep collocation method for delay Volterra integral equations. Comput. Methods Differ. Equ. 2016, 4, 205–216. [Google Scholar]
  26. Darania, P.; Pishbin, S. High-order collocation methods for nonlinear delay integral equation. J. Comput. Appl. Math. 2017, 326, 284–295. [Google Scholar] [CrossRef]
  27. Fazeli, S.; Hojjati, G. Numerical solution of Volterra integro-differential equations by superimplicit multistep collocation methods. Numer. Algorithms 2015, 68, 741–768. [Google Scholar] [CrossRef]
  28. Fazeli, S.; Hojjati, G.; Shahmorad, S. Multistep Hermite collocation methods for solving Volterra integral equations. Numer. Algorithms 2012, 60, 27–50. [Google Scholar] [CrossRef]
  29. Fazeli, S.; Hojjati, G.; Shahmorad, S. Super implicit multistep collocation methods for nonlinear Volterra integral equations. Math. Comput. Model. 2012, 55, 590–607. [Google Scholar] [CrossRef]
  30. Fazeli, S.; Hojjati, G.; Shahmorad, S. Multistep collocation and iterated multistep collocation methods for solving two-dimensional Volterra integral equations. J. Mod. Methods Numer. Math. 2015, 6, 1–21. [Google Scholar] [CrossRef]
  31. Ma, J.; Xiang, S. A collocation boundary value method for linear Volterra integral equations. J. Sci. Comput. 2017, 71, 1–20. [Google Scholar] [CrossRef]
  32. Sheng, C.; Wang, Z.; Guo, B. A multistep Legendre-Gauss spectral collocation method for nonlinear Volterra integral equations. SIAM J. Numer. Anal. 2014, 52, 1953–1980. [Google Scholar] [CrossRef]
  33. Lopez-Fernandez, M.; Lubich, C.; Schadle, A. Adaptive, fast, and oblivious convolution in evolution equations with memory. SIAM J. Sci. Comput. 2008, 30, 1015–1037. [Google Scholar] [CrossRef]
  34. Lopez-Fernandez, M.; Lubich, C.; Schadle, A. Fast and oblivious convolution quadrature. SIAM J. Sci. Comput. 2006, 28, 421–438. [Google Scholar]
  35. Crisci, M.R.; Russo, E.; Vecchio, A. Stability results for one-step discretized collocation methods in the numerical treatment of Volterra integral equations. Math. Comput. 1992, 58, 119–134. [Google Scholar] [CrossRef]
  36. Crisci, M.R.; Russo, E.; Jackiewicz, Z.; Vecchio, A. Global stability of exact collocation methods for Volterra integro-differential equations. Atti Sem. Mat. Fis. Univ. Modena 1991, 39, 527–536. [Google Scholar]
  37. Crisci, M.R.; Russo, E.; Vecchio, A. Stability of Collocation Methods for Volterra Integro-Differential Equations. J. Integral Equ. Appl. 1992, 4, 491–507. [Google Scholar] [CrossRef]
  38. Cardone, A.; Messina, E.; Vecchio, A. An adaptive method for Volterra–Fredholm integral equations on the half line. J. Comput. Appl. Math. 2009, 228, 538–547. [Google Scholar] [CrossRef]
  39. Conte, D.; D’Ambrosio, R.; Paternoster, B. GPU acceleration of waveform relaxation methods for large differential systems. Numer. Algorithms 2016, 71, 293–310. [Google Scholar] [CrossRef]
  40. Conte, D.; Paternoster, B. Parallel methods for weakly singular Volterra Integral Equations on GPUs. Appl. Numer. Math. 2017, 114, 30–37. [Google Scholar] [CrossRef]
  41. Cardone, A.; Jackiewicz, Z.; Sandu, A.; Zhang, H. Extrapolated Implicit-Explicit Runge-Kutta Methods. Math. Model. Anal. 2014, 19, 18–43. [Google Scholar] [CrossRef]
  42. Conte, D.; D’Ambrosio, R.; Jackiewicz, Z.; Paternoster, B. A practical approach for the derivation of two-step Runge-Kutta methods. Math. Model. Anal. 2012, 17, 65–77. [Google Scholar] [CrossRef]
  43. Conte, D.; D’Ambrosio, R.; Jackiewicz, Z.; Paternoster, B. Numerical search for algebraically stable two-step continuous Runge-Kutta methods. J. Comput. Appl. Math. 2013, 239, 304–321. [Google Scholar] [CrossRef]
  44. Cardone, A.; Conte, D.; Paternoster, B. Two-step collocation methods for fractional differential equations. Discret. Contin. Dyn. Syst. Ser. B 2017, 22, 1–17. [Google Scholar] [CrossRef]
  45. Burrage, K.; Cardone, A.; D’Ambrosio, R.; Paternoster, B. Numerical solution of time fractional diffusion systems. Appl. Numer. Math. 2017, 116, 82–94. [Google Scholar] [CrossRef]
  46. Conte, D.; Paternoster, B. Modified Gauss-Laguerre exponential fitting based formulae. J. Sci. Comput. 2016, 69, 227–243. [Google Scholar] [CrossRef]
  47. Ixaru, L.G.; Paternoster, B. A Gauss quadrature rule for oscillatory integrands. Comput. Phys. Commun. 2001, 133, 177–188. [Google Scholar] [CrossRef]
  48. Cardone, A.; D’Ambrosio, R.; Paternoster, B. High order exponentially fitted methods for Volterra integral equations with periodic solution. Appl. Numer. Math. 2017, 114, 18–29. [Google Scholar] [CrossRef]
  49. Cardone, A.; Ixaru, L.G.; Paternoster, B.; Santomauro, G. Ef-Gaussian direct quadrature methods for Volterra integral equations with periodic solution. Math. Comput. Simul. 2015, 110, 125–143. [Google Scholar] [CrossRef]
  50. Cardone, A.; D’Ambrosio, R.; Paternoster, B. Exponentially fitted IMEX methods for advection—Diffusion problems. J. Comput. Appl. Math. 2017, 316, 100–108. [Google Scholar] [CrossRef]
  51. D’Ambrosio, R.; Moccaldi, M.; Paternoster, B. Adapted numerical methods for advection-reaction-diffusion problems generating periodic wavefronts. Comput. Math. Appl. 2017, 74, 1029–1042. [Google Scholar] [CrossRef]
  52. D’Ambrosio, R.; Paternoster, B. Numerical solution of reaction–diffusion systems of λω type by trigonometrically fitted methods. J. Comput. Appl. Math. 2016, 294, 436–445. [Google Scholar] [CrossRef]
  53. D’Ambrosio, R.; De Martino, G.; Paternoster, B. General Nystrom methods in Nordsieck form: Error analysis. J. Comput. Appl. Math. 2016, 292, 694–702. [Google Scholar] [CrossRef]
  54. Butcher, J.; D’Ambrosio, R. Partitioned general linear methods for separable Hamiltonian problems. Appl. Numer. Math. 2017, 117, 69–86. [Google Scholar] [CrossRef]
  55. Muftahov, I.; Tynda, A.; Sidorov, D. Numeric solution of Volterra integral equations of the first kind with discontinuous kernels. J. Comput. Appl. Math. 2017, 313, 119–128. [Google Scholar] [CrossRef]
  56. Sidorov, D.N. On parametric families of solutions of Volterra integral equations of the first kind with piecewise smooth kernel. Differ. Equ. 2013, 49, 210–216. [Google Scholar] [CrossRef]
  57. Boykov, I.V.; Tynda, A.N. Numerical methods of optimal accuracy for weakly singular Volterra integral equations. Ann. Funct. Anal. 2015, 6, 114–133. [Google Scholar] [CrossRef] [Green Version]
  58. Castro, L.P.; Ramos, A. Hyers–Ulam–Rassias stability for a class of nonlinear Volterra integral equations. Banach J. Math. Anal. 2009, 3, 36–43. [Google Scholar] [CrossRef]
  59. Brzdek, J.; Eghbali, N. On approximate solutions of some delayed fractional differential equations. Appl. Math. Lett. 2016, 54, 31–35. [Google Scholar] [CrossRef]
  60. Bahyrycz, A.; Brzdek, J.; Lesniak, Z. On approximate solutions of the generalized Volterra integral equation. Nonlinear Anal. Real World Appl. 2014, 20, 59–66. [Google Scholar] [CrossRef]
  61. Gachpazan, M.; Baghani, O. Hyers-Ulam stability of nonlinear integral equation. Fixed Point Theory Appl. 2010, 2010, 927640. [Google Scholar] [CrossRef]
  62. Jung, S.M. A fixed point approach to the stability of a Volterra integral equation. Fixed Point Theory Appl. 2007, 2007, 057064. [Google Scholar] [CrossRef]
  63. Morales, J.R.; Rojas, E.M. Hyers–Ulam and Hyers–Ulam–Rassias stability of nonlinear integral equations with delay. Int. J. Nonlinear Anal. Appl. 2011, 2, 1–7. [Google Scholar]
  64. Kythy, P.K.; Puri, P. Computational Methods for Linear Integral Equations; Birkhauser: Boston, MA, USA, 2002. [Google Scholar]
  65. Muftahov, I.R.; Sidorov, D.N.; Sidorov, N.A. Lavrentiev regularization of integral equations of the first kind in the space of continuous functions. Izvestiya Irkutskogo Gosudarstvennogo Universiteta 2016, 15, 62–77. [Google Scholar]
  66. Muftahov, I.; Sidorov, D.; Zhukov, A.; Panasetsky, D.; Foley, A.; Li, Y.; Tynda, A. Application of Volterra Equations to Solve Unit Commitment Problem of Optimised Energy Storage and Generation. arXiv, 2016; arXiv:1608.05221. [Google Scholar]
Figure 1. Number of correct significant digits with respect to the number of mesh points. (a) Problem (37); (b) Problem (38).
Figure 1. Number of correct significant digits with respect to the number of mesh points. (a) Problem (37); (b) Problem (38).
Axioms 07 00045 g001
Figure 2. Error of two-step methods TS3 ( Axioms 07 00045 i001), TS3b ( Axioms 07 00045 i002) and TS5 ( Axioms 07 00045 i003), and slopes of order 3 (dashed line) and of order 5 (dash-dot line), applied to problem (61) (left) and on problem (62) (right).
Figure 2. Error of two-step methods TS3 ( Axioms 07 00045 i001), TS3b ( Axioms 07 00045 i002) and TS5 ( Axioms 07 00045 i003), and slopes of order 3 (dashed line) and of order 5 (dash-dot line), applied to problem (61) (left) and on problem (62) (right).
Axioms 07 00045 g002

Share and Cite

MDPI and ACS Style

Cardone, A.; Conte, D.; D’Ambrosio, R.; Paternoster, B. Collocation Methods for Volterra Integral and Integro-Differential Equations: A Review. Axioms 2018, 7, 45. https://doi.org/10.3390/axioms7030045

AMA Style

Cardone A, Conte D, D’Ambrosio R, Paternoster B. Collocation Methods for Volterra Integral and Integro-Differential Equations: A Review. Axioms. 2018; 7(3):45. https://doi.org/10.3390/axioms7030045

Chicago/Turabian Style

Cardone, Angelamaria, Dajana Conte, Raffaele D’Ambrosio, and Beatrice Paternoster. 2018. "Collocation Methods for Volterra Integral and Integro-Differential Equations: A Review" Axioms 7, no. 3: 45. https://doi.org/10.3390/axioms7030045

APA Style

Cardone, A., Conte, D., D’Ambrosio, R., & Paternoster, B. (2018). Collocation Methods for Volterra Integral and Integro-Differential Equations: A Review. Axioms, 7(3), 45. https://doi.org/10.3390/axioms7030045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop