A Novel Vieta–Fibonacci Projection Method for Solving a System of Fractional Integrodifferential Equations

: In this paper, a new approach for numerically solving the system of fractional integrodiffer-ential equations is devised. To approximate the issue, we employ Vieta–Fibonacci polynomials as basis functions and derive the projection method for Caputo fractional order for the ﬁrst time. An efﬁcient transformation reduces the problem to a system of two independent equations. Solving two algebraic equations yields an approximate solution to the problem. The proposed method’s efﬁciency and accuracy are validated. We demonstrate the existence of the solution to the approximate problem and conduct an error analysis. Numerical tests reinforce the interpretations of the theory.

A lot of issues in the fields of engineering, mathematics, science, and related disciplines have recently been expressed in fractional integrodifferential equations.
The study in [15] examined whether a set of nonlinear Hadamard fractional differential equations has any positive solutions.These equations have non-negative nonlinear components and are defined on an infinite interval.Additionally, they are sensitive to nonlocal coupled boundary conditions that include Hadamard fractional derivatives and Riemann-Stieltjes integrals.The Guo-Krasnosel'skii fixed-point theorem and the Leggett-Williams fixed-point theorem were both used by the authors to establish the key theorems.
In [16] The authors established the existence and uniqueness of solutions for a system of Caputo fractional differential equations with sequential derivatives, integral terms, two positive factors, and universal coupled Riemann-Stieltjes integral boundary conditions.Their results' proofs are founded on the Banach fixed-point theorem and the Leray-Schau alternative.The objective of [17] is to establish a discontinuous Galerkin method with a one-sided flux for a singularly perturbed regular Volterra integrodifferential equation.In [18], the authors proposed adequate requirements for the estimated controllability of a class of fractional differential systems in Banach space with a fixed delay.The existence of a moderate solution is examined using the fixed-point theorem.Using generalized Gronwall's inequality, Cauchy sequence, and functional analysis fundamentals, controllability results are derived.
In [19], the authors established the existence and uniqueness of results for a system of coupled differential equations involving both left Caputo and right Riemann-Liouville fractional derivatives and mixed fractional integrals, supplemented with nonlocal coupled fractional integral boundary conditions using standard fixed-point theorems.The objective of [20] is to devise an approximate spectral method for the nonlinear time-fractional partial integrodifferential equation with a weakly singular kernel.This strategy is predicated on establishing a new Hilbert space that satisfies the initial and boundary conditions.The new spectral collocation method is applied to acquire a precise numerical approximation using new basis functions based on shifted first-kind Chebyshev polynomials.The paper [21] investigated an effective method for solving fractional integrodifferential equations with various scientific applications.The proposed method is founded on the orthonormal polynomial of Legendre and the least squares method.The authors of [22] presented, analyzed, and devised spectral collocation procedures for a particular class of nonlinear singular Lane-Emden equations with generalized Caputo derivatives that manifest in the study of astronomical objects.Under the supposition that the exact solution is an element in L 2 , the proposed solution is approximated as a truncated series of normalized shifted Jacobi polynomials.The spectral collocation method is used as a solver to obtain the unknown expansion coefficients.
The purpose of [23] is to develop a new operational matrix method based on Vieta-Fibonacci polynomials and validate the proposed method using it to find an approximation to the solution of the fractional-order nonlinear reaction-advection-diffusion equation of the form subject to ( The study in question [24] examined the following specific category of fractional optimal control problem: which was addressed by the utilization of an approximate methodology relying on fractional shifted Vieta-Fibonacci functions.The operational matrix of fractional integral for this method is derived.In [25], the numerical solution of the following fractional-order stochastic integrodifferential equations is obtained using an operational matrix approach based on shifted Vieta-Fibonacci polynomials.Here, B is the Brownian motion.
The authors of [26] used a Vieta-Fibonacci collocation method to find a numerical solution to fractional delay integrodifferential equations with weakly singular kernels using the collocation method and resultant operational matrices linked to the Vieta-Fibonacci polynomials.Krasnoselskii's fixed-point theorem shows that the solution to this type of fractional delay singular integrodifferential equation exists and is unique.
Recently, many models, such as those in control interpreting, elasticity, signal analysis, and unusual propagation, have been formulated regarding fractional derivatives.These considerations lead us to propose a numerical method for solving a system of Caputo fractional integrodifferential equations.The reason for conducting our research on the Caputo derivative is because this kind of fractional operator, one of the most widely used fractional derivatives, in greater depth and takes into account the realistic initial conditions employed in physics.The following are some of the overarching goals of this work:

•
Examining a system of fractional integrodifferential equations for the first time; • The utilization of Vieta-Fibonacci polynomials was initially employed to numerically obtain the solution a system of fractional integrodifferential equations; • Developing a projection method for addressing a system of fractional integrodifferential equations using Vieta-Fibonacci polynomials; • Establishing a new error analysis in a Vieta-Fibonacci weighted space.
In this study, a new way of solving the system of fractional integrodifferential equations numerically is proposed.We use Vieta-Fibonacci polynomials as basis functions and derive for the first time the projection method for Caputo fractional order to obtain the approximation solution.The problem is reduced to a set of two separate equations by an excellent transformation.By solving two algebraic problems, we obtain the approximate solution to the problem.The suggested method is proven to work and be accurate.We show that there is a solution to the approximate problem.Numerical tests back up how the theory is interpreted.

Preliminaries
In this section, we begin by reviewing some of the basic terms and theoretical concepts used in fractional theory.
The Vieta-Fibonacci polynomials in H are used in this paper to provide a projection method for solving a coupled system of fractional integrodifferential equations of the type: where δ and λ are the unknown solutions and σ, ς and ψ(., .)are given real smooth functions.
The function ψ(., .) is called the kernel.We transform the above-coupled system into a system of two separable fractional integrodifferential equations, which is subsequently solved using the current method.To this end, as Mennouni suggests in [29,30], investigate the following transformation: Lemma 1. Problem ( 16) can be expressed in the following form: Proof.We have, By putting them into (16), we obtain Equations ( 20) and ( 21) are added together, and then ( 20) is subtracted from (21) to yield (18).
Denoting by A the integral operator, i.e., and setting System (18) can be represented this way in operator form: (23) Recall that the operator A is compact from H into itself.Also, and Moreover, J η 0 + : H → D is compact.

Vieta-Fibonacci 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-Fibonacci polynomials using recurrence relations and an analytical formula.
Vieta-Fibonacci polynomials F n of degree n ∈ N are defined as follows: The following iterative formula can be used to generate polynomial F n : Also, the explicit power series formula shown below can be used to calculate F n : In addition, F n are orthogonal polynomials with respect to the integral shown below: denote the corresponding orthogonal polynomials on [0, 1].Additionally, F T k are produced using the recurrence formula shown below: with We note that Moreover, where Letting Now, we introduce the first seven terms of F N j : Let π F n be the chain of bounded finite rank orthogonal projections described by Denote by • the corresponding norm on H. Thus, Let H n represent the space covered by the first n-orthonormal shifted Vieta-Fibonacci polynomials.Recall that

Development of the Method
Clearly that Λ n , ∆ n ∈ D ∩ H n+1 .Therefore, the system is approached by Alternatively, We suppose that −1 is not eigenvalue of We obtain 2n + 2 unknowns β n,k and γ n,k by solving the following two separate linear systems, As a result, two separate linear systems are produced, where, for (56)

Convergence Analysis
We now show how the current method converges.To that end, consider L 2 ([0, 1], R) to be the weighted space and .to be its norm.
Denote with I the identity operator.We recall that there exists B > 0 such that Since for n large enough and are uniformly bounded with respect to n. and (59) Proof.In fact, Moreover, In addition, Hence, Letting we obtain ≤ Υ B 3 (n + 1) 3 + B 4 (n + 1) 3 , for some constants B 3 , B 4 > 0.

Numerical Example
This section establishes numerical experiments to illustrate the results stated in the preceding one.The Maple programming language was used in these numerical evaluations.
Example 1.We study the fractional integrodifferential system (16) in this example, which has the exact solution as follows: and ψ(κ, ρ) = ρ − κ. So, and Also, We perform some numerical tests to demonstrate the effectiveness of this example.For example, for n = 5, unknowns β 5,0 • • • β 5,5 are described as follows: The approximate solution Λ 5 is offered by The achieved approximate solutions and exacts one for Example 1 are shown in Figure 1.The numerical outcomes for the second case using the current approximation are shown in Figure 2.

Conclusions
In this research, we generalized a projection method to solve a system fractional integrodifferential via the shifted Vieta-Fibonacci polynomials.The current method reduces the provided problem into two systems of algebraic equations.Approximate solutions to the given problem are obtained by solving the acquired systems.This integrodifferential system has obvious relevance to problems in mathematical research, especially those involving interactions in physics.The fractional operator has been found to affect the growth of numerical results significantly.This method can be used to study and find solutions to various fractional integrodifferential and integral problems.This method has the potential to be employed as a future project for solving the system of fractional partial differential equations outlined below under some conditions:

Figure 2 .
Figure 2. Comparison of exact solutions Λ and ∆ and approximate ones Ω n and ∆ n , respectively, for n = 10.