Next Article in Journal
Ant-Balanced Multiple Traveling Salesmen: ACO-BmTSP
Next Article in Special Issue
Improved Gradient Descent Iterations for Solving Systems of Nonlinear Equations
Previous Article in Journal
Special Issue on Algorithms for PID Controllers 2021
Previous Article in Special Issue
On the Numerical Treatment of the Temporal Discontinuity Arising from a Time-Varying Point Mass Attachment on a Waveguide
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Symbolic Method for Solving a Class of Convolution-Type Volterra–Fredholm–Hammerstein Integro-Differential Equations under Nonlocal Boundary Conditions

by
Efthimios Providas
* and
Ioannis Nestorios Parasidis
Department of Environmental Sciences, Gaiopolis Campus, University of Thessaly, 415 00 Larissa, Greece
*
Author to whom correspondence should be addressed.
Algorithms 2023, 16(1), 36; https://doi.org/10.3390/a16010036
Submission received: 16 December 2022 / Revised: 4 January 2023 / Accepted: 4 January 2023 / Published: 7 January 2023
(This article belongs to the Special Issue Computational Methods and Optimization for Numerical Analysis)

Abstract

:
Integro-differential equations involving Volterra and Fredholm operators (VFIDEs) are used to model many phenomena in science and engineering. Nonlocal boundary conditions are more effective, and in some cases necessary, because they are more accurate measurements of the true state than classical (local) initial and boundary conditions. Closed-form solutions are always desirable, not only because they are more efficient, but also because they can be valuable benchmarks for validating approximate and numerical procedures. This paper presents a direct operator method for solving, in closed form, a class of Volterra–Fredholm–Hammerstein-type integro-differential equations under nonlocal boundary conditions when the inverse operator of the associated Volterra integro-differential operator exists and can be found explicitly. A technique for constructing inverse operators of convolution-type Volterra integro-differential operators (VIDEs) under multipoint and integral conditions is provided. The proposed methods are suitable for integration into any computer algebra system. Several linear and nonlinear examples are solved to demonstrate the effectiveness of the method.

1. Introduction

Volterra–Fredholm integro-differential equations (VFIDEs) are encountered in many different branches of science and engineering [1,2,3,4]. As they are extremely difficult and in some cases even impossible to solve in exact closed form, many different techniques have been developed to construct approximate and numerical solutions. Some of these methods are mentioned below.
Two well-known, traditional techniques are the series solution method [1,5] and the variational iteration method [1]. A powerful and widely used method for solving linear and nonlinear integro-differential equations including VFIDEs is the Adomian decomposition method (ADM) and its variant, known as modified ADM [6,7]. Methods employing Taylor polynomial expansions were developed and applied to find the approximate solution of high-order linear VFIDEs under mixed conditions [8] and high-order nonlinear VFIDEs [9,10]. Collocation methods based on Haar wavelets and Legendre wavelets for the numerical solution of nonlinear VFIDEs were presented in [3,11], respectively. An approach using hybrid Legendre polynomials and block–pulse functions was also reported in [12]. The authors of [13] applied the polynomial least squares method to construct approximate analytical solutions for a very general class of nonlinear VFIDEs. A moving least squares (MLS) method was developed in [14]. The Tau method was also used for the numerical solution of the general form of linear VFIDEs [15]. In [16], the reproducing kernel Hilbert space method was employed to numerically solve first-order periodic VFIDEs. A method based on triangular functions and the relative operational matrix to determine the numerical solution of specific nonlinear VFIDEs was proposed in [17]. Fixed-point techniques for the approximate solution of linear and first-order nonlinear VFIDEs were discussed in [18,19], respectively. A quadrature and iterative technique for solving a nonlinear VFIDE with weakly singular kernels was developed in [20]. Numerical techniques for solving fractional VFIDEs were proposed in [21,22,23] using second Chebyshev wavelets (SCWs), fractional-order Bernoulli functions (FBFs) and Lucas wavelets (LWs), respectively. The numerical solution of multipoint boundary value problems for linear VFIDEs was studied using Legendre polynomials [24,25], and Bernstein polynomials [26,27]. A semi-analytical numerical method for solving multipoint boundary value problems for linear VFIDEs of the neutral type with linear functional arguments was presented in [28]. Finally, for other related types of problems, such as ordinary differential equations under multipoint and integral conditions and nonlinear Volterra equations with loads, see [29,30].
This paper deals with the exact closed-form solution of a class of Volterra–Fredholm integro-differential equations of the Hammerstein type. It is emphasized that due to the complexity involved in solving integro-differential equations (IDEs) and in contrast to numerical methods, exact analytical techniques cannot have a general purpose but are designed and applied to a specific class of IDEs. Methods for solving explicitly different types of Volterra and Fredholm IDEs are discussed in detail in [31]. A method for the exact solution of systems of first-order linear Fredholm IDEs under multipoint and integral conditions was presented in [32]. The exact solution of linear and nonlinear convolution-type Volterra–Fredholm integral equations (VFIEs) was efficiently treated in [33]. In [4], the existence and uniqueness criteria and the closed-form solution of the convolution-type linear VFIDEs of the following form were investigated:
i = 0 n a i u ( n i ) ( x ) + i = 0 n 0 x k i ( x t ) u ( n i ) ( t ) d t j = 1 m 0 b k ¯ j ( x , t ) i = 0 n a i u ( n i ) ( t ) d t = f ( x ) ,
subject to the boundary conditions
Φ i ( u ) = 0 , i = 1 , 2 , , n ,
where a i , i = 0 , 1 , , n ( a 0 0 ) and b > 0 are given real constants; k i ( x ) , i = 0 , 1 , , n , k ¯ j ( x , t ) , j = 1 , 2 , , m , and f ( x ) are given continuous functions; u ( x ) is the unknown function assumed n times continuously differentiable; and u ( i ) ( x ) = d i u d x i and Φ i , i = 1 , 2 , , n , are given linear bounded functionals. This technique was used to solve the second-order VFIDEs encountered in the bending analysis of elastic Euler–Bernoulli beams in Eringen’s two-phase nonlocal integral model in the case of a distributed transverse load and boundary conditions for a simply supported beam, cantilever beam, clamped pinned beam and clamped beam [4,34].
In this article, we extend the technique to Volterra–Fredholm integro-differential equation in the Hammerstein form
i = 0 n a i u ( n i ) ( x ) i = 0 n a x k i ( x , t ) u ( n i ) ( t ) d t j = 1 m a b q j ( x , t ) φ j t , u ( t ) , u ( t ) , , u ( n ) ( t ) d t = f ( x ) ,
under the general nonlocal boundary conditions
Φ i ( u ) j = 1 l ν i j Θ j ( u ) = 0 , i = 1 , 2 , , n ,
where a i , i = 0 , 1 , , n , ( a 0 0 ), a and b are given real constants; k i ( x , t ) , i = 0 , 1 , , n , q j ( x , t ) , φ j ( · ) , j = 1 , 2 , , m , and f ( x ) are given continuous functions; u ( x ) is the unknown function assumed n times continuously differentiable with u ( i ) ( x ) = d i u d x i ; Φ i , i = 1 , 2 , , n , are some special, linear bounded functionals; Θ j , j = 1 , 2 , , l , are given general linear bounded functionals involving values at fixed points and definite integrals of u ( i 1 ) ( x ) , i = 1 , 2 , , n , in the interval [ a , b ] ; and ν i j , i = 1 , 2 , , n , j = 1 , 2 , , l , are real constants.
The outline of the paper is as follows: In Section 2, a direct operator method for solving nonlocal boundary value problems for nonlinear VFIDEs is presented. In Section 3, existence and uniqueness criteria and a solution formula for linear VFIDEs are given. Section 4 deals with the inversion of linear Volterra integro-differential operators (VIDEs) of convolution type and the solution of the corresponding VFIDEs. The proposed methods are tested by solving several examples in Section 5. Finally, conclusions and future directions are discussed in Section 6.

2. Formulation and Solution of Nonlinear VFIDEs

First, we formulate the nonlocal boundary problem (3), (4) in a convenient operator form.
Let X = C [ a , b ] , a , b R , and let the operators A , K : X X defined as
A u = i = 0 n a i u ( n i ) ( x ) ,
K u = i = 0 n a x k i ( x , t ) u ( n i ) ( t ) d t ,
be an nth-order linear differential operator and a linear Volterra-type integral operator, respectively, where n N , a 0 0 , a i , i = 1 , 2 , , n , are real constants; u ( x ) C n [ a , b ] ; and kernel functions k i ( x , t ) X × X , i = 0 , 1 , , n .
Furthermore, let the nonlinear Fredholm.type integral operator K ¯ : X X be defined as
K ¯ u = j = 1 m a b q j ( x , t ) φ j t , u ( t ) , u ( t ) , , u ( n ) ( t ) d t ,
where q j ( x , t ) , φ j ( · ) , j = 1 , 2 , , m , are continuous functions. We assume that kernels q j ( x , t ) are degenerate, i.e., without loss of generality of the form
q j ( x , t ) = g j ( x ) h j ( t ) , j = 1 , 2 , , m ,
where g j ( x ) , h j ( t ) X . In this case, Equation (7) becomes
K ¯ u = j = 1 m g j ( x ) a b h j ( t ) φ j t , u ( t ) , u ( t ) , , u ( n ) ( t ) d t .
By defining the row vector of functions
g = g 1 g 2 g m , g j = g j ( x ) , j = 1 , 2 , , m ,
and the column vector of functionals
Ψ ( u ) = Ψ 1 ( u ) Ψ 2 ( u ) Ψ m ( u ) , Ψ j ( u ) = a b h j ( t ) φ j t , u ( t ) , u ( t ) , , u ( n ) ( t ) d t ,
where Ψ j , j = 1 , 2 , , m , are nonlinear functionals defined on X, Equation (8) is written compactly as
K ¯ u = j = 1 m g j Ψ j ( u ) = g Ψ ( u ) .
Let the column vectors of linear bounded functionals
Φ ( u ) = Φ 1 ( u ) Φ 2 ( u ) Φ n ( u ) , Θ ( u ) = Θ 1 ( u ) Θ 2 ( u ) Θ l ( u ) ,
and the m × constant matrix
N = ν 11 ν 12 ν 1 l ν 21 ν 22 ν 2 l ν n 1 ν n 2 ν n l ,
where ν i j , i = 1 , 2 , , n , j = 1 , 2 , , l , are real constants, and write the nonlocal boundary conditions in the matrix form
Φ ( u ) N Θ ( u ) = 0 ,
where 0 denotes the zero column vector.
By employing (5), (6) and (11) define the linear Volterra-type integro-differential operator V : X X ,
V u = A u K u , D ( V ) = u : u C n [ a , b ] , Φ ( u ) N Θ ( u ) = 0 ,
and then, by means of (10), the nonlinear Volterra–Fredholm-type integro-differential operator T : X X ,
T u = V u g Ψ ( u ) , D ( T ) = D ( V ) .
Thus, the boundary value problem (3), (4) can be written in the operator form
T u = f , f X .
We now state the following theorem.
Theorem 1.
If linear operator V is bijective on X and its inverse is denoted by V 1 , then the exact solution to the problem T u = f for all f X is given by
u = V 1 f + V 1 g c * ,
for every vector c * = Ψ ( u ) that is a solution of the nonlinear algebraic (transcendental) system
c Ψ V 1 f + V 1 g c = 0 .
Proof. 
By multiplying the nonlinear equation
T u = V u g Ψ ( u ) = f ,
by inverse operator V 1 , we obtain
u = V 1 f + V 1 g Ψ ( u ) .
By applying vector Ψ on both sides, we have
Ψ ( u ) = Ψ V 1 f + V 1 g Ψ ( u ) .
We set c = Ψ ( u ) and write
c = Ψ V 1 f + V 1 g c ,
which is a nonlinear algebraic (transcendental) system with vector c unknown. Let c * be a solution of this system. By substituting c * into (15), we obtain the solution in (13). □

3. Solution of Linear VFIDEs

In the case in which φ j t , u ( t ) , u ( t ) , , u ( n ) ( t ) is a linear function in u ( t ) , u ( t ) , , u ( n ) ( t ) , namely,
φ j t , u ( t ) , u ( t ) , , u ( n ) ( t ) = i = 0 n φ j i ( t ) u ( n i ) ( t ) , j = 1 , 2 , , m ,
where φ j i ( t ) X , j = 1 , 2 , , m , i = 1 , 2 , , n , Equation (8) is simplified as
K ¯ u = j = 1 m g j ( x ) a b h j ( t ) i = 0 n φ j i ( t ) u ( n i ) ( t ) d t ,
and Equation (9) is written as
Ψ ( u ) = Ψ 1 ( u ) Ψ 2 ( u ) Ψ m ( u ) , Ψ j ( u ) = a b h j ( t ) i = 0 n φ j i ( t ) u ( n i ) ( t ) d t ,
where now Ψ j , j = 1 , 2 , , m , are linear bounded functionals defined on C n [ a , b ] (see [35] (Theorem 3, p. 480)).
We introduce the m × m matrix
Ψ ( g ) = Ψ 1 ( g 1 ) Ψ 1 ( g m ) Ψ m ( g 1 ) Ψ m ( g m ) ,
where element Ψ i ( g j ) is the value of functional Ψ i on element g j . It is noted that for a m × k constant matrix C ,
Ψ ( g C ) = Ψ ( g ) C .
Define the linear Volterra–Fredholm integro-differential operator S : X X as
S u = V u g Ψ ( u ) , D ( S ) = D ( V ) = u : u C n [ a , b ] , Φ ( u ) N Θ ( u ) = 0 ,
and write the linear boundary value problem (3), (4) and (16) in the symbolic form
S u = f , f X .
We now give the following theorem, where we use the notation I m to denote the m × m identity matrix.
Theorem 2.
If linear operator V is bijective on X, then linear operator S is injective if and only if
det W ¯ = det I m Ψ V 1 g 0 ,
and the unique solution of the boundary value problem
S u = f , for   all f X ,
is given by
u = V 1 f + V 1 g W ¯ 1 Ψ V 1 f .
Proof. 
Suppose that det W ¯ 0 and u ker S . Then, S u = V u g Ψ ( u ) = 0 ; hence,
u = V 1 g Ψ ( u ) .
By applying vector Ψ on both sides and using (17), we have
Ψ ( u ) = Ψ V 1 g Ψ ( u ) = Ψ V 1 g Ψ ( u ) ,
and thus
I m Ψ V 1 g Ψ ( u ) = W ¯ Ψ ( u ) = 0 ,
from where it follows that Ψ ( u ) = 0 . Then, from (20), it is implied that u = 0 , i.e., ker S = { 0 } ; therefore, operator S is injective. Conversely, we assume that S is injective, and we will show that det W ¯ 0 , or equivalently, we suppose that det W ¯ = 0 , and we will show that S is not injective. Let c be a nonzero constant vector c such that W ¯ c = 0 , and let element v = V 1 g c X , which is different from zero; otherwise,
W ¯ c = I m Ψ V 1 g c = c Ψ V 1 g c = c Ψ V 1 g c = c Ψ v = c = 0 ,
which is a contradiction. Then,
S v = V v g Ψ ( v ) = g c g Ψ ( V 1 g c ) = g c g Ψ ( V 1 g ) c = g I m Ψ V 1 g c = 0 ,
which means that ker S { 0 } ; therefore, S is not injective.
To find the unique solution of the boundary value problem (18), we multiply (18) by operator V 1 to obtain
u = V 1 f + V 1 g Ψ ( u ) .
By applying vector Ψ as above and solving with respect to Ψ ( u ) , we obtain
Ψ ( u ) = I m Ψ V 1 g 1 Ψ V 1 f = W ¯ 1 Ψ V 1 f .
By substituting (22) into (21), we obtain
u = V 1 f + V 1 g W ¯ 1 Ψ V 1 f ,
which is the solution in (19). □

4. Volterra Integro-Differential Operators of Convolution Type

In the previous two sections, we assume that inverse linear operator V 1 is known a priori. However, finding the inverse of a Volterra integro-differential operator is not always an easy task.
In this section, we look at the special case where X = C [ 0 , b ] , the Volterra integral operator in (6) is of convolution type, namely,
K u = i = 0 n 0 x k i ( x t ) u ( n i ) ( t ) d t ,
and the vector of functionals Φ ( u ) in (11) is as follows:
Φ ( u ) = c o l Φ 1 ( u ) , Φ 2 ( u ) , , Φ n ( u ) = c o l u ( 0 ) , u ( 0 ) , , u ( n 1 ) ( 0 ) .
In this case, the inverse operator V 1 of the linear Volterra integro-differential operator
V u = A u K u , D ( V ) = u : u C n [ 0 , b ] , Φ ( u ) N Θ ( u ) = 0
can be constructed analytically using the Laplace transform for certain classes of functions as it is shown below.
First, we recall the following definitions and properties about the Laplace transform: For a function f ( x ) C [ 0 , ) of exponential order, that is, there exist real constants γ and M such that | f ( x ) | M e γ x for all x 0 , the Laplace transform of f defined by
F ( s ) = L { f } = 0 e s x f ( x ) d x ,
exists for all s > γ . Laplace transform operator L is linear and injective. For any two functions f ( x ) and g ( x ) of exponential order and any real constants a and b,
L { a f + b g } = a L { f } + b L { g } ,
and
f ( x ) = L 1 { F } .
The convolution of f and g is defined as
( f g ) ( x ) = 0 x f ( x t ) g ( t ) d t ,
and holds
L { f g } = L { f } L { g } ) .
For a function f ( x ) C n [ 0 , ) and its derivatives f ( i ) ( x ) C [ 0 , ) , i = 1 , 2 , , n , of exponential order, the Laplace transform of the nth derivative f n ( x ) exists and is given by
L { f ( n ) } = s n L { f } = 1 n s n f ( 1 ) ( 0 ) ,
see, for example, [36] and some new developments in [37].
We prove the following lemma.
Lemma 1.
Let u ( i ) ( x ) , k i ( x ) , i = 0 , 1 , , n , and f ( x ) be continuous functions of exponential order. Then, the equation
A u K u = f ,
is equivalent to
u = p ^ Φ ( u ) + f ^ ,
where Φ ( u ) is the functional defined in (23), and
F ( s ) = L { f } , K i ( s ) = L { k i } , i = 0 , 1 , , n , Q ( s ) = 1 i = 0 n a i K i ( s ) s n i , P ( s ) = i = 0 n a i K i ( s ) s n i , = 1 , 2 , , n , f ^ = L 1 { Q ( s ) F ( s ) } , p ^ = p ^ 1 p ^ 2 p ^ n , p ^ = L 1 { Q ( s ) P ( s ) } , = 1 , 2 , , n .
Proof. 
By applying the Laplace transform operator to (26), we obtain
L { A u K u } = L { A u } L { K u } = i = 0 n a i L { u ( n i ) } i = 0 n L 0 x k i ( x t ) u ( n i ) ( t ) d t = i = 0 n a i L { u ( n i ) } i = 0 n L { k i } L { u ( n i ) } = i = 0 n a i L { k i } L { u ( n i ) } = L { f } .
From (25), we have
L { u ( n i ) } = s n i L { u } = 1 n i s n i u ( 1 ) ( 0 ) ,
and after substituting into (29), we obtain
L { A u K u } = i = 0 n a i L { k i } s n i L { u } = 1 n i s n i u ( 1 ) ( 0 ) = i = 0 n a i L { k i } s n i L { u } i = 0 n a i L { k i } = 1 n i s n i u ( 1 ) ( 0 ) = i = 0 n a i L { k i } s n i L { u } = 1 n i = 0 n a i L { k i } s n i u ( 1 ) ( 0 ) = L { f } ,
or
i = 0 n a i K i ( s ) s n i U ( s ) = 1 n i = 0 n a i K i ( s ) s n i u ( 1 ) ( 0 ) = F ( s ) ,
where U ( s ) = L { u } and F ( s ) = L { f } . By solving for U and using relations (28), we obtain
U ( s ) = = 1 n Q ( s ) P ( s ) u ( 1 ) ( 0 ) + Q ( s ) F ( s ) .
By applying then the inverse operator, we obtain
L 1 { U ( s ) } = = 1 n L 1 { Q ( s ) P ( s ) } u ( 1 ) ( 0 ) + L 1 { Q ( s ) F ( s ) } ,
and thus
u = = 1 n p ^ l u ( 1 ) ( 0 ) + f ^ = p ^ Φ ( u ) + f ^ ,
which is Equation (27). □
From Lemma 1, the next two results follow.
Lemma 2.
The linear Volterra integro-differential operator V 0 : X X defined as
V 0 u = A u K u , D ( V 0 ) = u : u C n [ 0 , b ] , Φ ( u ) = 0
is injective, and the unique solution of the initial value problem V 0 u = f , for any f X , is u = f ^ .
Lemma 3.
The Volterra integro-differential operator V 1 : X X defined as
V 1 u = A u K u , D ( V 1 ) = u : u C n [ 0 , b ] , Φ ( u ) = c ,
where c is a vector of arbitrary real constants, is injective, and the unique solution of initial value problem V 1 u = f , for any f X , is u = p ^ c + f ^ .
We now state the following theorem for finding the inverse ( V 1 ) of linear operator V in (24).
Theorem 3.
Let u ( i ) ( x ) , k i ( x ) , i = 0 , 1 , , n , and f ( x ) be continuous functions of exponential order. Then, the linear Volterra integro-differential operator
V u = A u K u , D ( V ) = u : u C n [ 0 , b ] , Φ ( u ) N Θ ( u ) = 0 ,
is injective if
det W ^ = det [ I n N Θ ( p ^ ) ] 0 ,
and the unique solution of the nonlocal boundary value problem V u = f , for any f X , is given by
u = p ^ W ^ 1 N Θ ( f ^ ) + f ^ .
Proof. 
From Equation (30) and Lemma 1, we have that A u K u = f is equivalent to
u = p ^ Φ ( u ) + f ^ ,
where Φ ( u ) is to be determined.
By applying vector Θ on (33), we obtain
Θ ( u ) = Θ p ^ Φ ( u ) + Θ ( f ^ ) = Θ ( p ^ ) Φ ( u ) + Θ ( f ^ ) .
By substituting Θ ( u ) into
Φ ( u ) N Θ ( u ) = 0 ,
after some algebra, we obtain
I n N Θ ( p ^ ) Φ ( u ) = N Θ ( f ^ ) .
If condition (31) holds true, then Equation (34) can be solved uniquely with respect to Φ ( u ) to obtain
Φ ( u ) = I n N Θ ( p ^ ) 1 N Θ ( f ^ ) .
The substitution of (35) into (33) yields
u = p ^ I n N Θ ( p ^ ) 1 N Θ ( f ^ ) + f ^ ,
which is Equation (32). □
Finally, for completeness, we mention that Theorem 3 can be extended to problems with nonhomogeneous boundary conditions. For this, let operator V ˜ : X X be defined as
V ˜ u = A u K u , D ( V ˜ ) = u : u C n [ 0 , b ] , Φ ( u ) N Θ ( u ) = b .
where K is as usual the convolution-type integral operator; the linear bounded functionals Φ ( u ) = c o l ( u ( 0 ) , u ( 0 ) , , u ( n 1 ) ( 0 ) ) ; and b is a column vector of n real constants. It is emphasized that V ˜ is not a linear operator since domain D ( V ˜ ) is a nonlinear set. With the same arguments as above, the following theorem can be proved.
Theorem 4.
Under the same conditions as in Theorem 3, the Volterra integro-differential operator V ˜ is injective, and the unique solution of the nonlocal boundary value problem V ˜ u = f , for any f X , is given by
u = p ^ W ^ 1 N Θ ( f ^ ) + b + f ^ .

5. Examples

In this section, we solve four examples representative of the presented theory to demonstrate the ease of implementation and effectiveness of the method. The first example concerns a Volterra integro-differential equation (VIDE) under various types of nonhomogeneous boundary conditions. The remaining examples are one linear and two nonlinear boundary value problems for VFIDEs under nonlocal boundary conditions.

5.1. Example 1

Consider the first-order linear Volterra integro-differential equation
u ( x ) + a 1 u ( x ) 0 x e ( x t ) u ( t ) d t = f ( x ) , x [ 0 , 1 ] ,
where a 1 R , f ( x ) is the input or load function assumed to be a continuous function of exponential order, subject to the following:
(i)
Initial condition
u ( 0 ) = b 1 ,
(ii)
Two-point boundary condition
u ( 0 ) ν 11 u ( x 1 ) = b 1 ,
(iii)
Multipoint and integral boundary condition
u ( 0 ) ν 11 u ( x 1 ) ν 12 u ( x 2 ) ν 13 ξ 1 ξ 2 u ( t ) d t = b 1 ,
where ν 1 l , l = 1 , 2 , 3 , and b 1 are real constants, 0 x 1 < x 2 1 and 0 ξ 1 < ξ 2 1 .
We take n = 1 and
A u = u ( x ) + a 1 u ( x ) , K u = 0 x k 1 ( x t ) u ( t ) d t , k 1 ( x ) = e x , V u = A u K u , D ( V ) = u : u C 1 [ 0 , 1 ] , Φ ( u ) N Θ ( u ) = b ,
where
(i)
Φ ( u ) = u ( 0 ) , N = 0 , b = b 1 ,
(ii)
Φ ( u ) = u ( 0 ) , N = ν 11 , Θ ( u ) = u ( x 1 ) , b = b 1 ,
(iii)
Φ ( u ) = u ( 0 ) , N = ν 11 ν 12 ν 13 , Θ ( u ) = u ( x 1 ) u ( x 2 ) ξ 1 ξ 2 u ( t ) d t , b = b 1 .
We note that k 1 is a continuous function of exponential order; therefore, we can apply Theorem 4. Thus, the solution of the initial value problem (37), (38) for f ( x ) = 0 and b 1 = 1 turns out to be
u ( x ) = e a 1 + 1 2 x cosh ( c x ) + 1 a 1 2 c sinh ( c x ) ,
where c = 1 2 a 1 2 2 a 1 + 5 . This problem was studied in [38,39].
The solution of the two-point boundary value problem (37), (39) for a 1 = 2 , f ( x ) = 2 e x + x 2 + 3 x 5 , ν 11 = 1 , b 1 = 0 , and x 1 = 1 is
u ( x ) = x 2 x 1 .
Finally, the solution of the nonlocal boundary value problem (37), (40) for
a 1 = 1 , f ( x ) = ( x + 2 ) e x , ν 11 = 1 , ν 12 = e 1 / 2 , ν 13 = 1 , b 1 = 1 + e 1 4 e 3 4 ,
and the points x 1 = 0.25 , x 2 = 0.5 , ξ 1 = 0.75 and ξ 2 = 1 is
u ( x ) = e x .

5.2. Example 2

Solve the second-order linear VFIDE
u ( x ) + a 1 u ( x ) + a 2 u ( x ) 0 x cos ( x t ) u ( t ) d t 0 x sin ( x t ) u ( t ) d t 0 1 ( x 2 + x t ) u ( t ) d t = f ( x ) , x [ 0 , 1 ] ,
subject to homogeneous boundary conditions
u ( 0 ) = 0 , u ( 1 ) = 0 .
First, we note that n = 2 , m = 2 and = 2 and define
A u = u ( x ) + a 1 u ( x ) + a 2 u ( x ) , K u = 0 x k 1 ( x t ) u ( t ) d t + 0 x k 2 ( x t ) u ( t ) d t , k 1 ( x ) = cos x , k 2 ( x ) = sin x , K ¯ u = x 2 0 1 u ( t ) d t + x 0 1 t u ( t ) d t , Φ ( u ) = u ( 0 ) u ( 0 ) , N = 0 0 1 1 , Θ ( u ) = u ( 0 ) u ( 1 ) .
Next, we set the vectors
g = g 1 ( x ) g 2 ( x ) = x 2 x , Ψ ( u ) = 0 1 h 1 ( t ) u ( t ) d t 0 1 h 2 ( t ) u ( t ) d t = 0 1 u ( t ) d t 0 1 t u ( t ) d t ,
and define the linear operators
V u = A u K u , D ( V ) = u : u C 2 [ 0 , 1 ] , Φ ( u ) N Θ ( u ) = 0 , S u = V u g Ψ ( u ) , D ( S ) = D ( V ) .
We observe that k 1 , k 2 , g 1 and g 2 are continuous functions of exponential order. Therefore, we use Theorems 2 and 3 to construct the unique solution of the boundary value problem (41), (42).
Specifically, for a 1 = 1 , a 2 = 1 and f ( x ) = x 2 + x , we obtain the unique solution
u ( x ) = 12 745 e 4611 238 e x ( 25 e + 5 ) x 3 ( 71 e 81 ) x 2 ( 142 e 162 ) x 238 .
For f ( x ) = e x , we obtain
u ( x ) = 1 745 e 4611 ( 745 e 4611 ) x + 2111 e 836 e x + 300 e 2 + 340 e + 1180 x 3 + 852 e 2 + 2406 e 672 x 2 + 1704 e 2 + 4812 e 1344 x 2111 e + 836 .

5.3. Example 3

Let the second-order nonlinear VFIDE of Hammerstein type
u ( x ) 0 x cos ( x t ) u ( t ) d t 120 x 0 1 u ( t ) u ( t ) d t = f ( x ) , x [ 0 , 1 ] ,
subject to multipoint boundary conditions
u ( 0 ) = 2 u ( 1 2 ) 3 u ( 1 ) , u ( 0 ) = u ( 1 ) .
We start by setting n = 2 , m = 1 and = 3 , and
A u = u ( x ) , K u = 0 x cos ( x t ) u ( t ) d t , K ¯ u = 120 x 0 1 u ( t ) u ( t ) d t , Φ ( u ) = u ( 0 ) u ( 0 ) , N = 2 3 0 0 0 1 , Θ ( u ) = u ( 1 2 ) u ( 1 ) u ( 1 ) .
Next, we define the vectors
g = g 1 ( x ) = 120 x , Ψ ( u ) = 0 1 u ( t ) u ( t ) d t ,
and the operators
V u = A u K u , D ( V ) = u : u C 2 [ 0 , 1 ] , Φ ( u ) N Θ ( u ) = 0 , T u = V u g Ψ ( u ) , D ( T ) = D ( V ) .
Theorems 1 and 3 are applicable. Let f ( x ) = x 2 . It follows that (14) is a nonlinear algebraic quadratic equation in c that has the following two real solutions:
c 1 , 2 = 2351 ± 24 8922 149 , 520 .
Consequently, we obtain the following two solutions of the nonlinear boundary value problem (44), (45):
u 1 , 2 ( x ) = 72 ± 8922 598 , 080 96 x 5 + 1504 x 3 2496 x + 3 .

5.4. Example 4

As a last example, we consider the first-order nonlinear VFIDE of Hammerstein type
u ( x ) 0 x cos ( t x ) u ( t ) d t 2 0 1 ( x + t ) u 2 ( t ) d t = f ( x ) , x [ 0 , 1 ] ,
subject to the boundary condition
u ( 0 ) = u ( 1 ) .
It is easily seen that n = 1 , m = 2 and = 1 , and
A u = u ( x ) , K u = 0 x cos ( t x ) u ( t ) d t , K ¯ u = 2 0 1 ( x + t ) u 2 ( t ) d t , Φ ( u ) = u ( 0 ) , N = 1 , Θ ( u ) = u ( 1 ) .
Define the vectors
g = g 1 ( x ) g 2 ( x ) = 2 x 1 , Ψ ( u ) = 0 1 h 1 ( t ) u 2 ( t ) d t 0 1 h 2 ( t ) u 2 ( t ) d t = 0 1 u 2 ( t ) d t 0 1 2 t u 2 ( t ) d t ,
and then the operators
V u = A u K u , D ( V ) = u : u C 1 [ 0 , 1 ] , Φ ( u ) N Θ ( u ) = 0 , T u = V u g Ψ ( u ) , D ( T ) = D ( V ) .
The solution of T u = f follows using Theorems 1 and 3. Let
f ( x ) = 7 sin x 4 cos x 331 105 x + 125 28 .
It turns out that the nonlinear algebraic system in (14) consists of two quadratic equations with unknown vector c . It possesses two real and two complex roots, of which only the real root
c = c 1 c 2 = 331 210 43 28
is retrieved in exact form. In this case, the solution of the nonlinear boundary value problem (46), (47) is
u ( x ) = x 3 3 x 2 + 2 x + 1 .
We programmed the proposed methods into the free software Maxima computer algebra system (CAS) running on Microsoft Windows. We used the available built-in functions laplace, ilt and solve to compute the Laplace transform, the inverse Laplace transform and the solution of linear and nonlinear systems of algebraic equations, respectively.

6. Discussion and Conclusions

In general, there are not many procedures available in the literature for solving boundary value problems with nonlocal boundary conditions in closed form.
In this paper, a method is presented to solve in closed form a class of linear and nonlinear VFIDEs under homogeneous nonlocal boundary conditions. The method requires the existence of the explicit form of the inverse operator of the associated linear Volterra integro-differential operator. A procedure for constructing this inverse operator in the case of convolution-type VIDEs under homogeneous and nonhomogeneous nonlocal boundary conditions is also given. The methods are implemented in a computer algebra system, and several examples are solved.
The main advantage of the technique is that it is straightforward, does not require approximations or numerical procedures and provides the exact solution in closed form. Moreover, in the nonlinear case, unlike the numerical solutions, it can give all possible solutions that can be calculated analytically. It was proved to be easy to use and very effective. Its main drawback is that its applicability is limited to a class of continuous functions of exponential order for which the direct and inverse Laplace transform exist. It relies on the available facilities for the inverse Laplace transform and explicit integration in the computer algebra software used.
Finally, the solution of VFIDEs subject to nonhomogeneous nonlocal boundary conditions is a more challenging task that needs a different treatment that could be addressed in a separate article in the future.

Author Contributions

Conceptualization, E.P.; formal analysis, E.P. and I.N.P.; writing—original draft preparation, E.P.; writing—review and editing, I.N.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wazwaz, A.M. Volterra-Fredholm Integro-Differential Equations. In Linear and Nonlinear Integral Equations; Springer: Berlin/Heidelberg, Germany, 2011; pp. 285–309. [Google Scholar] [CrossRef]
  2. Rahman, M.; Jackiewicz, Z.; Welfert, B.D. Stochastic approximations of perturbed Fredholm Volterra integro-differential equation arising in mathematical neurosciences. Appl. Math. Comput. 2007, 186, 1173–1182. [Google Scholar] [CrossRef]
  3. Amin, R.; Nazir, S.; García-Magariño, I. A collocation method for numerical solution of nonlinear delay integro-differential equations for wireless sensor network and internet of things. Sensors 2020, 20, 1962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Providas, E. On the exact solution of nonlocal Euler–Bernoulli beam equations via a direct approach for Volterra-Fredholm integro-differential equations. AppliedMath 2022, 2, 269–283. [Google Scholar] [CrossRef]
  5. Turkyilmazoglu, M. rigorous power series expansion technique High-order nonlinear Volterra–Fredholm-Hammerstein integro-differential equations and their effective computation. Appl. Math. Comput. 2014, 247, 410–416. [Google Scholar] [CrossRef]
  6. Cherruault, Y. A reliable method for obtaining approximate solutions of linear and nonlinear Volterra-Fredholm integro-differential equations. Kybernetes 2005, 34, 1034–1048. [Google Scholar] [CrossRef]
  7. Hamoud, A.A.; Ghadle, K.P. Existence and uniqueness of the solution for Volterra-Fredholm integro-differential equations. J. Sib. Fed. Univ. Math. Phys. 2018, 11, 692–701. [Google Scholar]
  8. Yalçinbaş, S.; Sezer, M. The approximate solution of high-order linear Volterra–Fredholm integro-differential equations in terms of Taylor polynomials. Appl. Math. Comput. 2000, 112, 291–308. [Google Scholar] [CrossRef]
  9. Maleknejad, K.; Mahmoudi, Y. Taylor polynomial solution of high-order nonlinear Volterra–Fredholm integro-differential equations. Appl. Math. Comput. 2003, 145, 641–653. [Google Scholar] [CrossRef]
  10. Darania, P.; Ivaz, K. Numerical solution of nonlinear Volterra–Fredholm integro-differential equations. Comput. Math. Appl. 2008, 56, 2197–2209. [Google Scholar] [CrossRef] [Green Version]
  11. Sohaib, M.; Haq, S. An efficient wavelet-based method for numerical solution of nonlinear integral and integro-differential equations. Math. Methods Appl. Sci. 2020, 1–15. [Google Scholar] [CrossRef]
  12. Maleknejad, K.; Basirat, B.; Hashemizadeh, E. Hybrid Legendre polynomials and Block-Pulse functions approach for nonlinear Volterra–Fredholm integro-differential equations. Comput. Math. Appl. 2011, 61, 2821–2828. [Google Scholar] [CrossRef] [Green Version]
  13. Căruntu, B.; Paşca, M.S. Approximate Solutions for a Class of Nonlinear Fredholm and Volterra Integro-Differential Equations Using the Polynomial Least Squares Method. Mathematics 2021, 9, 2692. [Google Scholar] [CrossRef]
  14. Asgari, M.; Mesforush, A.; Nazemi, A. The numerical method for solving Volterra–Fredholm integro-differential equations of the second kind based on the meshless method. Asian-Eur. J. Math. 2022, 15, 2250002. [Google Scholar] [CrossRef]
  15. Shahmorad, S. Numerical solution of the general form linear Fredholm–Volterra integro-differential equations by the Tau method with an error estimation. Appl. Math. Comput. 2005, 167, 1418–1429. [Google Scholar] [CrossRef]
  16. Momani, S.; Arqub, O.A.; Hayat, T.; Al-Sulami, H. A computational method for solving periodic boundary value problems for integro-differential equations of Fredholm–Volterra type. Appl. Math. Comput. 2014, 240, 229–239. [Google Scholar] [CrossRef]
  17. Babolian, E.; Masouri, Z.; Hatamzadeh-Varmazyar, S. Numerical solution of nonlinear Volterra–Fredholm integro-differential equations via direct method using triangular functions. Comput. Math. Appl. 2009, 58, 239–247. [Google Scholar] [CrossRef] [Green Version]
  18. Berenguer, M.I.; Gámez, D.; López Linares, A.J. Fixed-Point Iterative Algorithm for the Linear Fredholm-Volterra Integro-Differential Equation. J. Appl. Math. 2012, 2012, 370894. [Google Scholar] [CrossRef] [Green Version]
  19. Berenguer, M.I.; Gámez, D.; López Linares, A.J. Fixed point techniques and Schauder bases to approximate the solution of the first order nonlinear mixed Fredholm–Volterra integro-differential equation. J. Comput. Appl. Math. 2013, 252, 52–61. [Google Scholar] [CrossRef]
  20. Touati, S.; Lemita, S.; Ghiat, M.; Aissaoui, M.Z. Solving a nonlinear Volterra-Fredholm integro-differential equation with weakly singular kernels. Fasc. Math. 2019, 62, 155–168. [Google Scholar] [CrossRef]
  21. Wang, Y.; Zhu, L. SCW method for solving the fractional integro-differential equations with a weakly singular kernel. Appl. Math. Comput. 2016, 275, 72–80. [Google Scholar] [CrossRef]
  22. Rahimkhani, P.; Ordokhani, Y.; Babolian, E. Fractional-order Bernoulli functions and their applications in solving fractional Fredholem–Volterra integro-differential equations. Appl. Numer. Math. 2017, 122, 66–81. [Google Scholar] [CrossRef]
  23. Dehestani, H.; Ordokhani, Y.; Razzaghi, M. Combination of Lucas wavelets with Legendre–Gauss quadrature for fractional Fredholm–Volterra integro-differential equations. J. Comput. Appl. Math. 2021, 382, 113070. [Google Scholar] [CrossRef]
  24. Yüzbaşı, S.; Sezer, M.; Kemancı, B. Numerical solutions of integro-differential equations and application of a population model with an improved Legendre method. Appl. Math. Model. 2013, 37, 2086–2101. [Google Scholar] [CrossRef]
  25. Rohaninasab, N.; Maleknejad, K.; Ezzati, R. Numerical solution of high-order Volterra–Fredholm integro-differential equations by using Legendre collocation method. Appl. Math. Comput. 2018, 328, 171–188. [Google Scholar] [CrossRef]
  26. Acar, N.I.; Daşcıoğlu, A. A projection method for linear Fredholm–Volterra integro-differential equations. J. Taibah Univ. Sci. 2019, 13, 644–650. [Google Scholar] [CrossRef] [Green Version]
  27. Hesameddini, E.; Shahbazi, M. Solving multipoint problems with linear Volterra–Fredholm integro-differential equations of the neutral type using Bernstein polynomials method. Appl. Numer. Math. 2019, 136, 122–138. [Google Scholar] [CrossRef]
  28. Reutskiy, S.Y. The backward substitution method for multipoint problems with linear Volterra–Fredholm integro-differential equations of the neutral type. J. Comput. Appl. Math. 2016, 296, 724–738. [Google Scholar] [CrossRef]
  29. Aida-zade, K.R.; Abdullaev, V.M. On the numerical solution of loaded systems of ordinary differential equations with nonseparated multipoint and integral conditions. Numer. Analys. Appl. 2014, 7, 1–14. [Google Scholar] [CrossRef]
  30. Sidorov, N.A.; Sidorov, D.N. Nonlinear Volterra Equations with Loads and Bifurcation Parameters: Existence Theorems and Construction of Solutions. Diff. Equat. 2021, 57, 1640–1651. [Google Scholar] [CrossRef]
  31. Polyanin, A.D.; Manzhirov, A.V. Handbook of Integral Equations, 2nd ed.; Chapman and Hall/CRC: New York, NY, USA, 2008. [Google Scholar] [CrossRef]
  32. Baiburin, M.M.; Providas, E. Exact Solution to Systems of Linear First-Order Integro-Differential Equations with Multipoint and Integral Conditions. In Mathematical Analysis and Applications; Springer Optimization and Its Applications; Rassias, T.M., Pardalos, P.M., Eds.; Springer: Cham, Switzerland, 2019; Volume 154, pp. 1–16. [Google Scholar] [CrossRef]
  33. Providas, E. An Algorithm for the Closed-Form Solution of Certain Classes of Volterra–Fredholm Integral Equations of Convolution Type. Algorithms 2022, 15, 203. [Google Scholar] [CrossRef]
  34. Providas, E. Closed-Form Solution of the Bending Two-Phase Integral Model of Euler-Bernoulli Nanobeams. Algorithms 2022, 15, 151. [Google Scholar] [CrossRef]
  35. Parasidis, I.N.; Providas, E. Extension Operator Method for the Exact Solution of Integro-Differential Equations. In Contributions in Mathematics and Engineering; Pardalos, P., Rassias, T., Eds.; Springer: Cham, Switzerland, 2016; pp. 473–496. [Google Scholar] [CrossRef]
  36. Kreyszig, E. Advanced Engineering Mathematics, 7th ed.; John Wiley & Sons. Inv.: Singapore, 1993. [Google Scholar]
  37. Rezaei, H.; Jung, S.M.; Rassias, T.M. Laplace transform and Hyers–Ulam stability of linear differential equations. J. Math. Anal. Appl. 2013, 403, 244–251. [Google Scholar] [CrossRef]
  38. Ma, J.; Brunner, H. A posteriori error estimates of discontinuous Galerkin methods for non-standard Volterra integro-differential equations. IMA J. Numer. Anal. 2006, 26, 78–95. [Google Scholar] [CrossRef]
  39. Hale, N. An ultraspherical spectral method for linear Fredholm and Volterra integro-differential equations of convolution type. IMA J. Numer. Anal. 2019, 39, 1727–1746. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Providas, E.; Parasidis, I.N. A Symbolic Method for Solving a Class of Convolution-Type Volterra–Fredholm–Hammerstein Integro-Differential Equations under Nonlocal Boundary Conditions. Algorithms 2023, 16, 36. https://doi.org/10.3390/a16010036

AMA Style

Providas E, Parasidis IN. A Symbolic Method for Solving a Class of Convolution-Type Volterra–Fredholm–Hammerstein Integro-Differential Equations under Nonlocal Boundary Conditions. Algorithms. 2023; 16(1):36. https://doi.org/10.3390/a16010036

Chicago/Turabian Style

Providas, Efthimios, and Ioannis Nestorios Parasidis. 2023. "A Symbolic Method for Solving a Class of Convolution-Type Volterra–Fredholm–Hammerstein Integro-Differential Equations under Nonlocal Boundary Conditions" Algorithms 16, no. 1: 36. https://doi.org/10.3390/a16010036

APA Style

Providas, E., & Parasidis, I. N. (2023). A Symbolic Method for Solving a Class of Convolution-Type Volterra–Fredholm–Hammerstein Integro-Differential Equations under Nonlocal Boundary Conditions. Algorithms, 16(1), 36. https://doi.org/10.3390/a16010036

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