Next Article in Journal
Design, Hardware Implementation on FPGA and Performance Analysis of Three Chaos-Based Stream Ciphers
Previous Article in Journal
Certain Sharp Coefficient Results on a Subclass of Starlike Functions Defined by the Quotient of Analytic Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Computational Technique for Solving Three-Dimensional Mixed Volterra–Fredholm Integral Equations

1
Department of Mathematics, Faculty of Science, Zagazig University, Zagazig P.O. Box 44519, Egypt
2
Department of Mathematics & Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(2), 196; https://doi.org/10.3390/fractalfract7020196
Submission received: 10 January 2023 / Revised: 1 February 2023 / Accepted: 10 February 2023 / Published: 16 February 2023

Abstract

:
In this article, a novel and efficient approach based on Lucas polynomials is introduced for solving three-dimensional mixed Volterra–Fredholm integral equations for the two types (3D-MVFIEK2). This method transforms the 3D-MVFIEK2 into a system of linear algebraic equations. The error evaluation for the suggested scheme is discussed. This technique is implemented in four examples to illustrate the efficiency and fulfillment of the approach. Examples of numerical solutions to both linear and nonlinear integral equations were used. The Lucas polynomial method and other approaches were contrasted. A collection of tables and figures is used to present the numerical results. We observe that the exact solution differs from the numerical solution if the exact solution is an exponential or trigonometric function, while the numerical solution is the same when the exact solution is a polynomial. The Maple 18 program produced all of the results.

1. Introduction

The multi-dimensional integral equations (MDIE) provide considerable flexibility in order to implement a wide range of problems and relationships as well as resolve boundary issues for differential equations (DE). In three-dimensional electromagnetic modeling, Hursan and Zhdanov [1] pioneered the Contraction integral equation method. Pachpatte presented numerous MDIE applications in [2]. Cheng investigated quantum mechanics with thermal radiation effect in [3]. Chew et al used integral equation methods for electromagnetic and Elastic Waves in [4]. Integral equations with two or more variables are challenging to resolve analytically. Thus, in such instances, it is essential to utilize numerical procedures to approach the solution.
The numerical solution of MDIE has been the subject of research. To illustrate, Guoqiang et al in [5], calculated the numeric solution of two-dimensional nonlinear Volterra integral equations (VIE) via employment and repeated employment techniques. Bazm, in [6], used Bernoulli polynomials to obtain the solution of a class of nonlinear two-dimensional integral equations. The numeric solution of two-dimensional nonlinear Fredholm integral equations (FIE), in the most general type of kernels, is based on Bernstien polynomials introduced by Maleknejad et al. in [7]. In [8], a three-dimensional differential converts technique is implemented for overcoming nonlinear three-dimensional VIE. Basseem, in [9], implemented the degenerate kernel technique for solving three-dimensional nonlinear integral equations of the two types. Kazemi et al., in [10], applied the Haar wavelet procedure for solving the three-dimensional nonlinear FIE. Using the 3D-block-pulse functions developed by Manafian and Bolghar in [11], the numerical solution of nonlinear three-dimensional Volterra integro-differential equations is presented.
Many works have been concentrated on the expansion of a better, higher, and more powerful technique for solving Volterra–Fredholm integral equations (VFIE), such as Maleknejad et al. in [12], applying the Bernstein’s approximation for solving the three-dimensional VFIE of the initial and two types. Mirzaee and Hadadiyan, in [13], applied a computing technique for overcoming the nonlinear composite VFIE. Mirzaee et al., in [14], presented triangular functions for solving the nonlinear mixed VFIE. Moreover, Mirzaee et al., in [15], implemented the amendment block-pulse functions to clarify the three-dimensional VFIE.
In this work, the approximate solution of 3D-MVFIEK2 have been obtained using three-dimensional Lucas polynomials. These polynomials are non-orthogonal. This polynomial was used by several scientists to solve fractional differential equations [16,17]. Cetin et al., in [18], utilized the Lucas polynomial procedure to research a system of higher-order DE. Baykus and Sezer, in [19], implemented a hybrid Taylor–Lucas employment procedure for the delay DE. Nadir, in [20], procured a solution of the differential-integro equation utilization of the Lucas series. Haq and Ali, in [21], analyze the approach to solving two-dimensional Sobolev equality utilization mixed Lucas and Fibonacci polynomials. Nisar et al., in [22], applied Lucas polynomials for solving 1D and 2D advection–diffusion–reaction equations. The properties of Lucas polynomials are presented in [23]. In [24], Chelyshkov polynomials approach the first class of two-dimensional nonlinear VIE. In [25], Lucas polynomials are used to approximate the solution of Cauchy integral equations. Modified Lucas polynomials were used by Youssri et al. [26] to numerically solve boundary value issues. The primary benefit of the current approach is that the issue at hand is converted into a set of algebraic equations that can be quickly resolved using computer programming.
The paper is set up as shows: Section 2 presents various Lucas polynomial properties. Lucas polynomials are used to solve 3D-MVFIEK2 in Section 3. Error analysis is discussed in Section 4. Numerical outcomes and comparisons with other techniques are presented in Section 5. The conclusion of the paper is presented in Section 6.
In this article, we shall contemplate the numerical solution of a size of 3D-MVFIEK2 in the next section, to create:
w ( x , y , z ) = g ( x , y , z ) + λ 0 1 0 1 0 1 K ( x , y , z , s , t , r ) w ( s , t , r ) d s d t d r +
+ λ 0 1 0 1 0 1 k 2 ( x , y , z , s , t , r ) w ( s , t , r ) d s d t d r ,
where
K ( x , y , z , s , t , r ) = k 1 ( x , y , z , s , t , r ) ( s , t , r ) ( x , y , z ) , 0 , ( s , t , r ) > ( x , y , z ) ,
x , y , z [ 0 , 1 ) , λ R , and w is the unknown function, g , k 1 and k 2 have analytic functions on D and D × D severally, where D = [ 0 , 1 ) × [ 0 , 1 ) × [ 0 , 1 ) .
Equation (1a) could be expressed as follows:
w ( x , y , z ) = g ( x , y , z ) + λ 0 1 0 1 0 1 H ( x , y , z , s , t , r ) w ( s , t , r ) d s d t d r ,
where
H ( x , y , z , s , t , r ) = K ( x , y , z , s , t , r ) + k 2 ( x , y , z , s , t , r ) .

2. Properties of Lucas Polynomials

Definition 1.
([23]) The Lucas polynomials have specified via the recurrence relationship
L m + 2 ( ζ ) = ζ L m + 1 ( ζ ) + L m ( ζ ) , m 0 ,
with L 0 ( ζ ) = 2 and L 1 ( ζ ) = ζ .
These polynomials possess the following expression, to establish:
L m ( ζ ) = ( ζ + ζ 2 + 4 ) m + ( ζ ζ 2 + 4 ) m 2 m ,
furthermore, they possess the subsequent explicit power form ofrepresentation
L m ( ζ ) = k = 0 m 2 n n k m k k ζ m 2 k , m 1 ,
where
m 2 = m 2 , m e v e n , m 1 2 , m o d d .
The initially smaller Lucas polynomials have expressed lower, and the collection of these factors is in Table 1.
Definition 2.
(Rodrigues’s formula).The Lucas polynomials L m ( ζ ) shall be obtained over the Rodrigue polynomials, andthe formula is specified via
L m ( ζ ) = 2 m ! ( 2 m ) ! ( ζ 2 + 4 ) 1 / 2 d m d ζ m { ( ζ 2 + 4 ) m 1 2 } .
Proposition 1.
The producing function for Lucas polynomials is determined as [23]
G ( ζ , t ) = m = 0 L m ( ζ ) t m = 1 + t 2 1 t 2 ζ t = 1 + ζ t + ( ζ 2 + 2 ) t 2 + ( ζ 3 + 3 ζ ) t 3 + .
Approximation function
Let w ( x ) be square integrable function on ( 0 , 1 ) , and we assume it shall be displayed in terms of the Lucas polynomials, as shows:
w ( x ) w N ( x ) = i = 0 M c i L i ( x ) ,
where L i ( x ) is the Lucas polynomials and c i , i = 0 , 1 , , M are the unknown coefficients. The Lucas series (2e) can be expressed as a matrix, as shows:
w ( x ) w M ( x ) = L ( x ) C ,
where L ( x ) = [ L 0 ( x ) L 1 ( x ) L M ( x ) ] is the ( M + 1 ) × ( M + 1 ) matrix of Lucas polynomials, C = [ c 0 c 1 c M ] T is the ( M + 1 ) × 1 unknown coefficients in a vector.
Then, by utilization of the Lucas polynomials L n ( x ) , specified via (2c), we write the matrix form L ( x ) , as follows:
L ( x ) = χ ( x ) Q T ,
where χ ( x ) = [ 1 x x 2 x M ] ,
If M is odd
Q = 2 0 0 0 0 1 1 1 0 0 0 2 1 1 1 0 2 2 2 0 0 0 3 2 2 1 0 0 ( m 1 ) ( m 1 ) / 2 ( m 1 ) / 2 ( m 1 ) / 2 0 ( m 1 ) ( m + 1 ) / 2 ( m + 1 ) / 2 ( m 3 ) / 2 0 0 ( m ) ( m + 1 ) / 2 ( m + 1 ) / 2 ( m 1 ) / 2 0 m m m 0 ,
if M even
Q = 2 0 0 0 0 1 1 1 0 0 0 2 1 1 1 0 2 2 2 0 0 0 3 2 2 1 0 0 0 ( m 1 ) m / 2 m / 2 ( m 2 ) / 2 0 0 m m / 2 m / 2 m / 2 0 m ( m + 2 ) / 2 ( m + 2 ) / 2 ( m 2 ) / 2 m m m 0 .
By the matrix (2f) and (2g), it follows that
w M ( x ) = χ ( x ) Q T C .
The function of two-dimensional continuous w ( x , y ) shall be expressed using Lucas polynomials in the following manner:
w ( x , y ) w M ( x , y ) = k = 0 M m = 0 M c k m L k ( x ) L m ( y ) = k = 0 M m = 0 M c k m L k m ( x , y ) = C L ( x , y ) ,
where C = [ c 00 c 01 c 0 M c M 0 c M M ] T is the order vector with an unknown coefficient for Lucas ( M + 1 ) 2 × 1 and L ( x , y ) = [ L 00 ( x , y ) L 0 M ( x , y ) L M 0 ( x , y ) L M M ( x , y ) ] is the square matrix of order ( M + 1 ) 2 .
Assume that the following Lucas polynomials are to be used to define a three-dimensional continuous function w ( x , y , z ) :
w ( x , y , z ) w M ( x , y , z ) = k = 0 M m = 0 M n = 0 M c k m n L k ( x ) L m ( y ) L m ( z )
= k = 0 M m = 0 M m = 0 M c k m n L k m n ( x , y , z ) = C ¯ L ( x , y , z ) ,
where C ¯ = [ c 000 c 00 M c 0 M M c M M M ] T is the order vector with an unknown coefficient for Lucas ( M + 1 ) 3 × 1 and L ( x , y , z ) = [ L 000 ( x , y , z ) L 00 M ( x , y , z ) L 0 M M ( x , y , z ) L M M M ( x , y , z ) ] is the square matrix of the order ( M + 1 ) 3 .
Let g ( x , y , z ) L 2 ( D ) be any function; we must calculate it using a truncated Lucas series:
g ( x , y , z ) k = 0 M m = 0 M m = 0 M g k m n L k ( x ) L m ( y ) L m ( z ) = k = 0 M m = 0 M m = 0 M g k m n L k m n ( x , y , z ) = L T ( x , y , z ) G ,
where G is an ( M + 1 ) 3 × 1 vector given by G = [ g 000 g 00 M g 0 M M g M M M ] T .
Assume k ( x , y , z , s , t , r ) L 2 ( D × D ) . It shall correspondingly expanded with respect to Lucas polynomials from:
k ( x , y , z , s , t , r ) L T ( x , y , z ) K L ( s , t , r ) ,
where K is the ( M + 1 ) 3 × ( M + 1 ) 3 Lucas coefficients.

3. Technique of Solution

In this part, we use the Lucas polynomials technique to solve Equation (1c), which corresponds to Equation (1a).
Approximating the functions w , g , and H with respect to three-dimensional Lucas polynomials, as specified in part 2, we find
w ( x , y , z ) L T ( x , y , z ) C ¯ ,
g ( x , y , z ) L T ( x , y , z ) G ,
H ( x , y , z , s , t , r ) L T ( x , y , z ) H L ( s , t , r ) .
Substituting from (3a) into (1c), we procure
L T ( x , y , z ) C ¯ L T ( x , y , z ) G + λ 0 1 0 1 0 1 L T ( x , y , z ) H L ( s , t , r ) L T ( s , t , r ) C ¯ d s d t d r .
Thus,
L T ( x , y , z ) C ¯ = L T ( x , y , z ) G + λ L T ( x , y , z ) H 0 1 0 1 0 1 L ( s , t , r ) L T ( s , t , r ) C ¯ d s d t d r .
By using
L ( x , y , z ) L T ( x , y , z ) C ¯ = C ¯ ¯ L ( x , y , z ) ,
where C ¯ ¯ = d i a g ( C ¯ ) is a ( M + 1 ) 3 × ( M + 1 ) 3 matrix in a diagonal.
It is clear that
0 1 0 1 0 1 L ( s , t , r ) d s d t d r = B 1 .
Substituting from (3c) and (3d) into (3b), we procure
L T ( x , y , z ) C ¯ = L T ( x , y , z ) G + λ L T ( x , y , z ) H C ¯ ¯ B 1 ,
therefore
L T ( x , y , z ) C ¯ = L T ( x , y , z ) G + λ L T ( x , y , z ) Y ,
where Y = H C ¯ ¯ B 1 .
By inserting the collocation point
( x i , y j , z k ) = ( i M , j M , k M ) , i , j , k = 0 , 1 , , M ,
we obtain
L T ( x i , y j , z k ) C ¯ = L T ( x i , y j , z k ) G + λ L T ( x i , y j , z k ) Y .
By simplifying (3g), we obtain
C ¯ = G + λ Y .
Next, when resolving the top linear system, we may discover the factor and we procure the estimated solution.

4. Error Analysis

To generalize Lucas polynomials of the fractional DE, Abd-Elhameed and Youssri researched the convergence and error analysis in their paper from [16]. This section covers the equation’s error analysis (1c).
Assume that ( C ( D ) , . ) has the Banach space of every continuous function on D with a norm of
w ( x , y , z ) = max ( x , y , z ) D w ( x , y , z ) .
Furthermore, we denote the error by
e M ( x , y , z ) = w M ( x , y , z ) w ( x , y , z ) ,
where w and w M are the analytical and approximate solutions of the three-dimensional integral Equation (1c) severally.
Theorem 1.
Suppose that g ( x , y , z ) in Equation (1c) is bounded ( x , y , z ) D and H ( x , y , z , s , t , r ) M 1 ( x , y , z , s , t , r ) D × D . Subsequently, Equation (1c) has a unique solution whenever 0 < α < 1 , λ M 1 = α .
Proof. 
Let
e M ( x , y , z ) = w M ( x , y , z ) w ( x , y , z ) = max ( x , y , z ) D w M ( x , y , z ) w ( x , y , z )
= max ( x , y , z ) D λ 0 1 0 1 0 1 H ( x , y , z , s , t , r ) w M ( s , t , r ) d s d t d r λ 0 1 0 1 0 1 H ( x , y , z , s , t , r ) w ( s , t , r ) d s d t d r
= max ( x , y , z ) D λ 0 1 0 1 0 1 H ( x , y , z , s , t , r ) ( w M ( s , t , r ) w ( s , t , r ) d s d t d r
w M ( x , y , z ) w ( x , y , z ) max ( x , y , z ) D λ 0 1 0 1 0 1 H ( x , y , z , s , t , r ) d s d t d r .
Subsequently,
e M ( x , y , z ) e M ( x , y , z ) ( λ M 1 ) α e M ( x , y , z ) ,
therefore
( 1 α ) e M ( x , y , z ) 0 .
Thus, if 0 < α < 1 , we have e M ( x , y , z ) 0 as M , and this completes the proof of the theorem.
For an error estimation of the approximate solution of Equation (1c), we consider
r M ( x , y , z ) + g ( x , y , z ) = w M ( x , y , z ) λ 0 1 0 1 0 1 H ( x , y , z , s , t , r ) w M ( x , y , z ) d s d t d r ,
where r M is the perturbation function.
By deducting Equation (4c) from Equation (1c) we get
r M ( x , y , z ) e M ( x , y , z ) + λ M 1 e M ( x , y , z ) .
Therefore,
r M ( x , y , z ) ( 1 + λ M 1 ) e M ( x , y , z ) ,
then r M ( x , y , z ) is bounded. □

5. Numerical Examples

In this part, four numerical illustrations are introduced to explain the suggested technique. All computations have been formed by using Maple 18.
Example 1.
Let the next three-dimensional FIE of the second type be:
w ( x , y , z ) = g ( x , y , z ) + 0.01 0 1 0 1 0 1 ( e ( x y z ) s 2 t r 2 ) w p ( s , t , r ) d s d t d r , ( x , y , z ) D , p = 1 , 2 ,
with the analytical solution w ( x , y , z ) = x 2 y 2 z 2 .
Case (I) if p = 1 , we obtain the following linear integral equation (LIE):
w ( x , y , z ) = x 2 y 2 z 2 1 10000 e ( x y z ) + 0.01 0 1 0 1 0 1 ( e ( x y z ) s 2 t r 2 ) w ( s , t , r ) d s d t d r , ( x , y , z ) D .
Case (II) if p = 2 , we obtain the following nonlinear integral equation (NIE):
w ( x , y , z ) = x 2 y 2 z 2 1 29400 e ( x y z ) + 0.01 0 1 0 1 0 1 ( e ( x y z ) s 2 t r 2 ) w 2 ( s , t , r ) d s d t d r , ( x , y , z ) D .
Implementation Lucas polynomial technique for Equation (5b) with M = 2 and via utilization of the employment points (3f). We receive a set of 27 linear algebraic equations with an equivalent number of unknowns. We reverse this scheme and obtain the sum calculated, as follows:
c 0 , 0 , 0 = 1 , c 0 , 0 , 1 = 0 , c 0 , 0 , 2 = 1 , c 0 , 1 , 0 = 0 , c 0 , 1 , 1 = 0 , c 0 , 1 , 2 = 0 , c 0 , 2 , 0 = 1 ,
c 0 , 2 , 1 = 0 , c 0 , 2 , 2 = 1 , c 1 , 0 , 0 = 0 , c 1 , 0 , 1 = 0 , c 1 , 0 , 2 = 0 , c 1 , 1 , 0 = 0 , c 1 , 1 , 1 = 0 ,
c 1 , 1 , 2 = 0 , c 1 , 2 , 0 = 0 , c 1 , 2 , 1 = 0 , c 1 , 2 , 2 = 0 , c 2 , 0 , 0 = 1 , c 2 , 0 , 1 = 0 , c 2 , 0 , 2 = 1 ,
c 2 , 1 , 0 = 0 , c 2 , 1 , 1 = 0 , c 2 , 1 , 2 = 0 , c 2 , 2 , 0 = 1 , a 2 , 2 , 1 = 0 , c 2 , 2 , 2 = 1
and, we receive the approximate resolution that is the same as the accurate resolution.
Similarly, applying the Lucas polynomial method for Equation (5c), we obtain the approximate solution in the form of
w 2 ( x , y , x ) = x 2 y 2 z 2 1.011599486 × 10 10 x y z 3.102400059 × 10 10 x 2 y z 2 3.102400059 × 10 10 x 2 y 2 z
+ 1.917200016 × 10 10 x y z 2 + 1.917200016 × 10 10 x y 2 z 3.102400059 × 10 10 x y 2 z 2 + 1.917200016 × 10 10 x 2 y z .
This is a comparison between the exact solution, approximate solution and absolute error presented in Table 2. Figure 1 represented absolute error of example 1 when z = 0.5 .
Example 2.
Let the next three-dimensional IE of the second type be [9,10]:
w ( x , y , z ) = g ( x , y , z ) + 0.01 0 1 0 1 0 1 ( x 2 s r z t 2 + x 2 s y r z ) w p ( s , t , r ) d s d t d r , ( x , y , z ) D , p = 1 , 2 ,
with the exact solution w ( x , y , z ) = x 2 y 2 z .
Comparing Equation (5d) with Equation (1c), we note that
H ( x , y , z , s , t , r ) = k 1 ( x , y , z , s , t , r ) + k 2 ( x , y , z , s , t , r ) = x 2 s r z t 2 + x 2 s y r z ,
k 1 ( x , y , z , s , t , r ) = x 2 s r z t 2 ( s , t , r ) ( x , y , z ) , 0 , ( s , t , r ) > ( x , y , z ) ,
k 2 ( x , y , z , s , t , r ) = x 2 s y r z .
Case (I) if p = 1 , we obtain the following LIE
w ( x , y , z ) = x 2 y 2 z 1 6000 x 2 z 1 3600 x 2 y z + 0.01 0 1 0 1 0 1 ( x 2 s r z t 2 + x 2 s y r z ) w ( s , t , r ) d s d t d r .
Case (II) if p = 2 , we obtain the following NIE
w ( x , y , z ) = x 2 y 2 z 1 16800 x 2 z 1 12000 x 2 y z + 0.01 0 1 0 1 0 1 ( x 2 s r z t 2 + x 2 s y r z ) w 2 ( s , t , r ) d s d t d r .
Similarly, as in example 1, and applying the Lucas polynomial method for Equation (5e) with M = 2 via utilization, the employment points are (3f). We receive a set of 27 linear algebraic equations with an equivalent number of unknowns. We reverse this scheme and obtain the sum calculated, as follows:
c 0 , 0 , 0 = 0 , c 0 , 0 , 1 = 1 , c 0 , 0 , 2 = 0 , c 0 , 1 , 0 = 0 , c 0 , 1 , 1 = 0 , c 0 , 1 , 2 = 0 , c 0 , 2 , 0 = 0 ,
c 0 , 2 , 1 = 1 , c 0 , 2 , 2 = 0 , c 1 , 0 , 0 = 0 , c 1 , 0 , 1 = 0 , c 1 , 0 , 2 = 0 , c 1 , 1 , 0 = 0 , c 1 , 1 , 1 = 0 ,
c 1 , 1 , 2 = 0 , c 1 , 2 , 0 = 0 , c 1 , 2 , 1 = 0 , c 1 , 2 , 2 = 0 , c 2 , 0 , 0 = 0 , c 2 , 0 , 1 = 1 , c 2 , 0 , 2 = 0 ,
c 2 , 1 , 0 = 0 , c 2 , 1 , 1 = 0 , c 2 , 1 , 2 = 0 , c 2 , 2 , 0 = 0 , a 2 , 2 , 1 = 1 , c 2 , 2 , 2 = 0 ,
and we receive the approximate resolution, which is the same as the accurate resolution.
Similarly, when applying the Lucas polynomial method for Equation (5f), we obtain the approximate solution in the form
w 2 ( x , y , x ) = 1.62042 × 10 9 x y z + 1.80048 × 10 9 x y z 2 + 2.52052 × 10 9 x y 2 z 2.80064 × 10 9 x y 2 z 2
+ 2.482404895 × 10 9 x 2 y z 2.80056 × 10 9 x 2 y z 2 + 4.0008 × 10 9 x 2 y 2 z 2 + 0.9999999965 x 2 y 2 z .
Comparisons between the absolute errors of the presented technique and the Haar wavelet technique [10] are shown in Table 3. Figure 2 represented the absolute error of example 2 when z = 0.5 .
Example 3.
Let the next three-dimensional FIE of the second type be:
w ( x , y , z ) = g ( x , y , z ) + 0.01 0 1 0 1 0 1 ( e r x 2 y 2 s t ) w p ( s , t , r ) d s d t d r , ( x , y , z ) D , p = 1 , 2 ,
with analytical solution w ( x , y , z ) = x y e z .
Case (I) if p = 1 we obtain the following LIE:
w ( x , y , z ) = x y e z ( 1 1800 y 2 x 2 1 1800 e 2 x 2 y 2 ) + 0.01 0 1 0 1 0 1 ( e r x 2 y 2 s t ) w ( s , t , r ) d s d t d r .
Case (II) if p = 2 , we obtain the following NIE:
w ( x , y , z ) = x y e z ( 1 4800 x 2 y 2 1 4800 x 2 y 2 e 3 ) + 0.01 0 1 0 1 0 1 ( e r x 2 y 2 s t ) w 2 ( s , t , r ) d s d t d r .
Similarly as examples 1 and 2, we obtain the approximate solutions of Equations (5h) and (5i), respectively, in the forms
w 2 ( x , y , z ) = 1.3 × 10 16 x + 1.2 × 10 16 y 0.9417567917 x y z + 0.3096362366 x y z 2 9.1 × 10 9 x y 2 z
+ 5.290847763 × 10 9 x y 2 z 2 1.08 × 10 8 x 2 y z + 6.616953330 × 10 9 x 2 y z 2 + 8 × 10 9 x 2 y 2 z 3.7500021 × 10 9 x 2 y 2 z 2
+ 0.9999999962 x y + 4.281458270 × 10 9 x y 2 + 4.571508262 × 10 9 x 2 y + 7.348829758 × 10 7 x 2 y 2 + 4.0 × 10 16 y 2 ,
and
w 2 ( x , y , z ) = 1.28 × 10 9 x y 2 z 2 + 1.04 × 10 9 x 2 y z 1.28 × 10 9 x 2 y z 2 + 9.6 × 10 10 x 2 y 2 z 3.2 × 10 10 x 2 y 2 z 2
0.9417568054 x y z + 0.3096362465 x y z 2 + 1.040000000 × 10 9 x y 2 z + 1.000000001 x y 3.96292120 × 10 10 x y 2
3.96292120 × 10 10 x 2 y + 9.724822738 × 10 7 x 2 y 2 .
Absolute errors of example 3 is shown in Table 4 and Figure 3 and Figure 4.
Example 4.
Let the next three-dimensional Fredholm integral equation of the second type be:
w ( x , y , z ) = x y z e x 2 y 2 z 2 + 1 16 ( 1 e 1 ) 3 x y z + 1 2 0 1 0 1 0 1 x y z w ( s , t , r ) d s d t d r ,
with exact solution w ( x , y , z ) = x y z e x 2 y 2 z 2 .
The numeric outcomes of example 4 are displayed in Table 5 according to the suggested procedure, and Figure 5 represented the Absolute error when z = 0.5 .

6. Conclusions

A numeric technique depending on the Lucas polynomial is proposed in this document for the approximate solution of 3D-MVFIEK2. As the numeric outcome has demonstrated, by comparing the present method with other numerical techniques, it is clear that the results of the Lucas polynomial method are better than the results of the Haar wavelet technique [10], and the CPT time used for solving the examples is very small when using the Maple program.

Future Work

We can develop this method for solving the singular 3D-VFIEK2 with some modifications. Moreover, we can use the finite difference method for solving different kinds of integral equations.

Author Contributions

All authors (A.M.S.M., A.S.N., K.M.H., D.S.M.) contributed to the analysis, formulated the approach, implemented the study, and wrote the document. The final document is acknowledged by all contributors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

In this novel copy, no instructional reports are produced or reach easily with the running assessment information.

Acknowledgments

The researchers would like to sincerely thank M. A. Abdou for his insightful counsel and scholarly recommendations during the course of the study.

Conflicts of Interest

Regarding the original copy, beginning, or appearance of this document, there are no conflict of interest.

References

  1. Hursan, G.; Zhdanov, M.S. Contraction integral equation method in three-dimensional electromagnetic modeling. Radio Sci. 2002, 6, 1–13. [Google Scholar] [CrossRef] [Green Version]
  2. Pachpatte, B.G. Multidimensional Integral Equations and Inequalities; Springer Science, Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  3. Cheng, Z. Quantum effects of thermal radiation in a Kerr nonlinear blackbody. J. Optical Soc. Amer. B 2002, 19, 1692–1705. [Google Scholar] [CrossRef]
  4. Chew, W.C.; Tong, M.S.; Hu, B. Integral Equation Methods for Electromagnetic and Elastic Waves. Synth. Lect. Comput. Elect. 2008, 3, 1–241. [Google Scholar]
  5. Guoqiang, H.; Hayami, K.; Sugihara, K.; Jiong, W. Extrapolation method of iterated collocation solution for two-dimensional nonlinear Volterra integral equations. App. Math. Comput. 2009, 112, 70–76. [Google Scholar] [CrossRef]
  6. Bazm, S. Numerical solution of a class of nonlinear two-dimensional integral equations using bernoulli polynomials. Sahand Commun. Math. Anal. (SCMA) 2016, 3, 37–51. [Google Scholar]
  7. Babaaghaie, A.; Maleknejad, K. A new approach for numerical solution of two-dimensional nonlinear Fredholm integral equations in the most general kind of kernel, based on Bernstein polynomials and its convergence analysis. J. Comput. Appl. Math. 2018, 344, 482–494. [Google Scholar] [CrossRef]
  8. Ziqan, A.; Armiti, S.; Suwan, I. Solving three-dimensional Volterra integral equation by the reduced differential transform method. Int. J. Of Applied Math. Res. 2016, 5, 103–106. [Google Scholar] [CrossRef] [Green Version]
  9. Basseem, M. Degenerate Kernel Method for Three Dimension Nonlinear Integral Equations of the Second Kind. Univers. J. Integral Equations 2015, 3, 61–66. [Google Scholar]
  10. Kazemia, M.; Torkashvand, V.; Ezzati, R. A new method for solving three-dimensional nonlinear Fredholm integral equations by Haar wavelet. Int. J. Nonlinear Anal. Appl. 2021, 12, 115–133. [Google Scholar]
  11. Manafian, J.; Bolghar, P. Numerical solutions of nonlinear 3-dimensional Volterra integral-differential equations with 3D-block-pulse functions. Math. Meth. Appl. Sci. 2018, 41, 4867–4876. [Google Scholar] [CrossRef]
  12. Maleknejad, K.; Rashidinia, J.; Eftekhari, T. Numerical solution of three-dimensional Volterra-Fredholm integral equations of the first and second kinds based on Bernsteins approximation. Appl. Math. Comput. 2018, 339, 272–285. [Google Scholar]
  13. Mirzaee, F.; Hadadiyan, E. A computational method for nonlinear mixed Volterra-Fredholm integral equations. CJMS 2013, 2, 113–123. [Google Scholar]
  14. Mirzaee, F.; Hadadiyan, E. Three-dimensional triangular functions and their applications for solving nonlinear mixed Volterra-Fredholm integral equations. Alex. Eng. J. 2016, 55, 2943–2952. [Google Scholar] [CrossRef] [Green Version]
  15. Mirzaee, F.; Hadadiyan, E. Applying the modified block-pulse functions to solve the three-dimensional Volterra-Fredholm integral equations. Appl. Math. Comput. 2015, 265, 759–767. [Google Scholar] [CrossRef]
  16. Abd-Elhameed, W.M.; Youssri, Y.H. Generalized Lucas polynomial sequence approach for fractional differential equations. Nonlinear Dyn. 2017, 89, 1341–1355. [Google Scholar] [CrossRef]
  17. Abd-elhameed, W.M.; Youssri, Y.H. Spectral solutions for fractional differential equations via a novel lucas operational matrix of fractional derivatives. Rom. J. Phys. 2016, 61, 795–813. [Google Scholar]
  18. Cetin, M.; Sezer, N.; Gler, C. Lucas Polynomial Approach for System of High-Order Linear Differential Equations and Residual Error Estimation. Math. Eng. 2015, 2015, 625984. [Google Scholar]
  19. Baykus, N.M. Sezer: Hybrid Taylor-Lucas collocation method for numerical solution of high-order pantograph type delay differential equations with variables delays. Appl. Math. Inf. Sci. 2017, 11, 1795–1801. [Google Scholar] [CrossRef]
  20. Nadir, M. Lucas polynomials for solving linear integral equations. Theor. Appl. Comput. Sci. 2017, 11, 13–19. [Google Scholar]
  21. Haq, S.; Ali, I. Approximate solution of two.dimensional Sobolev equation using a mixed Lucas and Fibonacci polynomials. Eng. Comput. 2021, 3, 2059–2068. [Google Scholar]
  22. Ali, I.; Haq, S.; Nisar, K.S.; Arifeen, S. Numerical study of 1D and 2D advection-diffusion-reaction equations using Lucas and Fibonacci polynomials. Arab. J. Math. 2021, 10, 513–526. [Google Scholar] [CrossRef]
  23. Koshy, T. Fibonacci and Lucas Numbers with Applications; Wiley-Interscience Publication: New York, NY, USA, 2001. [Google Scholar]
  24. Mahdy, A.M.S.; Shokry, D.; Lotfy, K. Chelyshkov polynomials strategy for solving 2-dimensional nonlinear Volterra integral equations of the first kind. Comput. Appl. Math. 2022, 41, 1–13. [Google Scholar] [CrossRef]
  25. Mahdy, A.M.S.; Mohamed, D.S. Approximate solution of Cauchy integral equations by using Lucas polynomials. Comput. Appl. Math. 2022, 41, 1–20. [Google Scholar] [CrossRef]
  26. Youssri, Y.H.; Sayed, S.H.M.; Mohamed, A.S.; Aboeldahab, M.E.; Abd-Elhameed, W.M. Modified Lucas polynomials for the numerical treatment of second-order boundary value problems. Comput. Differ. Equations 2023, 11, 12–31. [Google Scholar]
Figure 1. Absolute error of example 1, M = 2 and z = 0.5.
Figure 1. Absolute error of example 1, M = 2 and z = 0.5.
Fractalfract 07 00196 g001
Figure 2. Absolute error of example 2, M = 2, z = 0.5.
Figure 2. Absolute error of example 2, M = 2, z = 0.5.
Fractalfract 07 00196 g002
Figure 3. Absolute error of example 3, case I, M = 2 and z = 0.5.
Figure 3. Absolute error of example 3, case I, M = 2 and z = 0.5.
Fractalfract 07 00196 g003
Figure 4. Absolute error of example 3 case II, M = 2 and z = 0.5.
Figure 4. Absolute error of example 3 case II, M = 2 and z = 0.5.
Fractalfract 07 00196 g004
Figure 5. Absolute error of example 4, M = 2 and z = 0.5.
Figure 5. Absolute error of example 4, M = 2 and z = 0.5.
Fractalfract 07 00196 g005
Table 1. The initially smaller Lucas polynomials are the array of these factors.
Table 1. The initially smaller Lucas polynomials are the array of these factors.
Polynomials’ Lucas Coefficients’ Array 
L 1 ( ζ ) = ζ 1
L 2 ( ζ ) = ζ 2 + 2 12
L 3 ( ζ ) = ζ 3 + 3 ζ 13
L 4 ( ζ ) = ζ 4 + 4 ζ 2 + 2 142  
L 5 ( ζ ) = ζ 5 + 5 ζ 3 + 5 ζ 155   
L 6 ( ζ ) = ζ 6 + 6 ζ 4 + 9 ζ 2 + 2 1692
L 7 ( ζ ) = ζ 7 + 7 ζ 5 + 14 ζ 3 + 7 ζ 17147
L 8 ( ζ ) = ζ 8 + 8 ζ 6 + 20 ζ 4 + 16 ζ 2 + 2 1820162
L 9 ( ζ ) = ζ 9 + 9 ζ 7 + 27 ζ 5 + 30 ζ 3 + 9 ζ 1927309
L 10 ( ζ ) = ζ 10 + 10 ζ 8 + 35 ζ 6 + 50 ζ 4 + 25 ζ 2 + 2 1103550252
Table 2. Numerical outcomes of example 1, case II.
Table 2. Numerical outcomes of example 1, case II.
( x , y , z ) Exact SolutionApproximate SolutionAbsolute Error
( 0 , 0 , 0 ) 000
( 0.1 , 0.1 , 0.1 ) 0.000001 9.999999471 × 10 7 5.295114830 × 10 14
( 0.2 , 0.2 , 0.2 ) 0.000064 0.00006399999982 1.868539870 × 10 13
( 0.3 , 0.3 , 0.3 ) 0.000729 0.0007289999997 3.341722165 × 10 13
( 0.4 , 0.4 , 0.4 ) 0.004096 0.004095999998 1.280713568 × 10 12
( 0.5 , 0.5 , 0.5 ) 0.015625 0.01562499999 5.782493832 × 10 12
( 0.6 , 0.6 , 0.6 ) 0.046656 0.04665599996 1.968260082 × 10 11
( 0.7 , 0.7 , 0.7 ) 0.117649 0.1176490000 5.302805863 × 10 11
( 0.8 , 0.8 , 0.8 ) 0.262144 0.2621439999 1.211866911 × 10 10
( 0.9 , 0.9 , 0.9 ) 0.531441 0.5314409998 2.459639868 × 10 10
( 1 , 1 , 1 ) 1 0.9999999996 4.567199615 × 10 10
Table 3. Numerical outcomes of example 2, case II.
Table 3. Numerical outcomes of example 2, case II.
( x , y , z ) Abs Error of Our Technique, M = 2Abs Error of [10], M = 4Abs Error of [10], M = 8
( 0 , 0 , 0 ) 000
( 0.1 , 0.1 , 0.1 ) 1.027090710 × 10 12 3.513425711 × 10 9 8.93185232 × 10 10
( 0.2 , 0.2 , 0.2 ) 4.734244968 × 10 12 3.075653759 × 10 8 7.81583453 × 10 9
( 0.3 , 0.3 , 0.3 ) 7.843093150 × 10 12 1.127441345 × 10 7 2.86408818 × 10 8
( 0.4 , 0.4 , 0.4 ) 6.348725880 × 10 12 2.884384112 × 10 7 7.32523191 × 10 8
( 0.5 , 0.5 , 0.5 ) 7.60305950 × 10 13 6.047489575 × 10 7 1.53545196 × 10 7
( 0.6 , 0.6 , 0.6 ) 1.066256719 × 10 11 1.116532759 × 10 6 2.83425621 × 10 7
( 0.7 , 0.7 , 0.7 ) 1.874489071 × 10 11 1.886594200 × 10 6 4.78810760 × 10 7
( 0.8 , 0.8 , 0.8 ) 2.352410432 × 10 11 2.985685054 × 10 6 7.57628837 × 10 7
( 0.9 , 0.9 , 0.9 ) 3.444933645 × 10 11 4.492504497 × 10 6 1.13981913 × 10 6
( 1 , 1 , 1 ) 8.2584895 × 10 11 6.493699097 × 10 6 1.64733199 × 10 6
Table 4. Numeric outcomes of example 3.
Table 4. Numeric outcomes of example 3.
( x , y , z ) = ( 2 , 2 , 2 ) Abs. Error for Case I, M = 2Error for Case II, M = 2
= 1 4.641 × 10 8 6.08 × 10 8
= 2 1.9520440 × 10 4 3.1952140 × 10 4
= 3 7.221208467 × 10 5 7.221218 × 10 5
= 4 1.147130303 × 10 5 1.1471321 × 10 5
= 5 1.600822047 × 10 6 1.6008264 × 10 6
= 6 2.109874999 × 10 7 2.109887 × 10 7
Table 5. Numeric outcomes of example 4.
Table 5. Numeric outcomes of example 4.
( x , y , z ) = ( 2 , 2 , 2 ) Abs. Error for  M = 2
= 1 9.8458149 × 10 5
= 2 1.957533070 × 10 3
= 3 6.463230447 × 10 4
= 4 1.190393926 × 10 4
= 5 1.772988091 × 10 5
= 6 2.409198170 × 10 6
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

Mahdy, A.M.S.; Nagdy, A.S.; Hashem, K.M.; Mohamed, D.S. A Computational Technique for Solving Three-Dimensional Mixed Volterra–Fredholm Integral Equations. Fractal Fract. 2023, 7, 196. https://doi.org/10.3390/fractalfract7020196

AMA Style

Mahdy AMS, Nagdy AS, Hashem KM, Mohamed DS. A Computational Technique for Solving Three-Dimensional Mixed Volterra–Fredholm Integral Equations. Fractal and Fractional. 2023; 7(2):196. https://doi.org/10.3390/fractalfract7020196

Chicago/Turabian Style

Mahdy, Amr M. S., Abbas S. Nagdy, Khaled M. Hashem, and Doaa Sh. Mohamed. 2023. "A Computational Technique for Solving Three-Dimensional Mixed Volterra–Fredholm Integral Equations" Fractal and Fractional 7, no. 2: 196. https://doi.org/10.3390/fractalfract7020196

APA Style

Mahdy, A. M. S., Nagdy, A. S., Hashem, K. M., & Mohamed, D. S. (2023). A Computational Technique for Solving Three-Dimensional Mixed Volterra–Fredholm Integral Equations. Fractal and Fractional, 7(2), 196. https://doi.org/10.3390/fractalfract7020196

Article Metrics

Back to TopTop