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

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, ∫ l 0 ∫ h 0 P(l, h, r, u, x(r, u))drdu, x(l, h) )) × G ( l, h, ∫ a 0 ∫ a 0 Q (l, h, r, u, x(r, u)) drdu, 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.


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 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., 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: 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).

Preliminaries
Let R denote the set of real numbers, and write R + = [0, ∞) .Suppose Ē is a real Banach space with the norm ., and let X( = φ) ⊆ Ē.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 Ē centred at a and with radius b is denoted by B(a, b).In addition, we use the symbol M Ē to denote the family of all non-empty and bounded subsets of Ē and use N Ē 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: 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 xy ∈ A is defined with the properties: ) α(xy) = (α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: xy ≤ 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: ) is said to be a measure of noncompactness in Ē if: (i) for all X ∈ M Ē, we have that µ(X) = 0 implies that X is precompact.
where convX is the convex closure of set X.
(vi) µ (λX 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 ∞ ∈ kerµ.
For given subsets X, Y in a Banach algebra E, the product XY defined by: XY = {xy : 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: 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 µ(TY) ≤ 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 Ē.Let S : Z → Z be a continuous mapping.Suppose that there is a constant k ∈ [0, 1) such that: 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: 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 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 xy ∈ E, where: Since t, s are arbitrary, therefore (xy)z = x(yz).
Similarly, it can be shown that: x(y + z) = xy + xz, (x + y)z = xz + yz and: Therefore E is an algebra.
The space E is also a normed space with the norm: For all x, y ∈ E and l, h ∈ [0, a], and so: sup 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: Therefore, for all t, s ∈ [0, a], we get: 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]: which yields: sup Thus: 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: Further, we define: 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).

Main Result
In this section, we study the existence of solutions of the integral Equation (2).We consider the following assumptions: (1) The functions where l, h ∈ [0, a] and x, x 1 , x 2 , y, y 1 , y 2 ∈ R. Furthermore, let: where K ≥ 0. Proof.Let us consider the operators F and Ĝ defined on E by: ( Fx)(l, h) = U(l, h, x(l, h)) + F l, h, l 0 h 0 P(l, h, r, u, x(r, u))drdu, x(l, h) and: From Assumptions (1)-( 3), we get that F and Ĝ map C(I × I) into itself.Furthermore, let us define another operator T on C(I × I) as follows: It is obvious that T maps C(I × I) into itself. Let: and: Let x ∈ C(I × I) be fixed and l, h ∈ I.We get: Let α = 2k, β = L + M.Then, we have: for x ∈ C(I × I).From (3), we have that the operator T maps B d ⊂ C(I × I) into B d , where: , where: and: Furthermore, we have: Let > 0 be fixed and x(l, h), y(l, h) ∈ B d such that: x − y ≤ , (l, h ∈ I).
Since P is continuous, so it is uniformly continuous on the compact set Thus, F is continuous on B d .Similarly, one can prove that Ĝ 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 , and: Furthermore: (P(l 2 , h 2 , r, u, x(r, u)) − P(l 1 , h Therefore: This gives: 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: Similarly, we can show that: 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: 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 This completes the proof.

An Illustrative Example
We construct the following example to illustrate the obtained result in the previous section.

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: A modified homotopy perturbation for the above problem can be introduced as follows: 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: 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: and: Since in ( 8) f (l, h) = 0, therefore f 1 (l, h) = f 2 (l, h) = 0. From ( 9)-( 14), we have: Now, we use Adomain polynomials for simplicity for the nonlinear terms: where the Adomain polynomials are given by: 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).

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].
are continuous, and there exist nonnegative constants L, M such that: