Modeling dynamics of fast food and obesity for evaluating the peer pressure effect and workout impact

In recent years, chronic diseases, such as high blood pressure, diabetes, heart attack, and cancer, have increased around the world. Obesity is a common factor that makes individuals susceptible to these diseases. One reason for excessive weight gain is the frequent consumption of fast food. This study examined the impact that fast food has on obesity by analyzing the influence of peer pressure on fast food consumption and the role that exercise plays in weight gain. Two mathematical models were proposed to study each factor separately. The models demonstrated the dynamics of individuals transitioning from one weight class to another according to their body mass index (BMI). The models were governed by systems of nonlinear ordinary differential equations. A qualitative approach for analyzing the systems was utilized, and the numerical simulations supported the qualitative results. The results of this study may contribute to the development of strategies for controlling obesity. In particular, the models showed the importance of resisting peer pressure that drives individuals to eat fast food meals, as well as maximizing the role that quitters can play in convincing obese individuals to stop consuming fast food meals. In addition to the two previous strategies, the research revealed that physical activity has a significant role in reducing weight.


Introduction
Obesity is a condition defined by increased levels of fat in the human body. Individuals are classified according to their weight using the global body mass index (BMI). This measurement is calculated by dividing the body mass (weight) over the square of the body height. If the output is less than 18.5, then the person is classified as underweight. The normal and ideal weight for a person occurs if the BMI is between 18.5 and 24.9. If the BMI is between 25 and 30, then the person is overweight and is more likely to become obese. Once the BMI score is more than 30, then the person is considered obese.
Many factors may contribute to weight gain, such as environment, culture, metabolism, and genes. However, the primary cause is an excessive intake of calories over the daily expenditure of energy. Unfortunately, the nature of modern life has influenced individuals to be less physically active, and there has been a dramatic increase in fast food restaurants. Since fast food is cheap and readily available, individuals often make it their first choice. Although it has a rich and satisfying taste, fast food is high in calories and contains excessive fat and sugar. Studies have shown that frequent visits to fast food restaurants affect weight gain [5,25,29,30].
The Organization for Economic Co-operation and Development (OECD) published a report [26] identifying a significant increase in obesity rates in countries around the world. The United States was at the top of the OECD list, where obesity increased from 14% to approximately 34%, which is about one third of the population. Also, the Ministry of Health in Saudi Arabia conducted a statistical survey on obesity in 2013 [22]. These statistics showed that obesity reached 24% among males and 34% among females in Saudi Arabia. In recent years, more chronic diseases, such as high blood pressure, diabetes, heart attack, and cancer, have spread around the world. One primary factor to being susceptible to these diseases is obesity [7,11,23,27]. As a result, obesity plays a vital role in reducing the average life of individuals.
It is highly likely that social interaction influences individual's habits. Researchers have specifically been examining the relationship between peer pressure and obesity. A study conducted on European children [16] reported a link between peer influences and obesity, with the largest effects on children in Spain, Italy, and Cyprus. The same connection was found in [24] for Chinese children and adolescents. One of the first studies to measure peer pressure found that adolescents with higher BMIs have a significant impact on their peers and that females are influenced more by their peers than males [32]. Similar evidence stated in [21] supports the claim that peer pressure significantly impacts the weight of adolescent females more than that of adolescent males. In addition, they found that same-gender peers have a stronger influence among females. Peer pressure in relation to weight-related behavior has also been investigated, especially in terms of how peer pressure among adolescents affects the frequency of eating fast food [4]. The findings in [15] showed a relation between fast food consumption among adolescents and their peers. However, there is minimal difference whether the effect of peers is used or ignored as a factor in the statistical analysis. Along with reducing the intake of calories, physical activity is an essential way to lose excess weight [8,13,33]. Hankinson et al. [17] conducted a longitudinal study over 20 years to evaluate the impact of high levels of physical activity on BMI and waist circumference. They found that young adults who maintained high levels of physical activity (at least 30 minutes a day) throughout their lives gained fewer kilograms during the transition to middle age. Similarly, Donnelly et al. [10] studied the effect of exercise on sedentary, overweight, and obese young adults. Without any change in their diet, a program consisting of 16 months of exercise produced a significant reduction in BMI for males and prevented weight gain in females.
Mathematical models have been used to study and analyze many social phenomena and epidemic diseases. They are capable of predicting global dynamic behavior and shedding light on factors with the most significant influence. Lately, researchers have employed models to investigate the role that fast food plays in relation to obesity. In 2004, Evangelista et al. [14] studied the role of peer pressure in fast food consumption, as well as its effect on individual weight in the United States. They found that a reduction in peer pressure to eat fast food could reduce obesity rates. Jodar et al. [20] developed an obesity model for children 3 to 5 years old, and they studied the effect of the bakery, fried meals, and soft drinks (BFS) on childhood obesity in the Spanish region of Valencia. They presented numerical simulations showing an increasing trend in childhood obesity from the previous years. Arenas et al. [6] extended the same model to the whole population, but with a non-autonomous term. They modified the transmission rate due to social pressure on BFS consumption as a time interval parameter. They obtained sufficient conditions to guarantee the existence of a periodic positive solution. Finally, Santonga et al. [31] investigated the effect of public health campaigns on reducing excess weight in Valencia. Using the numerical simulations, they concluded that preventive health campaigns to stop people from increasing their food intake were more effective than treatment; however, the combination of both was more efficient.
Most previous studies have examined the obesity models numerically due to their complexity. However, qualitative analysis of mathematical models [9,12,18,28] presents more general results and allows for better understanding of the critical parameters that play a significant role in the dynamics of the models, see for example [1][2][3]19].
In this work, we investigate the proposed mathematical models using both qualitative and numerical analysis. We build two models: model I studies the effect of peer pressure on obesity and fast food consumption; model II is an extension to model I, where workout factor is included. In both models, we find the equilibrium points and their existing conditions in relation to threshold numbers. Also, we explore the local and global conditions for the stability of equilibria using the linearization method and Lyapunov functions. Furthermore, we illustrate numerical simulations to support the qualitative results as well as analyzing the critical parameters of the models.
The outline of the paper is organized as follows. In Sect. 2, we present statistical analysis, which was carried out by the authors, regarding the peer pressure effects on fast food consumption and the impact of physical activity on obesity. In Sect. 3, we formulate model I by demonstrating the transmission mechanism, and then we investigate the model qualitatively and numerically. In Sect. 4, we add the workout factor to model I and examine model II by analyzing the role of exercise in reducing weight. Finally, we give a brief conclusion in Sect. 5.

Statistical survey
A cross-sectional online survey was carried out in Jeddah, Saudi Arabia from 26 April to 28 July 2016. Our aim was to explore the extent to which peer pressure influences individuals to eat at fast food restaurants in Jeddah. Also, we investigated the impact of consuming fast food meals and the role of physical activity on the weight of individuals.
Participants were recruited through technological social networks. The survey consisted of an online questionnaire that participants were asked to complete. The number of participants who completed the questionnaire was n = 1229. The questionnaire was divided into two sections. The first section contained demographic questions about age, gender, BMI, and living area. The second section requested the participant' behavior towards consuming fast food and engaging in physical activity.
We analyzed the data obtained from the study by using SPSS (descriptive analysis) and Excel. The results are presented in this section.
The age distribution of the participants was 138 adolescents (11%), 838 young adults (68%), 175 middle-aged adults (14%), and 78 older adults (6%). Thus, the majority of the participants were young adults. Also, most of the participants, at a total of 1008 (82%), were female. The participants were classified according to their weight using BMI as follows: 501 individuals (41%) were normal-weight, 354 individuals (29%) were overweight, 277 individuals (23%) were obese, and 97 individuals (8%) were underweight. In terms of the living areas of the participants, the majority (45%) lived in north Jeddah, followed by 24% in the center of Jeddah, 17% in the southeast of Jeddah, and 14% in the southwest of Jeddah (see Table 1).
The investigation of the behavior of participants regarding fast food consumption revealed that most of them (86%) ate fast food meals and the majority (93%) had been eating at fast food restaurants for more than one year. The frequency of consumption varied as follows: 378 (36%) ate fast food meals 1-2 times a month; 301 (29%) ate fast food meals once a week; 292 (28%) ate fast food meals 2 to 3 times a week, and only (8%) ate fast food meals four or more times a week. Additionally, more than half of the participants (64%) felt that consuming fast food meals increased their weight. In terms of the main objective that the study aimed to investigate, we found that 711 participants (58%) reported that peer pressure from family and friends was one of the reasons they ate fast food meals. In terms of physical activity, almost half of the participants (45%) performed exercises, with most of them (38%) spending 2 to 3 hours a week working out. Furthermore, 65% of the participants felt that physical activity was effective in reducing their weight (see Table 2).
The analysis in Table 3 shows that 85% of overweight and 88% of obese participants consumed fast food. Also, 66% of obese and 65% of overweight participants ate fast food meals at least once a week. This suggests that there is a strong relationship between weight gain and fast food. In addition, we found that peer pressure was an influential factor that drove 70% of obese and 71% of overweight participants to consume fast food. As for physical activity, 56% of overweight participants who consumed fast food also worked out, compared to 37% of obese participants.

Formulation of the model
We build a mathematical model to study how peer pressure may affect the temptation that individuals face to eat at fast food restaurants and consequently gain weight. We consider a homogeneously mixed constant population at time t denoted by P(t). The population is divided into four separate classes of individuals: normal-weight individuals, N(t); overweight individuals who eat fast food meals, S(t); obese individuals who eat fast food meals, O(t); and quitters who no longer eat fast food meals, Q(t). We assume that normal-weight individuals enter compartment S with rate β, as a result of peer pressure from overweight individuals who eat at fast food restaurants frequently. Also, overweight individuals enter compartment O with rate γ due to peer pressure from obese individuals who eat at fast food restaurants normally. On the other hand, quitters who no longer eat fast food may have positively influenced overweight and obese individuals by convincing them to refrain from consuming fast food to reduce their weight and susceptibility to harmful diseases. Accordingly, overweight and obese individuals move to compartment Q with rates ε 1 and  Fig. 1, where all parameters belong to the interval (0, 1]. Based on the above assumptions, the model is governed by the following nonlinear system of ordinary differential equations: where Note that P (t) = 0, which means that the total population is constant. By rescaling the system with the following substitutions: with x + y + z + q = 1. We consider the domain of the system to be It can be easily proved that is a positively invariant set, i.e., every solution of system (2) starting with initial conditions in remains there for all t > 0.

Equilibrium points and basic reproduction number
The free equilibrium point of system (2) is E 1,0 = (1, 0, 0, 0), which represents a community with normal-weight individuals. Following epidemiological models, we obtain the basic reproduction number R 1,0 of the model by using the free equilibrium point and utilizing the method of next generation matrix [9]. In our model, R 1,0 represents the number of secondary cases of normal-weight individuals who start to gain weight due to consuming fast food meals under peer pressure from a single overweight individual who eats fast food frequently. Now, we define F and V as follows: The Jacobian matrices of F and V at the free equilibrium point E 1,0 are Thus, the next generation matrix is Hence, the basic reproduction number of the system is the spectral radius of matrix To find the equilibrium points of system (2), we solve the equations: x = 0, y = 0, z = 0, and q = 0, hence, the equilibrium points of the model are: ). The positive equilibrium point exists if y 1,4 , z 1,4 , q 1,4 > 0. Set

Stability analysis
We present in this section the study of local and global stability for the equilibrium points using the linearization method and Lyapunov function [14,19,24].
In the same manner, we can prove the local stability of E 1,2 and E 1,3 . Therefore, we only state the following theorems.

Global stability.
Here, we introduce a Lyapunov function [23] to prove the global stability of the positive equilibrium point E 1,4 of system (2). Proof Consider the Lyapunov function: Since L (x, y, z, q) ≤ 0, thus E 1,4 is globally stable in .

Numerical analysis
In this section, we solve the model numerically and present simulations of the mathematical model to show its agreement with previous qualitative results. We considered the per capita peer pressure rate β to be the product of two probabilities: the probability of a normal-weight individual meeting an overweight individual per week and the probability of them eating a fast food meal together. In the same manner, the remaining per capita peer pressure rates γ , ε 1 , and ε 2 were estimated. The numerical simulations of the equilibrium points of model I are illustrated in Fig. 2 with initial conditions from the data in Tables 1 and 2. (i) Free equilibrium E 1,0 . The equilibrium point satisfies the stability condition R 1,0 < 1. Therefore, we choose the following parameters: μ = 0.02, β = 0.01, γ = 0.01, ε 1 = 0.01, ε 2 = 0.01. The numerical simulations in Fig. 2(a) show that the long-term behavior of the solution curves of the system approaches the free equilibrium E 1,0 . Hence the system is locally asymptotically stable about E 1,0 for the above set of parameters. In particular, we see the overweight, obese, and quitters compartments decrease to the value zero, and the normalweight compartment increases to the value 1. This would mean, theoretically, that all individuals in the society eventually achieve normal-weight and refrain from consuming fast food.
(ii) Overweight equilibrium E 1,1 . Figure 2(b) illustrates the changes in the four classes over time under the following parameter settings: μ = 0.02, β = 0.05, γ = 0.01, ε 1 = 0.01, ε 2 = 0.01. The parameters satisfy the qualitative conditions of the overweight equilibrium: R 1,0 > 1, R 1,1 < 1, and R 1,2 < 1. The numerical simulations agree with the qualitative results since all solution curves of the system tend to the overweight equilibrium E 1,1 . There is an increase in the number of normal-weight individuals as a result of the small number of overweight individuals. However, as time passes, the overweight individuals influence normal-weight individuals to eat fast food. With greater consumption of fast food, the number of normal-weights begins to decline to a balance level after reaching a certain peak. Consequently, the number of overweight individuals rises with increasing β to an equilibrium level that is higher than the normal-weight level. In addition, the effect of obese and quitter individuals, in this case, is small in the society, and their numbers in time tend to zero.
(iii) Free quitters equilibrium E 1,2 . In the same manner, Fig. 2(c) numerically indicates that the free quitters equilibrium E 1,2 is locally asymptotically stable. Under the qualitative conditions of E 1,2 : 1 < R 1,1 < R 1,3 , we set the parameters to be: μ = 0.02, β = 0.1, γ = 0.1, ε 1 = 0.01, ε 2 = 0.01. The variation trend in the simulation is different from the previous cases. The number of obese individuals increases to a certain level as a result of their influence on overweight individuals to consume more fast food meals. Consequently, a decline in the number of overweight individuals leads to a rise in the number of normal-weight individuals. This indicates that the impact of overweight individuals on normal-weight individuals over time decreases, since the rise of obese individuals made it easy for overweight individuals to find friends to engage in eating at fast food restaurants with them. The effect of quitters, in this case, is still small and their numbers in time tend to zero.
(iv) Free obese equilibrium E 1,3 . The numerical simulations in Fig. 2(d) exhibit the local asymptotic stability of the free obese equilibrium E 1,3 . We set the parameters to be: μ = 0.02, β = 0.1, γ = 0.1, ε 1 = 0.1, ε 2 = 0.1, which satisfy the conditions of E 1,3 : R 1,2 > 1 and R 1,2 > R 1,3 . In this case, the simulation shows that the normal-weight and overweight classes at first behave in the same manner as in the previous case. However, as the number of quitters starts to rise, the number of overweight and obese individuals decreases below the number of quitters. This reflects the influence that quitters had on overweight and obese individuals to stop eating fast food meals; for this reason, the number of obese individuals eventually tends toward zero.

Parameter analysis
The variation in the classes in model I is affected by the influential forces of peer pressure from overweight, obese, and quitter individuals. In this section, we examine the effect of peer pressure parameters, particularly on the obese class.
(i) Parameters β, γ . We investigate the effects on obesity as a result of the changes in the overweight and obese peer pressure rates. The numerical simulations in Fig. 3 demonstrate the time variation of obesity when the overweight peer pressure rate β is compared with the obese peer pressure rate γ . Note that all other parameters are set to be the following values: μ = 0.02, ε 1 = 0.01, ε 2 = 0.01. The simulations show that obesity has higher values when β > γ compared to when β < γ . In addition, the level of obesity is greater when β = γ = 0.2 than when β = γ = 0.1. Under these circumstances, we find that obesity is high when β is large and vice versa. In other words, the comparison of the overweight and obese peer pressure rates indicates that the overweight peer pressure rate plays a more significant role in confining obesity. (ii) Parameters ε 1 , ε 2 . We explore how the levels of obesity are impacted by alterations to the peer pressure rate of quitter on overweight (ε 1 ) and the pressure rate of quitter on obese (ε 2 ). To begin, we set obesity at the highest level (from the previous case), that is, when β = 0.2 and γ = 0.1. The numerical simulations in Fig. 4 illustrate the time variation of obesity in four cases: ε 1 = ε 2 = 0.01, ε 1 < ε 2 , ε 1 = ε 2 = 0.05, and ε 1 > ε 2 . From these simulations, we find that the level of obesity declines when ε 1 < ε 2 . Furthermore, obesity decreases to a lower level when ε 1 = ε 2 = 0.05. It appears as if the quitters' peer pressure rates had a major impact in helping to control the fast food eating habits of obese individuals.
In conclusion, the analysis designated two key parameters: β and ε 2 as control strategies for obesity. The first strategy is to minimize the peer pressure from overweight individuals rather from obese individuals. The second strategy is to maximize the role of quitters in convincing obese individuals to stop consuming fast food meals.

Model II: workout effect 4.1 Formulation of the model
We modify model I to take into account the individual workout factor. We assume that exercising causes overweight individuals to enter compartment N with rate α 1 and obese individuals to move to compartment S with rate α 2 . Individuals in class Q may lose weight by avoiding fast food meals as well as performing exercises. After losing weight, some quitters may be susceptible to returning to the fast food life style and enter compartment N with rate α 3 . Our aim is to determine the effect of exercise in reducing weight among overweight and obese individuals who frequently eat at fast food restaurants and to inspect the fallback effect of quitters on obesity. The dynamics of the model are presented in Fig. 5. Again, all parameters belong to the interval (0, 1].
The model is governed by the following nonlinear system: Figure 5 The dynamics of model II Note that P (t) = 0, The total population here is again constant. By rescaling the system with the same substitution as in model I: with x + y + z + q = 1. Consider the domain of the system to be = (x, y, z, q) : Again, is a positively invariant set.

Equilibrium points and basic reproduction number
The free equilibrium point of system (4) is E 2,0 = (1, 0, 0, 0), that is, all the individuals in the society have normal-weight. By using the next generation matrix [9], we obtain the basic reproduction number R 2,0 of model II: R 2,0 = β α 1 +μ . As in model I, R 2,0 represents the number of normal-weight individuals affected by the peer pressure of a single overweight individual to start gaining weight due to consuming fast food.

Stability analysis
In this section, we present a local stability analysis of all equilibrium points of system (4) using the linearization method [28] and Routh-Hurwitz criteria [12].
Local stability.
In the same manner, we can prove the local stability of E 2,2 and E 2,3 . Therefore, we state the following theorems.

Numerical analysis
The numerical simulations of model II are presented in this section. The peer pressure parameters β, γ , ε 1 , and ε 2 are defined in the same manner as model I. We set the following assumptions for the workout parameters. To begin, we suppose that an overweight individual who frequently eats at fast food restaurants gains 0.25 Kg per week. Also, an overweight individual who walks at a fast pace for seven hours per week burns 3010 cal which is approximately equal to 0.39 Kg. Therefore, an overweight individual who eats fast food meals and workout burns 0.14 Kg per week. We assume that overweight individuals move to the normal-weight class if they burn 10 Kg. This process takes approximately 72 weeks. Thus, the overweight exercise rate α 1 is equal to 1/72 ≈ 0.0138 per week. Likewise, we assume that an obese individual gains 0.25 Kg from fast food meals but burns 0.44 Kg per week from exercise. Therefore, the obese exercise rate α 2 is equal to 1/105 ≈ 0.0095 per week, when assuming that obese individuals have to lose 20 Kg to move to the overweight class. Furthermore, since the quitters class includes both overweight and obese individuals who do not eat fast food meals, we assume that the average calories burnt in this class from exercise is equal to 0.415 Kg per week. Thus, the quitters exercise rate is equal to 1/48 ≈ 0.0208 per week, when assuming that individuals have to lose 20 Kg to be in the normal-weight class. However, if we presume that 50% of the quitters, who have lost weight from exercising, enter a state of delusion and become vulnerable to fast food meals, then the rate α 3 is equal to 0.0104 per week. From the previous assumptions, we set the workout parameters to be α 1 = α 2 = α 3 = 0.01. On the other hand, the peer pressure parameters are chosen arbitrarily to verify the existence and local stability conditions. The initial conditions in the simulations are chosen from the data in Tables 1 and 2. (i) Free equilibrium point E 2,0 . The parameters were chosen to satisfy the existence and local stability conditions of E 2,0 , that is, R 2,0 < 1. We set the parameters to be as follows: μ = 0.02, β = 0.01, γ = 0.01, ε 1 = 0.01, ε 2 = 0.01, α 1 = 0.01, α 2 = 0.01, α 3 = 0.01. The numerical simulations in Fig. 6(a) show that the long-term behavior of the solution curves of the system approaches E 2,0 . In fact, we notice that the size of normal-weight class increases to its maximum value of one, while the overweight, obese, and quitter classes tend to zero. Thus, eventually, all individuals in the society have normal-weight. When comparing Figs. 6(a) and 2(a), the effect of exercise is demonstrated in lessening the time the system takes to reach the free equilibrium point.
(ii) Overweight equilibrium point E 2,1 . The numerical simulations in Fig. 6(b) exhibit the local stability of the overweight equilibrium E 2,1 with the following choices of the parameters: μ = 0.02, β = 0.05, γ = 0.01, ε 1 = 0.01, ε 2 = 0.01, α 1 = 0.01, α 2 = 0.01, α 3 = 0.01. The parameters satisfy the conditions of E 2,1 : R 2,0 > 1, R 2,1 < 1, R 2,2 < 1. In this case, we note that the number of normal-weight and overweight individuals settles to a certain level in the long term, while the size of the obese and quitters groups tends to zero. When comparing the simulations in Fig. 6(b) with those in Fig. 2(b), we see that the size of the overweight class decreased due to the workout effect.
(iv) Free obese equilibrium E 2,3 . In the same manner, Fig. 6(d) numerically shows that the free obese equilibrium E 2,3 is locally asymptotically stable. Under the qualitative conditions of E 2,3 : R 2,2 > 1 and R 2,5 < 1, we set the parameters to be as follows: μ = 0.02, β = 0.1, γ = 0.1, ε 1 = 0.1, ε 2 = 0.1, α 1 = 0.01, α 2 = 0.01, α 3 = 0.01. Figure 6(d) shows that the population in the long term is divided into normal-weight, overweight, and quitter individuals, but no obese individuals. Furthermore, a comparison between Figs. 2(d) and 6(d) shows that the number of quitters has decreased due to the fallback of those who lost weight from exercise and returned to being vulnerable to fast food meals. Consequently, the size of the overweight class has increased owing to the effects of peer pressure.
(v) Positive equilibrium E 2,4 . In the last case, the parameters satisfy the existence conditions of E 2,4 . Therefore, we choose the following parameters: μ = 0.02, β = 0.1, γ = 0.1, ε 1 = 0.05, ε 2 = 0.05, α 1 = 0.01, α 2 = 0.01, α 3 = 0.005. The numerical simulation in Fig. 6(e) shows that the long-term behavior of the solution curves of the system approaches the positive equilibrium E 2,4 . Hence the system is locally asymptotically stable at E 2,4 for the above set of parameters. When comparing Figs. 6(e) and 2(e), we see that the size of quitters has decreased as a result of exercise and being susceptible to returning to fast food meals. Moreover, the size of overweight and obese individuals has increased due to the effect of exercise and peer pressure.

Parameter analysis
In this section, we investigate the effects of exercise rates and peer pressure rates on obesity. For this purpose, we present numerical simulations for low and high workout rates, as opposed to low and high peer pressure rates. A society with only peer pressure rates and no workout rates is presented in model I, whereas model II undertakes to show the impact of both peer pressure and workouts. We set the workout parameters to be as follows: α 1 = α 2 = α 3 = 0.01 (for low workout rates) and α 1 = α 2 = α 3 = 0.08 (for high workout rates). The peer pressure parameters are set according to the following cases: Case (1): We examine the case of low fast food peer pressure rates, that is, β = γ = 0.05, with low quitters' peer pressure rates of ε 1 = ε 2 = 0.01. Under these considerations, Figure 7 Time variation of obesity for low and high workout raters with (a) low peer pressure rates, (b) low fast food pressure rates and high quitters peer pressure rates, (c) hight fast food peer pressure rates and low quitters peer pressure rates, and (d) high peer pressure rates Fig. 7(a) demonstrates the time variation of obesity. The simulation shows that high workout rates fight obesity better than low workout rates. Also, high workout rates take less time to diminish obesity than low workout rates.
Case (2): We let the fast food peer pressure rates to be the same, that is, β = γ = 0.05, but with high quitters' peer pressure rates of ε 1 = ε 2 = 0.1. Following these assumptions, Fig. 7(b) illustrates how obesity eventually reaches zero in less time than in case (1). This reflects the significance of workouts as well as the role of quitters in the society.
Case (3): We adjust the fast food peer pressure rates to a higher value, that is, β = γ = 0.2, but with low quitters' peer pressure rates of ε 1 = ε 2 = 0.01. Accordingly, the numerical results in Fig. 7(c) exhibit the eventual disappearance of obesity with high workout rates. However, for low workout rates, obesity remains in the society. This indicates that the low workout rates with high fast food peer pressure rates do not help in minimizing obesity.
Case (4): We set all peer pressure rates to be high, that is, β = γ = 0.2 and ε 1 = ε 2 = 0.1. According to these settings, Fig. 7(d) displays the time variations of obesity. From the simulations, we find that obesity could be eliminated over time when quitters' peer pressure rates are high. The same applied to workout rates; however, it took longer for obesity to completely disappear with low workout rates than it did with high workout rates.
In summary, the analysis shows that exercise is a significant factor in battling obesity when fast food peer pressure rates are low. However, when these rates are high, then exer-cise does not play a major role in diminishing obesity unless workout rates are also high. Moreover, higher values of quitters' peer pressure rates support workout rates in decreasing obesity regardless of the values of fast food peer pressure rates. Furthermore, we find that increasing workout rates reduces the time it took for obesity to be eliminated from society. In conclusion, obesity may be limited by three control measures: first, reducing peer pressure from overweight and obese individuals to consume fast food meals; second, increasing the efforts of quitters in persuading others to avoid fast food meals; and third, maximizing society's awareness of the correlation between exercise and weight loss.

Conclusion
In this research, two mathematical models were built to examine the distribution of individuals between weight classes based on fast food consumption. The models were represented through a system of nonlinear ordinary differential equations. To understand the resulting behavior, a qualitative approach was utilized by finding the equilibrium points and analyzing their stability. Furthermore, numerical simulations were displayed to support the qualitative results and to analyze critical parameters that suggest how some strategies control obesity.
The first model examined the influence of peer pressure from family members and friends on the average intake of fast food meals. We divided the population into four classes: normal-weight individuals, overweight individuals who ate fast food meals, obese individuals who ate fast food meals, and quitters who no longer eat fast food meals. The results of the qualitative study showed that the system had five equilibrium points. In particular, the free equilibrium point (E 1,0 ) existed with no conditions and was locally stable when R 1,0 < 1. However, if R 1,0 > 1, then the overweight equilibrium (E 1,1 ) existed but was locally stable when R 1,1 < 1 and R 1,2 < 1. Meanwhile, if R 1,1 > 1, then the free quitters equilibrium (E 1,2 ) existed and was locally stable when R 1,1 < R 1,3 . Moreover, if R 1,2 > 1, then the free obese equilibrium (E 1,3 ) existed and was locally stable when R 1,2 > R 1,3 . As for the positive equilibrium (E 1,4 ), it existed when 0 < R 1,2 < R 1,3 < R 1,1 < R 1,0 ; however, it was globally stable with no conditions. Furthermore, we solved model I numerically and the simulations of the model agreed with the previous qualitative results. Finally, the analysis of the model' parameters suggested that reducing the peer pressure from overweight individuals inflicted on normal-weight individuals had the most significant impact on controlling obesity, along with maximizing the role of quitters in convincing obese individuals to stop consuming fast food meals.
The second model examined the effect of physical activity on reducing weight gain brought on by peer pressure to consume fast food meals. This model was built by adding the exercise factor to model I. The qualitative study of model II produced the same equilibrium points as of model I, but with different existence and local stability conditions. The free equilibrium point (E 2,0 ) existed with no conditions and was locally stable if R 2,0 < 1. However, when R 2,0 > 1, then the overweight equilibrium (E 2,1 ) existed and was locally stable if R 2,1 < 1 and R 2,2 < 1. On the other hand, when R 2,1 > 1, then the free quitters equilibrium (E 2,2 ) existed and was locally stable if R 2,4 < 1. Moreover, when R 2,2 > 1, then the free obese equilibrium (E 2,3 ) existed and was locally stable if R 2,5 < 1. Finally, the positive equilibrium (E 2,4 ) existed when R 2,3 > 1 and α 2 +μ γ < y 2,4 < α 3 +μ ε 1 , but it was locally stable if the eigenvalues of the Jacobian matrix satisfied the Routh-Hurwitz conditions. The research established agreement between the qualitative and numerical results. In addition, the analysis of the workout parameters revealed the effects of physical activity in reducing weight.
In conclusion, the results of the mathematical models were consistent with what was observed in the statistical survey (see Sect. 2) as well as in the literature. Thus, the models are suitable for future analysis. From the results, we suggest the following control strategies to combat the development of obesity: • Emphasizing through education the importance of fighting peer pressure that drives individuals to eat fast food meals; • Maximizing the effort of quitters in convincing obese individuals to refrain from consuming fast food meals; • Spreading awareness in the society of the significance of physical activity in preventing weight gain.