Next Article in Journal
Multi-Criteria Decision Making Approach for Hybrid Operation of Wind Farms
Next Article in Special Issue
Geometric Properties of Certain Classes of Analytic Functions Associated with a q-Integral Operator
Previous Article in Journal
Temporal-Spatial Nonlinear Filtering for Infrared Focal Plane Array Stripe Nonuniformity Correction
Previous Article in Special Issue
Generalized Mittag-Leffler Input Stability of the Fractional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence of Solution for Non-Linear Functional Integral Equations of Two Variables in Banach Algebra

1
Department of Mathematics and Statistics, University of Victoria, Victoria, BC V8W 3R4, Canada
2
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
3
Department of Mathematics, Rajiv Gandhi University, Rono Hills, Doimukh 791112, Arunachal Pradesh, India
4
Department of Mathematics, Gauhati University, Guwahati 781014, Assam, India
5
Operator Theory and Applications Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(5), 674; https://doi.org/10.3390/sym11050674
Submission received: 30 March 2019 / Revised: 2 May 2019 / Accepted: 5 May 2019 / Published: 16 May 2019
(This article belongs to the Special Issue Integral Transformations, Operational Calculus and Their Applications)

Abstract

:
The aim of this article is to establish the existence of the solution of non-linear functional integral equations x ( l , h ) = U ( l , h , x ( l , h ) ) + F l , h , 0 l 0 h P ( l , h , r , u , x ( r , u ) ) d r d u , x ( l , h ) × G l , h , 0 a 0 a Q l , h , r , u , x ( r , u ) d r d u , x ( l , h ) of two variables, which is of the form of two operators in the setting of Banach algebra C [ 0 , a ] × [ 0 , a ] , a > 0 . Our methodology relies upon the measure of noncompactness related to the fixed point hypothesis. We have used the measure of noncompactness on C [ 0 , a ] × [ 0 , a ] and a fixed point theorem, which is a generalization of Darbo’s fixed point theorem for the product of operators. We additionally illustrate our outcome with the help of an interesting example.

1. Introduction

Many real-life problems in which we go over the investigation of various branches of mathematical physics, for example, gas kinetic theory, radiation, and neutron transportation, can be depicted and demonstrated by methods of non-linear functional integral equations (for example, we refer to [1,2,3,4]). Banaś and Lecko [5] introduced the concept of fixed points of product operators in Banach algebra. Dhage [6,7] used the concept of the fixed point theorem to find the solution of functional integral equations in Banach algebra. Banaś and Olszowy [8] used the class of measures of noncompactness to obtain the existence of solutions of nonlinear integral equations in Banach algebra. Deepmala and Pathak [9] studied the existence of the solution of nonlinear functional integral equations of a single variable in Banach algebra C [ a , b ] of all real-valued continuous functions on the interval [ a , b ] equipped with the maximum norm.
Kuratowski [10], in the year 1930, first introduced the idea of the measure of noncompactness (denoted by “ α ”). For any bounded subset A of a metric space X ,
α A = inf δ > 0 : A j = 1 m A j , A j X , diam A j < δ ( j = 1 , , m ) , m N ,
where:
diam A j = sup d ( a 1 , a 2 ) : a 1 , a 2 A j .
Using this idea, Darbo [11] exhibited a fixed point theorem that plays a very significant role in the finding of existence theorems. In the recent past, there have been a few fruitful endeavours to apply the idea of the measure of noncompactness in the investigation of the existence of solutions for various kinds of differential and integral equations, for example one can refer to [12,13,14,15,16,17,18,19,20].
In many physical problems, we come across nonlinear integral equations. The fixed point theory plays a significant role to obtain the solutions of such equations. Deepmala and Pathak, in [9], studied the following nonlinear functional integral equation, which can be considered as a particular case of many nonlinear functional integral equations that are applicable in mechanics, physics, economics, etc.,
x ( t ) = u ( t , x ( t ) ) + f t , 0 t p ( t , s , x ( s ) ) d s , x ( α ( t ) ) × g t , 0 a q t , s , x ( s ) d s , x ( β ( t ) ) for t [ 0 , a ] .
The authors of [9] used the measure of noncompactness to obtain the existence of the solution of the integral Equation (1) in Banach algebra C [ 0 , a ] with the help of the fixed point theorem.
Motivated by the work of [9], in this article, we study the solvability of non-linear functional integral equations of two variables, which we come across in various branches of nonlinear analysis. We consider an integral equation in the following form:
x ( l , h ) = U ( l , h , x ( l , h ) ) + F l , h , 0 l 0 h P ( l , h , r , u , x ( r , u ) ) d r d u , x ( l , h ) × G l , h , 0 a 0 a Q l , h , r , u , x ( r , u ) d r d u , x ( l , h ) for l , h [ 0 , a ] .
The right-hand side of the above integral equation that we are considering is the product of two functional operators involving integral operators and applying a fixed point theorem, which is a generalization of Darbo’s fixed point theorem for the product of operators to check the existence of the solution of the integral equation in Banach algebra. It can be seen that Equation (2) is a generalization of Equation (1) in two variables. Here, we used a fixed point theorem associated with Darbo’s condition of the measure of noncompactness in Banach algebra of continuous functions in [ 0 , a ] × [ 0 , a ] to establish the solvability of Equation (2). Furthermore, we used the modified homotopy perturbation analytic method to find the solution of Equation (2).

2. Preliminaries

Let R denote the set of real numbers, and write R + = 0 , . Suppose E ¯ is a real Banach space with the norm . , and let X ( ϕ ) E ¯ . The closure and convex closure of X will be denoted by X ¯ and convX, respectively. The convex closure of a set X of points in the Euclidean plane or in a Euclidean space over the reals is the smallest convex set that contains X . A closed ball in E ¯ centred at a and with radius b is denoted by B ( a , b ) . In addition, we use the symbol M E ¯ to denote the family of all non-empty and bounded subsets of E ¯ and use N E ¯ to denote its subfamily consisting of all relatively compact sets.
Definition 1.
Let X be a linear space over R . A norm on X is a function from X to R + , commonly denoted | | . | | such that:
(N1) 
| | x | | 0 and | | x | | = 0 x = 0 ;
(N2) 
| | α x | | = | α | | | x | | ;
(N3) 
| | x + y | | | | x | | + | | y | | for all x , y X and α R .
The pair ( X , | | . | | ) is called a normed space. A complete normed space is called a Banach space.
Definition 2.
An algebra A is a vector space A over a field K such that for each ordered pair of elements x , y A , a unique product x y A is defined with the properties:
(A1) 
( x y ) z = x ( y z ) ,
(A2) 
x ( y + z ) = x y + x z ,
(A3) 
( x + y ) z = x z + y z ,
(A4) 
α ( x y ) = ( α x ) y = x ( α y ) for all x , y , z A and scalars α .
A normed algebra A is normed space, which is an algebra such that for all x , y A :
x y x y
and if A has an identity e, then e = 1 .
A Banach algebra is a normed algebra that is complete, considered as a normed space.
The notion of the measure of noncompactness due to Banaś and Goebel [21] is as follows:
Definition 3.
A function μ : M E ¯ 0 , is said to be a measure of noncompactness in E ¯ if:
(i) 
for all X M E ¯ , we have that μ ( X ) = 0 implies that X is precompact.
(ii) 
the family ker μ = X M E ¯ : μ X = 0 is non-empty, and ker μ N E ¯ .
(iii) 
X Z μ X μ Z .
(iv) 
μ X ¯ = μ X .
(v) 
μ c o n v X = μ X where c o n v X is the convex closure of set X .
(vi) 
μ λ X + 1 λ Z λ μ X + 1 λ μ Z for λ 0 , 1 .
(vii) 
if X n M E ¯ , X n = X ¯ n , X n + 1 X n for n = 1 , 2 , 3 , and lim n μ X n = 0 , then n = 1 X n ϕ .
The family ker μ is called the kernel of measure μ . Note that the intersection set X from the above condition (vii) is a member of the family ker μ . Since μ ( X ) μ ( X n ) for any n , we deduce μ ( X ) = 0 . Consequently, X k e r μ .
For given subsets X , Y in a Banach algebra E , the product X Y defined by:
X Y = { x y : x X , y Y } .
In [8], Banaś and Olszowy defined the measure of noncompactness μ on the Banach algebra E, which satisfies condition ( m ) if for arbitrary sets X , Y M E such that:
μ ( X Y ) | | X | | μ ( Y ) + | | y | | μ ( X ) .
Deepmala and Pathak [9] used this concept of measure of noncompactness and obtained the existence of the solution of Equation (1).
Definition 4
([21]). Let E be a Banach space. Consider a non-empty subset X of E and a continuous operator T : X E transforming the bounded subset of X to the bounded ones. We say that T satisfies the Darbo condition with a constant k with respect to measure μ provided μ ( T Y ) k μ ( Y ) for each Y M E such that Y X . If k < 1 , then T is called a contraction with respect to μ .
Remark 1.
The Darbo condition has many applications, particularly in fixed point theorems, which can be applied to check the existence of the solution of different types of integral, differential, and integro differential equations. The Darbo condition can be potentially applied to extend the linear space in the work of Shang [22]. The assumptions (1)–(4) of the next section have been utilized in the study of consensus problems (see [23,24]).
We recall the following important theorems:
Theorem 1
([11]). Assume that Z is a non-empty, closed, bounded, and convex subset of a Banach space E ¯ . Let S : Z Z be a continuous mapping. Suppose that there is a constant k 0 , 1 such that:
μ ( S M ) k μ ( M ) , M Z .
Then, S has a fixed point.
Theorem 2
([8]). Suppose that X is a non-empty, bounded, convex, and closed subset of a Banach algebra E , and the operators P and T transform continuously the set X into E such that P ( X ) and T ( X ) are bounded. Furthermore, suppose that the operator S = P . T transforms X into itself. If P and T satisfy on the set X the Darbo condition with respect to the measure of noncompactness μ with the constants k 1 and k 2 , respectively, then S satisfies on X the Darbo condition with constant P ( X ) k 2 + T ( X ) k 1 . Particularly, if:
P ( X ) k 2 + T ( X ) k 1 < 1 ,
then S is a contraction with respect to the measure of noncompactness μ and has at least one fixed point in X .
We consider the space E = C ( [ 0 , a ] × [ 0 , a ] ) , which consists of the set of real-valued continuous functions on [ 0 , a ] × [ 0 , a ] . It is obvious that E is the vector space over the field of scalars R with the following operations:
( x + y ) ( t , s ) = x ( t , s ) + y ( t , s )
and:
( α x ) ( t , s ) = α x ( t , s ) ,
where x , y E , α R and t , s [ 0 , a ] . Since x , y E , i.e., both x , y are real-valued continuous functions on [ 0 , a ] × [ 0 , a ] and the product of two real-valued continuous functions is also a real-valued continuous function, therefore x y E , where:
( x y ) ( t , s ) = x ( t , s ) y ( t , s ) , t , s [ 0 , a ] .
Let z E . For all t , s [ 0 , a ] ,
x y z ( t , s ) = x y ( t , s ) z ( t , s ) = x ( t , s ) y ( t , s ) z ( t , s ) = x ( t , s ) y z ( t , s ) = x y z ( t , s ) .
Since t , s are arbitrary, therefore ( x y ) z = x ( y z ) .
Similarly, it can be shown that:
x ( y + z ) = x y + x z ,
( x + y ) z = x z + y z
and:
α x y = α x y = x α y .
Therefore E is an algebra.
The space E is also a normed space with the norm:
x = sup x ( l , h ) : l , h [ 0 , a ] , a > 0 , x E .
For all x , y E and l , h [ 0 , a ] ,
( x y ) ( l , h ) = x ( l , h ) y ( l , h ) = x ( l , h ) y ( l , h )
and so:
sup l , h [ 0 , a ] ( x y ) ( l , h ) sup l , h [ 0 , a ] x ( l , h ) sup l , h [ 0 , a ] y ( l , h ) ,
i.e.,
x y x y .
Thus, E is a normed algebra.
Let x n ( t , s ) n = 1 be a Cauchy sequence in E where x n ( t , s ) R × R for all n N and t , s [ 0 , a ] . Then:
x n x m 0 ( n , m ) .
Therefore, for all t , s [ 0 , a ] , we get:
x n ( t , s ) x m ( t , s ) 0 ( n , m ) .
For fixed t , s [ 0 , a ] , the sequence x n ( t , s ) is a Cauchy sequence of real numbers, so it is a convergent sequence and converging to x 0 ( t , s ) E (say) as the limit of the continuous function is also continuous. Therefore, for all t , s [ 0 , a ] :
x n ( t , s ) x 0 ( t , s ) 0 ( n , m )
which yields:
sup t , s [ 0 , a ] x n ( t , s ) x 0 ( t , s ) 0 ( n , m )
Thus:
x n x 0 0 ( n )
which proves that E is complete normed space. Hence, we conclude that the space E has the Banach algebra structure.
Let X be a fixed non-empty and bounded subset of E = C ( [ 0 , a ] × [ 0 , a ] ) , and for x X and ϵ > 0 , the modulus of the continuity function (denoted by ω ( x , ϵ ) ) is given by the formula:
ω ( x , ϵ ) = sup x ( l , h ) x ( v , w ) : l , h , v , w [ 0 , a ] , l v ϵ , h w ϵ .
Further, we define:
ω ( X , ϵ ) = sup ω ( x , ϵ ) : x X , ω 0 ( X ) = lim ϵ 0 ω ( x , ϵ ) .
Similar to [5], it can be shown that the function ω 0 ( X ) is a regular measure of non-compactness in the space C ( [ 0 , a ] × [ 0 , a ] ) . Apart from this, it is easy to check that the measure ω 0 ( X ) satisfies condition ( m ) .

3. Main Result

In this section, we study the existence of solutions of the integral Equation (2). We consider the following assumptions:
(1)
The functions U : [ 0 , a ] × [ 0 , a ] × R R , F : [ 0 , a ] × [ 0 , a ] × R × R R , and G : [ 0 , a ] × [ 0 , a ] × R × R R are continuous, and there exist nonnegative constants L , M such that:
U ( l , h , 0 ) L , F ( l , h , M 1 , 0 ) M a n d F ( l , h , M 2 , 0 ) M ,
where M 1 , M 2 R .
(2)
Let A i : [ 0 , a ] × [ 0 , a ] R + ( i = 1 , 2 , 3 , 4 , 5 ) be continuous functions such that:
U ( l , h , x 1 ) U ( l , h , x 2 ) A 1 ( l , h ) x 1 x 2 ,
F ( l , h , y , x 1 ) F ( l , h , y , x 2 ) A 2 ( l , h ) x 1 x 2 ,
G ( l , h , y , x 1 ) G ( l , h , y , x 2 ) A 3 ( l , h ) x 1 x 2 ,
F ( l , h , y 1 , x ) F ( l , h , y 2 , x ) A 4 ( l , h ) y 1 y 2
and:
G ( l , h , y 1 , x ) G ( l , h , y 2 , x ) A 5 ( l , h ) y 1 y 2 ,
where l , h [ 0 , a ] and x , x 1 , x 2 , y , y 1 , y 2 R . Furthermore, let:
K = max { A i ( l , h ) : i = 1 , 2 , 3 , 4 , 5 ; l , h [ 0 , a ] } ,
where K 0 .
(3)
The functions P , Q are continuous functions from [ 0 , a ] × [ 0 , a ] × [ 0 , a ] × [ 0 , a ] × R to R .
(4)
Furthermore, 4 α β < 1 for α = 2 k , β = L + M .
Theorem 3.
Under the hypotheses (1)–(4), Equation (2) has at least one solution in E = C ( I × I ) , where I = [ 0 , a ] .
Proof. 
Let us consider the operators F ^ and G ^ defined on E by:
( F ^ x ) ( l , h ) = U ( l , h , x ( l , h ) ) + F l , h , 0 l 0 h P ( l , h , r , u , x ( r , u ) ) d r d u , x ( l , h )
and:
( G ^ x ) ( l , h ) = G l , h , 0 a 0 a Q l , h , r , u , x ( r , u ) d r d u , x ( l , h ) , where l , h [ 0 , a ] .
From Assumptions (1)–(3), we get that F ^ and G ^ map C ( I × I ) into itself. Furthermore, let us define another operator T ^ on C ( I × I ) as follows:
T ^ x = ( F ^ x ) ( G ^ x ) .
It is obvious that T ^ maps C ( I × I ) into itself.
Let:
I 1 ( x ) = 0 l 0 h P ( l , h , r , u , x ( r , u ) ) d r d u
and:
I 2 ( x ) = 0 a 0 a Q l , h , r , u , x ( r , u ) d r d u .
Let x C ( I × I ) be fixed and l , h I . We get:
( T ^ x ) ( l , h ) = ( F ^ x ) ( l , h ) . ( G ^ x ) ( l , h ) = U ( l , h , x ( l , h ) ) + F l , h , I 1 ( x ) , x ( l , h ) × G l , h , I 2 ( x ) , x ( l , h ) ( U ( l , h , x ( l , h ) ) U ( l , h , 0 ) + U ( l , h , 0 ) + | F l , h , I 1 ( x ) , x ( l , h ) F l , h , I 1 ( x ) , 0 | + F l , h , I 1 ( x ) , 0 ) × ( | G ( l , h , I 2 ( x ) , x ( l , h ) ) G ( l , h , I 2 ( x ) , 0 ) | + G ( l , h , I 2 ( x ) , 0 ) ) A 1 ( l , h ) x ( l , h ) + L + A 2 ( l , h ) x ( l , h ) + M × A 3 ( l , h ) x ( l , h ) + M 2 K x + L + M K x + M 2 K x + L + M 2 .
Let α = 2 k , β = L + M . Then, we have:
F ^ x α x + β ,
G ^ x α x + β
and:
T ^ x α x + β 2
for x C ( I × I ) .
From (3), we have that the operator T ^ maps B d C ( I × I ) into B d , where:
B d = x ( l , h ) I : x ( l , h ) d
for d 2 d d 1 , where:
d 1 = 1 2 α β 1 4 α β 2 α 2
and:
d 2 = 1 2 α β + 1 4 α β 2 α 2 .
Furthermore, we have:
F ^ B d α d + β
and:
G ^ B d α d + β .
Let ϵ > 0 be fixed and x ( l , h ) , y ( l , h ) B d such that:
x y ϵ , ( l , h I ) .
Then, we have:
( F ^ x ) ( l , h ) ( F ^ y ) ( l , h ) = | U l , h , x ( l , h ) + F l , h , I 1 ( x ) , x ( l , h ) U l , h , y ( l , h ) F l , h , I 1 ( y ) , y ( l , h ) | U l , h , x ( l , h ) U l , h , y ( l , h ) + F l , h , I 1 ( x ) , x ( l , h ) F l , h , I 1 ( x ) , y ( l , h ) + F l , h , I 1 ( x ) , y ( l , h ) F l , h , I 1 ( y ) , y ( l , h ) A 1 ( l , h ) x ( l , h ) y ( l , h ) + A 2 ( l , h ) x ( l , h ) y ( l , h ) + A 4 ( l , h ) I 1 ( x ) I 1 ( y ) 2 K x y + K 0 l 0 h P ( l , h , r , u , x ( r , u ) ) P ( l , h , r , u , y ( r , u ) ) d r d u 2 K x y + K a 2 ω ( P , ϵ ) ,
where:
ω ( P , ϵ ) = sup P ( l , h , r , u , x ( r , u ) ) P ( l , h , r , u , y ( r , u ) ) : l , h , r , u I , x , y [ d , d ] , x y < ϵ .
Since P is continuous, so it is uniformly continuous on the compact set I × I × I × I × [ d , d ] ; therefore:
ω ( P , ϵ ) 0 a s ϵ 0 .
Thus, F ^ is continuous on B d . Similarly, one can prove that G ^ is continuous on B d . Thus, we can conclude that T ^ is continuous on B d .
Let us consider a non-empty subset X of B d and x X . Then, for a fixed ϵ > 0 and l 1 , l 2 , h 1 , h 2 I such that l 1 l 2 , h 1 h 2 , l 1 l 2 ϵ , h 1 h 2 ϵ , one obtains:
( F ^ x ) ( l 2 , h 2 ) ( F ^ x ) ( l 1 , h 1 ) = | U l 2 , h 2 , x ( l 2 , h 2 ) + F l 2 , h 2 , 0 l 2 0 h 2 P ( l 2 , h 2 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) U l 1 , h 1 , x ( l 1 , h 1 ) F l 1 , h 1 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 1 , h 1 ) | U l 2 , h 2 , x ( l 2 , h 2 ) U l 2 , h 2 , x ( l 1 , h 1 ) + U l 2 , h 2 , x ( l 1 , h 1 ) U l 1 , h 1 , x ( l 1 , h 1 ) + | F l 2 , h 2 , 0 l 2 0 h 2 P ( l 2 , h 2 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) F l 2 , h 2 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) | + | F l 2 , h 2 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) F l 1 , h 1 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) | + | F l 1 , h 1 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) F l 1 , h 1 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 1 , h 1 ) | A 1 ( l , h ) x ( l 2 , h 2 ) x ( l 1 , h 1 ) + U l 2 , h 2 , x ( l 1 , h 1 ) U l 1 , h 1 , x ( l 1 , h 1 ) + A 4 ( l , h ) 0 l 2 0 h 2 P ( l 2 , h 2 , r , u , x ( r , u ) ) d r d u 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u + | F l 2 , h 2 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) F l 1 , h 1 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) | + A 2 ( l , h ) x ( l 2 , h 2 ) x ( l 1 , h 1 )
which yields:
( F ^ x ) ( l 2 , h 2 ) ( F ^ x ) ( l 1 , h 1 ) 2 K x ( l 2 , h 2 ) x ( l 1 , h 1 ) + U l 2 , h 2 , x ( l 1 , h 1 ) U l 1 , h 1 , x ( l 1 , h 1 ) + K | 0 l 2 0 h 2 P ( l 2 , h 2 , r , u , x ( r , u ) ) d r d u
0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u | + | F l 2 , h 2 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) F l 1 , h 1 , 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u , x ( l 2 , h 2 ) | .
Let:
ω ( U , ϵ ) = sup U l 2 , h 2 , x ( l 2 , h 2 ) U l 1 , h 1 , x ( l 1 , h 1 ) : l 1 , l 2 , h 1 , h 2 I , l 2 l 1 ϵ , h 2 h 1 ϵ , x [ d , d ] ,
ω ( P , ϵ ) = sup P l 2 , h 2 , r , u , x ( r , u ) P l 1 , h 1 , r , u , x ( r , u ) : l 1 , l 2 , h 1 , h 2 , r , u I , l 2 l 1 ϵ , h 2 h 1 ϵ , x [ d , d ] ,
k ¯ = sup P ( l , h , r , u , x ( r , u ) ) : l , h , r , u I , x [ d , d ]
and:
ω ( F , ϵ ) = sup F l 2 , h 2 , z , x ( l 2 , h 2 ) F l 1 , h 1 , z , x ( l 1 , h 1 ) : l 1 , l 2 , h 1 , h 2 I , l 2 l 1 ϵ , h 2 h 1 ϵ , x [ d , d ] , z [ k ¯ a 2 , k ¯ a 2 ] .
Furthermore:
0 l 2 0 h 2 P ( l 2 , h 2 , r , u , x ( r , u ) ) d r d u 0 l 1 0 h 1 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u 0 l 2 0 h 2 P ( l 2 , h 2 , r , u , x ( r , u ) ) P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u + l 1 l 2 h 1 h 2 P ( l 1 , h 1 , r , u , x ( r , u ) ) d r d u a 2 ω ( P , ϵ ) + k ¯ ϵ 2 .
Therefore:
( F ^ x ) ( l 2 , h 2 ) ( F ^ x ) ( l 1 , h 1 ) 2 K x ( l 2 , h 2 ) x ( l 1 , h 1 ) + ω ( U , ϵ ) + K a 2 ω ( P , ϵ ) + k ¯ ϵ 2 + ω ( F , ϵ ) .
This gives:
ω ( F ^ x , ϵ ) 2 K ω ( x , ϵ ) + ω ( U , ϵ ) + K a 2 ω ( P , ϵ ) + k ¯ a 2 + ω ( F , ϵ ) .
Since U and F are continuous on I × I × R and I × I × R × R , respectively, therefore we get:
ω ( U , ϵ ) 0 , ω ( P , ϵ ) 0 and ω ( F , ϵ ) 0 as ϵ 0 .
Thus:
ω 0 ( F ^ X ) 2 K ω 0 ( X ) .
Similarly, we can show that:
ω 0 ( G ^ X ) 2 K ω 0 ( X ) .
From (4)–(7) and Theorem 2 (for the details of this theorem, we refer to [8]), we get that T ^ satisfies the Darbo condition on B d with respect to measure ω 0 with constant:
2 K ( α d + β ) + 2 K ( α d + β ) = 4 K ( α d + β ) = 4 K ( α d 1 + β ) = 4 K α 1 2 α β 1 4 α β 2 α 2 + β = 2 K 1 1 4 α β α < 1 .
This implies that T ^ is a contraction operator on B d with respect to ω 0 . Thus, by Theorem 2, we have that T ^ has at least one fixed point in B d . Hence, Equation (2) has at least one solution in B d C ( [ 0 , a ] × [ 0 , a ] ) . This completes the proof. □

4. An Illustrative Example

We construct the following example to illustrate the obtained result in the previous section.
Example 1.
Consider the following integral equation:
x ( l , h ) = 1 6 cos l + h 2 + 1 9 0 l 0 h r u e l h 3 + x 2 ( r , u ) d r d u 1 8 0 1 0 1 l h 6 + x ( r , u ) d r d u
for l , h [ 0 , 1 ] = I . Here, we have:
U ( l , h , x ( l , h ) ) = 1 6 cos l + h 2 ,
F ( l , h , y , x ( l , h ) ) = y 9 ,
G ( l , h , y , x ( l , h ) ) = y 8 ,
P ( l , h , r , u , x ( r , u ) ) = r u e l h 3 + x 2 ( r , u ) ,
Q ( l , h , r , u , x ( r , u ) ) = l h 6 + x ( r , u )
and a = 1 ; x , y R .
It is obvious that all the functions U , F , G , P , and Q are continuous. We have:
U ( l , h , x 1 ) U ( l , h , x 2 ( l , h ) ) = 0 . x 1 ( l , h ) x 2 ( l , h ) ,
F ( l , h , y , x 1 ( l , h ) ) F ( l , h , y , x 2 ( l , h ) ) = 0 . x 1 ( l , h ) x 2 ( l , h ) ,
G ( l , h , y , x 1 ( l , h ) ) G ( l , h , y , x 2 ( l , h ) ) = 0 . x 1 ( l , h ) x 2 ( l , h ) ,
F ( l , h , y 1 , x ( l , h ) ) F ( l , h , y 2 , x ( l , h ) ) = 1 9 y 1 y 2 ,
and:
G ( l , h , y 1 , x ( l , h ) ) G ( l , h , y 2 , x ( l , h ) ) = 1 8 y 1 y 2 .
It follows that:
A 1 ( l , h ) = A 2 ( l , h ) = A 3 ( l , h ) = 0 , A 4 ( l , h ) = 1 9 a n d A 5 ( l , h ) = 1 8 .
Consequently, we get K = 1 8 .
Furthermore,
U ( l , h , 0 ) 1 6 ,
F ( l , h , y 1 , 0 ) 1 27 e 2 1 198 . 29
and:
G ( l , h , y 2 , 0 ) 1 48 .
Thus:
M = 1 48 , L = 1 6
and:
4 α β = 9 48 < 1 .
Hence, all the assumption from (1)–(4) are satisfied. Thus, by applying Theorem 3, we conclude that Equation (8) has at least one solution in the Banach algebra C ( [ 0 , 1 ] × [ 0 , 1 ] ) .

5. An Iterative Algorithm Created by a Coupled Semi-Analytic Method to Find the Solution of the Integral Equation

To find an approximation of solution for Equation (8), we make an iterative algorithm by a coupled method created by modified homotopy perturbation and the Adomian decomposition method in the case of two-dimensional functions. Applications of the modified homotopy perturbation method to solve nonlinear integral equations, nonlinear singular integral equations, and nonlinear differential equations can be seen in [25,26,27], respectively. The Adomian decomposition method to solve physical problems was used in [28] and also to solve integro-differential equations system in [29]. However, in this article, we introduce a modified homotopy perturbation method in terms of a function with two variables, and for simplification of nonlinear terms, we use the Adomian decomposition method in the suitable form; therefore, we make an effective algorithm by the above process. Equation (8) can be shown in a general form of the two-dimensional nonlinear problem:
A x ( l , h ) f ( l , h ) = 0
with ( l , h ) I × I , where A is a general nonlinear operator and f is a known analytic function. Similar to [26,27], we divide the general operator A into two nonlinear operators as M 1 and M 2 . Of course, M 1 or M 2 can be linear operators in the special case that also f is converted to f 1 and f 2 functions; in other words, we have:
M 1 x ( l , h ) f 1 ( l , h ) + M 2 x ( l , h ) f 2 ( l , h ) = 0 .
A modified homotopy perturbation for the above problem can be introduced as follows:
H u ( l , h ) , p = M 1 u ( l , h ) f 1 ( l , h ) + p M 2 u ( l , h ) f 2 ( l , h ) = 0 , p [ 0 , 1 ] ,
where p is an embedding parameter and u is an approximation of x . According to the variations of p = 0 to p = 1 , it can be observed that M 1 u ( l , h ) = f 1 ( l , h ) to A u ( l , h ) = f ( l , h ) . This implies that for p = 1 in (9), we get the solution of (5).
We consider the above solution as the series:
x ( l , h ) u ( l , h ) = k = 0 p k u k ( l , h )
and:
x ( l , h ) = lim p 1 u ( l , h ) .
To solve Equation (8), M 1 , M 2 and f can be defined as follows:
M 1 ( x ( l , h ) ) = x ( l , h ) ,
M 2 ( x ( l , h ) ) = 1 6 cos l + h 2 + 1 9 0 l 0 h r u e l h 3 + x 2 ( r , u ) d r d u 1 8 0 1 0 1 l h 6 + x ( r , u ) d r d u ,
and:
f ( l , h ) = f 1 ( l , h ) + f 2 ( l , h ) .
Since in (8) f ( l , h ) = 0 , therefore f 1 ( l , h ) = f 2 ( l , h ) = 0 . From (9)–(14), we have:
k = 0 p k u k ( l , h ) p 1 6 cos l + h 2 + 1 9 0 l 0 h r u e l h 3 + k = 0 p k u k ( r , u ) 2 d r d u 1 8 0 1 0 1 l h 6 + k = 0 p k u k ( r , u ) d r d u = 0 .
Now, we use Adomain polynomials for simplicity for the nonlinear terms:
1 6 cos l + h 2 + 1 9 0 l 0 h r u e l h 3 + k = 0 p k u k ( r , u ) 2 d r d u = k = 0 p k A k ( l , h )
and:
1 8 0 1 0 1 l h 6 + k = 0 p k u k ( r , u ) d r d u = k = 0 p k A ^ k ( l , h ) ,
where the Adomain polynomials are given by:
A k ( l , h ) = 1 k ! d k d p k 1 6 cos l + h 2 + 1 9 0 l 0 h r u e l h 3 + k = 0 p k u k ( r , u ) 2 d r d u p = 0
and:
A ^ k ( l , h ) = 1 k ! d k d p k 1 8 0 1 0 1 l h 6 + k = 0 p k u k ( r , u ) d r d u p = 0
Therefore, we have:
k = 0 p k u k ( l , h ) p k = 0 p k A k ( l , h ) k = 0 p k A ^ k ( l , h ) = 0 .
By rearranging the terms in powers in p of (15) and using modified homotopy perturbation (9), the coefficients of p powers must be equal to zero, so we obtain an iterative algorithm (Algorithm 1) to solve for the numerical solution of Equation (8).
Algorithm 1. Algorithm of calculating u k ( l , h )
u 0 ( l , h ) = 0 ,
u 1 ( l , h ) = A 0 ( l , h ) A ^ 0 ( l , h ) ,
u k ( l , h ) = i = 0 k 1 A i ( l , h ) A ^ k 1 i ( l , h ) , k = 2 , 3 , .
Calculating the sequence u 0 ( l , h ) , u 1 ( l , h ) , , we can obtain a closed form of the solution for (8) using the above algorithm.
We compute the Adomain polynomial for k = 0 ,
A 0 ( l , h ) = 1 6 cos l + h 2 + 1 9 0 l 0 h r u e l h 3 + u 0 2 ( r , u ) d r d u = 1 6 cos l + h 2 + l 2 h 2 e l h 108
and:
A ^ 0 ( l , h ) = 1 8 0 1 0 1 l h 6 + u 0 ( r , u ) d r d u = l h 48 .
Therefore, we obtain by the algorithm:
u 1 ( l , h ) = 1 6 cos l + h 2 + l 2 h 2 e l h 108 l h 48 .
We use (10) to approximate x ( l , h ) by a few term of u k ( l , h ) as follows:
x 1 ( l , h ) = u 0 ( l , h ) + u 1 ( l , h ) = 1 6 cos l + h 2 + l 2 h 2 e l h 108 l h 48 .

6. Conclusions

In our present investigation, we have established the existence of the solution of a functional integral equation of two variables, which is of the form of the product of two operators in the Banach algebra C ( [ 0 , a ] × [ 0 , a ] ) , a > 0 and illustrated our results with the help of an example. We also constructed an iteration algorithm to get the solution of Equation (8). Further, one can solve Equation (8) using different numerical, as well as analytical methods in the setting of Banach sequence spaces and Banach algebra. Moreover, due our existence theorem for Equation (8) of two variables, we therefore conclude that our existence result is more general than the one obtained earlier by Deepmala and Pathak [9].

Author Contributions

All authors contributed equally in this work.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chandrasekhar, S. Radiative Transfer; Oxford University Press: London, UK, 1950. [Google Scholar]
  2. Hu, S.; Khavanin, M.; Zhuang, W. Integral equations arising in the kinetic theory of grass. Appl. Anal. 1989, 34, 261–266. [Google Scholar] [CrossRef]
  3. Kelly, C.T. Approximation of solution of some quadratic integral equations in transport theory. J. Integral Equ. 1982, 4, 221–237. [Google Scholar]
  4. Zabrejko, P.P.; Koshelev, A.I.; Krasnosel’skii, M.A.; Mikhlin, S.G.; Rakovshchik, L.S.; Stetsenko, V.J. Integral Equations; Nauka: Moscow, Russia, 1968. [Google Scholar]
  5. Banaś, J.; Lecko, M. Fixed points of the product of operators in Banach algebra. Panamer Math. J. 2002, 12, 101–109. [Google Scholar]
  6. Dhage, B.C. A flxed point theorem in Banach algebras with applications to functional integral equations. Kyungpook Math. J. 2004, 44, 145–155. [Google Scholar]
  7. Dhage, B.C. On a flxed point theorem in Banach algebras with applications. Appl. Math. Lett. 2005, 18, 273–280. [Google Scholar] [CrossRef]
  8. Banaś, J.; Olszowy, L. On a class of measure of noncompactness in Banach algebras and their application to nonlinear integral equations. J. Anal. Appl. 2009, 28, 1–24. [Google Scholar] [CrossRef]
  9. Pathak, H.K. A study on some problems on existence of solutions for nonlinear functional-integral equations. Acta Math. Sci. 2013, 33, 1305–1313. [Google Scholar]
  10. Kuratowski, K. Sur les espaces complets. Fund. Math. 1930, 15, 301–309. [Google Scholar] [CrossRef] [Green Version]
  11. Darbo, G. Punti uniti in trasformazioni a codominio non compatto. Rend. Sem. Mater. Univ. Padova 1955, 24, 84–92. [Google Scholar]
  12. Alotaibi, A.; Mursaleen, M.; Mohiuddine, S.A. Application of measure of noncompactness to infinite system of linear equations in sequence spaces. Bull. Iran. Math. Soc. 2015, 41, 519–527. [Google Scholar]
  13. Aghajani, A.; Allahyari, R.; Mursaleen, M. A generalization of Darbo’s theorem with application to the solvability of systems of integral equations. J. Comput. Appl. Math. 2014, 260, 68–77. [Google Scholar] [CrossRef]
  14. Arab, R.; Allahyari, R.; Haghighi, A.S. Existence of solutions of infinite systems of integral equations in two variables via measure of noncompactness. Appl. Math. Comput. 2014, 246, 283–291. [Google Scholar] [CrossRef]
  15. Darwish, M.A. On quadratic integral equation of fractional orders. J. Math. Anal. Appl. 2005, 311, 112–119. [Google Scholar] [CrossRef] [Green Version]
  16. Mishra, L.N.; Sen, M.; Mohapatra, R.N. On existence theorems for some generalized nonlinear functional-integral equations with applications. Filomat 2017, 31, 2081–2091. [Google Scholar]
  17. Mohiuddine, S.A.; Srivastava, H.M.; Alotaibi, A. Application of measures of noncompactness to the infinite system of second-order differential equations in lp spaces. Adv. Differ. Equ. 2016, 2016, 317. [Google Scholar] [CrossRef]
  18. Mursaleen, M.; Mohiuddine, S.A. Applications of measures of noncompactness to the infinite system of differential equations in lp spaces. Nonlinear Anal. 2012, 75, 2111–2115. [Google Scholar] [CrossRef]
  19. Mursaleen, M.; Rizvi, S.M.H. Solvability of infinite systems of second order differential equations in c0 and 1 by Meir-Keeler condensing operators. Proc. Am. Math. Soc. 2016, 144, 4279–4289. [Google Scholar] [CrossRef]
  20. Srivastava, H.M.; Das, A.; Hazarika, B.; Mohiuddine, S.A. Existence of solutions of infinite systems of differential equations of general order with boundary conditions in the spaces c0 and l1 via the measure of noncompactness. Math. Meth. Appl. Sci. 2018, 41, 3558–3569. [Google Scholar] [CrossRef]
  21. Banaś, J.; Goebel, K. Measure of Noncompactness in Banach Spaces, Lecture Notes in Pure and Applied Mathematics; Marcel Dekker: New York, NY, USA, 1980; Volume 60. [Google Scholar]
  22. Shang, Y. Subspace confinement for switched linear systems. Forum Math. 2017, 29, 693–699. [Google Scholar] [CrossRef]
  23. Shang, Y. Fixed-time group consensus for multi-agent systems with non-linear dynamics and uncertainties. IET Control Theory Appl. 2018, 12, 395–404. [Google Scholar] [CrossRef]
  24. Shang, Y. Resilient consensus of switched multi-agent systems. Syst. Control Lett. 2018, 122, 12–18. [Google Scholar] [CrossRef]
  25. Glayeri, A.; Rabbani, M. New Technique in Semi-Analytic Method for Solving Non-Linear Differential Equations. Math. Sci. 2011, 5, 395–404. [Google Scholar]
  26. Rabbani, M. Modified homotopy method to solve non-linear integral equations. Int. J. Nonlinear Anal. Appl. 2015, 6, 133–136. [Google Scholar]
  27. Rabbani, M. New Homotopy Perturbation Method to Solve Non-Linear Problems. J. Math. Comput. Sci. 2013, 7, 272–275. [Google Scholar] [CrossRef]
  28. Adomian, G. Solving Frontier Problem of Physics: The Decomposition Method; Kluwer Academic Press: Dordrecht, The Netherlands, 1994. [Google Scholar]
  29. Rabbani, M.; Arab, R.; Hazarika, B. Solvability of nonlinear quadratic integral equation by using simulation type condensing operator and measure of noncompactness. Appl. Math. Comput. 2019, 349, 102–117. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Srivastava, H.M.; Das, A.; Hazarika, B.; Mohiuddine, S.A. Existence of Solution for Non-Linear Functional Integral Equations of Two Variables in Banach Algebra. Symmetry 2019, 11, 674. https://doi.org/10.3390/sym11050674

AMA Style

Srivastava HM, Das A, Hazarika B, Mohiuddine SA. Existence of Solution for Non-Linear Functional Integral Equations of Two Variables in Banach Algebra. Symmetry. 2019; 11(5):674. https://doi.org/10.3390/sym11050674

Chicago/Turabian Style

Srivastava, Hari M., Anupam Das, Bipan Hazarika, and S. A. Mohiuddine. 2019. "Existence of Solution for Non-Linear Functional Integral Equations of Two Variables in Banach Algebra" Symmetry 11, no. 5: 674. https://doi.org/10.3390/sym11050674

APA Style

Srivastava, H. M., Das, A., Hazarika, B., & Mohiuddine, S. A. (2019). Existence of Solution for Non-Linear Functional Integral Equations of Two Variables in Banach Algebra. Symmetry, 11(5), 674. https://doi.org/10.3390/sym11050674

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