Journal of Applied Mathematics and Physics
Vol.07 No.06(2019), Article ID:93293,15 pages
10.4236/jamp.2019.76088

Mathematical Model of the Dynamics of Rumor Propagation

Samuel Musa1, Mohammed Fori2

1Modibbo Adama University of Technology Yola, Yola, Adamawa State, Nigeria

2Department of Mathematics, College of Education Waka Biu, Waka-Biu, Borno State, Nigeria

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: April 28, 2019; Accepted: June 24, 2019; Published: June 27, 2019

ABSTRACT

In this study, we proposed a deterministic mathematical model that attempts to explain the propagation of a rumor using epidemiological models approach. The population is divided into four classes which consist of ignorant individuals, I ( t ) , spreaders targeting community through media, M ( t ) , spreaders targeting community through verbal communication, G ( t ) and stiflers, R ( t ) . We explored existence of the equilibria and analyzed its stability. It was established that rumour-free equilibrium E o is locally asymptotically stable if R o < 1 ; meaning rumor can seize spreading in a population, and unstable if R o > 1 leads to new rumor spreading in the population. Numerical simulations of the dynamic model are carried out on the system to confirm the analytical results. We see that the dynamics of rumor spreading show similar behavior to that found in the dynamics of infectious diseases except that the spread depends on the classes of spreader.

Keywords:

Rumor, Propagation, Rumor-Free Equilibrium, Rumor-Endemic Equilibrium

1. Introduction

The propagation of rumors is a typical form of social interaction that exerts a powerful influence on human affairs. Rumors represent unproven expositions about or interpretations of news, events, or problems that are of public interest. Because rumors are unconfirmed information, it is hard to determine whether they are true or false [1] . Traditionally, rumors are propagated by word of mouth [2] [3] [4] . Nowadays, with the emergence of the internet, rumors spread by instant messengers, emails or publishing blogging and social media that provide a faster speed of transmission [5] [6] . Once a rumor spreads, the truth is at risk of being distorted in the public sphere. Many case studies indicate that most rumors are false news, which exert negative impacts on the public and on social security [7] . Rumour is a generally circulated story, report or statement without facts to confirm its truth; it is like a virus, once it is transmitted to an individual, an outbreak of events happens in a short period of time. It is thus that spreading a rumor is like a contagious disease that does not need more than a few infected individuals to infect the entire population. If a rumor starts in a population, sooner or later everyone will know it [8] .

In recent years, online social media is growing rapidly. Social media provides a convenient communication scheme for people. Meanwhile, the scheme enables unreliable sources to spread large amounts of unverified information among people [9] [10] [11] [12] [13] . Rumors are thus possible to spread more quickly and widely through online social media compared to traditional offline social communities. The wide spread of misinformation may bring disorder to people especially when they are facing religion or political crises. This indicates that it is crucial for social media to identify misinformation in time so as to limit the spread of rumors [14] .

At the beginning, mathematical models for the rumors were considered merely speculative and imprecise, but for the fact that rumor spreading is now seen like the transmission of disease [15] , it leads to a deeper understanding of the future spread of rumor. Rumor transmitted diseases are particularly well suited for modelling because the route of contact is clearly defined: it is through human activities [16] [17] . In addition, there are so many questions that come to play when it comes to addressing the spread of rumor. Modeling is the only practical approach. The spreading of rumor is in many ways similar to the spreading of epidemic infection by the spreader or the infections to notify or infect the susceptible [18] [19] [20] [21] [22] . After notified or infected, that the spreader can become the stifler is similar as the infections can recover after sometimes [23] .

2. The Rumor Models

Rumor can be viewed as an infection of the mind [24] [25] [26] [27] . Therefore, the rumor propagation problem can be studied with many modeling techniques used in the study of epidemics. Some of these approaches include deterministic models, stochastic models and complex network [28] .

Different possible behaviors in the dynamics of rumor spreading were studied by [25] . In their study, they proposed a deterministic propagation rumor model and detected relevant conditions, derived from local stability analysis of the rumor free equilibrium and rumor-endemic equilibrium.

Rumor propagation through different types of mathematical models was investigated by [29] . In their study, they consider rumor propagation with truth-spreading and determine the threshold which governs the dynamics of the

system.

The dynamics of rumor model in complex heterogonous networks was analyzed by [6] . They first introduced a useful stochastic method that allows them to obtain meaningful time profiles for the quantities characterizing the propagation process.

3. Model Formulation

Following the classical assumption on the variables described in Table 1, we formulate a deterministic mathematical model to describe the propagation dynamics of rumors. In the model, the total population is divided into four epidemiological classes of ignorant I ( t ) , spreaders targeting community through media M ( t ) , spreaders targeting community through verbal communication G ( t ) and stiflers R ( t ) . In the formulation of the model, we wish to clearly state some essential assumptions that:

1) Recruitment into the ignorant compartment is in constant rate,

2) There is no movement between the classes of spreaders,

3) Ignorant individual become spreader after interaction with spreader through media or spreader through verbal communication,

Table 1. Variables and parameters of the model and their descriptions.

4) The probability of an ignorant to become a spreader through media or verbally are equal,

5) The spreading process are mutually exclusive event,

6) Spreaders through media and verbal communications become stifler at the rate α 1 and α 2 respectivelly,

7) Proportion of ignorant will interact with the spreaders at the rate σ and become stifler.

From the model diagram (Figure 1), we have ignorant class I ( t ) to be individuals who are not aware of the rumor. The model is structured such that ignorant individual can either be infected with rumor by individual in M ( t ) or G ( t ) . Note that, proportion σ of ignorant individual can easily progress to R ( t ) without spreading the rumor. The spreader class M ( t ) are individual who are aware of the rumor and are actively spreading it through media. The spreader class G ( t ) are individuals who have been informed about the rumor and they are actively spreading it orally to the community. The stifler R ( t ) , are individuals who are individuals who are aware of the rumor, spent time spreading it but are no longer spreading the rumor due to some reasons and they also discourage the spread of the rumor.

Thus, the differential equations for the deterministic model are as follow:

d I d t = π μ I θ 2 β 1 I M θ 2 β 2 I G ( M + G ) σ I (1.1)

d M d t = θ 2 β 1 I M μ M α 1 M (1.2)

d G d t = θ 2 β 2 I G μ G α 2 G (1.3)

d R d t = ( M + G ) σ I + α 1 M + α 2 G μ R (1.4)

N = I + M + G + R (1.5)

4. Basic Properties

Since the model (1.1) to (1.4) monitors human populations, all the variables and the associated parameters are non-negative at all time. It is important to show that the model variables of the model remain non negative for all non-negative initial conditions.

Lemma 1: The region D = { ( I , M , G , R ) R + 4 : I + M + G + R π μ } is positively invariant and attract all solutions in R + 4 .

Proof:

Adding all the equations from (1.1) to (1.4), gives the rate of change of the total human population

d N d t = π μ N

Figure 1. Model flow diagram.

Since d N d t = π μ N whenever N ( t ) > π μ , then d N d t < 0 , implying d N d t is bounded by π μ N .

Thus, a standard comparison theorem by [30] can be used to show that

N ( t ) N ( 0 ) e μ t + π μ ( 1 e μ t )

In particular, N ( t ) π μ if N ( 0 ) π μ . Thus, R is positively invariant (i.e. all solution in D remain in D for all time). Furthermore, if N ( t ) > π μ then either the solution enters R in finite time or N ( t ) approaches π μ and the spreader variables M and G approaches zero. Hence, D is attracting (i.e. all solution in R + 4 eventually approach, enter or stay in D). Therefore, the model is epidemiologically and mathematically well posed since all the variables remain non-negative for all t 0 . Hence it is sufficient to study the dynamics of the system (1.1) to (1.4) in D.

5. Analysis of the Model

In this section, we discuss the existence and uniqueness of Rumor Free Equilibrium (RFE) of the model and its analysis. The model Equations (1.1) to (1.4) has an RFE given by

E o = ( I * , M * , G * , R * ) = ( π μ , 0 , 0 , 0 ) .

The local stability of RFE given will be investigated using the next generation matrix method. We calculate the next generation matrix for the system of the question (1.1) to (1.4) by enumerating the number of ways that: 1) new spreaders arise 2) number of ways that individuals can move but only one way to create a spreader. So, let

F = rate of appearance of new spreaders into the compartment and,

V = rate of transfer into (out) of compartment

F = ( θ 2 β 1 I * 0 0 θ 2 β 2 I * )

V = ( μ + α 1 0 0 μ + α 2 )

Hence the NGM with large domain is two dimensional and is given by F V 1

F V 1 = ( θ β 1 π 2 μ ( μ + α 1 ) 0 0 θ β 2 π 2 μ ( μ + α 2 ) ) (1.6)

The dominant eigenvalue of (1.6) is equal to R o , therefore we evaluate the characteristic equation of (1.6) by using | ( F V 1 ) λ I | = 0 , which gives

R o = R o m + R o v

R o = π θ [ β 1 ( μ + α 2 ) + β 2 ( μ + α 1 ) ] 2 μ ( μ + α 1 ) ( μ + α 2 )

where

R o m = θ β 1 π 2 μ ( μ + α 1 )

R o v = θ β 2 π 2 μ ( μ + α 2 )

The jacobian of (1.1) to (1.4) at the equilibrium point E o = ( π μ , 0 , 0 , 0 ) is

= ( μ ( θ β 1 + σ μ ) π 2 μ ( θ β 2 + σ μ ) π 2 μ 0 0 π θ β 1 2 μ μ α 1 0 0 0 0 π θ β 2 2 μ μ α 2 0 0 α 1 + σ π μ α 2 + σ π μ μ ) (1.7)

Now we try to calculate the eignvalues of (1.7) by finding the characteristic equation using the formula | J E λ I | = 0

| ( μ λ ) ( θ β 1 + σ μ ) π 2 μ ( θ β 2 + σ μ ) π 2 μ 0 0 ( π θ β 1 2 μ μ α 1 λ ) 0 0 0 0 ( π θ β 2 2 μ μ α 2 λ ) 0 0 α 1 + σ π μ α 2 + σ π μ ( μ λ ) | = 0 (1.8)

Solving (1.8), we have

λ 1 = μ

λ 2 = ( R o M 1 ) ( μ + α 1 )

λ 3 = ( R o V 1 ) ( μ + α 2 )

λ 4 = μ

Theorem 1: The rumor-free equilibrium of the model equation (1.1) and (1.4) given by E o , is locally asymptotically stable if R o < 1 and unstable if R o > 1 . Thus this theorem 1 implies that for any given rumor in a population, it can be eliminated when R o < 1 .

Proof

Having λ 1 and λ 4 to be negative, we also see that λ 2 and λ 3 are both negative too when R o m < 1 and R o v < 1 . Since all the eigenvalues of (1.8) have negative real parts when R o m < 1 and R o v < 1 , we conclude that rumor-free equilibrium is locally asymptotically stable.

Stability Analysis of the Rumor Endemic Equilibrium

When rumor persists in a population (i.e. at least M 0 or G 0 ), the model question (1.1) to (1.4) has two equilibrium points denoted by

E 1 = ( I Λ , M Λ , G Λ , R Λ ) and E 2 = ( I Λ Λ , M Λ Λ , G Λ Λ , R Λ Λ )

called rumor verbal-endemic equilibrium and rumor media-endemic equilibrium point. For the existence and uniqueness of E 1 and E 2 their coordinate has to satisfy the following I Λ > 0 , M Λ > 0 G Λ > 0 , R Λ > 0 and I Λ Λ > 0 , M Λ Λ > 0 , G Λ Λ > 0 , R Λ Λ > 0 respectively.

Equating (1.1) to (1.4) to zero, we obtained

π μ I θ 2 β 1 I M θ 2 β 2 I G ( M + G ) = 0 (1.9)

θ 2 β 1 I M μ M α 1 M = 0 (1.10)

θ 2 β 2 I G μ G α 2 G = 0 (1.11)

( M + G ) σ I + α 1 M + α 2 G μ R = 0 (1.12)

From (1.10), we have

I Λ = 2 ( μ + α 1 ) θ β 1 (1.13)

M Λ = 0 (1.14)

Substituting (1.13) and (1.14) in (1.9), we have

G Λ = π θ β 1 2 μ ( μ + α 1 ) β 2 θ ( μ + α 1 ) + 2 σ ( μ + α 1 ) (1.15)

Using (1.13), (1.14) and (1.15) in (1.12), we obtained

R Λ = ( 2 σ ( μ + α 1 ) + α 1 θ β 1 θ β 1 μ ) ( π θ β 1 2 μ ( μ + α 1 ) β 2 θ ( μ + α 1 ) + 2 σ μ ( μ + α 1 ) ) (1.6)

Thus E 1 = ( I Λ , M Λ , G Λ , R Λ )

Local Stability of rumor verbal-endemic equilibrium

We used the Jacobian Stability approach to prove the stability of the rumor verbal-endemic equilibrium.

The Jacobian of (1.9) to (1.12) at the equilibrium point E 1 = ( I Λ , M Λ , G Λ , R Λ ) is

( μ θ 2 β 1 M Λ θ 2 β 2 G Λ ( M Λ + G Λ ) σ θ 2 β 1 I Λ σ I Λ θ 2 β 2 I Λ σ I Λ 0 θ 2 β 1 M Λ θ 2 β 1 I Λ μ α 1 0 0 θ 2 β 2 G Λ 0 θ 2 β 2 I Λ μ α 2 0 ( M Λ + G Λ ) σ α 1 + σ I Λ α 2 + σ I Λ μ ) (1.17)

If we evaluate (1.17) at E 1 and find the eigenvalue using characteristic equation | J ( E 1 ) λ I | = 0 , we will obtain

λ 1 = μ ( θ β 2 + 2 σ 2 ) ( π θ β 1 2 μ ( μ + α 1 ) θ β 2 ( μ + α 1 ) + 2 σ σ ( μ + α 1 ) )

λ 2 = 0

λ 3 = β 2 ( μ + α 1 ) β 1 μ α 2

λ 4 = μ

From the above questions, it is observed clearly that only λ 1 < 0 and λ 4 < 0 , therefore we conclude that rumor verbal-endemic equilibrium is unstable since some of the eigenvalue are greater than zero. Furthermore, once a rumor spreads, the truth is at risk of being distorted in the public sphere, therefore it must spread first before it decline naturally or by counter rumor.

Rumor media-endemic equilibrium

To obtain the rumor media endemic equilibrium we use the same approach as in rumor verbal-endemic equilibrium.

From (1.11), we have

I Λ Λ = 2 ( μ + α 1 ) θ β 2 (1.18)

G Λ Λ = 0 (1.19)

Substituting (1.18) and (1.19) in (1.9), we have

M Λ Λ = π θ β 2 2 μ ( μ + α 2 ) β 1 θ ( μ + α 2 ) + 2 σ ( μ + α 2 ) (1.20)

Using (1.18), (1.19) and (1.20) in (1.12), we obtained

R Λ Λ = ( 2 σ ( μ + α 1 ) + α 1 θ β 1 θ β 1 μ ) ( π θ β 1 2 μ ( μ + α 1 ) β 2 θ ( μ + α 1 ) + 2 σ μ ( μ + α 1 ) ) (1.21)

Thus E 2 = ( I Λ Λ , M Λ Λ , G Λ Λ , R Λ Λ )

Local Stability of Rumor Media-endemic Equilibrium

We also use the Jacobin stability approach to prove the stability of the rumor media endemic. The Jacobian of (1.9) to (1.12) at equilibrium point E 2 = ( I Λ Λ , M Λ Λ , G Λ Λ , R Λ Λ ) is

( μ θ 2 β 1 M Λ Λ θ 2 β 2 G Λ Λ ( M Λ Λ + G Λ Λ ) σ θ 2 β 1 I Λ Λ σ I Λ Λ θ 2 β 2 I Λ Λ σ I Λ Λ 0 θ 2 β 1 M Λ Λ θ 2 β 1 I Λ Λ μ α 1 0 0 θ 2 β 2 G Λ Λ 0 θ 2 β 2 I Λ Λ μ α 2 0 ( M Λ Λ + G Λ Λ ) σ α 1 + σ I Λ Λ α 2 + σ I Λ Λ μ ) (1.22)

If we evaluate (1.22) at E 2 and find the eigenvalue using characteristic equestion | J ( E 1 ) λ I | = 0 , we will obtain

λ 1 = μ ( θ β 1 + 2 σ 2 ) ( π θ β 2 2 ( μ + α 2 ) θ β 1 ( μ + α 2 ) + 2 σ ( μ + α 2 ) )

λ 2 = β 1 ( μ + α 2 ) β 2 μ α 1

λ 3 = 0

λ 4 = μ

Since some of the eigenvalues are not less than zero, we conclude that rumor through media-endemic equilibrium is unstable. Consequently, when rumor occurs in an ignorant population it must spread first before it declines either naturally or by counter rumor.

6. Numerical Solution and Discussion of Results

This section deals with the numerical studies of the developed model and discussion of the simulation result using estimated parameters and adopted some from other models.

Model Simulation

We performed some numerical experiments using ode45 function from MATLAB to study the behavior of the systems on the ignorant, spreaders through media (spreader 1), spreader through verbal communication (spreader 2) and stiflers population. The parameter values used in the simulation of this model are presented in Table 2.

The results of the simulation of the model (1.1) to (1.4) with parameter value in Table 2 are presented in the figures below: (Figures 2-7).

Experiment 1: Investigating the impact of contact rates between the ignorant and spreaders.

Experiment 2: Impact of rate at which spreaders become stiflers.

Figure 2. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when β 1 = 0.8 and β 2 = 0.5 .

Figure 3. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when β 1 = 0.8 and β 2 = 0.3 .

Figure 4. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when β 1 = 0.8 , β 2 = 0.3 .

Figure 5. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population α 1 = 0.122 and α 2 = 0.120 .

Figure 6. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when α 1 = 0.32 , α 2 = 0.29 .

Figure 7. Time evolution of ignorant, spreader 1, Spreader 2 and Stifler population when α 1 = 0 , α 2 = 0 .

Table 2. Estimated and adopted parameter values and variables used in the simulation of the model.

7. Discussion of Results

Here we discuss the results obtained in the analytical studies and the numerical experiments.

Analytical Results

The modified model consists of a 4-dimensional system of ordinary differential equations. The Rumour Free Equilibrium was established for the system (1.1) to (1.4). We also obtained the reproduction number of the two spreader comportments using next generation matrix method and we established the stability of the rumor Free Equilibrium and Endemic Equilibrium of the modified model using linearization methods. We observed that all the eigenvalues have negative real parts, implying that the rumor free equilibrium is locally asymptotically stable and unstable otherwise. For the endemic equilibrium point, it is observed that some of the eigenvalues of polynomial (1.17) and (1.22) are positive and we conclude that rumor-endemic equilibrium is unstable (Figure 5-7). That is if there exists a rumor in population, it has to spread through first before it decline naturally or suppress by another rumor.

It is observed that the impact of R o M is the result of the interaction of each parameters and it is easy to decrease R o M by regulating and controlling a few parameters so also R o V . The parameter with the most positive impact on R o M is the rate of contact between the spreaders targeting the community through media and the ignorant. Similarly the most positive impact on R o V is the rate of contact between the spreaders targeting the community through verbal communication and the ignorant. While the parameter that has moderate influence is the recruitment rate. In this study, we proposed two classes of spreaders and we observed that the propagation of rumor through media and verbal communication are mutually exclusive. Therefore, the appropriate increase in the propagation of rumor through media is conducive to the elimination of the propagation of rumor through verbal communication and vice versa (Figures 2 to 7). On the other hand, the propagation of rumor through media and verbal communication coexist for long time under certain condition (Figure 4). At this point, the improvement of the propagation of the rumor through media is also conducive to the propagation of rumor through verbal communication in a certain range (Figures 2 to 7).

8. Conclusions

Deterministic model on rumor propagation was analyzed to get insight into its dynamical features and to obtain associated spreaders thresholds. Some of the theoretical findings of the study are as follows:

1) rumor spreading has a life span which depends on the type of class of the spreader [see (Figures 2-4)].

2) the analysis carried out shows that the increase in stiflers population increased the rate at which spreaders convinced.

Contribution

The model can be used to control rumor propagation by varying the parameters in the model, not only that it can also be modified to manage some other social negative behaviors like smoking, drug abuse, corruption, prostitution, gang, etc.

Conflicts of Interest

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

Cite this paper

Musa, S. and Fori, M. (2019) Mathematical Model of the Dynamics of Rumor Propagation. Journal of Applied Mathematics and Physics, 7, 1289-1303. https://doi.org/10.4236/jamp.2019.76088

References

  1. 1. Zhao, L.J., Wang, X.L., Wang, J.J., Qiu, X.Y. and Xie, W.L. (2014) Rumor-Propagation Model with Consideration of Refutation Mechanism in Homogeneouse Social Networks. Discrete Dynamics in Nature and Society, 2014, Article ID: 659273. https://doi.org/10.1155/2014/659273

  2. 2. Di Fonzo, N. and Bordia, P. (2007) Rumor, Gossip and Urban Legends. Diogenes, 54, 19-35. https://doi.org/10.1177/0392192107073434

  3. 3. Dunbar, R. (1996) Grooming, Gossip, and the Evolution of Language. Faver and Faber, London.

  4. 4. Guerin, B. and Miyazaki, Y. (2006) Analyzing Rumors, Gossip, and Urban Legends through Their Conversational Properties. The Psychological Record, 56, 23-33. https://doi.org/10.1007/BF03395535

  5. 5. Fang, J., Edward, D., Yang, C. and Naren, R. (2013) Epidemiological Modelling of News and Rumors on Twitter. Proceedings of the 7th Workshop on Social Network Mining and Analysis, Chicago, Illinois, 11 August 2013, 1-9.

  6. 6. Moreno, Y., Necovee, M. and Pacheco, A.F. (2008) Dynamics of Rumor Spreading in Complex Networks. Unpublished.

  7. 7. Zhao, L., Wang J. and Huang, R. (2015) Immunization against the Spread of Rumors in Homogenous Networks. PLoS ONE, 10, e0124978. https://doi.org/10.1371/journal.pone.0124978

  8. 8. Musa, S. and Shamaki, T.A. (2017) Dynamics of Rumors Spread with Moral Etiquette Class. National Conference of NMS, University of Makurdi, Makurdi, Nigeria.

  9. 9. Auwal, A.M. and Farouk, A.U. (2016) Rumor Spreading Model: Vulnerability and its Psychological Effect in Study Salt and Water Ebola Virus Cure Rumor of August 2014 in Nigeria. Proceedings of 56th ISERD International Conference, Mecca, Saudi Arabia, 16-17 December 2016, 30-38.

  10. 10. Benjamin, D., Mahmoud, F. and Tobias, F. (2012) Why Rumors Spread Fast in social Networks. Unpublished, Germany, 1-10.

  11. 11. Deering, S. (1988) Multicast Routing in Internetworks and Extended LANs. SIGCOMM '88 Symposium Proceedings on Communications Architectures and Protocols, Stanford, CA, USA, 16-18 August 1988, 55-64. https://doi.org/10.1145/52324.52331

  12. 12. Nekovee, M., Moreno, Y., Bianconi, G. and Marsili, M. (2007) Theory of Rumor Spreading in a Complex Social Networks. Physica A: Statistical Mechanics and its Applications, 374, 457-470. https://doi.org/10.1016/j.physa.2006.07.017

  13. 13. Yusuf, I., Abdulrahman, S. and Adamu, G. (2016) Controlling the Spread of Corruption through Social Media: Mathematical Approach. IOSR Journal of Mathematics, 12, 75-81.

  14. 14. Zhao, L., Wang, J., Chen, Y., Wang, Q. and Cui, H.X. (2012) SIHR Rumor Spreading Model in Social Networks. Physica A: Statistical Mechanics and Its Application, 2444-2453. https://doi.org/10.1016/j.physa.2011.12.008

  15. 15. Kermack, W.O. and Mckendrick, A.G. (1927) A Contribution to the Mathematical Theory of Epidemics, Proceedings of the Royal Society A, 115, 700-721. https://doi.org/10.1098/rspa.1927.0118

  16. 16. Gabriel, C. and Carter, T. (1996) Information Transmission through Human Informants. Simulation, 5, 12-24.

  17. 17. Piqueira, J.R.C., Castana, M.C. and Monteiro, L.H.A. (2004) Modeling the Spreading of HIV in Homosexual Populations with Heterogeneous Preventive Attitude. Journal of Biological Systems, 12, 439-456. https://doi.org/10.1142/S0218339004001312

  18. 18. Anderson, R.M. and May, R. (1992) Infectious Diseases of Humans: Dynamics and Control of (Paper Backed) Oxford. Oxford University Press, New York.

  19. 19. Bailey, N.T.J. (1975) The Mathematical Theory of Infectious DISEASES and Its Applications. 2nd Edition, Hafner Press, New York.

  20. 20. Frauenthal, J.C. (1980) Mathematical Modelling in Epidemiology. Springer, Berlin. https://doi.org/10.1007/978-3-642-67795-3

  21. 21. Murray, J. (2002) Mathematical Biology. 3rd Edition, Springer, New York, USA.

  22. 22. Rodrigues, H.S. (2016) Application of SIR Epidemiological Model: New Trends. International Journal of Applied Mathematics and Informatics, 10, 92-97.

  23. 23. Zhao, X.X. and Wang, J.Z. (2013) Dynamical Model about Rumor Spreading with Medium. Discrete Dynamics in Nature and Society, 2019, Article ID: 586867. https://doi.org/10.1155/2013/586867

  24. 24. Allport, G.W. and Postman, L. (1947) The Psychology of Rumor. Journal of Clinical Psychology, 3, 247.

  25. 25. Rene, E. and Marco, O. (2017) A Deterministic Mathematical Model for the Spread of Two Rumors. 18 p.

  26. 26. Rosnow, R.L. (1991) The Psychology of Rumor. American Psychologist, 46, 484-496. https://doi.org/10.1037/0003-066X.46.5.484

  27. 27. Squires, H. and Tappender, P. (2011) Mathematical Modelling and Its Application to Social Care. Methods Review, 7, 10-22.

  28. 28. Daley, D. and Kendall, D. (1965) Stochastic Rumor. Journal of Apply Probability, 42-55. https://doi.org/10.1093/imamat/1.1.42

  29. 29. Liu, Y.J., Zeng, C.M. and Luo, Y.Q. (2018) Dynamics of a New Rumor Propagation Model with the Spread of Truth. Applied Mathematics, 9, 536-549. https://doi.org/10.4236/am.2018.95038

  30. 30. Lakshimikanthan, V., Leela, S. and Martynyuk, A.A. (1989) Stability of Analysis of Nonlinear System. Mareel Dekker, New York.

  31. 31. Roberto, J. and Piqueira, C. (2010) Rumour Propagation Model: An Equilibrium Study. Mathematical Problems in Engineering, 2010, Article ID: 631357. https://doi.org/10.1155/2010/631357