Fractional dynamic system simulating the growth of microbe.

There are different approaches that indicate the dynamic of the growth of microbe. In this research, we simulate the growth by utilizing the concept of fractional calculus. We investigate a fractional system of integro-differential equations, which covers the subtleties of the diffusion between infected and asymptomatic cases. The suggested system is applicable to distinguish the presentation of growth level of the infection and to approve if its mechanism is positively active. An optimal solution under simulation mapping assets is considered. The estimated numerical solution is indicated by employing the fractional Tutte polynomials. Our methodology is based on the Atangana-Baleanu calculus (ABC). We assess the recommended system by utilizing real data.

Introduction

Integro-differential dynamic system of equations simulates various states from science and engineering corresponding to the analysis, control, and optimization studies. The main model in this direction is the Wilson–Cowan system, which designs the dynamics of connections between populations of very inhibitory system in cells or neurons. It was developed by Hugh R. Wilson and Jack D. Cowan [1]. The system and its generalizations have been extensively utilized in forming neuronal or cell populations [2]. The system is significant traditionally because it utilizes phase plane approaches and mathematical solutions to designate the reactions of neuronal populations to motivations. The general system involves simple integro-differential equations, therefore, limit cycle performance (neural fluctuations) and stimulus-dependent suggested reactions are expected. The key results contain the solvability of multiple stable situations and hysteresis in the population’s reaction.

Coronavirus (COVID-19) has been an infectious virus molded by a recently exposed coronavirus. It has been recorded by the World Health Organization (WHO), it is a pandemic. The first WHO warning of dyed-in-the-wool cases of COVID-19 indicated on January 2020 with 282 cases (see [3, 4]). There is an increasing number of research works that develop the growth of the COVID-19 infection by using an ordinary dynamic system [5] and fractal-fractional dynamic system [6]. Utilizing the recent information from European and African countries, Atangana and Araz presented numerous statistical analyses [7, 8]. Musa et al. [9] introduced a nonlinear 4D-system of ordinary differential equations describing COVID-19. Atangana [10] formulated a numerical design using the Newton polynomial. Other strategies can be located in efforts by Memon et al. [11]. Newly, numerous mathematical simulations have been indicated to realize the coronavirus infection. Supreme of these representations are based on classical integer-order derivative or classical fractional differential operators, which cannot get the vanishing memory and boundary performance found in numerous biological phenomena. Consequently, we investigate the coronavirus disease in this study by discovering the dynamics of COVID-19 infection utilizing the fractional Caputo derivative.

The terminal coronavirus continues to blow out across the globe, and mathematical models can be utilized to display suspected, recovered, and deceased coronavirus patients, as well as how many persons have been tested or even vaccinated. Consequently, mathematical and statistical solutions of the infected human beings overall can decrease the risk of future COVID-19 spread. In this study, we aim to generalize the Wilson–Cowan system (WCS) utilizing the concept of fractional calculus to study the growth of COVID-19 population. This investigation includes a dynamic term, which is the exponential law to discover and realize the graph of the growth. The solvability of the system is indicated by using the optimal point theorem of simulation function. Other behaviors are indicated such as the approximated solvability using the fractional Tutte polynomials.

Preparations

This section deals with some concepts and the properties of these concepts.

ABC-definition

In recent decades, numerous physical issues have been exposed using the fractional calculus. The essential explanations for employing fractional calculus are that various measures, structures, and inequities display capability to remember the past or nonlocal possessions. The basic outlook and appearances of fractional calculus and fractional differential equations are recognized in various reviews. Most researches focus on the derivatives, which include kernels. For instance, the main difference between the Caputo operator, the Caputo–Fabrizio operator [12], and others is that the Caputo operator is communicated by giving a power law, the Caputo–Fabrizio operator is adapted by utilizing an exponential growth act. The Atangana–Baleanu operator is presented by signifying the extended Mittag-Leffler function [13].

Definition 2.1

Let , be the Atangana–Baleanu operator of order μ of a function χ fulfilling where indicates a normalization function, Ξ represents the Mittag-Leffler function. Corresponding to , the ABC integral is realized by

Example 2.1

Consider the function , then the ABC integral becomes

In our study, since we focus on the approximated solutions, we assume that for all . Applications of this calculus in COVID-19 can be located in [14-16].

Approximate point theorem

We deal with the resulting notion of cyclic ∁-condensing operator. Let ð be a measure of noncompactness on a Banach space and A, B be nonempty and convex subsets of (see [17]).

Definition 2.2

A function is called simulation if and only if , , where and for and .

Definition 2.3

Let be an operator. It is f-cyclic condensing (CC) if, for every nonempty, bounded, convex, and closed subset ,

Lemma 2.2

A relatively nonexpansive cyclic f-condensing mapping recognizes a best proximity point.

Results

We let be the accumulation number of infected people, which represents the sum of the number of the standard infected persons and that of the asymptomatic transmission ones : . Take into account that includes people who have previously been diseased. Therefore, there are rate functions joining χ and ϒ. The following integro-differential system indicates the generalization of WCS. We suggest the generalization by using the ABC formula of fractional calculus as follows: where the variables are defined in different intervals: , , , and .

To study the solvability of system (1), we formulate the following assumptions: Here, we introduce our theorem for the solvability of system (1). We define an operator as follows:

All the functions are continuous in such that , , , and χ, ϒ are inside the nonempty, bounded, closed, and convex sets and respectively.

For a sup.norm, we suppose that , , so that . In addition, for all and , we suppose that there occurs a positive constant fulfilling

For any , , there exists a positive function which is upper semi-continuous and achieves and

Theorem 3.1

Consider system (1) satisfying hypotheses (A1)–(A3). Then it has an optimal solution in , whenever

Proof

Consider the operator and . We aim to show that is a cyclic operator. Let , then we get where . By letting , where and , we have Thus, . In the same manner, we conform that, for , this indicates that and hence . We conclude that is cyclic. The above conclusion shows that the set is bounded in ⊺2, and the set is bounded in ⊺1.

Recall that indicates an optimum outcome of system (1) if and only if . Therefore, we have to prove this fact. Next, we aim to prove that is equicontinuous in ⊺2. For ς and , we have which implies that is equicontinuous in ⊺2. In a similar manner, we confirm that is equicontinuous in ⊺1. As a consequence and via the Arzela–Ascoli theorem, we point that the pair is relatively compact. Now, we have to show that is relatively nonexpansive.

For , we inform that But ϵ is an arbitrary constant, thus when , we have the inequality This indicates that is relatively nonexpansive.

We proceed to show that is f-condensing. Assume that is a nonempty, bounded, closed, and convex set such that Thus, we obtain By putting , then we arrive at Hence, the necessary requirements of Lemma 2.2 are achieved. Thus, the operator has the best proximity point and thus system (1) has an optimal solution.

This completes the proof. □

Numerical structures

In this subsection, we introduce some numerical systems pointing to utilize Theorem 3.1. The chief assumption in Theorem 3.1 is . This inequality is very informal to fulfill likening with other existence theorems attaining (A1–A3). Theorem 3.1 designates that the system attractive formula (1) admits an optimal solution. This kind of result is very significant in dynamic and control systems. By this result, one can investigate the stability, oscillatory solution, and other performances of the resolution.

Example 3.2

By using Mathematica 11.2, the solution is realized by the integral where c is a constant and W represents the product log function. As an application of Theorem 3.1, we assume that and , where For instance, , we attain ; thus, by Theorem 3.1, system (3) admits an optimal solution converging to a limit cycle. In another case, suppose that , then . This implies that system (3) admits optimal solution converging to a limit cycle. Similarly, for and . Figure 1 represents various cases considering the value of ρ.

Figure 1

Results of (3) with various values of ρ. From the upper left, the solution is for , while the lower left is the optimal outcome for . The upper and lower right graphs represent the solutions for and respectively. We indicate that the cyclic result for these values is based on the fact that is cyclic

Example 3.3

Consider the following system: where the value of ρ achieves For instance, when , we have an optimal solution with the initial condition . Furthermore, it is unstable cyclic because it indicates a portrait unstable limit cycle (see Fig. 2, the upper graphs). When , the system admits an optimal solution with a portrait unstable limit cycle (see Fig. 2, the lower graphs).

Figure 2

Solutions of (4) for different cases, depending on the value of ρ. We recognize that the cyclic solution for these cases is based on the fact that is cyclic

Remark 3.1

The upper value of ρ plays an important role in describing the behavior of the optimal solution. It represents that when the upper value of ρ is equal to 0.33, then we have a stable solution (see Example 3.2). If the value is greater than 0.33, we have an unstable optimal solution (see Example 3.3). We shall take into account this fact when we establish the connections of Tutte polynomials in the next section.

Approximated solution

In this section, we aim to use a type of polynomial to approximate the solution of (1). In this place, we suggest to use the Tutte polynomial because the formality of this polynomial involves two variables as follows (see [18]): By using the construction of Example 2.1, we have the fractional Tutte polynomial as follows: To determine the upper bound (approximated value) of , based on Theorem 2.2, we shall consider that these weights satisfy the upper bound of ρ, which is given by the formula Note that providing that . This value is approximated with the upper bound of ρ in Example 3.2. By suggesting the solution of system (3) in terms of fractional Tutte polynomials, we have

We suppose that χ and ϒ have the same gathering of roots as the original polynomials. That is, these polynomials can be recognized in a Grobner basis (GB). For linear functions in any number of variables, GB is analogous to Gaussian elimination. For example, if , then (see Fig. 3). For nonlinear cases, we have the following results.

Figure 3

GB for the linear case. From the left: the graph represents 3D solution, contour plot, and integral curve. Similarly, for the case for all real coefficients. The second row is GB of fractional order , which indicates the set

Example 3.4

with root (see Fig. 4).

Figure 4

GB for the nonlinear case. From the left: the graph represents , τς, and respectively. The second row represents their fractional case when

Example 3.5

with the real root (see Fig. 5).

Figure 5

From the left: the graph represents and respectively. The second row represents the fractional case with

Example 3.6

with two real roots and (see Fig. 6).

Figure 6

GB for the nonlinear case with the fractional case when

Example 3.7

with one real root (see Fig. 7).

Figure 7

GB for the nonlinear case in Example 3.7

Application

In this subsection, we shall utilize live data regarding COVID-19, which were recorded in May. Figures 8 and 9 indicate live data which were recognized in May for Brazil and USA. By using the approximated solution, we see that the data converge to the result of different BG values where the value of τ indicates the starting situation. Experimentally, we indicate that the good result appears, when the value of fractional order and the interval of convergence solution is .

Figure 8

Approximate solution by using GB in May, where the data was given by . For Brazil, we select , while for USA is . The fractional order is and respectively.

Figure 9

Approximate solution by using GB in May, where the data was given by . For Brazil, we select , while for USA is . The fractional order is and respectively.

Conclusion

By using the fractional calculus, type ABC, we have generalized WCS. We have got two kinds of solutions, the first one is the optimal solution (see Theorem 3.1) using the concept of simulation function and the second is approximated solution using the GB set of polynomials of two variables . The optimal solution brings the stability, oscillation, and periodicity. The second solution is validated for discrete data. In this investigation, we generalized the most popular graph polynomials called the Tutte polynomial and a variety of carefully related graph polynomials such as the harmonic, movement, reliability, and shelling polynomials. We also used the Tutte polynomial to demonstrate how graph polynomials may be both dedicated and generalized, and how they can convert information relevant to medical applications. We concluded with a brief conversation of computational complexity deliberations. Different studies are presented using fractional calculus together with statistics and polynomials categories, which can be located in [19-25].