Solvability and stability of a fractional dynamical system of the growth of COVID-19 with approximate solution by fractional Chebyshev polynomials

Lately, many studies were offered to introduce the population dynamics of COVID-19. In this investigation, we extend different physical conditions of the growth by employing fractional calculus. We study a system of coupled differential equations, which describes the dynamics of the infection spreading between infected and asymptomatic styles. The healthy population properties are measured due to the social meeting. The result is associated with a macroscopic law for the population. This dynamic system is appropriate to describe the performance of growth rate of the infection and to verify if its control is appropriately employed. A unique solution, under self-mapping possessions, is investigated. Approximate solutions are presented by utilizing fractional integral of Chebyshev polynomials. Our methodology is based on the Atangana–Baleanu calculus, which provides various activity results in the simulation. We tested the suggested system by using live data. We found positive action in the graphs.


Introduction
Coronavirus disease COVID-19 is an infectious disease caused by a newly discovered coronavirus. It has been diffusing quickly in the world and the World Health Organization (WHO) characterized it as a pandemic. The first WHO indication of dyed-in-the-wool situations of COVID-19 was labeled on January 21, 2020 with 282 recognized cases, which is developed with the most present certificate on March 18, 2020, which extends to 191,127 concluded cases (see [1,2]). Normal growth approaches have been tested to describe the time development of the COVID-19 infection [3]. Fundamentally, by applying the system d dt where ϕ denotes the sum of infected people and the spreading phase, the increasing number of asymptomatic infected individuals was labeled. Recently, an extensive presentation of the fractal-fractional dynamic system of COVID-19 spread was introduced by Atangana [4]. The current investigation concerns the fractional dynamic system of the growth laws by exploiting the idea of fractional calculus. This idea contains an important term, which is the exponential law to find and accept the graph of the growth. The existence and uniqueness consequences are deliberated in the application of the fixed-point theory of self-mappings. Other properties are observed, such as the approximate solvability using the fractional Chebyshev polynomials.

Fractional dynamic system (FDS)
In this section, we construct the dynamic system of coupled equations. Before that, we need the following preliminaries about the conformable calculus.

Atangana-Baleanu calculus (ABC)
In recent decades, several physical complications have been displayed when employing the fractional calculus. The fundamental clarifications for applying fractional derivative illustrations are that various arrangements, constructions, and inequalities show the ability to remember past, or nonlocal properties, which cannot be stimulated using normal order derivatives. The fundamental ideas and applications of fractional calculus and fractional differential equations can now be found in many surveys. While most of the principal studies were based on the procedure of the Riemann-Liouville fractional order derivative, or the Caputo fractional order derivative, it has been observed that these derivatives have the property that their kernels have a singularity at the end of the interval of interest. The essential differences among the arbitrary derivatives are their unlike kernels which can be selected to fit the requirements of different applications. For example, the central variations between the Caputo fractional derivative, the Caputo-Fabrizio derivative [5], and others are that the Caputo calculus is expressed by employing a power law, while the Caputo-Fabrizio derivative is characterized by using an exponential decay. Atangana-Baleanu operator is introduced by suggesting the generalized Mittag-Leffler function [6].
Definition 1 (Fractional differential operator) A differential operator ν , ν ∈ (0, 1) is called a fractional Atangana-Baleanu derivative of order ν of a function ϕ if and only if ν can be written as where Ξ indicates the Mittag-Leffler function. The fractional integral is formulated by

Construction of FDS
In the construction of FDS, we denote by Ξ (t) the increasing overall number of infected persons, which is the sum of the number of the increasing recognized infected individuals ϕ(t) and of the asymptomatic transmissions ψ(t), that is, Ξ (t) = ϕ(t) + ψ(t). Regarding the statistics of ϕ(t), the number of transitions on and off, and cured people are convoluted, because they have been previously ill. Therefore, there are rate functions linking ϕ and ψ. We formulate the coupled-FDS as follows: where α, α 1 and β, β 1 are the connection rate continuous functions of ψ in μ ϕ(t) and ϕ in ν ψ(t), respectively. They describe the damping properties in line for the control energy.

Results
In this section, we proceed to discuss the solution existence and uniqueness for system (1). Moreover, we investigate the controller solution from different views.

Stability of solution
In this section, we deal with the stability of the unique solution via fixed point theorem. System (1) can be expressed by the general system satisfying the following hypotheses: (A1) Assume that X : [0, T] × R × R → R is a nondecreasing continuously differentiable function with X(0, 0, 0) = 0 and nonvanishing in a compact interval (0, T]. Furthermore, there is a positive constant κ such that (A2) Assume that A : [0, T] × R × R → R is a nondecreasing continuously differentiable function with A(0, 0, 0) = 0 and nonvanishing in a compact interval (0, T]. In addition, assume that there exists a positive constant K such that We aim to establish the existence and uniqueness of solution to system (2) using selfmapping fixed point theorem [7]. Put ℵ = R and define an operator P : R × R → R × R as follows: Since (ϕ, ψ) ∈ R × R, P is a self-mapping.

Lemma 3.2 Let the functions B
: R 3 → R + be defined as follows: Then the function B ∈ R forms a metric.
This metric forms the maximum measurement between the three cases of growth of COVID-19. Note that this metric can be extended to include other cases in dynamic systems.

Theorem 3.3
Suppose that the dynamic system (2) satisfies hypotheses (A1) and (A2). If the positive constants κ and K are such that , T < ∞, then P has a unique fixed point in the ball B r , where r ≤ 1.
Proof In view of the assumption on κ, and the definition of the metric in Lemma 3.2, we have This proves the boundedness of the operator P 1 in the unit ball B r1 of radius 0 < r 1 < 1. Similarly for P 2 , which is bounded in the ball B r2 , 0 < r 2 < 1. Combining the above conclusions, we obtain that the operator P = (P 1 , We proceed to investigate other properties of operator P 1 . Let t, τ ∈ (0, T) be such that if t > τ then ϕ(t) > ϕ(τ ) (increasing function). A simple calculation implies that Thus, P 1 is equicontinuous on B r1 . Similarly for P 2 , Thus, the integral operator P is equicontinuous on B r .
Hence, this implies that inequality (3) holds. Similarly, for P 2 , which implies that the integral operator P has a unique fixed point lying in B r = (B r1 , B r2 ), r ≤ 1.

(t), α(t)}, max{β 1 (t), β(t)}).
Proof Define an operator Q : R × R → R × R as follows: This yields the boundedness of the operator Q 1 in the unit ball B α max . Similarly for Q 2 , which is bounded in the ball B β max . Combining the above conclusions, we obtain that the operator Q = (Q 1 , Q 2 ) is bounded in B r = (B α max , B β max ). We proceed to investigate other properties of operator Q 1 . Let t, τ ∈ (0, T) be such that if t > τ then ϕ(t) > ϕ(τ ) (increasing function). A simple calculation implies that Thus, Q 1 is equicontinuous on B α max . Similarly for Q 2 , Thus, the integral operator Q is equicontinuous on B r = (B α max , B β max ).

Approximate solvability
In this section, we consider a generalization for the Chebyshev polynomials of the first type by using the ABC operator. We shall present two cases. The first one uses constant connections, and is called the symmetric solution. This case represents the setting when both ϕ and ψ have the same number of infected and cured. While the second case considers the connections as functions with respect to t.

Symmetric solvability with constant connections
Define the expanded formula of the solution by where T n (t) indicates the Chebyshev polynomials of the first kind such that T 0 (t) = 1, T n+1 (t) = 2tT n (t) -T n-1 (t), n ≥ 1. Chebyshev polynomials are of unlimited significance in various parts of mathematics, mainly approximation theory. The integrals of Chebyshev polynomials are (see [8]) The solution of (1) is given by the following construction: Now, by using the definition of ϕ in (7) and the integral formula of T n in (8), we have By symmetry, we obtain By the assumption t ≤ t ν < T ν ≤ ν ≤ 1 (see Theorem 3.4), we have that the asymptotic behavior of the Chebyshev polynomials is T n (t) ∼ 1, ∀n, t → 1.
Thus, the finite case of (11) becomes Consequently, by the convexity of the functions ϕ(t) = ψ(t) which are majored by 1 1-t , t ∈ (0, 1), we have the following construction: , , . . . Thus, the approximate symmetric solution can be seen as follows: ,

Symmetric solvability with functional connections
In this case, we have the following power series of the solution of (1): The approximate solvability of (1) can be presented in the next result.
China statistics in March has steady circulation; consequently, we propose (ϕ 1 , ψ 1 ) = ( 2t+2 π , 2t+2 π ). For USA data, the chart shows high rising confirmed cases, therefore we apply exponential connections ϕ 1 (t)(exp( ϕ 1 (t)) -1) similarly for ψ(t). Note that the data are shown in March   Figure 2 The dynamic evolution of system (21) when ν = 0.5, with the approximate solution by fractional Chebyshev polynomials (ϕ 2 , ψ 2 ) = ( 2(t+2t 2 ) π , 2(t+2t 2 ) π ) for Russia in March and April, respectively. In March the number of infections was per person, while in April it was per K (thousand). The connection coefficient in March is C 0.5,2 = 0.79, while in April it is C 0.5,2 = 0.099 but it is still approximated by (ϕ 2 , ψ 2 ). We confirmed that the approximate solution by fractional Chebyshev polynomials fits the future expectation of the number of infections.
We added also the case of Brazil. The picture is similar also for Brazil data, which indicates huge changes from March to April (see Fig. 3).

Figure 3
The dynamic evolution of system (21) when ν = 0.5, for Brazil in March and April, respectively. In March the number of infections was per person, therefore, we used (ϕ 3 , ψ 3 ). In April the data was per K = 1000, therefore, we find that (ϕ 2 , ψ 2 ) is a suitable solution

Conclusion
Based on the above, we conclude that the fractional-fractal dynamic system based on the Atangana-Baleanu fractional operator indicates flexibility and accuracy of introducing approximate solutions by fractional Chebyshev polynomials. For our future work, we aim to use the same calculus (ABC) to generalize different polynomials to get an optimal solution.