Abstract

This paper is devoted to the model of Lassa hemorrhagic fever (LHF) disease in pregnant women. This disease is a biocidal fever and epidemic. LHF disease in pregnant women has negative impacts that were initially appeared in Africa. In the present study, we find an approximate solution to the fractional-order model that describes the fatal LHF disease. Laplace transforms coupled with the Adomian decomposition method (ADM) are applied. In addition, the fractional-order LHF model is numerically simulated in terms of a varied fractional order. Furthermore, a fractional order optimal control for the LHF model is studied.

1. Introduction

Infectious diseases remain a major threat to human health and welfare throughout history. Their spread is affected by several factors such as mode of transmission, infectious agent, infectious periods, incubation period, susceptibility, and resistance [1]. Some devastating infectious diseases like Lassa fever are endemic in many parts of the world and continue to emerge.

Lassa fever, also known as Lassa hemorrhagic fever (LHF), is an animal-borne, or zoonotic, acute viral hemorrhagic illness caused by Lassa virus, a member of the family Arenaviridae. The reservoir is the rats Mastomys natalensis, M. erythroleucus, and Hylomyscus pamfi. Food contaminated with rodent urine, saliva, or feces is the main cause of human cases. Human-to-human transmission can occur via exposure to the blood, feces, urine, saliva, or vomitus of infected patients [2, 3]. The incubation period is ca. 2-21 days. LHF is basically endemic in some countries of West Africa, and outbreaks have occurred in Benin, Ghana, Guinea, Liberia, Nigeria, Sierra Leone, and Togo. Cases have been imported to Germany, United Kingdom, United States, and Sweden. This disease with high mortality rates (e.g., 80% among pregnant women) is likely affecting between 100,000 and 300,000 people every year, and it kills around 5,000 people, almost all in West Africa alone [4, 5].

Fractional Calculus (FC) is a fruitful field of mathematical research with numerous applications in engineering [6], nanotechnology [7], optics [8], human diseases [9], and chaos soliton theory [10]. Application of FC is also adopted in biology, heat transfer, system identification, genetic algorithms, traffic systems, telecommunications, physics as well as finance, and economics [1120].

Herein, we are interested in the paradigm showing the nature of LHF, since it is lethal and transmissible [1, 2, 21, 22]. Using mathematical modeling to probe this fatal illness may provide insight, answers, and guidance useful for controlling its spread, and the originating applications. The given model is studied through Caputo fractional operator. In this model, the Caputo fractional operator is used since it permits both initial and boundary conditions to involve in modeling problems. The Caputo fractional operator can be applied on sufficiently differentiable functions only. One of the advantages of the fractional derivative of Caputo compared with the fractional derivative 3 of Riemann-Liouville is that the fractional derivative of Caputo for the constant is equal to zero. We develop an approximate solution to the fractional-order model describing the LHF affecting pregnant women using the Laplace transform coupled with the ADM. Numerical simulations are also presented. A fractional optimal control for the SIRD model is ultimately offered. For a good survey on the approximating of the differential equations, see [2326]. Some recent research fractional optimal controls can be found in [2729]. Atangana [21] studied the LHF using beta differential operator, where the model was proposed as follows in (1). Table 1 shows the meaning of the parameters in the model (1). where .

Also, the model has been studied by Goyal et al. [22]. The -homotopy analysis transform method (q-HATM) has been applied for solving the LHF model [30]. The LHF model via Atangana-Baleanu fractional derivative has been studied in [31]. The SIR model has been studied through Caputo-Fabrizio fractional operator in [32].

In this paper, we reconsider the system (1) by the following fractional sense of the system (1) where is the Caputo fractional derivative of order

We organize the rest of our paper as follows. In Section 2, we provide some basic mathematical concepts that will be necessary for our work. Also, the description of the proposed model is presented. We then follow by presenting the series solution for the proposed model in Section 3. We provide the approximate solution construction in Section 4. The numerical solution of the proposed model, having arbitrary order, is given in Section 5. The optimal control case is studied in Section 6. In Section 7, a concise conclusion of our paper results is presented.

2. Preliminaries

For studying the proposed model (2), we present in this section the basic definitions concerned with fractional calculus. For a good survey on the basic definitions and several applications on fractional calculus, see [6, 7, 33].

2.1. Some Important Definitions

Definition 1. The Riemann–Liouville integral is represented by where is the known gamma function and is the order of this fractional integral.

Definition 2. The Caputo-type fractional derivative of order is defined as

Definition 3. The Laplace transform of Caputo fractional derivative of order [33] is given by

2.2. Mathematical Formulation and Description for the Proposed Model

The meaning of the parameters is shown in Table 1. Now, the description and the mathematical formulation of the model can be shown as follows.

With time we can express the rate of change in susceptible population as follows

With time we can express the rate of change in the population of the infected pregnant women as follows

Here, is the number of pregnant women excluded from the group of susceptible.

With time the rate of change in recovery population is written as

Finally, the rate of change of dying population with time is given as

Hence, the mathematical model for the LHF is given by model (2), with initial conditions

Also Here, and are sufficiently differentiable functions. Model (2) has not yet been studied by using the Laplace transform coupled with the ADM. Also, we introduce numerical simulations and fractional optimal controls for the SIRD model.

The interaction between different compartments of LHF model (2) can be viewed via the following signal flow graph shown in Figure 1.

2.3. Calculating the Basic Reproduction Number

Basic reproduction number is the expected production of newly infected females, in a totally susceptible women, by meeting a typical infective woman. The local stability and instability of the disease-free equilibrium (DFE) depend on the value of . Also, it defines the local stability threshold for the DFE. Moreover, it is a very important tool for controlling the disease, and it is an essential epidemiological criteria of disease. When the DFE is locally asymptotically stable; a small number of infections into the population may cause it to evolve into an endemic prevalence. From another point of view, when the DFE is locally unstable; a sufficiently small number of infected women will produce an outbreak. Here, is estimated from the paradigm of nonlinear FDE’s (2) via the next-generation matrix method [34].

Since the Jacobian matrix is

At the disease-free equilibrium point and then by putting the model in the form [34] where

then

Then, the basic reproduction number of the given paradigm is . If then, infection cannot develop because an infected woman infects less than one new woman during their infectious period on average. Conversely, if the disease can enter the population because on average each infected woman will infect more than one new woman.

3. Series Solution for the Proposed Model

Applying Laplace transformation of (2) gives which gives where .

By dividing the th equation of (18) by , we get

Taking the inverse of Laplace transformation on (19) gives

Now, we will find the series solution by putting

Decomposing the terms gives

The polynomial is known as “Adomian polynomial” with the form

By substituting from (21) and (22) into (20), we get

Note that in (29).

4. Approximate Solution

4.1. Construction of the Approximate Solution

Hereafter, we will discuss an approximate solution for the proposed model. Taking the inverse of Laplace transformation for (24), (25), (26), and (27) gives

Since then, if we take the first four terms from the above series and take into consideration the equations from (30) to (46), we get

After the compensation by the values of we obtain

Now, by taking the first three terms of , , , and and putting , , , , , , , , = , in the proposed model yield (see [30]).

Figure 2 exhibits the dynamical behavior of and against various fractional orders ().

4.2. Convergence and Error Analysis

It is clear that the previous solution is a series. This solution converges rapidly and uniformly to the exact solution. We will use classical techniques to prove the convergence of the series (21), and hence, we can get the sufficient conditions for the convergence of the used method.

Theorem 4. Suppose is the Banach space and is a contractive nonlinear operator with for all Then, based on Banach contraction principle, there is a unique point such that where By applying the ADM, the series in (21) can be written as follows: and suppose that where then, we obtain

Proof. For the first part by the mathematical induction if we obtain Suppose the result is true at hence, We have This shows that .☐

For the second part (b) .

Since and So, we obtain

5. Fractional-Order Numerical Simulations of the SIRD System

First, we recall the basics of the applied numerical technique that have been given for the numerical simulation of fractional IVPs with Caputo derivatives (2). The technique is an extension of the well-known Adams-Bashforth-Moulton (ABM) integrator which is common for the numerical simulation of 1st order differential equations [35]. This technique comes from the idea that the IVP is equal to the equation of Volterra integral. The fractional ABM technique is fully presented by the following formulas (all other states can be found same as ). Assuming that is the domain of the solution and , , : where

A complete study of the stability of the method has been produced in [36].

In this part, we have found the numerical simulation of the proposed paradigm via the predictor-corrector PECE method of ABM [35] coded in MATLAB. The behavior of the obtained solution is shown by taking suitable values of the parameters as in [30]: , , , , , , , , = , and . In Figure 3, the behavior of and populations versus time with different fractional derivative orders is shown. In Figure 4, the behavior of and populations versus time with different fractional derivative orders is shown. In Figure 5, the dynamics of the populations in and 3D space with different fractional derivative orders are shown. In Figure 6, the dynamics of the populations in and 3D-space with different fractional derivative orders are shown. In Figure 7, the dynamics of the populations in , and 3D space with different fractional derivative orders are shown. In Figure 8, the dynamics of the populations in , and 3D space with different fractional derivative orders are shown. In Figure 9, the Lyapunov exponents concerning the fraction derivative order are displayed. From Figure 9, we can conclude that the system is unstable in the transient period since LE1 is positive. But as shown LE1 loses its positivity after a short period which indicates that the system settles to its fixed point. In Figure 10, we show the 3D dynamics of the infected population concerning both time and infection rate at fractional derivative order 0.95. In Figure 11, we show the 3D dynamics of the susceptible population concerning both time and infection rate at fractional derivative order 0.95. There is a quick convergence between the obtained and the exact solution. In complex models, it is well known that a small variation in physical behavior stimulates several new results to understand and analyze nature in a systematic manner.

6. Fractional Optimal Control (FOCP) for SIRD Paradigm

In this position, we discuss the optimal control of the fractional-order SIRD paradigm:

According to references [3740], consider the model (70) inside The proposed control functions can be formulated as follows: and refers to the end of the time period, provides the rate of susceptibility per year, and provides the rate of curing from infection population to recovered population per year. with is the control functions. Here, the function space norm is defined as . This concept of the norm considered to be equal to taking the limit as of the norm.

The following is the definition of the objective function (the control variables are squared)

Such that , and clarify, respectively, the quantity of infectious, a per year rate of susceptible, and a per year rate of progress from infective to recovered.

Finding the optimal controls is the basic goal in FOCPs, in order to minimize the predefined objective function: where subjected to the constraint where

The initial conditions are as follows:

To construct the FOCP, consider the selected objective function defined in [4144]:

Such that

In equation (77), (78), we write the Hamiltonian of the cost (objective) function (76) and the control of the fractional-order SIRD paradigm (70): then

From (76) and (78), following the method in [37, 39, 4345], the FOPC’s necessary and sufficient conditions are formulated as additionally, where represent the Lagrange multipliers. In equations (79) and (80), the necessary conditions using FOPC Hamiltonian are represented.

The main result is recorded in Theorem 5.

Theorem 5. If are optimal controls related to the controlled states , then, the adjoint variables exist that satisfy the following conditions: (i)Adjoint conditions (costate conditions)Considering the related conditions represented in the recorded theorem and using equations (82), see [37, 39, 43, 44], we obtain the following four conditions that can be formulated as (ii)Transversality conditions (iii)Conditions of optimality moreover, by applying equations (80), the control functions are given as follows:

Proof. The system of costate conditions (83)–(86) is produced from equation (81) such that the Hamiltonian is presented as Furthermore, the conditions in equation (82) are guaranteed, and the optimal control given by equations (91)–(92) can be produced via equation (80).☐

Substituting in (70), the state system is formulated as follows:

For more details about FOC, see [37, 39, 40, 44].

The existence of the control system (83)–(86) is to be proved here, the same can be found in [3743] as follows:

Let

Let .

We have, at , , , , where , and are constants greater than zero. Hence, each of the four functions , and agrees with the Lipschitz conditions. Then, each of the four kernels is continuous respecting the four features, then. See for instance [3743] for more information about the existence and uniqueness.

7. Conclusion

In this work, we have provided the approximate solution constructions for the fractional-order LHF model. In addition, the numerical simulations of the fractional-order LHF model have been experimented with many arbitrary orders. Also, an optimal control case has been investigated. The outcomes acquired by utilizing Laplace transform combined with the ADM are very compatible with results existing in the previous research. Further, the applied method for finding the numerical solution of the proposed model is introduced without considering any perturbation, discretization, or transformations. The current examination enlightens the considered nonlinear manner that rely upon the time history and the time moment by applying the idea of fractional order differentiation. Moreover, the obtained results show that Laplace transform combined with the ADM is amazingly deliberate, more powerful, and exceptionally precise, and which can be applied to break down the different classes of nonlinear marvels that emerge in science and innovation. For future work, we suggest generalizing the studied LHF model to show the effect of the infection from other sources as animals, men, and children.

Data Availability

The data used to support the findings of this study are available from the corresponding author on request.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Conceptualization was done by M.H.; methodology was done by M.H, A.E., and A.M.; software was done by M.H, A.E., and A.M.; validation was done by M.H, A.E., S.U., A.G., and A.M.; formal analysis was done by M.H, A.E., S.U., A.G., and A.M.; investigation was done by M.H, A.E., S.U., A.G., and A.M.; resources was done by M.H, A.E., S.U., A.G., and A.M.; data curation was done by M.H, A.E., S.U., A.G., and A.M.; writing—original draft preparation was done by M.H. and A.E.; writing—review and editing was done by M.H. and A.E.; visualization was done by M.H. and A.E.; supervision was done by M.H.; project administration was done by M.H.; funding acquisition was done by A.G. and S.U. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

The authors acknowledge the support and funding of Research Center for Advanced Material Science (RCAMS) at King Khalid University through Grant No. RCAMS/KKU/008-21.