Next Article in Journal
Statistical Validation of the “ECODIES” Questionnaire to Measure the Digital Competence of Colombian High School Students in the Subject of Mathematics
Next Article in Special Issue
From Wavelet Analysis to Fractional Calculus: A Review
Previous Article in Journal
Design of Finite Time Reduced Order H∞ Controller for Linear Discrete Time Systems
Previous Article in Special Issue
Parameter Identification for Lithium-Ion Battery Based on Hybrid Genetic–Fractional Beetle Swarm Optimization Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Projection Method for a System of Fractional Cauchy Integro-Differential Equations via Vieta–Lucas Polynomials

by
Abdelkader Moumen
1,† and
Abdelaziz Mennouni
2,*,†
1
Department of Mathematics, Faculty of Science, University of Ha’il, Ha’il 55425, Saudi Arabia
2
Department of Mathematics, LTM, University of Batna 2, Mostefa Ben Boulaïd, Fesdis, Batna 05078, Algeria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 16 November 2022 / Revised: 14 December 2022 / Accepted: 17 December 2022 / Published: 22 December 2022

Abstract

:
This work presents a projection method based on Vieta–Lucas polynomials and an effective approach to solve a Cauchy-type fractional integro-differential equation system. The suggested established model overcomes two linear equation systems. We prove the existence of the problem’s approximate solution and conduct an error analysis in a weighted space. The theoretical results are numerically supported.

1. Introduction

In recent years, numerous issues in mathematics, engineering, physics, and allied fields have been formulated in integral equations, notably, singular integral equations. Singular integral equations with a Cauchy kernel are an important class of such equations.
Fractional-order differential equations have been utilized extensively in quantum mechanics, astrobiology, medical science, chemical engineering, robust control, engineering, physiology, and hydrodynamics to describe various events.
The purpose of [1] is to investigate the dynamical evolution of symmetric oscillator with a fractional Caputo operator. The dynamical properties of the considered model such as equilibria and its stability are also examined. The existence, results, and uniqueness of proposed model solutions are considered using methods from fixed point theory. The authors of [2] decoupled a Lotka–Volterra model to explore the critical typical form coefficients of bifurcations for one-parameter and two-parameter bifurcations using a newly disclosed nonstandard finite difference method. In [3], a neural network strategy for solving the spatiotemporal fractional advection–diffusion equation with a nonlinear source term was described. Utilizing shifted Legendre orthogonal polynomials with variable coefficients, the network is created. The loss function of a neural network can be determined theoretically based on the features of unstable fractional derivatives. Multiple and generic discrete-time planar bifurcations were examined in [4]; using bifurcation theory and numerical continuation approaches, the Hindmarsh–Rose oscillator was analyzed for the study of three types of one-parameter bifurcation and five types of two-parameter bifurcation. Complex dynamics of the Kopel model with non-symmetric response among oligopolies were described in [5]. The studies indicated that a fixed point in a non-symmetric model may undergo fold, transcriptional, pitchfork, and Neimark–Sacker bifurcation under specific parameters.
The most valuable features for simulating functional problems and mathematical modeling are projection methods. These methods are efficient techniques for numerically solving integral and integro-differential equations. Over the past two decades, practical approaches for solving compact operator equations utilizing the Galerkin and Kulkarni approaches have been established. These two strategies inspired [6] to solve the following bounded equation. In [7], the authors recently proved how to solve fuzzy integro-differential equations with weak singularities using airfoil polynomials. The authors of [8] presented a novel projection approach based on Legendre polynomials for evaluating integro-differential equations of the Cauchy type. In addition, they investigated, in [9], a projection approximation for solving integro-differential problems of the Cauchy kind using first-order airfoil polynomials.
The orthogonality property of various significant polynomials, such as Vieta–Lucas polynomials, is used to approach the solution of different functional equations. These polynomials are crucial to solving functional equations.
The recurrence interaction of Vieta–Pell and Vieta–Pell–Lucas polynomials was introduced in [10]. The authors defined the associated sequences and obtained the Binet form and generating functions of Vieta–Pell and Vieta–Pell–Lucas polynomials. They also described some differentiation rules as well as finite summation formulas.
The authors of [11] developed a collocation approach based on a novel family of orthogonal functions for numerically treating a class of second-order singular multi-pantograph delay differential equations. The Vieta–Lucas functions were examined as differential bases, and the maximum norm errors were calculated.
The variable-order fractional form of the coupled nonlinear Ginzburg–Landau equations was formulated using the non-singular variable-order fractional derivative in the Heydari–Hosseininia concept in [12]. A numerical scheme based on shifted Vieta–Lucas polynomials was used to solve this system.
The authors of [13] proposed a fractional model of non-Newtonian Casson and Williamson boundary layer flow in fluid flow that takes the heat flux and slip velocity into account. The governing nonlinear system of PDEs is transformed into a nonlinear set of coupled ODEs, which is then solved using Vieta–Lucas polynomials, which are used to implement the spectral collocation method.
In [14], a new approximation technique combined with Vieta–Lucas orthogonal polynomials is used to solve the advection–dispersion equation, a fractional-order mathematical physics model.
In the last two decades, numerous results for solving compact operator equations using the Galerkin and Kulkarni approaches have been established. These two approaches are used to approach the solution of the following bounded equation, as cited in [6]:
x A x = y ,
where A is a bounded linear operator, y is a know function, and x is a unknown function. The author examined an approximation of the following equation:
x n G A n x n G = Π n y ,
with approximate solution x n and approximate operator A n . The author originally established approximate operator A n as a Galerkin approximation.
A n G : = Π n A Π n .
Here, G stands for Galerkin, and ( Π n ) n 1 is a sequence of bounded projections each one of finite rank, that is,
Π n 2 = Π n also Π n * = Π n .
Second, he utilized Kulkarni’s approximation in the following manner:
A n K : = Π n A + A Π n Π n A Π n .
More recently, the authors of [15] introduced a novel projection approach for the following system of classical Cauchy integro-differential equations on L 2 ( [ 0 , 1 ] , C ) using shifted Legendre polynomials.
φ ( τ ) + 0 1 ω ( ς ) ς τ d ς = h ( τ ) , 0 τ 1 , ω ( τ ) + 0 1 φ ( ς ) ς τ d ς = k ( τ ) , 0 τ 1 , φ ( 0 ) = 0 , ω ( 0 ) = 0 .
This work presents a projection method for solving the following system of fractional Cauchy integro-differential equations on L 2 ( [ 0 , 1 ] , C ) :
c D 0 + η φ ( τ ) + 0 1 ω ( ς ) ς τ d ς = h ( τ ) , 0 τ 1 , c D 0 + η ω ( τ ) + 0 1 φ ( ς ) ς τ d ς = k ( τ ) , 0 τ 1 , k = 0 r 1 τ k k ! φ ( k ) ( 0 ) = 0 , k = 0 r 1 τ k k ! ω ( k ) ( 0 ) = 0 , r 1 < η < r .
We turn the problem into a system of two separate equations that looks like this:
c D 0 + η Ω + K Ω = H , c D 0 + η ϕ K ϕ = K .
Our aim is employ Vieta–Lucas polynomials and present an approximation scheme to solve the above system. The solution is found for two different linear equation systems. The existence of a solution to the approximation equation is demonstrated, and an investigation of error analysis is presented. The theory is illustrated with numerical examples.
The remaining sections of this research paper are described in the following: The following section discusses some of the fundamental terms and theoretical concepts used in fractional theory. Section 3 discusses the system of fractional Cauchy integral equations. Section 4 discusses some key aspects of Vieta–Lucas polynomials and the development of the method. Section 5 improves the convergence of the approximate solution and estimates the error analysis. Section 6 explores some numerical examples.

2. Preliminaries

In this section, we begin by reviewing some of the basic terms and theoretical concepts used in fractional theory.
Let Γ indicate the fundamental Euler Gamma function in the analysis of fractional differential equations.
Definition 1.
The left-sided Riemann–Liouville fractional integral of order η > 0 of an integrable function ω : ( 0 , ) R is described by
J 0 + η ω ( τ ) : = 1 Γ ( η ) 0 τ ω ( ς ) ( τ ς ) 1 η d ς , τ > 0 .
Remark 1.
The above integral can be represented in convolution form as follows:
J 0 + η ω ( τ ) = ( ψ η 🟉 ω ) ( τ ) ,
where
ψ η ( τ ) : = τ η 1 / Γ ( η ) , τ > 0 , 0 , τ 0 .
Definition 2.
The left-sided Riemann–Liouville fractional derivative of order η > 0 of a continuous function ω : ( 0 , ) R is defined by
D 0 + η ω ( τ ) : = 1 Γ ( 1 η ) d d τ 0 τ ω ( ς ) ( τ ς ) η d ς .
Remark 2.
For η > 0 , we have
D 0 + η ω ( τ ) = d d τ J 0 + 1 η ω ( τ ) ,
J 0 + η D 0 + η ( ω ( τ ) ω ( 0 ) ) = ω ( τ ) ω ( 0 ) .
Definition 3.
For an absolutely continuous function ω, the Caputo fractional derivative of order η > 0 is defined by
c D 0 + η ω ( τ ) = J 0 + 1 η d d τ ω ( τ ) = 1 Γ ( 1 η ) 0 τ ( τ ς ) η ω ( ς ) d ς .
Remark 3.
For a continuous function ω, the relationship between the Caputo and Riemann–Liouville fractional derivatives is provided by
c D 0 + η ω ( τ ) = D 0 + η ( ω ( τ ) ω ( 0 ) ) , η > 0 .
In addition,
c D 0 + η ω ( τ ) = D 0 + η ω ( τ ) k = 0 r 1 τ k k ! ω ( k ) ( 0 ) , τ > 0 , r 1 < η < r .

3. System of Fractional Cauchy Integro-Differential Equations

Let H : = L 2 ( [ 0 , 1 ] , C ) be the space of complex-valued Lebesgue square integrable functions on [ 0 , 1 ] . In this study, a new projection technique for solving a system of fractional Cauchy integro-differential equations using Vieta–Lucas polynomials in H is presented.
Consider the system of fractional Cauchy integro-differential equations of the following form:
c D 0 + η φ ( τ ) + 0 1 ω ( ς ) ς τ d ς = h ( τ ) , 0 τ 1 , c D 0 + η ω ( τ ) + 0 1 φ ( ς ) ς τ d ς = k ( τ ) , 0 τ 1 , k = 0 r 1 τ k k ! φ ( k ) ( 0 ) = 0 , k = 0 r 1 τ k k ! ω ( k ) ( 0 ) = 0 , r 1 < η < r .
The integral of each equation denotes the Cauchy principal value as follows:
0 1 ω ( ς ) ς τ d ς = lim ε 0 0 τ ε ω ( ς ) ς τ d ς + τ + ε 1 ω ( ς ) ς τ d ς , 0 1 φ ( ς ) ς τ d ς = lim ε 0 0 τ ε φ ( ς ) ς τ d ς + τ + ε 1 φ ( ς ) ς τ d ς .
As in [16], let us consider the following transform:
ϕ : = ω φ , H : = k + h ;
Ω : = ω + φ , K : = k h .
Lemma 1.
Problem (1) can be expressed in the following form:
c D 0 + η Ω ( τ ) + 0 1 Ω ( ς ) ς τ d ς = H ( τ ) , 0 τ 1 , c D 0 + η ϕ ( τ ) 0 1 ϕ ( ς ) ς τ d ς = K ( τ ) , 0 τ 1 , k = 0 r 1 τ k k ! Ω ( k ) ( 0 ) = 0 , k = 0 r 1 τ k k ! ϕ ( k ) ( 0 ) = 0 , r 1 < η < r .
Proof. 
In fact,
ω = Ω + ϕ 2 , h = H K 2 ;
φ = Ω ϕ 2 , k = H + K 2 .
By substituting them into (1), we obtain
c D 0 + η Ω ϕ ( τ ) + 0 1 Ω + ϕ ( ς ) ς τ d ς = H K ( τ ) ,
c D 0 + η Ω + ϕ ( τ ) + 0 1 Ω ϕ ( ς ) ς τ d ς = H + K ( τ ) .
We derive (1) by adding Equations (1) and (2) together and then subtracting (1) from (2). □
System (1) can be rewritten in operator form as follows:
c D 0 + η Ω ( τ ) + K Ω ( τ ) = H ( τ ) , c D 0 + η ϕ ( τ ) K ϕ ( τ ) = K ( τ ) ,
where K is the Cauchy integral operator, i.e.,
K φ ( τ ) : = 0 1 φ ( ς ) ς τ d ς , 0 τ 1 .
Following [8], operator K is bounded from H into itself.
Letting
D : = ω H : ω ( k ) H , k = 0 r 1 τ k k ! ω ( k ) ( 0 ) = 0 , r 1 < η < r .
It is well known that
J 0 + η c D 0 + η ω ( τ ) = ω ( τ ) k = 0 r 1 τ k k ! ω ( k ) ( 0 ) .
So,
J 0 + η c D 0 + η ω ( τ ) = ω ( τ ) , for all ω D .
In addition, J 0 + η : H D is compact.

4. Vieta–Lucas Polynomials

In this section, we look at a class of orthogonal polynomials. These polynomials can be used to construct a new family of orthogonal polynomials known as Vieta–Lucas polynomials using recurrence relations and an analytical formula.
Vieta–Lucas polynomials V n of degree n N are defined as follows:
V n ( τ ) = 2 cos ( n θ ) , θ = cos 1 τ 2 , θ [ 0 , π ] for all τ 2 .
The following iterative formula can be used to generate polynomial V n :
V n ( τ ) = τ V n 1 ( τ ) V n 2 ( τ ) , n = 2 , 3 , , V 0 ( τ ) = 2 , V 1 ( τ ) = τ .
Also, the explicit power series formula shown below can be used to calculate V n :
V n ( τ ) = i = 0 n 2 ( 1 ) i n Γ ( n i ) Γ ( i + 1 ) Γ ( n + 1 2 i ) τ n 2 i , n = { 2 , 3 , } .
In addition, V n are orthogonal polynomials with respect to the integral shown below:
V m , V n = 2 2 1 4 τ 2 V m ( τ ) V n ( τ ) d τ = 0 , m n 0 , 4 π , m = n = 0 , 2 π , m = n 0 .
Let
V k S ( τ ) = V k ( 4 τ 2 ) = V 2 k ( 2 τ ) , k = 0 , 2 , 3 ,
denote the corresponding normalized sequence. Additionally, V k S are produced using the recurrence formula shown below:
V k + 1 S ( τ ) = ( 4 τ 2 ) V k S ( τ ) V k 1 S ( τ ) , k = 1 , 2 , ,
with
V 0 S ( τ ) = 2 , V 1 S ( τ ) = 4 τ 2 .
We note that
V p S ( τ ) = 2 p j = 0 p ( 1 ) j 4 p j Γ ( 2 p j ) Γ ( j + 1 ) Γ ( 2 p 2 j + 1 ) τ p j , p = { 2 , 3 , }
and
V k S , V j S ω = 0 1 V k S ( τ ) V j S ( τ ) ω ( τ ) d τ = 0 , k j 0 , 4 π , k = j = 0 , 2 π , k = j 0 ,
where
ω ( τ ) = 1 τ τ 2 .
Letting
Δ j S : = V j S 2 π , j = 0 , V j S 2 π , j 0 .
Now, we introduce the first six terms of Δ j S :
Δ 0 S ( s ) = 1 π , Δ 1 S ( s ) = 2 2 s 1 π , Δ 2 S ( s ) = 2 8 s 2 8 s + 1 π , Δ 3 S ( s ) = 2 32 s 3 48 s 2 + 18 s 1 π , Δ 4 S ( s ) = 2 128 s 4 256 s 3 + 160 s 2 32 s + 1 π , Δ 5 S ( s ) = 2 512 s 5 1280 s 4 + 1120 s 3 400 s 2 + 50 s 1 π .
Let π n S be the chain of bounded finite rank orthogonal projections described by
π n S ψ : = j = 0 n 1 ψ , Δ j S ω Δ j S , where ψ , Δ j S ω : = 0 1 ω ( σ ) ψ ( σ ) Δ j S ( σ ) d σ .
Denote with · the corresponding norm on H . Thus,
lim n π n S ϑ ϑ = 0 , for all ϑ H .
Let H n represent the space covered by the first n-shifted Vieta–Lucas polynomials. It is obvious that J 0 + η ( H n ) = H n + 1 .
We note that Ω n , ϕ n D H n + 1 . Thus, the system
Ω + J 0 + η K Ω = J 0 + η H , ϕ J 0 + η K ϕ = J 0 + η K
is approximated by
Ω n + J 0 + η π n S K Ω n = J 0 + η π n S H , ϕ n J 0 + η π n S K ϕ n = J 0 + η π n S K .
We assume that 1 and 1 are not eigenvalues of J 0 + η K . Thus, both operators I + J 0 + η K and I J 0 + η K are invertible.
We recall that J 0 + η is compact and
lim n J 0 + η π n S K J 0 + η K J 0 + η K = 0 , lim n J 0 + η π n S K J 0 + η K J 0 + η π n S K = 0 .
Writing
Ω n = j = 0 n a n , j Δ j S , ϕ n = j = 0 n b n , j Δ j S .
We obtain 2 n + 2 unknowns a n , j and b n , j by solving the two separate linear systems,
j = 0 n a n , j D 0 + η Δ j S + π n S K Δ j S = π n S J 0 + η H , with j = 0 n a n , j Δ j S ( 0 ) = 0 , j = 0 n b n , j D 0 + η Δ j S π n S K Δ j S = π n S J 0 + η K , with j = 0 n b n , j Δ j S ( 0 ) = 0 .
As a result, two separate linear systems are produced,
j = 0 n M n ( i , j ) a n , j = E n , i , i = 0 , , n , j = 0 n M ^ n ( i , j ) b n , j = E ^ n , i , i = 0 , , n .
where, for i = 0 , , n 1 and j = 0 , , n ,
M n ( i , j ) : = 0 1 c D 0 + η Δ j S ( s ) Δ i S ( s ) ω ( s ) d s + 0 1 0 1 Δ j S ( σ ) σ s d σ Δ i S ( s ) ω ( s ) d s , M ^ n ( i , j ) : = 0 1 c D 0 + η Δ j S ( s ) Δ i S ( s ) ω ( s ) d s 0 1 0 1 Δ j S ( σ ) σ s d σ Δ i S ( s ) ω ( s ) d s , M n ( n , j ) : = Δ j S ( 0 ) , M ^ n ( n , j ) : = Δ j S ( 0 ) , E n ( i ) : = 0 1 H ( s ) c D 0 + η Δ i S ( s ) ω ( s ) d s , E ^ n ( i ) : = 0 1 K ( s ) c D 0 + η Δ i S ( s ) ω ( s ) d s , E n ( n ) : = 0 , E ^ n ( n ) : = 0 .

5. Convergence Analysis

We now show how the current method converges. To that end, consider L ϖ 2 ( [ 0 , 1 ] , C ) to be the weighted space and . ϖ to be its norm.
Denote with I the identity operator. We recall that there exists M > 0 such that
( I π n S ) ψ ϖ M 12 n 3 / 2 , for all ψ L ϖ 2 ( [ 0 , 1 ] , C ) .
Since J 0 + η K is compact, according to [8], operators ( I + J 0 + η π n S K ) 1 and ( I J 0 + η π n S K ) 1 exist for n large enough and are uniformly bounded with respect to n.
Theorem 1.
Assume that k , h , K ω , K φ L ϖ 2 ( [ 0 , 1 ] , C ) . Then, there exist M 1 , M 2 > 0 such that
Ω n Ω ϖ M 1 12 n 3 / 2 ,
and
ϕ n ϕ ϖ M 2 12 n 3 / 2 .
Proof. 
In fact,
Ω n + J 0 + η π n S K Ω n = J 0 + η π n S H , ϕ n J 0 + η π n S K ϕ n = J 0 + η π n S K .
Moreover,
Ω Ω n = I + J 0 + η K 1 J 0 + η H I + J 0 + η π n S K 1 J 0 + η π n S H + I + J 0 + η π n S K 1 J 0 + η H I + J 0 + η π n S K 1 J 0 + η H = I + J 0 + η π n S K 1 J 0 + η I π n S K Ω + I π n S H .
In addition,
ϕ ϕ n = I J 0 + η K 1 J 0 + η K I J 0 + η π n S K 1 J 0 + η π n S K + I J 0 + η π n S K 1 J 0 + η K I J 0 + η π n S K 1 J 0 + η K = I J 0 + η π n S K 1 J 0 + η I π n S K ϕ + I π n S K .
Hence,
Ω Ω n = I + J 0 + η π n S K 1 J 0 + η I π n S K ( ω + φ ) + I π n S ( k + h ) ,
and thus
ϕ n ϕ = I J 0 + η π n S K 1 J 0 + η I π n S K ( ω φ ) + I π n S ( k h ) .
Letting
δ : = I + J 0 + η π n S K 1 J 0 + η , γ : = I J 0 + η π n S K 1 J 0 + η ,
we obtain
Ω n Ω ϖ δ I π n S K ( ω + φ ) ϖ + I π n S ( k + h ) ϖ δ I π n S K ( ω + φ ) ϖ + I π n S ( k + h ) ϖ δ M 3 12 n 3 / 2 + M 4 12 n 3 / 2 ,
for some constants M 3 , M 4 > 0 . Moreover,
ϕ n ϕ ϖ = I J 0 + η π n S K 1 J 0 + η I π n S K ϕ + J 0 + η I π n S K ϖ γ I π n S K ( ω φ ) ϖ + I π n S ( k h ) ϖ γ M 5 12 n 3 / 2 + M 6 12 n 3 / 2 ,
for some constants M 5 , M 6 > 0 .
Letting
M 1 : = δ max M 3 , M 4 a n d M 2 : = γ max M 5 , M 6 ,
we obtain the desired results. □

6. Numerical Example

In this section, numerical experiments are established to illustrate the results presented in the previous section. In these numerical evaluations, the Maple programming language was implemented.
Example 1.
We study the fractional integro-differential system (1) in this example, which has the exact solution as follows:
φ ( τ ) = 1 2 τ 3 , ω ( τ ) = 1 2 τ 3 2 τ 2 , η = 1 3 .
In this case, we obtain
Ω ( τ ) = τ 3 τ 2 , ϕ ( τ ) = τ 2 .
H ( τ ) = 1 120 Γ 2 3 [ 243 τ 8 / 3 + 120 Γ 2 3 τ 3 ln ( τ ) + 216 τ 5 / 3 120 Γ 2 3 τ 3 ln ( τ + 1 ) 120 Γ 2 3 τ 2 ln ( τ ) + 120 Γ 2 3 τ 2 ln ( τ + 1 ) 120 Γ 2 3 τ 2 + 60 τ + 20 Γ 2 3 ]
and
K ( τ ) = 1 10 Γ 2 3 [ 18 τ 5 / 3 10 Γ 2 3 τ 2 ln ( 1 τ ) + 216 τ 5 / 3 + 10 Γ 2 3 τ 2 ln ( τ ) 10 Γ 2 3 τ 5 Γ 2 3 ] .
To show the efficiency of this example, we shall offer some numerical testes. For example, for n = 7 , unknowns a 7 , 0 a 7 , 7 are as follows:
a 7 , 0 = 0.78334 × 10 1 , a 7 , 1 = 0.39165 × 10 1 , a 7 , 2 = 0.78330 × 10 1 , a 7 , 3 = 0.39165 × 10 1 , a 7 , 4 = 5.6631 × 10 7 , a 7 , 5 = 0.17120 × 10 5 , a 7 , 6 = 3.2531 × 10 7 , a 7 , 7 = 0.14416 × 10 5 .
Approximate solution Ω 7 is given by
Ω 7 ( τ ) = 0.16755 × 10 3 τ + 0.98870 τ 3 + 0.30814 × 10 1 τ 4 0.44449 × 10 1 τ 5 + 0.32446 × 10 1 τ 6 0.94226 × 10 2 τ 7 0.99779 τ 2 0.21845 × 10 5 .
In addition,
b 7 , 0 = 0.46998 , b 7 , 1 = 0.62666 , b 7 , 2 = 0.15668 , b 7 , 3 = 0.30641 × 10 5 , b 7 , 4 = 0.28020 × 10 5 , b 7 , 5 = 0.44393 × 10 5 , b 7 , 6 = 0.20660 × 10 5 , b 7 , 7 = 0.25847 × 10 5 .
Approximate solution ϕ 7 is provided by
ϕ 7 ( τ ) = 0.71009 × 10 3 τ + 0.29927 × 10 1 τ 3 1.0066 τ 2 0.11933 × 10 4 0.71647 × 10 1 τ 4 + 0.93245 × 10 1 τ 5 0.62505 × 10 1 τ 6 + 0.16894 × 10 1 τ 7 .
Table 1 shows the numerical findings obtained for Example 1 using the method we suggest.
Figure 1 depicts the graphs of the achieved absolute errors for Example 1. Here, we note that the τ symbol in the figures refers to the τ -axis labels.
Example 2.
Various types of fractional integro-differential and integral systems can be studied and solved with the proposed method. As a second example for testing, we consider the system of fractional logarithmic integro-differential equations of the following form:
c D 0 + η φ ( τ ) + 0 1 ω ( ς ) ln ς τ d ς = h ( τ ) , 0 τ 1 , c D 0 + η ω ( τ ) + 0 1 φ ( ς ) ln ς τ d ς = k ( τ ) , 0 τ 1 , k = 0 r 1 τ k k ! φ ( k ) ( 0 ) = 0 , k = 0 r 1 τ k k ! ω ( k ) ( 0 ) = 0 , r 1 < η < r ,
which has the following exact solution:
φ ( τ ) = 1 2 τ 3 + τ 4 , ω ( τ ) = 1 2 τ 4 τ 3 , η = 1 2 .
Here, we have
Ω ( τ ) = τ 3 , ϕ ( τ ) = τ 4 .
We provide some numerical tests to demonstrate the efficacy of this illustration. If n = 5 , then the values of a 5 , 0 a 5 , 5 are as shown below:
a 5 , 0 = 0.39168 , a 5 , 1 = 0.58757 , a 5 , 2 = 0.23503 , a 5 , 3 = 0.39169 × 10 1 , a 5 , 4 = 0.15704 × 10 5 , a 5 , 5 = 9.5589 × 10 7 .
Approximate solution Ω 5 is given by
Ω 5 ( τ ) = 0.15159 × 10 3 τ 0.30425 × 10 3 τ 2 0.99947 τ 3 + 0.39048 × 10 3 τ 5 + 0.68731 × 10 5 .
In addition,
b 5 , 0 = 0.34276 , b 5 , 1 = 0.54834 , b 5 , 2 = 0.27418 , b 5 , 3 = 0.78336 × 10 1 , b 5 , 4 = 0.97933 × 10 2 , b 5 , 5 = 5.6341 × 10 7 .
Approximate solution ϕ 5 is provided by
ϕ 5 ( τ ) = 0.79785 × 10 5 τ + 0.99955 τ 4 + 0.27127 × 10 3 τ 3 0.23936 × 10 4 τ 2 + 0.45188 × 10 4 + 0.23015 × 10 4 τ 5 .
Figure 2 illustrates the numerical results produced for Example 2 using the indicated approach.
Example 3.
Now, we consider the above system of fractional logarithmic integro-differential equations, with the same exact solutions.
φ ( τ ) = 1 2 τ 3 + τ 4 , ω ( τ ) = 1 2 τ 4 τ 3 .
In Table 2, we examine the influence of fractional order η on the approximate solutions.

7. Conclusions

The application of projection methods is extended in this study so that it can be applied to a fractional Cauchy singular integro-differential system. Shifted Vieta–Lucas polynomials are a class of orthogonal polynomials that serve as the foundation of this approach. The significance of this fractional Cauchy singular integro-differential system is evident in mathematical sciences issues, particularly in physics interactions. It has been observed that the fractional operator has a substantial impact on the development of the numerical results. Numerous kinds of fractional integro-differential and integral systems can be investigated with this method, and their solutions can be found.

Author Contributions

Conceptualization, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Methodology, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Investigation, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Resources, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Data curation, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Writing—original draft preparation, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Writing—review and editing, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Supervision, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni); Funding acquisition, A.M. (Abdelkader Moumen) and A.M. (Abdelaziz Mennouni). All authors have read and agreed to the published version of the manuscript.

Funding

The research of the first author was funded by King Khalid University Researchers Supporting project number RGP.1/350/43, King Khalid University, Abha, Saudi Arabia.

Data Availability Statement

Not applicable.

Acknowledgments

The first author extends his appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through Small Groups (RGP.1/350/43).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, C.; ur Rahman, M.; Baleanu, D. On fractional-order symmetric oscillator with offset-boosting control. Nonlinear Anal. Model. Control. 2022, 27, 994–1008. [Google Scholar] [CrossRef]
  2. Eskandari, Z.; Avazzadeh, Z.; Khoshsiar Ghaziani, R.; Li, B. Dynamics and bifurcations of a discrete-time Lotka–Volterra model using nonstandard finite difference discretization method. Math Meth. Appl. Sci. 2022. [Google Scholar] [CrossRef]
  3. Qu, H.D.; Liu, X.; Lu, X.; ur Rahman, M.; She, Z.H. Neural network method for solving nonlinear fractional advection-diffusion equation with spatiotemporal variable-order. Chaos Solitons Fractals 2022, 156, 111856. [Google Scholar] [CrossRef]
  4. Li, B.; Liang, H.; He, Q. Multiple and generic bifurcation analysis of a discrete Hindmarsh-Rose model. Chaos Solitons Fractals 2021, 146, 110856. [Google Scholar] [CrossRef]
  5. Li, B.; Liang, H.; Shi, L.; He, Q. Complex dynamics of Kopel model with nonsymmetric response between oligopolists. Chaos Solitons Fractals 2022, 156, 111860. [Google Scholar] [CrossRef]
  6. Mennouni, A. Two projection methods for skew-hermitian operator equations. Math. Comput. Model. 2012, 55, 1649–1654. [Google Scholar] [CrossRef]
  7. Araour, M.; Mennouni, A. A New Procedures for Solving Two Classes of Fuzzy Singular Integro-Differential Equations: Airfoil Collocation Methods. Int. J. Appl. Comput. Math. 2022, 8, 35. [Google Scholar] [CrossRef]
  8. Mennouni, A. A projection method for solving Cauchy singular integro-differential equations. Appl. Math. Lett. 2012, 25, 986–989. [Google Scholar] [CrossRef] [Green Version]
  9. Mennouni, A. Airfoil polynomials for solving integro-differential equations with logarithmic kernel. Appl. Math. Comput. 2012, 218, 11947–11951. [Google Scholar] [CrossRef]
  10. Tasci, D.; Yalcin, F. Vieta-Pell and Vieta-Pell-Lucas polynomials. Adv. Differ. Equ. 2013, 2013, 224. [Google Scholar] [CrossRef]
  11. Izadi, M.; Yüzbaşı, Ş.; Ansari, K.J. Application of Vieta–Lucas Series to Solve a Class of Multi-Pantograph Delay Differential Equations with Singularity. Symmetry 2021, 13, 2370. [Google Scholar] [CrossRef]
  12. Heydari, M.H.; Avazzadeh, Z.; Razzaghi, M. Vieta-Lucas polynomials for the coupled nonlinear variable-order fractional Ginzburg-Landau equations. Appl. Numer. Math. 2021, 165, 442–458. [Google Scholar] [CrossRef]
  13. Adel, M.; Assiri, T.A.; Khader, M.M.; Osman, M.S. Numerical simulation by using the spectral collocation optimization method associated with Vieta-Lucas polynomials for a fractional model of non-Newtonian fluid. Results Phys. 2022, 41, 105927. [Google Scholar] [CrossRef]
  14. Agarwal, P.; El-Sayed, A.A. Vieta–Lucas polynomials for solving a fractional-order mathematical physics model. Adv. Differ. Equ. 2020, 626, 2020. [Google Scholar] [CrossRef]
  15. Althubiti, S.; Mennouni, A. A Novel Projection Method for Cauchy-Type Systems of Singular Integro-Differential Equations. Mathematics 2022, 10, 2694. [Google Scholar] [CrossRef]
  16. Mennouni, A. A new efficient strategy for solving the system of Cauchy integral equations via two projection methods. Transylv. J. Math. Mech. 2022, 14, 63–71. [Google Scholar]
Figure 1. Comparison of exact solutions Ω and ϕ and approximate solutions Ω n and ϕ n , respectively, for the first example, where n = 7 .
Figure 1. Comparison of exact solutions Ω and ϕ and approximate solutions Ω n and ϕ n , respectively, for the first example, where n = 7 .
Mathematics 11 00032 g001
Figure 2. Comparison of exact solutions Ω and ϕ and approximate solutions Ω n and ϕ n , respectively, for the second example, where n = 5 .
Figure 2. Comparison of exact solutions Ω and ϕ and approximate solutions Ω n and ϕ n , respectively, for the second example, where n = 5 .
Mathematics 11 00032 g002
Table 1. Numerical results for Example 1.
Table 1. Numerical results for Example 1.
                  n                                  Ω Ω n ϖ                                    ϕ ϕ n ϖ                   
34.1428 × 10 5 1.5520 × 10 4
58.6994 × 10 6 3.5924 × 10 5
79.6542 × 10 8 8.2544 × 10 8
137.2543 × 10 13 7.2547 × 10 12
179.2541 × 10 15 3.8531 × 10 14
Table 2. Influence of fractional order η on the approximate solutions.
Table 2. Influence of fractional order η on the approximate solutions.
                  η                                  Ω Ω n ϖ                                    ϕ ϕ n ϖ                   
0.31.0548 × 10 5 3.1473 × 10 5
0.41.2668 × 10 4 1.0739 × 10 4
0.67.9977 × 10 6 3.2081 × 10 5
0.78.4737 × 10 5 4.7114 × 10 5
0.89.8905 × 10 5 9.2993 × 10 5
1.39.1313 × 10 5 9.7317 × 10 5
1.68.4393 × 10 5 2.3766 × 10 4
2.45.1952 × 10 5 7.0123 × 10 4
3.74.0924 × 10 5 1.4281 × 10 3
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

Moumen, A.; Mennouni, A. A New Projection Method for a System of Fractional Cauchy Integro-Differential Equations via Vieta–Lucas Polynomials. Mathematics 2023, 11, 32. https://0-doi-org.brum.beds.ac.uk/10.3390/math11010032

AMA Style

Moumen A, Mennouni A. A New Projection Method for a System of Fractional Cauchy Integro-Differential Equations via Vieta–Lucas Polynomials. Mathematics. 2023; 11(1):32. https://0-doi-org.brum.beds.ac.uk/10.3390/math11010032

Chicago/Turabian Style

Moumen, Abdelkader, and Abdelaziz Mennouni. 2023. "A New Projection Method for a System of Fractional Cauchy Integro-Differential Equations via Vieta–Lucas Polynomials" Mathematics 11, no. 1: 32. https://0-doi-org.brum.beds.ac.uk/10.3390/math11010032

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