# Transmission Dynamics of Fractional Order Brucellosis Model Using Caputo–Fabrizio Operator

Authors: Olumuyiwa James Peter
Journal: International Journal of Differential Equations
 (2020)

---

## Abstract

In this paper, a noninteger order Brucellosis model is developed by employing the Caputo–Fabrizio noninteger order operator. The approach of noninteger order calculus is quite new for such a biological phenomenon. We establish the existence, uniqueness, and stability conditions for the model via the fixed-point theory. The initial approachable approximate solutions are derived for the proposed model by applying the iterative Laplace transform technique. Finally, numerical simulations are conducted for the analytical results to visualize the effect of various parameters that govern the dynamics of infection, and the results are presented using plots.


---

## Body

## 1. Introduction

Brucellosis is a contagious zoonotic infection that is attributable to diverse species of Brucella. Although the disease is well contained in Australia and the UK, the annual global incidence of Brucellosis is estimated above 500 000 [1, 2]. Brucellosis is endemic in the Caribbean, South America, Mediterranean Europe, Africa, Asia, and in the Middle East and poses considerable economic consequences in some of the regions despite the enforcement of firm veterinary hygiene interventions [3]. Typically, animal Brucellosis is attributed to either direct or indirect contacts with the Brucellosis bacteria [4]. While direct contact is from contact with the infected secretions, droppings, or vertical transmission from mother-to-offspring during delivery, indirect contact occurs through exposure to contaminated environment or objects. It may also occur through respiratory or digestive tracks especially during feeding [5]. Human Brucellosis, on the contrary, is triggered by contamination from infectious animals or their products. Orchitis, epididymitis, placentitis, and abortions are the commonest consequences of Brucellosis [6]. The primary hosts for Brucellosis are sheep, goats, cattle, buffalo, swine, horses, dogs, and man [1]. Control strategies for Brucella within domestic animals involve detection, inoculation, and elimination of the affected livestock [7]. However, rapid control of Brucellosis is achieved by the combination of the interventions.The quantification of intricate and complex biological systems in terms of mathematical structures has posed a serious challenge to contemporary scientists worldwide. A good number of models with soliton solutions that are designed in terms of integer-order derivatives had been extensively studied in the last two decades [8–12]. Some of these models in the literature include Tzizeica-type equation [13], Fokas–Lenells equation [14], complex Ginzburg–Landau equation [15], differential-difference equation [16], nonlinear oscillator with discontinuities [17], and Bernstein wavelet equations [18]. The aforementioned models play significant roles in advancing suitable control procedures in a biological system [19].Scientific and technical models of Brucellosis dynamics can be used to interpret the primary outcomes and condense the discoveries to the best advantage of protection at various stages; still, a number of those models in the form of nonfractional derivative were thoroughly investigated [20–23]. Researchers in mathematics have been engaged in the last decade in the formulation of mathematical models of Brucellosis to investigate the stability and dynamical behaviours of the disease [24, 25]. The primary model of these kinds, which examine the indirect transmission mechanism of Brucellosis with the influence of immunology, was jointly investigated by Zhang et al. [26]. In [1], Li et al. proposed a model that quantified the spread of Brucellosis from sheep to man and among sheep. Unlike in [1], Brucellosis models that described the transmission of Brucellosis in the cattle population and that characterized stillbirths in dairy cattle are designed and analyzed by Tumwiine and Robert [27] and Nie et al. [6], respectively.As time goes by, advancing mathematical models using noninteger order became a significant area of study since the evolution-related realities and evidence are represented more effectively in terms of noninteger calculus [28, 29]. The elaborate investigations of such models play an exemplary role through advancing implementations of various real-life problem. A large number of logical and numerical procedures and schemes are proposed to define arbitrary and real solutions with noninteger operators for differential equations have been playing an exceptional task in the advancement of implementations of numerous problems. A good number of logical and numerical procedures and schemes are proposed to figure out indefinite and actual solutions of the differential equations with noninteger operators [30–35]. However, up till now, there is no work which has been designed to analyze the dynamics of Brucellosis with the fractional-order calculus.There are four compartments into which the Brucellosis model in [26] is separated, denoted as susceptible dairy cows S(t), the exposed dairy cow E(t), infectious dairy cows I(t), and the Brucella contaminated environment B(t). The sum of the population at any time t is represented by N = S + E + I + B. In [26], Zhang et al. proposed a Brucellosis model and analyzed the qualitative characteristics of the suggested model which is denoted by(1)S=θ+pS+E+I−αSI−dS−βBS,E=αSI+βBS−dE−δE,I=δE−dI−aI,B=nE+I−ω+chB,where θ denotes the yearly recruitment rate of dairy cows, the yearly birth rate of cow is p, the annual culling rate is d, the treatment outcome rate is δ, the rate of yearly elimination for the positive dairy cows is a, the yearly number of Brucella is n, the yearly rate of natural death of Brucella is ω, the rate of sterilization in a disinfection is c, and the yearly number of disinfections is h, and α is the indirect transmission rate and β is the pathogen density-dependent component.Motivated by the previously mentioned literature, we have applied Caputo–Fabrizio fractional orderϕ, where ϕ∈0,1 to explore dynamics of the Brucellosis model discussed in Zhang et al. [26]:(2)DtϕCFSt=θ+pS+E+I−αSI−dS−βBS,DtϕCFEt=αSI+βBS−dE−δE,DtϕCFIt=δE−dI−aI,DtϕCFBt=nE+I−ω+chB,with initial conditions S0=b1,E0=b2,I0=b3,andB0=b4.The typical model in [26] is derived for ϕ=1. As the Caputo–Fabrizio fractional order quantifies the previously mentioned phenomenon accurately. The numerical results are obtained by employing the Laplace transform method. Arbitrary values are allocated to the initial conditions and parameters to validate our results.The remainder of this script is arranged as follows. In Section2, some fundamental remarks and definitions associated with noninteger calculus are stated. In Section 3, existence and uniqueness conditions are derived by using the fixed-point theory. In Section 4, stability results are obtained via iterative Laplace transform technique. In Section 5, the proposed method is employed to noninteger order Brucellosis model and the results of the simulations are illustrated graphically. Section 6 is devoted to the discussions of numerical results. Finally, the conclusion of the work is stated in Section 7.


## 2. Preliminaries

Definition 1.
(see [36]). Let u1∈H′0,a,a>0,0<ϕ<1, then the Caputo fractional operator is defined as(3)DtϕCFu1t=2−ϕMϕ21−ϕ∫0texp−ϕt−S1−ϕu1′tdt,where Mϕ is the normalization function which depends on ϕ, such that Mϕ=M1=1.Definition 2.
(see [37]). The Caputo–Fabrizio fractional integral operator of order 0<ϕ<1 is defined as(4)JtϕCFu1t=2−ϕMϕ21−ϕu1t+2ϕ2−ϕMϕ∫0tu1tdt.Definition 3.
(see [36]). The Laplace transformation for the Caputo fractional operator of order 0<ϕ≤1 and n∈M is expressed as(5)LJtn+ϕCFu1tρ=11−ϕLun+1Lexp−ϕ1−ϕt=ρn+1Lu1−ρnu10−ρn−1u1′0…u1n0ρ+ϕ1−ρ.
Generally, we have(6)LDtϕCFutρ=ρLu1tρ+ϕ1−ρ,n=0,LDtϕ+1CFu1tρ=ρ2Lu1t−ρu10−u1′0ρ+ϕ1−ρ,n=1.


## 3. Existence and Uniqueness

This section deals with the existence and uniqueness of the fractional Brucellosis model using fixed-point theory [38–45]. By applying the fractional integral operator (2) into (4), we have(7)St=S0+21−ϕ2ϕMϕθ+pSt+Et+It−αStIt−dSt−βBtSt+21−ϕ2ϕMϕ∫0tθ−pSε+Eε+Iε−αSεIt−dSε−dSε−βBεSε,Et=E0+21−ϕ2ϕMϕαStIt+βBtSt−dEt−δEt+21−ϕ2ϕMϕ∫0tαSεIε+βBεSε−dEε−δEεdε,It=I0+21−ϕ2ϕMϕδEt−dIt−aIt+21−ϕ2ϕMϕ∫0tδEε−dIε−aIεdε,Bt=B0+21−ϕ2ϕMϕnEt+It−ω+cBt+21−ϕ2ϕMϕ∫0tnEε+Iε−ω+chBεdε.For clarity, we consider the following kernels:(8)ϕ1=θ+pS+E+I−αSI−dS−βBS,ϕ2=αSI+βBS−dE−δE,ϕ3=δE−dI−aI,ϕ4=nE+I−ω+chB.Theorem 1.
The kernelsϕ1, ϕ2, ϕ3, and ϕ4 in (8) satisfy the Lipschitz condition, provided the following inequality holds:(9)0<Y1,Y2,Y3,Y4<1.Proof.
LetS1 and S2 for the kernel given ϕ1, E1 and E2 for the kernel ϕ2, I1 and I2 for the kernel ϕ3, and B1 and B2 for the kernel ϕ4 be the respective functions associated with the following:(10)ϕ1t,S1t−ϕ1t,S2t=−p+αI+d+βBS1t−S2t,ϕ2t,E1t−ϕ2t,E2t=−d+δE1t−E2t,ϕ3t,I1t−ϕ1t,I2t=−d+aI1t−I2t,ϕ4t,B1t−ϕ4t,B2t=−ω+chB1t−B2t.
By applying the properties of the norm on (10), respectively, we have(11)ϕ1t,S1t−ϕ1t,S2t=−p+αI+d+βBS1t−S2t≤αI+βBS1t−S2t≤a1S1t−S2t,where a1=αi+βb,i=maxt∈ΩIt,b=maxt∈ΩBt, are bounded functions. Subsequently, we have(12)ϕ2t,E1t−ϕ1t,E2t=≤a2E1t−E2t,ϕ3t,I1t−ϕ3t,I2t=≤a3I1t−I2t,ϕ4t,B1t−ϕ4t,B2t=≤a4B1t−B2t,where a2=d+δ,a3=d+a,anda4=ω+ch.
We construct the following recursive formula:(13)Snt=21−ϕ2−ϕMϕϕ1t,Sn−1t+2ϕ2−ϕMϕ∫0tϕ1ε,Sn−1εdε,Ent=21−ϕ2−ϕMϕϕ2t,En−1t+2ϕ2−ϕMϕ∫0tϕ1ε,En−1εdε,Int=21−ϕ2−ϕMϕϕ3t,In−1t+2ϕ2−ϕMϕ∫0tϕ1ε,In−1εdε,Bnt=21−ϕ2−ϕMϕϕ4t,Bn−1t+2ϕ2−ϕMϕ∫0tϕ1ε,Bn−1εdε.
Furthermore, by applying triangular inequality, we have(14)Dnt=Snt−Sn−1t≤21−ϕ2−ϕMϕϕ1t,Sn−1t−ϕ1t,Sn−2t+2ϕ2−ϕMϕ∫0tϕ1ε,Sn−1ε−ϕ1ε,Sn−2εdε,Fnt=Ent−En−1t≤21−ϕ2−ϕMϕϕ2t,En−1t−ϕ2t,En−2t+2ϕ2−ϕMϕ∫0tϕ2ε,En−1ε−ϕ2ε,En−2εdε,Int=Int−In−1t≤21−ϕ2−ϕMϕϕ3t,In−1t−ϕ3t,In−2t+2ϕ2−ϕMϕ∫0tϕ3ε,In−1ε−ϕ1ε,In−2εdε,Bnt=Bnt−Bn−1t≤21−ϕ2−ϕMϕϕ4t,Bn−1t−ϕ4t,Bn−2t+2ϕ2−ϕMϕ∫0tϕ4ε,Bn−1ε−ϕ4ε,Bn−2εdε.
Snt=∑m=0∞Dmt, Int=∑m=0∞Fmt, Int=∑m=0∞Gmt, and Bnt=∑m=0∞Hmt.
The kernelsϕ1, ϕ2, ϕ3, and ϕ4 which satisfy the Lipschitz condition are therefore obtained:(15)Dnt=Snt−Sn−1t≤21−ϕ2−ϕMϕa1Sn−1t−Sn−2t+2ϕ2−ϕMϕa1∫0tSn−1ε−Sn−2εdε,Fnt=Ent−En−1t≤21−ϕ2−ϕMϕa2En−1t−En−2t+2ϕ2−ϕMϕa2∫0tEn−1ε−En−2εdε,Gnt=Int−In−1t≤21−ϕ2−ϕMϕa3In−1t−In−2t+2ϕ2−ϕMϕa3∫0tIn−1ε−In−2εdε,Hnt=Bnt−Bn−1t≤21−ϕ2−ϕMϕa4Bn−1t−Bn−2t+2ϕ2−ϕMϕa4∫0tBn−1ε−Bn−2εdε.
This proves the result.Theorem 2.
TheSEIB Brucellosis model involving the fractional operator in (2) has a solution (existence of the solution).Proof.
Following the result in (15) and with respect to the recursive expression, we get the following:(16)Dnt≤S0+21−ϕ2−ϕMϕa1n+2ϕ2−ϕMϕa1tn,Fnt≤E0+21−ϕ2−ϕMϕa2n+2ϕ2−ϕMϕa2tn,Gnt≤I0+21−ϕ2−ϕMϕa3n+2ϕ2−ϕMϕa3tn,Hnt≤B0+21−ϕ2−ϕMϕa4n+2ϕ2−ϕMϕa4tn.
Hence, (16) exists, additionally, we show that (16) is the solution of (2) with the following assumption:(17)St=Snt−Δ1nt,Et=Ent−Δ2nt,It=Int−Δ3nt,Bt=Bnt−Δ4nt,where the remaining terms of the solution are Δ1n,Δ2n,Δ3n, and Δ4n. Hence, we have(18)Snt−Sn−1t=21−ϕ2−ϕMϕϕ1t,St−Δ1nt+2ϕ2−ϕMϕ∫0tϕ1ε,Sε−Δ1nεdε,Ent−En−1t=21−ϕ2−ϕMϕϕ2t,Et−Δ2nt+2ϕ2−ϕMϕ∫0tϕ2ε,Eε−Δ2nεdε,Int−In−1t=21−ϕ2−ϕMϕϕ3t,It−Δ3nt+2ϕ2−ϕMϕ∫0tϕ3ε,Iε−Δ3nεdε,Bnt−Bn−1t=21−ϕ2−ϕMϕϕ4t,Bt−Δ4nt+2ϕ2−ϕMϕ∫0tϕ4ε,Bε−Δ4nεdε.
We make use of norm and Lipschitz condition to obtain(19)St−S0−21−ϕ2−ϕMϕϕ1t,St−Δ1nt−2ϕ2−ϕMϕ∫0tϕ1ε,Sε−Δ1nεdε≤Δ1n1+21−ϕ2−ϕMϕa1+2ϕ2−ϕMϕa1t.
Usinglimitn⟶∞ in (19), we have(20)Δ1nt⟶0.
Therefore,(21)St=S0+21−ϕ2−ϕMϕϕ1t,St+2ϕ2−ϕMϕ∫0tϕ1ε,Sεdε.
In the same way, asn⟶∞, we obtain the following.
Using the same approach, in (21), we have(22)Δ2nt⟶0,Δ3nt⟶0,Δ4nt⟶0.
Similarly, from (21), there exists a solution of (2).Theorem 3.
TheSEIB Brucellosis model involving the fractional operator in (2) has a unique solution.Proof.
Proof We adopt another solution (2) say S∗,E∗,I∗, and B∗; then, we have(23)St−S∗t=21−ϕ2−ϕMϕϕ1t,St−ϕ1t,S∗t+2ϕ2−ϕMϕ∫0tϕ1ε,Sε−ϕ1t,S∗εdε,Et−E∗t=21−ϕ2−ϕMϕϕ2t,Et−ϕ2t,E∗t+2ϕ2−ϕMϕ∫0tϕ2ε,Eε−ϕ2t,E∗εdε,It−I∗t=21−ϕ2−ϕMϕϕ3t,It−ϕ3t,I∗t+2ϕ2−ϕMϕ∫0tϕ3ε,Iε−ϕ3t,I∗εdε,Bt−B∗t=21−ϕ2−ϕMϕϕ4t,Bt−ϕ4t,B∗t+2ϕ2−ϕMϕ∫0tϕ4ε,Bε−ϕ4t,B∗εdε.
By applying the norm on (23),(24)St−S∗t≤21−ϕ2−ϕMϕϕ1t,St−ϕ1t,S∗t+2ϕ2−ϕMϕ∫0ttϕ1ε,Sε−ϕ1t,S∗εdε,Et−E∗t≤21−ϕ2−ϕMϕϕ2t,Et−ϕ2t,E∗t+2ϕ2−ϕMϕ∫0tϕ2ε,Eε−ϕ2t,E∗εdε,It−I∗t≤21−ϕ2−ϕMϕϕ3t,It−ϕ3t,I∗t+2ϕ2−ϕMϕ∫0tϕ3ε,Iε−ϕ3t,I∗εdε,Bt−B∗t≤21−ϕ2−ϕMϕϕ4t,Bt−ϕ4t,B∗t+2ϕ2−ϕMϕ∫0tϕ4ε,Bε−ϕ4t,B∗εdε.
By considering Theorems1 and 2, we obtain the following results:(25)St−S∗t≤a121−ϕ2−ϕMϕSt−S∗t+a12ϕ2−ϕMϕtSt−S∗t,Et−E∗t≤a221−ϕ2−ϕMϕEt−E∗t+a22ϕ2−ϕMϕtEt−E∗t,It−I∗t≤a321−ϕ2−ϕMϕIt−I∗t+a32ϕ2−ϕMϕtIt−I∗t,Bt−B∗t≤a421−ϕ2−ϕMϕBt−B∗t+a42ϕ2−ϕMϕtBt−B∗t.
The solution functions in (19) satisfy the following inequalities:(26)St−S∗t≤1−2a12−ϕMϕ1−ϕ−tϕ≤0,Et−E∗t≤1−2a22−ϕMϕ1−ϕ−tϕ≤0,It−I∗t≤1−2a32−ϕMϕ1−ϕ−tϕ≤0,Bt−B∗t≤1−2a42−ϕMϕ1−ϕ−tϕ≤0.
Going by the last equation, we conclude that(27)St=S∗t,Et=E∗t,It=I∗t,Bt=B∗t.


## 4. Stability

This section deals with the application of iterative Laplace transformation technique on the fractional Brucellosis model, obtaining the stability condition for the approximate solution.


### 4.1. Iterative Laplace Transformation Technique

We consider the Brucellosis model in (2) with the initial conditions in (3). Applying the Laplace transformation to both sides of equation (2), we have(28)ρLSt−S0ρ+ϕ1−ρ=Lθ−pS+E+I−αSI−dS−βBS,ρLEt−E0ρ+ϕ1−ρ=LαSI+βBS−dE−δE,ρLIt−I0ρ+ϕ1−ρ=LδE−dI−aI,ρLBt−B0ρ+ϕ1−ρ=LnE+I−ω+chB.By rearranging the following, we have(29)LSt=S0ρ+ρ+ϕ1−ρρLθ−pS+E+I−αSI−dS−βBS,LEt=E0ρ+ρ+ϕ1−ρρLαSI+βBS−dE−δE,LIt=I0ρ+ρ+ϕ1−ρρLδE−dI−aI,LBt=B0ρ+ρ+ϕ1−ρρLnE+I−ω+chB.The inverse Laplace transform of (29) gives(30)St=S0+L−1ρ+ϕ1−ρρLθ−pS+E+I−αSI−dS−βBS,Et=E0+L−1ρ+ϕ1−ρρLαSI+βBS−dE−δE,It=I0+L−1ρ+ϕ1−ρρLδE−dI−aI,Bt=B0+L−1ρ+ϕ1−ρρLnE+1−ω+chB.We obtain the following infinite series solution for theSEIB Brucellosis model in (2):(31)St=S0+S1+S2+S3+⋯,Et=E0+E1+E2+E3+⋯,It=I0+I1+I2+I3+⋯,Bt=B0+B1+B2+B3+⋯.


## 4.1. Iterative Laplace Transformation Technique

We consider the Brucellosis model in (2) with the initial conditions in (3). Applying the Laplace transformation to both sides of equation (2), we have(28)ρLSt−S0ρ+ϕ1−ρ=Lθ−pS+E+I−αSI−dS−βBS,ρLEt−E0ρ+ϕ1−ρ=LαSI+βBS−dE−δE,ρLIt−I0ρ+ϕ1−ρ=LδE−dI−aI,ρLBt−B0ρ+ϕ1−ρ=LnE+I−ω+chB.By rearranging the following, we have(29)LSt=S0ρ+ρ+ϕ1−ρρLθ−pS+E+I−αSI−dS−βBS,LEt=E0ρ+ρ+ϕ1−ρρLαSI+βBS−dE−δE,LIt=I0ρ+ρ+ϕ1−ρρLδE−dI−aI,LBt=B0ρ+ρ+ϕ1−ρρLnE+I−ω+chB.The inverse Laplace transform of (29) gives(30)St=S0+L−1ρ+ϕ1−ρρLθ−pS+E+I−αSI−dS−βBS,Et=E0+L−1ρ+ϕ1−ρρLαSI+βBS−dE−δE,It=I0+L−1ρ+ϕ1−ρρLδE−dI−aI,Bt=B0+L−1ρ+ϕ1−ρρLnE+1−ω+chB.We obtain the following infinite series solution for theSEIB Brucellosis model in (2):(31)St=S0+S1+S2+S3+⋯,Et=E0+E1+E2+E3+⋯,It=I0+I1+I2+I3+⋯,Bt=B0+B1+B2+B3+⋯.


## 5. Description

In this section, we analyze the numerical simulations of the Caputo–Fabrizio Brucellosis model (2) with the following initial conditions S0=a1=60,E0=a2=40,I0=a3=20, and B0=a4=10 for specific fractional-order value ϕ∈0,1. The parameters used and their values are given in Table 1.Table 1
Details of parameters [26, 46].
ParameterValueDescriptionω3.0Annual natural mortality of Brucellad2.0Annual natural elimination ratea0.8Annual elimination rate for the positive cowsδ4.0Clinical outcome ratep1.0The annual birth rate of dairy cowsα0.2Indirect transmission rateθ2The annual introduction number of dairy cowsn0.2The annual quantity of Brucellac0.52The rate of sterilization in a disinfectionh0.5The yearly number of disinfections


## 6. Numerical Results and Discussion

Figures1–4 display the numerical results for the susceptible cows S(t), exposed cows E(t), infectious cows I(t), and the Brucella contaminated environment B(t) obtained for different values of (ϕ = 1, 0.9, 0.8, 0.7) with different initial conditions S0=a1=60,E0=a2=40,I0=a3=20, and B0=a4=10 by utilizing the iterative Laplace transform technique of noninteger order Brucellosis model which confirms that the iterative Laplace transform method is promising and can predict the behaviour of the variables accurately within the region of consideration. The integer order has a little impact on the dynamics of the Brucellosis model. This negligible effect has been taken care of in our proposed fractional-order model and our simulations have shown that a change in the values of the order influences the transmission dynamics of the Brucellosis model. Figure 1 shows that susceptible cows S(t) increase steadily in the first six months when ϕ = 0.8 but falls afterwards. However, the population of the susceptible cows falls continuously at other various fractional values of ϕ. The plot in Figure 2 of exposed cows E(t) demonstrates that the population of the exposed cow rises rapidly in the first three months when ϕ = 0.8 but drops afterwards. The population of the exposed cows also falls continuously for other fractional values of ϕ in Figure 2. The populations of the infectious cows I(t) and the Brucella contaminated environment B(t) depicted in Figures 3 and 4, respectively, also follow similar pattern and behaviour as in Figures 1 and 2. It can be inferred that, as the fractional-order rises from ϕ = 0.5 to ϕ = 1, the population of cows infected with Brucellosis asymptotically tends to zero.Figure 1
Approximate solution of susceptible cows.Figure 2
Approximate solution of exposed cows.Figure 3
Approximate solution of infectious cows.Figure 4
Approximate solution of concentration of bacteria in the environment.Figures5–8 shows different behaviours of the population of susceptible cows S(t), exposed cows E(t), infectious cows I(t), and the Brucella contaminated environment B(t), where the required simulations have been carried out for t = 12 months while we set the value of the fractional-order parameter ϕ=1. Furthermore, it is observed that the numerical outcome of the simulation is a function of the noninteger derivative ϕ. This shows that the special noninteger operator such as the Caputo–Fabrizio operator is efficient and plays significant roles in inaccurate predictions. Besides, the hybrid feature of Caputo–Fabrizio is strong to capture the complex dynamics and offer reliable predictions.Figure 5
Approximate solution of susceptible cows forϕ=1.Figure 6
Approximate solution of exposed cows forϕ=1.Figure 7
Approximate solution of infectious cows forϕ=1.Figure 8
Approximate solution of concentration of bacteria forϕ=1.


## 7. Conclusion

In this paper, we have analyzed the Caputo–Fabrizio noninteger order Brucellosis model. Using fixed-point theorem, the steady results for the existence and uniqueness of solutions and the stability of the proposed noninteger order Brucellosis model have been derived. The iterative Laplace transform technique is applied to conduct the simulations by using a set of values whose sources are from the literature for the model parameters. The effects of various values for the order of the Caputo–Fabrizio noninteger order derivative Brucellosis model on the transmission dynamics of the disease are revealed by the results of the simulations. The results of the simulations, therefore, confirm the validity and efficiency of the Caputo–Fabrizio noninteger derivative in quantifying dynamics of Brucellosis disease and other related kind of problems.



---

*Source: 2791380-2020-08-14.xml*