Abstract

Khat is a green leaf and greenish plant where its branches and leaves are chewed to discharge liquid having active chemicals that change the user’s mood. The purpose of this article is to develop and analyze a mathematical model that can be used to understand the dynamics of chewing Khat. The proposed model monitors the dynamics of five compartments, namely, a group of people who do not chew Khat, designated as ; a group of people who are surrounded by Khat chewers but do not chew at present and may chew Khat in the future, denoted this as ; a group of people who chew Khat, which is represented in ; a group of people contains individuals who consumed Khat quite temporarily for social, spiritual, and recreational purposes, and we describe this group in ; and a group of people those who constantly chew Khat, and they are denoted by . We determined the Khat chewing generation number using the next-generation matrix method, and we have examined the biological meaningfulness, mathematical wellposedness, and stability of both Khat chewing-free and Khat chewing-present equilibrium points of the model analytically. Numerical simulations were presented by solving our dynamical system using Matlabode45 to check the analytical results by considering parameter estimations. The results of this study show that, for , the Khat chewing-free equilibrium point is stable, and it is unstable for , and the Khat chewing-present equilibrium point is stable if , and it is unstable if . The stability of both equilibrium points implies that, for a high rate of conversion from non-Khat chewer to exposed groups , the inflow of an insignificant number of Khat chewers to the community produces a significant number of Khat chewers , and if the return back from Khat chewing to the exposed group because of socio-economic, environmental, and religious influences grows exponentially, the inflow of an insignificant number of Khat chewers to the community produces an insignificant number of Khat chewers. It is found that increasing the rate of conversion from non-Khat chewer to exposed groups makes the disease eradication more challenging. We, therefore, strongly urge religious leaders, social committee leaders, elders, and health experts to teach their followers to reduce their Khat-chewing habits.

1. Introduction

Khat has different names in different countries. For example, chat in Ethiopia, Qat in Yemen [1], Miraa in Kenya, and Jaad in Somalia [2]. Khat is also known as the “Leaf of Allah,” [3] and others express it as a‟ flower of paradise” [4]. Khat (Catha Edulis) is a wild green leaf and greenish plant that grows commercially in East Africa and the Arabian Peninsula whose branches and leaves are chewed to release liquids having the active chemicals cathinone and cathine, which alter the user’s mood [4, 5]. Chewing Khat has a long history dating back to the early 14th century and is commonly used for social and spiritual purposes [6]. Although initially Muslim prayers used Khat for ceremonies, it is also widely used by Christian followers for recreational purposes [7]. Ethiopia is currently a growing producer of Khat [810]. Land for Khat production from 2001/2002 to 2014/2015 increased by 160%, and new areas of the country, which had not been used before, cultivated rapidly [11]. It plays a significant role in the social life of many people in most parts of Ethiopia and is economically beneficial to producers.

The practice of chewing Khat in the Horn of Africa and the Arabian Peninsula has been a popular practice for centuries. Particularly in Ethiopia, the number of Khat users has increased dramatically in recent years, with an average distribution of 15.3% [12, 13]. In the new era, for example, in Harar, even children under the age of five harvest Khat [2]. However, the average age of Khat users in Dessie, Ethiopia, is 19.23 years [14]. Regular users of this plant will find comfort, care, and a feeling of joy. Postchewing effects are often accompanied by sleeplessness, irritability, and poor concentration [4]. Excessive consumption of Khat can lead to serious family, health, and financial problems [15, 16]. Unfortunately, this practice has spread in the Western countries [17, 18] and has become a global problem [19].

Chewing Khat is now a public health problem, including heart failure and high blood pressure. At present, Khat chewing is a means for disseminating COVID 19 for the reason that, during harvesting, transporting, distributing, selling, and selecting a package, of course, the communal and individual act of Khat is exposed to human contact, and most of the consumers are chewing unwashed (contaminated) leaves of Khat in public places [20].

The prevalence of chewing Khat is high among high school and college students due to peer pressure, family chewing habits, and living conditions. Many researchers say that a complete ban on the production and trade of Khat could harm Ethiopia’s economy in the short term [21]; some say the Ethiopian government has imposed two taxes at the local and federal levels. According to the review, chewing habits have affected the majority of productive age groups, confirmed by many students [7].

According to Dhiraj Kumar Das et al. [22], media awareness is a very effective tool to generate preventive control measures of communicable diseases. To account for the effect of media awareness, there are two approaches: updating the disease transmission rate to accumulate the significant fall in transmission due to preventive measures and by incorporating a mass media compartment to represent the public interaction with mass media.

Mathematical modeling plays an important role in understanding the complexities of communicable diseases such as TB, ebola, rabies, influenza, West Nile virus, zika, measles, pertussis, and HIV/AIDS. The primary purpose of studying communicable disease modeling is to design better health measures to control and ultimately to exterminate the disease. In the last few decades, both mathematicians and epidemiologists have developed many mathematical models to analyze the transmission procedure of several infectious diseases. To determine appropriate public health measures, a threshold quantity, called the basic reproduction number , plays a significant role. In general, if , the disease will die out in time whereas if , the disease will persist and thus transmit through the individuals. Such type of characterization is familiar in most of the epidemiological models. [2325].

As stated by Samui et al. [26], mathematical modeling can aid in understanding the following: how transmissible the disease is, when the infectivity becomes high during the course of an epidemic, how acute the disease is, and how effective interventions has been and ought to be.

Makembo, Karandzha, and David [5] developed a mathematical model of Khat (Miraa) chewing, by dividing the entire population into three parts: the number of susceptible persons , the number of light users , and the number of addicts , and they used the SIS compartment model.

We proposed here a deterministic ordinary differential equation model that can represent the dynamics of Khat-chewing communicable diseases, and its control mechanisms in public health and socio-economic crises of a country.

It may be useful to assess the effectiveness of control strategies and to study the Khat-chewing methods used to predict Khat-chewing patterns. We stratified the total human population into five compartments, namely, a group of people who do not chew Khat, designated by , a group of people who are surrounded by Khat chewers but do not chew at present; they may chew Khat in the future, and we denote this as , a group of people who chew Khat, which we represent in , a group of people containing individuals who consumed Khat quite temporarily for social, spiritual, and recreational purposes, and we describe this group in , and a group of people those who constantly chew Khat, and they are denoted by , and the general population is denoted by and is expressed as . We also note that people who chew Khat may return to vulnerable groups due to social, economic, environmental, and religious factors with a rate and constant and diverse population. In this study, we used a compartment model similar to SEIR [2229]. In this deterministic model, we have considered a system of first-order nonlinear ordinary differential equations, which is derived and analyzed analytically and analyzed only numerically using Matlab ode45. From our dynamical system, we got the free equilibrium point of Khat chewing, that is, the point where Khat chewing is not practiced in a given population, and Khat chewing presents equilibrium points and the Khat-chewing generation number, denoted by . The generation number is intended to determine whether the chewing habits of the population have disappeared or increased.

The simulations of our model were performed using parameters with their estimated values. Through simulations, the model is analyzed to obtained different situations that produce interesting results among specific classes. The article is organized as follows; in the second section, we formulate a Khat-chewing model, by classifying a constant and heterogeneous entire population into five groups; in the third section, we analyzed our model, checked its biological completeness and mathematically correct significance, determined equilibrium points and the Khat-chewing generation number, and tested existence and stability of equilibrium points of the dynamical system. In section four, we have presented numerical simulations that support our results. Finally, we conclude the manuscript with a conclusion and recommendation in the fifth section.

2. Model Formulation for Khat Chewing

Our model contains five first-order nonlinear ordinary differential equations. The parameters we used in our study are all positive and ranges between and are listed and described in Table 1, and our models’ diagram representation is shown in Figure 1.

We symbolized the entire population using a flow chart as follows.

The mathematical representation of the above flow chart is described as follows:

Subject to the initial conditions, , , , , and .

3. Model Analysis

3.1. Positivity of Model Solutions

Theorem 1. If initial conditions of (1) are all positive, then it has a positive solution.

Proof. Suppose , , , , and , then(i) Hence, the solution of is positive.(ii) Hence, the solution of is positive.(iii)Hence, the solution of is positive.
Using a similar method, solutions of and are positive.
Since all solutions of the dynamical system are positive, the feasible solution containing all solutions of (1) has the form:

3.2. Boundedness of Solutions

To say that the dynamical system (1) is biologically complete and mathematically correct, all its solutions must be bounded.

Theorem 2. For all , a region exists in which all solutions of (1) are contained and bounded.

Proof. The rate of change of total population size is defined as follow:By substituting all values of left-hand sides of (1) in the above rate of change of total solution, we obtained ; assuming that disease-induced death rate , initial population size , and , for all ; hence, all solutions of (1) are bounded.

3.3. Existence and Stability of Equilibrium Points
3.3.1. Scaling

To reduce the complexity of our model, we resize the fraction by scaling it as follows. Assuming that , , , , and , our dynamical system (1) is reduced to the following:subject to the initial conditions, , , , , and .

Since the compartment does not appear in the first four equations of the model (5), the dynamics of (5) is the same as the following one:

3.3.2. Equilibrium Points

Khat chewing-free equilibrium points: Khat chewing-free equilibrium point of (1) is a balance point that exists when Khat chewing is not practiced in a given population, that is, if there are no people chewing Khat. We determine Khatchewing-free equilibrium point of (6) by assuming and substituting it into the right side of (6), and finally, equate the result to zero.

Hence, and .

Khat chewing-present equilibrium points: Khat chewing-present equilibrium point of (1) is the equilibrium point that exists when Khat chewing is practiced in a given population, that is, if there are people who chew Khat. We determine Khat chewing-present equilibrium point of (6) by assuming at least one of the variables, .

Hence,

3.3.3. Khat-Chewing Generation Number

The Khat-chewing generation number is used to determine secondary Khat chewers, i.e., one Khat chewer individual would create another chewer over the duration of the Khat-chewing period, as long as everyone else is exposed to Khat chewers. The Khat-chewing generation number plays an important role in controlling the significant increase in the number of Khat chewers. In a more general way, the Khat-chewing generation number can be stated as the number of new infections (new Khat chewers) created by a typical infective population at a disease-free equilibrium. < 1 determines, on average, an infected population creates less than one new infected population during the course of its infective period, and the infection can die out but for  = 1, an infected individual can infect on an average 1 person, that is, the spread of the disease is stable. In the reverse way, > 1 determines each infected population creates, on average, more than one new infection, and the disease can spread over the population [26, 29, 30]. To determine the Khat-chewing generation number, we used the production matrix method described in [26, 2933].

Afterwards, we determine the square matrix , which is the emergence rate of new Khat-chewing devices (new infections), and which is the transmission rate of individuals. The order of both matrices is , where, is the number of infected compartments. The new infected compartment and transfer compartments are and , respectively, and we get a new subsystem.

Then, , and .

Because is nonsingular, is also nonsingular. is called the next-generation matrix of our model, and its maximum positive Eigen value is known as the Khat-chewing generation number:

The eigenvalues are and .

The maximum positive eigenvalue of which is the Khat-chewing generation number is

If , then , implies no Khat chewers or the number of Khat chewers reduced to zero [26, 29, 30].

If , then , which indicates that the insignificant inflow of Khat chewers into the community will produce a significant number of Khat chewers [26, 29, 30].

Stability analysis of the Khat chewing-free equilibrium point (analytically): the Jacobian matrix of (6) is defined as follows:

The Jacobian matrix of (6) at Khat chewing-free equilibrium point is as follows:

Theorem 3. The Khat chewing-free equilibrium point of (1) is stable if , and it is unstable if .

Proof. Proving local stability of leads us to proving local stability of .
The characteristics equation is determined using :, whereThe eigenvalues are , and we can find roots of using the Routh–Herwitz stability criterion.

Definition 1. The Routh–Herwitz stability criterion [34] states the following:(a)If there is no sign change in the first column of the Routh–Herwitz table, the system is stable.(b)The number of roots of a polynomial lying in the right half-plane is equal to the number of sign changes in the first column. Consequently, the system is unstable if the poles are on the right.If , then there are no sign changes in the first column of the Routh–Herwitz table. Using the above definition, the system is stable; this implies that (the Jacobian matrix has a negative eigenvalue if ), and is stable. Hence, is stable.
If at least one of are positive, there is a sign change in the first column of the Routh–Herwitz table. Then, the system is unstable; this implies that , and is unstable. Hence, is unstable.
Stability analysis of the Khat chewing-present equilibrium point (analytically): the Jacobian matrix of (6) is defined as follows:The Jacobian matrix of (6) at Khat chewing-present equilibrium point iswhere

Theorem 4. The Khat chewing-present equilibrium point of (1) is stable if and unstable if .

Proof. Proving stability of leads us to proving local stability of .
The characteristics equation of the Jacobian matrix of (6) at is, whereIf , or all are negative, then there are no sign changes in the first column of the Routh–Herwitz table; hence, the system is stable; this implies that , and is unstable. Therefore, is unstable.
If at least one of are negative or positive, there are sign changes in the first column of the Routh–Herwitz table. Then, the system is unstable; this implies that , and is stable. Hence, is stable.

4. Model Numerical Simulation

4.1. Parameter Estimation and Stability Analysis of Equilibrium Points Numerically

Table 2 shows the estimated and recommended parameter values and their corresponding Khat-chewing generation numbers. Using the results, we analyzed our model numerically.

4.2. Stability Analysis of Khat Chewing-Free Equilibrium Point

We analyzed the stability of Khat chewing-free equilibrium point using estimated values of column 1 and 4 of Table 2 as follows.For column 1,The characteristics equation of is determined using ; then, its characteristic’s equation isThe eigenvalues are . Since all Eigen values are negative (no sign change in the first column of the Routh–Herwitz table), then the system is stable which implies; then, Khat chewing-free equilibrium point is stable. Hence, Khat chewing-free equilibrium point is also stable. Therefore, the number of Khat chewers is reduced to zero. For column 4,

Its characteristic equation is , and its roots are . There is a sign change in the first column of the Routh–Herwitz table or the roots have no same signs; then, the system is unstable which implies ; then, Khat chewing-free equilibrium point is unstable. Hence, Khat chewing-free equilibrium point is also unstable. is unstable or , which indicates an insignificant inflow of Khat chewers into the community will produce a large number of Khat chewers.

4.3. Stability Analysis of Khat Chewing-Present Equilibrium Point

We analyzed the stability of Khat chewing-present equilibrium point , using estimated values of column 2 of Table 2 as follows. where

Its characteristics equation is :

There is no sign change in the first column of the Routh–Herwitz table or the roots have the same sign; then, the system is stable; this implies that , and is unstable. Hence, is unstable.

4.4. Numerical Simulation

In this section, we cherry-pick which parameter affects the model in the event of an uncontrolled spread of Khat chewers in the community and which parameter affects the model when the number of Khat chewers in the community decreases. For initial conditions , , , , and and using column 1 of Table 2, we obtained Figure 2 which illustrates that the number of non-Khat chewers, the number of exposed, number of Khat chewers, and the number of temporary Khat chewers reduced or approaches to zero as time rises, and the number of permanent Khat chewers increases as time rises.

For initial conditions , , , , and and using column 1 of Table 2, we obtained Figure 3 which demonstrates that the Khat-chewing generation number increases as the value of increases. Hence, if , , and the number of Khat chewers reduces to zero, and if , , andthe inflow of an insignificant number of Khat chewers into the community produces a significant number of Khat chewers.

For initial conditions , , , , and and using column 1 of Table 2, we obtained Figure 4 which demonstrates that the Khat-chewing generation number decreases as the value of increases; then, the spreading of Khat chewers declines as increases.

For initial conditions , , , , and and using column 4 of Table 2, we obtained Figure 5 which represents that the Khat-chewing generation number increases as the value of increases. Therefore, for , , andthe inflow of an insignificant number of Khat chewers into the community produces a significant number of Khat chewers.

For initial conditions , , , , and and using column 4 of Table 2, we obtained Figure 6 which symbolizes that the Khat-chewing generation number increases as the value of increases in column 4 of Table 2 for . From this graph we conclude that the inflow of a insignificant number of Khat chewers into the community produces a significant number of Khat chewers even the rate of return back from Khat chewing to exposed because of socio-economic, environmental, and religious factors increases. Thus, to turn down the spreading of Khat chewers to the community, the rate of return from Khat chewing to an exposed group because of socio-economic, environmental and religious influences must grow exponentially.

5. Conclusion and Recommendation

In this study, we examined the return from Khat chewing to the exposed group due to socio-economic, environmental, and religious factors, as well as the rate of conversion from Khat chewing to the exposed group of consequences in the Khat chewing dynamics’ model. We have derived the Khat chewing generations’ number and Khat chewing-free and Khat chewing-present equilibrium points. For , the Khat chewing-free equilibrium point is stable, and it is unstable for , and the Khat chewing-present equilibrium point is stable if , and it is unstable if .

The stability of the two equilibrium points implies that the influx of the insignificant number of Khat chewers into the community produces a significant number of Khat chewers if the conversion rate of nonchewers to vulnerable groups is high, and the inflow of Khat chewers to the community produces an insignificant number of Khat chewers if they return back from Khat chewing to the exposed group because socio-economic, environmental, and religious influences raise rapidly. We recommend to the community to stop or minimize chewing Khat habits because the practice of chewing Khat leads to economic crises, social crises (turn into isolated), and health crises, and since media awareness are a very effective tool to generate precautionary control measures of communicable diseases [22], so both governmental and nongovernmental health office stakeholders must alert their people about the impact of Khat chewing through media. In addition, to future researchers, we recommend them to extend the model by grouping the entire population based on sex, religion, and geographic locations.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.