# A Mathematical Modeling Analysis of Racism and Corruption Codynamics with Numerical Simulation as Infectious Diseases

Authors: Belela Samuel Kotola; Shewafera Wondimagegnhu Teklu
Journal: Computational and Mathematical Methods in Medicine
 (2022)

---

## Abstract

Racism and corruption are mind infections which affect almost all public and governmental sectors. However, we cannot find enough published literatures on mathematical model analyses of racism and corruption coexistence. In this study, we have contemplated the dynamics of racism and corruption coexistence in communities, using deterministic compartmental model to analyze and suggest proper control strategies to stakeholders. We used qualitative and comprehensive mathematical methods and analyzed both the racism model in the absence of corruption and the corruption model in the absence of racism. We have computed basic reproduction numbers by applying the next generation matrix method. The developed model has a disease-free equilibrium point that is locally asymptotically stable whenever the reproduction number is less than one. Additionally, we have done sensitivity analysis to observe the effect of the parameters on the incidence and transmission of the mind infections that deduce the transmission rates of both the racism and corruption are highly sensitive. The numerical simulation we have simulated showed that the endemic equilibrium point of racism and corruption coexistence model is locally asymptotically stable whenmaxRr,Rc>1, the effects of parameters on the basic reproduction numbers, and the effect of parameter on the infectious groups. Finally, the stakeholders must focus on minimizing the transmission rates and increasing the recovery (removed) rate for both racism and corruption action which can be considered prevention and controlling strategies.


---

## Body

## 1. Introduction

Three decades ago, David Mason, a sociologist, who defined the term “institutional racism” would remain a political catchphrase devoid of analytical rigor [1]. Marx and Engels used the term “race” to refer to a wide range of human collectives, including ethnic groups based on skin color, nations, and even social classes [2]. Nowadays, racism is defined as behaviors rooted in beliefs about the innate inferiority of others, scholars working in the Allport tradition argued that racism was more typically expressed as a perception that certain racial groups did not abide by norms of hard work and patriotism and this newer symbolic form of racism has been conceptualized as a set of attitudes acquired through socialization, as perceptions that people teach to each other through interpersonal interaction and learn through education, mass media, religious institutions, and other important sources of communication [3]. Racism is when someone treats another person unfairly because of their membership in a group or because of their opinions about that group’s members. It is also when someone has strong negative feelings against another person because of their race [4, 5]. Many of us think that racism is an act of abuse or harassment. However, it does not need to involve brutal or overawing activities. Taking the racial nickname, considering situations, and when people may be ignored by members or participants due to their clan are also racism [4, 6, 7]. Racism in today’s society breeds horrible diseases in communities including serious dangers of lingering flaws, confounding of corporate structures, and a decline in the social coherence of present associations [8–10]. Violence is the acyclic act that contains an accumulation of stress and acute violence that can be the heart of racism [5, 6, 11]. Many literatures have been done by scholars on the expansion of racism in the community [5–14].Corruption is an unpleasant act for communities in general; however, it does not trouble and upset everyone on an equal level. Although corruption is harmful to society as a whole, it frequently has a greater negative impact on existing marginalized groups and is primarily practiced in developing nations [15–17]. The illegal gratitude and abuse of public office for private gain and embezzlement of public funds can be considered corruption [18, 19]. Indeed, it is an infectious activity and dishonest behavior in the body of the society, which seems to be cancer to an economic, social, and political renaissance in the country [20–22]. This illegal act can be conducted by a person or institution entrusted with a position of authority often to acquire inappropriate benefits [21]. The willingness to act of corruption can be offered from either the receptor side or the provider side [15]. The act of corruption is found in every sector of society and can affect a small group of people (petty corruption) or affect part or the entire government (grand corruption) [18, 21]. Even if most countries including Ethiopia have anticorruption policies as well as measures that are being made to eradicate corruption, it still remains a worldwide problem among the communities [15, 23]. People who are unable to receive services they are entitled to without using bribes or contacts have practical difficulties and irritation in their daily lives as a result [24, 25]. Different studies on corruption stated in [15, 18, 22, 26–32] clarify the impact of corruption, economically, socially, and politically among the communities.Today in our globe, ethnic and racial diversity is increasing rapidly [12, 33]. Those who are already impoverished of opportunities due to racism have a great chance of being exacerbated by corruption [4]. Through our review process, we have adhered that there is a tough kinship between racism and corruption. However, in the current research study, there are only a few studies [34–36] that exhibit and examine the problems due to corruption and racism coexistence in the communities. A solid grasp of mathematics for communities in nations is essential for the advancement of science, technology, and economic growth. This is because mathematics skills are very widely essential in understanding other disciplines including social sciences, engineering, sciences, arts, and outspread to all areas of science, technology as well as business enterprises, and hence, it has been becoming a key in all sciences [37]. Li et al. [38] formulated and analyzed a nonlinear dynamical analysis and optimal control strategies for a new rumor-spreading model with comprehensive interventions. They calculated the basic reproduction number with important biological significance, and the stability of equilibriums is proved. Applying the optimal control theory, the expression of optimal control pairs is obtained. In the simulation part, they examined the optimal control under 11 control strategies and through the data analysis of incremental cost-effectiveness ratio and infection averted ratio of all control strategies and provide a flexible control strategy for the security management department. Teklu and Terefe [37] formulated and analyzed a mathematical model on the dynamics of university students’ with animosity towards mathematics with optimal control theory. They have shown that the animosity-free equilibrium point is local and global stability when the basic reproduction number is less than unity, and the animosity-dominance equilibrium points local and global stability whenever the basic reproduction number is greater than unity. They carried out numerical simulations, and from the result, they recommend that prevention and treatment control measures are the best strategies to minimize and possibly to eradicate the animosity-infection throughout the community.One of the phenomena that can be presented by the mathematical model is the impact and expansion of racism-corruption coexistence among the communities. In this study, we have reviewed a literature done by other scholars to examine the spread and transmission of racism and corruption single existence as infectious diseases with a mathematical modelling approach as stated in [8, 9, 12, 15, 18, 21–23, 27, 28, 30, 31, 39–45], some applied modeling approach for social media addiction [20], and some used modeling approach for violence [6, 7, 10, 11]. However, to the best of our knowledge, no one has developed and analyzed a mathematical model on racism-corruption coexistence among individuals in a given society. Consequently, our newly proposed study contemplates the dynamics of racism and corruption coexistence in communities, using a deterministic compartmental model to analyses and suggests proper control strategies to stakeholders. Therefore, in this newly proposed racism-corruption coexistence model, we are motivated and interested to examine this connection by constructing a mathematical model of racism-corruption coexistence contagion with controlling strategies. The structure of the rest of this study is organized as follows. In Section 2, we describe and formulate the compartmental mathematical model of racism-corruption coexistence. Section 3 is dedicated to examining model analysis including the equilibrium points, basic reproduction numbers, and the stability analysis of the submodels and the main model. Section 4 presents sensitivity analysis and numerical simulations. Finally, we carried out discussions and conclusions in Sections 5 and 6, respectively.


## 2. Mathematical Model Formulation

### 2.1. Model Descriptions and Assumptions

We have assumed that all the parameters used in this mode are nonnegative. The recruitment rate entering into the susceptible class is from birth and immigration. Moreover, we have considered that the susceptible individuals are equally likely to be corrupt and/or racialist and the corrupt and/or racialist individual compels susceptible individuals into corruption and/or racism practice(s) as they effectively interact. Upon being recovered, individuals become either susceptible or honest from the act of corruption and/or racism. Using the above basic assumptions and descriptions, we have divided the total populationN into seven distinct classes. These classes are those individuals who are susceptible to corruption or racism St, those who are corrupted Ct, those who stopped corruption R1t, those who are racist Rt, those who stopped racism R2t, and those who are both corrupted and racist C1t and those who stopped both corruption and racism at the same timeR3t.The susceptible individuals become corrupted with standard incidence rate given by(1)λCt=βNCt+ω1C1t,where ω1≥1 is the modification parameter that increases infectivity and β is the corruption transmission rate. Moreover, we have used the racism mass action incidence rate given by
(2)λRt=αRt+ω2C1t,where ω2≥1 is the modification parameters that increase infectivity and α is the racism transmission rate.Using the model assumptions and descriptions stated above, the flow chart of the racism and corruption dynamics is given byUsing Figure1, the corresponding dynamical system of coexistence transmission dynamics is given by
(3)dSdt=Λ+θ2R2+θ3R3+θ1R1−λR+λC+μS,dCdt=λCS+σ1C1–γ1+ξλR+μC,dRdt=λRS+σ2C1−ηλC+γ2+μR,dC1dt=ηλCR+ξλRC−σ2+γ3+σ1+μC1,dR1dt=γ1C−θ1+μR1,dR2dt=γ2R−θ2+μR2,dR3dt=γ3C1−θ3+μR3.Figure 1
Flow chart of the transmission dynamics whereλC and λR are given in (1) and (2), respectively.


## 2.1. Model Descriptions and Assumptions

We have assumed that all the parameters used in this mode are nonnegative. The recruitment rate entering into the susceptible class is from birth and immigration. Moreover, we have considered that the susceptible individuals are equally likely to be corrupt and/or racialist and the corrupt and/or racialist individual compels susceptible individuals into corruption and/or racism practice(s) as they effectively interact. Upon being recovered, individuals become either susceptible or honest from the act of corruption and/or racism. Using the above basic assumptions and descriptions, we have divided the total populationN into seven distinct classes. These classes are those individuals who are susceptible to corruption or racism St, those who are corrupted Ct, those who stopped corruption R1t, those who are racist Rt, those who stopped racism R2t, and those who are both corrupted and racist C1t and those who stopped both corruption and racism at the same timeR3t.The susceptible individuals become corrupted with standard incidence rate given by(1)λCt=βNCt+ω1C1t,where ω1≥1 is the modification parameter that increases infectivity and β is the corruption transmission rate. Moreover, we have used the racism mass action incidence rate given by
(2)λRt=αRt+ω2C1t,where ω2≥1 is the modification parameters that increase infectivity and α is the racism transmission rate.Using the model assumptions and descriptions stated above, the flow chart of the racism and corruption dynamics is given byUsing Figure1, the corresponding dynamical system of coexistence transmission dynamics is given by
(3)dSdt=Λ+θ2R2+θ3R3+θ1R1−λR+λC+μS,dCdt=λCS+σ1C1–γ1+ξλR+μC,dRdt=λRS+σ2C1−ηλC+γ2+μR,dC1dt=ηλCR+ξλRC−σ2+γ3+σ1+μC1,dR1dt=γ1C−θ1+μR1,dR2dt=γ2R−θ2+μR2,dR3dt=γ3C1−θ3+μR3.Figure 1
Flow chart of the transmission dynamics whereλC and λR are given in (1) and (2), respectively.


## 3. Qualitative Analysis of the Model (3)

Before we analyze the racism-corruption coexistence model (3), we need to gain some background about the racism submodel and corruption submodel expansion dynamics.


### 3.1. Racism Submodel Analysis

We have derived the mathematical model of racism in the absence of corruption from the full racism and corruption coexistence model by makingC=C1=,R1=R3=0, so that we do have the dynamical system
(4)dSdt=Λ+θ2R2−λR+μS,dRdt=λRS−γ2+μR,dR2dt=γ2R−θ2+μR2.Theorem 1 (Positivity).
The solutionsSt, Rt, and R2t of the racism dynamical system (4) are nonnegative for all timet>0.Proof.
Let us defineτ=supt>0:St>0,Rt>0andR2t>0.
Since all,St, Rt, and R2t are continuous so that we can say τ>0. If τ=+∞, then positivity holds. Nevertheless, if 0<τ<+∞, thenSt=0 or t=0R2t=0.
From the first equation of the racism model, we do havedS/dt=Λ+θ2R2−λR+μS. Then, after applying the integrating factor method with some mathematical calculations, we have obtained Sτ=gS0+g∫0τexp∫λR+μdtπλ+θ2R2dt>0, where g=exp−μτ+∫0τλR+μ>0,S0>0,R2t>0. Moreover using the definition of τ, the solution Sτ>0 so that Sτ≠0. Using the same procedure all, the solutions of the dynamical system are nonnegative.


#### 3.1.1. Racism-Free Equilibrium Point of the Submodel

Racism-free equilibrium point of the racism model in the absences of corruption is obtained by making the right-hand side of equation is equal to zero providing that the racist class is equal to zero asR=0 which gives following result. Er0=SO,RO,R20=Λ/μ,0,0.


#### 3.1.2. Reproduction Number of Racism Model in the Absences of Corruption

The reproduction number is the average number of people that become racist because of the entry of one racial person into a completely susceptible population in the absence of intervention. Moreover, reproduction number utilizes to determining the effect of the control measures and to understand the capability of the corruption to disseminate in the entire community when the control strategies are applied [18].The reproduction number of racism in the absence of corruption model denoted byRr which is manipulated by Van den Driesch, Pauline, and James Warmouth next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiErO/∂xj∂νiErO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another, and Er0 is the disease-free equilibrium point ErO=SO,ROR20=Λ/μ,0,0The general transmission matrixFix and the transition matrix Vix are given by
(5)Fix=λRS00andVix=γ2+μRθ2+μR2−γ2RμS−Λ−θ2R2.Then, we have obtained(6)F=Λαμ000,V=γ1+μ0−γ2θ2+μ,(7)FV−1=Λαγ2+μμ000.Thus the eigenvalues ofFV−1 are 0,αΛ/γ2+μμ.The reproduction number of racism in the absence of corruption model is given by(8)Rr=Λαγ2+μμ.


#### 3.1.3. Local Stability of Racism-Free Equilibrium Point

Theorem 2.
The racism-free equilibrium pointEr0=SO,ROR20=Λ/μ,0,0 of the system 4 is locally asymptotically stable if the reproduction number Rr<1, and it is unstable ifRr>1.Proof.
The Jacobean matrix at racism-free equilibrium point isJEr0 of the model 4 given by
(9)JEr0=−μ0θ20αΛμ−γ2+μ00γ2−θ2+μ⟹−μ−λ0θ20αΛμ−γ2+μ−λ00γ2−θ2+μ−λ=0.
By using Software Wolfram Mathematica, we have the eigenvaluesλ1=−μ, λ2=αΛ/μ−γ2+μ, and λ3=−θ2+μ. But λ2=αΛ/μ−γ2+μ can have the form λ2=αΛ/μαΛ/μγ2+μ−1=αΛ/μRr−1. Moreover, λ2=αΛ/μRr−1<0 if and only if Rr<1; hence, all the eigenvalues are negative which implies the racism-free equilibrium point is locally asymptotically stable if and only if Rr<1; otherwise, it is unstable.


#### 3.1.4. Global Stability of Racism-Free Equilibrium Point

Theorem 3.
The racism-free equilibrium is globally asymptotically stable ifRr<1.Proof.
To prove the global asymptotic stability (GAS) of the racism-free equilibrium point, we have used the method of Lyapunov functions.
We defined a Lyapunov functionl1 such that; l1=aR, wherea=μ/αΛ.
(10)⇒dl1dt=λRμSαΛ−R.
But we do haveλR=αR and N=S+R+R1=λ/μ.
(11)⇒dl1dt≤Rr−1Rsodl1dt<0,ifRr<1.
Moreover,dl1/dt=0 if R=0 or R=1. From this fact, we do have λ/μ,0,0 is the only singleton set in S,R,R1∈Ω1:dl1/dt=0. Therefore, by the principle of LaSalle (1976), the racism-free equilibrium point is globally asymptotically stable ifRr<1.


#### 3.1.5. Existence of Endemic Equilibrium Point of Racism Model in the Absences of Corruption

It is mandatory to be sure about number of endemic equilibrium of the model before investigating the global asymptotic stability of the disease-free equilibrium point (DFE). The endemic equilibrium point of the dynamical system of (4) is solved by making right side of the system equal to zero providing that R≠0. Suppose the endemic equilibrium point of the model is denoted byEr∗=S∗,R∗R2∗.The corresponding force of infection isλRt=αRt, and we have derived the following:
(12)λR=αΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR⇒λRλR+μθ2+μγ2+μ−θ2γ2λR=αΛλRθ2+μ.⇒λR=0 or λR+μθ2+μγ2+μ−θ2γ2λR=αΛθ2+μ after simplification and rearrangement of the terms; we have λR=0 or λR=γ2+μθ2+μ/γ2+θ2+μαΛ/μγ2+μ−1.
(13)⇒λR=0orλR=γ2+μθ2+μγ2+θ2+μRr−1⇒λR>0,ifRr>1.Therefore, there is unique endemic equilibrium point for the racism model in the absence of corruption given byEr∗=S∗,R∗,R2∗ exist when Rr>1 where
(14)S∗=Λθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λR,R∗=ΛλRθ2+μλR+μθ2+μγ2+μS−θ2γ2λR,R2∗=Λγ2γ2+μλR+μθ2+μγ2+μ−θ2γ2λR.Theorem 4.
The racism model in the absence of corruption has a unique endemic equilibrium point wheneverRr>1.Theorem 5.
The endemic equilibrium pointEr∗=S∗,R∗R2∗ is locally asymptotically stable if the Rr>1, otherwise unstable. To deduce the local stability of the endemic equilibrium point, we use the method of Jacobian matrix and Routh Hurwitz stability criteria. The corresponding Jacobian matrix of the dynamical system at the endemic equilibrium point Er∗=S∗,R∗R2∗is
(15)JEr∗=−αR∗+μ−αS∗θ2αR∗αS∗−γ2+μ00γ2−θ2+μ.Then, the characteristic equation of the above Jacobian matrix is given by(16)–αR∗+μ−λ−αS∗θ2αR∗αS∗−γ2+μ−λ00γ2−θ2+μ−λ=0,A−λBθ2CD−λ00EF−λ=0,where A=–αΛλRθ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR+μ, B=−αΛθ2+μγ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR, C=αΛλRθ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR, D=αΛθ2+μγ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR−γ2+μ, E=γ2, and F=−θ2+μ.
(17)⇒A−λD−λF−λ−BCF−λ+θ2CE=0,⇒λ3−A+D+Fλ2+AF+DF−BC+ADλ+ADF+ECθ2−BCF=0,⇒a0λ3+a1λ2+a2λ+a3=0,where a0=1,a1=−A+D+F, a2=AF+DF−BC+AD, and a3=ADF+ECθ2−BCF.Buta0=1>0 and a1=−θ2γ2λR/λR+μθ2+μγ2+μ−θ2γ2λR1−αΛ/γ2+μμ.
(18)a1=−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−Rr>0.Following the same algebraic manipulation, all the coefficients of the characteristic’s polynomial are positives wheneverRr>1. Now, we can determine the local stability of endemic equilibrium point by applying the Routh-Hurwitz criteria ona0λ3+a1λ2+a2λ+a3=0.
(19)λ3λ2λ1λ0a0a2a1a3b10c10,where
(20)b1=−1a1a0a2a1a3=−1a1a0a3−a1a2=b1=−1a1a0a3−a1a2,b1=−1a1a3−a1a2,⇒b1=λR+μθ2+μγ2+μ+θ2γ2λRθ2γ2λR1−RrαΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μαΛθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λRθ2+μ−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−RrαΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μ>0.⇒b1>0 if Rr>1.In the same procedure,(21)c1=−1b1a1a3b1b1=−1b1b1a1−a3b1=a3−a1,c1=αΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μαΛθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λRθ2+μ−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−Rr.⇒c1>0 ifRr>1.We have observed that the first column of the Routh-Hurwitz array has no sign change; thus, the endemic equilibrium point of the dynamical system is locally asymptotically stable forRr>1.


### 3.2. Mathematical Analysis of the Corruption Model in the Absences of Racism

The mathematical model of corruption in the absence of racism is obtained from the full racism and corruption coexistence model (3) by making R=R2=C1=0 so that we do have the dynamical system.
(22)dSdt=Λ+θ1R1−λC+μS,dCdt=λCS–γ1+μC,dR1dt=γ1C−θ1+μR1.Theorem 6 (Positivity of the submodel solutions).
The solutionsSt, Ct, and R1t of the dynamical system of corruption model (22) are nonnegative for all timet>0.Proof.
Let us defineτ=supt>0:St>0,Ct>0,andR1t>0.
Since all,St, Ct, and R1t are continuous so that we can say τ>0. If τ=+∞, then positivity holds. Nevertheless, if 0<τ<+∞, then St=0 or t=0 andR1t=0.
From the first equation of the racism model, we do havedS/dt=Λ+θ1R1−λC+μS.
Then, applying the integrating factor method with some mathematical calculations, we have obtainedSτ=fS0+f∫0τexp∫λc+μdtΛ+θ1R1dt>0, where f=exp−μτ+∫0τλc+μ>0,S0>0,R1t>0. Moreover, using the definition of τ, the solution Sτ>0 so that Sτ≠0. Using the same procedure, all the solutions of the dynamical system are nonnegative.


#### 3.2.1. Corruption-Free Equilibrium Point

Corruption-free equilibrium point of the corruption model in the absences of racism is obtained by making the right-hand side of equation equal to zero providing that the corrupted class is equal to zero asC=0 which gives resultEcO=SO,CO,R1O=Λ/μ,0,0.


#### 3.2.2. Reproduction Number of Corruption Model in the Absences of Racism

The reproduction number of corruption in the absence of racism model denoted byRc which is manipulated by next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiEcO/∂xj∂νiEcO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another where 1≤i,j≤m, m is the number of infected compartments, and Ec0 is the disease-free equilibrium point EO=SO,COR1O=Λ/μ,0,0.The general transition matrixFix and the transmission matrix Vix are given by
(23)Fix=λCS00andVix=γ1+μRθ1+μR1−γ1RμS−Λ−θ1R1.Then, we have obtained(24)F=β000,V=γ1+μ0−γ1θ1+μ,(25)FV−1=βγ1+μ000.Thus, the eigenvalues ofFV−1 are 0,β/γ1+μ.Therefore, the reproduction number of the corruption model in the absence of racism is given byRC=β/γ1+μ.


#### 3.2.3. Local Stability of Corruption-Free Equilibrium Point

Theorem 7.
The racism-free equilibrium pointECO=SO,COR1O=Λ/μ,0,0 of the system (22) is locally asymptotically stable if the reproduction number Rc<1, and it is unstable ifRc>1.Proof.
The Jacobean matrix at corruption-free equilibrium point isJEcO of the model (3) is given by
(26)JErO=−μ0θ10β−γ1+μ00γ1−θ1+μ,⟹−μ−λ0θ10β−γ1+μ−λ00γ1−θ1+μ−λ=0.
Using Software Wolfram Mathematica, we have obtained the eigenvaluesλ1=−μ, λ2=β−γ1+μ, and λ3=θ1+μ.
Butλ2=β−γ1+μ can have the form λ2=γ2+μβ/γ2+μ−1=γ2+μRc−1. Moreover, λ2=γ2+μRc−1<0 if and only if Rc<1; hence, all the eigenvalues are negative which implies the disease-free equilibrium point is locally asymptotically stable if and only if Rc<1; otherwise, it is unstable.


#### 3.2.4. Global Stability of Corruption-Free Equilibrium Point

Theorem 8.
The corruption-free equilibrium is globally asymptotically stable ifRc<1.Proof.
To prove the global asymptotic stability (GAS) of the corruption-free equilibrium point, we have used the method of Lyapunov functions.
We defined a Lyapunov functionl3 such that l3=bC whereb=μ/Λ, ⇒dl3/dt=λCμS/Λ−C.
But we do haveλC=βC/N and N=S+C+R2=λ/μ, ⇒dl2/dt≤Rr−1C, and thus dl2/dt<0 if Rc<1.
Moreover,dl2/dt=0 if C=0 or Rc=1. From this fact, we do have λ/μ,0,0 is the only singleton set in S,C,R2∈Ω:dl2/dt=0. Therefore, by the principle of LaSalle (1976), racism-free equilibrium point is globally asymptotically stable ifRC<1.


#### 3.2.5. Existence and Uniqueness of Endemic Equilibrium Point of Corruption Model in the Absences of Racism

It is crucial to be sure about the number of endemic equilibrium of the model before investigating the global asymptotic stability of the DFE. The endemic equilibrium point of the dynamical system of (22) is solved by making right side of the system equal to zero providing that C≠0. Suppose the endemic equilibrium point of the model is denoted by EC∗=S∗,C∗R1∗. The corresponding force of infection is λCt=β/NCt, and we have derived the following.
(27)λC=βCγ1+μ/λCCγ1/θ1+μ+C+Cγ1/θ1+μ,⇒γ1γ1+μθ1+μ+γ1λCθ1+μ+λC=β.After some algebraic simplification and rearrangement of the terms, we have(28)λC=γ1γ1+μθ1+μβγ2+μ−1,⇒λC=γ1γ1+μθ1+μRC−1.Therefore, there is unique endemic equilibrium point for corruption model in the absence of racism given byEc∗=S∗,C∗,R1∗, where
(29)S∗=Λθ1+μγ1+μλC+μθ1+μγ1+μ−θ1γ1λC,C∗=ΛλCθ1+μλC+μθ1+μγ1+μ−θ1γ1λC,R1∗=Λγ1γ1+μλC+μθ1+μγ1+μ−θ1γ1λC.


#### 3.2.6. Local Stability of Endemic Equilibrium Point of the Corruption Model in the Absence of Racism

Theorem 8.
The endemic equilibrium pointEc∗=S∗,C∗R1∗ is locally asymptotically stable if the RC>1, otherwise unstable.Proof.
See the Appendix.


### 3.3. The Racism and Corruption Coexistence Model Analysis

#### 3.3.1. Basic Properties of the Coexistence Model (3)

The mathematical modeling is the representation of real world phenomena that can be demonstrated by dealing with different quantitative and qualitative attributes. In this newly extended model, we have represented human populations, which cannot be negative. Therefore, we need to show that all the state variables in our model are always nonnegative as well as the solutions of the dynamical system remains positive with positive initial conditions in the bounded region given byΩ=S,C,R,C1,R1,R2,R3∈ℝ+7,N≤Λ/μ.Theorem 9 (Positivity of the model solutions).
The solutionsSt, Ct, Rt, C1t,R1t, R2t,and R3t of the racism and corruption coexistence model (22) are nonnegative for all timet>0.Proof.
By defining,τ=supt>0:St>0,Ct>0,Rt>0,C1t>0,R1t>0,R2t>0,andR3t>0.
AllSt, Ct, Rt, C1t,R1t, R2t,and R3t are continuous so that we can deduce that τ>0. If τ=+∞, then positivity holds. However, if 0<τ<+∞, then St=0 or Ct=0 or Rt=0 or C1t=0 or R1t=0 or R2t=0 orR3t=0.
From the first equation of the racism and corruption coexistence model, we do havedS/dt=Λ+θ2R2+θ3R3+θ1R1−λR+λC+μS, and applying the integrating factor method and after some calculations, we have obtained Sτ=hS0+h∫0τexp∫λR+λC+μdtΛ+θ2R2+θ3R3+θ1R1dt>0, where h=exp−μτ+∫0τλR+λC+μ>0, S0>0, R2t>0, R3t>0 and R1>0. Finally, using the definition of τ, the solution Sτ>0 so that Sτ≠0.Following the same procedure all the solutions of the dynamical system are nonnegative.Theorem 10 (Boundedness of the model solutions).
The regionΩ=S,C,R,C1,R1,R2,R3∈ℝ+7,N≤Λ/μ is bounded inℝ+7.Proof.
Since all the state variables are nonnegative in the absence of infections, we have obtaineddN/dt≤Λ−μN. By incorporating standard comparison theorem, we have obtained ∫dN/Λ−μN≤∫dt and integrating both sides gives −1/μlnΛ−μN≤t+c, where c is some constant, and after some mathematical calculations and simplifications, we have obtained 0≤Nt≤Λ/μ. This result implies all the possible solutions of the given dynamical system with positive initial conditions are bounded.


#### 3.3.2. Racism and Corruption-Free Equilibrium Point of the Model (3)

The racism and corruption-free equilibrium point of the model is obtained by making the right-hand side of the system (3) is equal to zero providing that all the infected classes are zero as C=R=C1=0 given byE0=S0,C0,R0,C10,R10,R20,R30=Λ/μ,0,0,0,0,0,0.


#### 3.3.3. Basic Reproduction Number of the Coexistence Model (3)

The reproduction number of racism and corruption coexistence model denoted byRcr the Van den Driesch, Pauline, and James Warmouth next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiEO/∂xj∂νiEO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another where 1≤i,j≤m, m is the number of infected compartments, and Ech0 is the disease-free equilibrium point EO=S0,C0,R0,C10,R10R20,R30=Λ/μ,0,0,0,0,0,0.The general transmission matrixFix and the transition matrix Vix are given by
(30)Fix=λCSλRS0000andVix=γ1+ξλR+μC−σ1C1ηλC+γ2+μR−σ2C1σ2+γ3+σ1+μC1−ηλCR−ξλRCθ1+μR1−γ1Cθ2+μR2−γ2Rθ3+μR3−γ3C1.Then, after some calculations, we have obtained(31)F=β0βω10ΛαμΛαω2μ000,V=γ1+μ0−σ10γ2+μ−σ200σ2+γ3+σ1+μ,FV−1=βA0−Bβω1AE0ΛαCμ−DΛαλω2CEμ000,whereA=γ1+μ, B=−σ1, C=γ2+μ, D=−σ2, and E=σ2+γ3+σ1+μ.The eigenvalues ofFV−1 are 0,αΛ/γ2+μμ,β/γ1+μ. Thus, the reproduction number of racism and corruption coexistence denoted by Rrc and ismaxΛα/γ2+μμ,β/γ1+μ.That meansRrc=maxRr,Rc=maxΛα/γ2+μμ,β/γ1+μ.


#### 3.3.4. Local Stability of Coexistence-Free Equilibrium Point of the Model (3)

Theorem 11.
The racism and corruption free equilibrium pointEO=SO,CO,RO,C1O,R1O,R2O,R3O=Λ/μ,0,0,0,0,0,0 of the model is locally asymptotically stable if the reproduction number Rrc<1, and it is unstable ifRrc>1.Proof.
The Jacobean matrixJX of the model (3) is given by
(32)JX=−λR+λC+μ−βNS−αS−ω1βNS−ω1αSθ1θ2θ3λCβNS–γ1+ξλR+μ−ξαCσ1000λR−βNηRαS−ηλC+γ2+μσ20000βηRN+ξλRηλC+ξαC−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ,⟹JEO=−μ−β−απλμ−ω1β+ω1απλμθ1θ2θ30β–γ1+μ0σ100000απλμ−γ2+μσ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ.
Using Wolfram Mathematica, we have obtained the eigenvalues ofJEO as λ1=−μ, λ2=β–γ1+μ, λ3=απλ/μ−γ2+μ, λ4=−σ2+γ3+σ1+μ, λ5=−θ1+μ, λ6=−θ2+μ,λ7=−θ3+μ.
Butλ2 and λ3 can be rewritten as follows: λ2=β–γ1+μ=γ1+μβ/γ1+μ−1 and λ3=αΛ/μ−γ2+μ=γ2+μΛα/γ2+μμ−1.
Hence, all the eigenvalues are negative ifRrc<1. Therefore, the racism and corruption coexistence free equilibrium point is locally asymptotically stable if and only if Rrc<1; otherwise, it is unstable.


#### 3.3.5. Existence of Racism-Corruption Coexistence Equilibrium Point

The racism-corruption coexistence endemic equilibrium point of the full model (3) is denoted by E∗=S∗,C∗,R∗,C1∗,R1∗,R2∗,R3∗ which occurs when the mind infection persist in the community, and we computed by making the right hand side of the model as zero and obtained as
(33)S∗=Λ+θ1R1∗+θ2R2∗+θ3R3∗λR∗+λC∗+μ,C∗=λC∗S∗+σ1C1∗γ1+ξλR∗+μ,C1∗=ηλC∗R∗+ξλR∗C∗σ2+γ3+σ1+μ,R∗=λR∗S∗+σ2C1∗ηλC∗+γ2+μ,R1∗=γ1C∗θ1+μ,R2∗=γ2R∗θ2+μ,R3∗=γ3C1∗θ3+μ.From the analysis of the corruption only submodel (22) and the racism only submodel (4), we have shown that there is no endemic equilibrium point if Rc<1 and Rr<1, respectively, implying that there is no endemic equilibrium point if Rrc<1 for the coexistence model (3); in other words, the racism-corruption coexistence free equilibrium point is globally stable if Rrc<1.The summary of the racism-corruption mind infection persistence equilibrium points: The explicit computation of the mind infection persistence equilibrium point of the coinfection model (3) in terms of model parameters is difficult analytically since the system is highly nonlinear; however, the model (3) endemic equilibriums corresponds to
(1)
E1∗=S∗,0,R∗,0,R1∗,0,0,if Rr>1 is the corruption-free (racism persistence) equilibrium point. The analysis of the equilibrium E1∗ is similar to the endemic equilibrium Er∗ in the model (2)(2)
E2∗=S∗,C∗,0,0,0,R2∗,0, if Rc>1 is the racism-free (corruption persistence) equilibrium point. The analysis of the equilibrium E2∗ is similar to the endemic equilibrium Ec∗ in equations (3)(3)
Erc∗=S∗,C∗,R∗,C1∗,R1∗,R2∗,R3∗ is the racism-corruption coexistence persistence equilibrium point


#### 3.3.6. Bifurcation Analysis of the Racism-Corruption Coexistence Model

In this section, we apply the center manifold theory given by Theorem2 of Castillo-Chavez and Song [47] to ascertain the local asymptotic stability of the endemic equilibrium due to the convolution of the first approach (eigenvalues of the Jacobian). To make use of the center manifold theory, the following change of variables is made by symbolizing: S=x1, C=x2,R=x3, C1=x4, and R1=x5,R2=x6, and R3=x7 such that N=x1+x2+x3+x4+x5+x6+x7. Furthermore, by using vector notation X=x1,x2,x3,x4,x5,,x6,x7T, the Racism-Corruption coexistence model (3) can be written in the form dX/dt=FX with F=f1,f2,f3,f4,f5,f6,f7T, as follows:
(34)dx1dt=f1=Λ+θ1x5+θ2x6+θ3x7−λR+λC+μx1,(35)dx2dt=f2=λCx1+σ1x4–γ1+ξλR+μx2,(36)dx4dt=ηλCx3+ξλRx2−σ2+γ3+σ1+μx4,(37)dx5dt=γ1x2−θ1+μx5,(38)dx6dt=γ2x3−θ2+μx6,withλR=αx3+ω2x4 and λC=β/Nx2+ω1x4.Then, the method entails evaluating the Jacobian of the system (34) at the DFE point ECR0, denoted by JECR0, and this gives us
(39)JECR0=−μ−β−αx10−βω1+αω2x10θ1θ2θ30β−γ1+μ0βω1+σ100000αx10−γ2+μαω2x10+σ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200θ2+μ0000γ300θ3+μ.ConsiderRCR=1, and suppose that β=β∗is chosen as a bifurcation parameter.FromRCR=1 as RCR=β/γ1+μ=1 and solving for β, we have obtainedβ=β∗=γ1+μ.
(40)Jβ∗=−μ−β∗−αx10−β∗ω1+αω2x10θ1θ2θ30β∗−γ1+μ0βω1+σ100000αx10−γ2+μαω2x10+σ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ.After some steps of the calculation, we have obtained the eigenvalues ofJβ∗ as λ1=−μ, λ2=0 or λ3=αx10−γ2+μ=αΛ/μ−γ2+μ=γ2+μαΛ/μγ2+μ−1=γ2+μRR−1<0 if RR<1 or λ4=−σ2+γ3+σ1+μ or λ5=−θ1+μ or λ6=−θ2+μ, and λ7=−θ3+μ.It follows that the JacobianJECR0 of equation (34) at the disease-free equilibrium with β=β∗, denoted by Jβ∗, has a simple zero eigenvalue with all the remaining eigenvalues have negative real part. Hence, the Centre Manifold theory given by Theorem 2 of Castillo-Chavez and Song [47] can be used to analyze the dynamics of the model (3). In particular, it will be used to show that the model (34) undergoes backward bifurcation atRCR=1.In eigenvectors ofJβ∗, for the case RCR=1, it can be shown that the Jacobian of the system (34) at β=β∗ (denoted by Jβ∗) has right eigenvectors associated with the zero eigenvalue given by u=u1,u2,u3,u4,u5,u6,u7T as
(41)−μ−β∗−αx10M1θ1θ2θ30β∗−γ1+μ0β+σ100000M2αω2x10+σ2000000M30000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μu1u2u3u4u5u6u7=0000000,where M1=−β∗ω1−αω2x10, M2=αx10−γ2+μ, and M3=−σ2+γ3+σ1+μ.Then solving equation (41) the right eigenvectors associated with the zero eigenvalue are given by
(42)u1=−θ1+μβ∗u2+θ1γ1u2,μθ1+μ<0,u2=u2>0,u3=0,u4=0,u5=γ1θ1+μu2>0,u6=0,u7=0.Similarly, the left eigenvector associated with the zero eigenvalues atβ=β∗ given by v=v1,v2,v3,v4,v5,v6,v7T as
(43)v1v2v3v4v5v6v7T∗−μβ∗αx10M1θ1θ2θ30β∗−γ1+μ0β+σ100000M2αω2x10+σ2000000M30000γ100−θ1+μ0000γ200θ2+μ0000γ300θ3+μ=00000000.Then, solving equation (43) the left eigenvectors associated with the zero eigenvalue are given by
(44)v1=v3=v4=v5=v6=v7=0andv2=v2>0.After some steps of calculations, the bifurcation coefficientsa and b are obtained as
(45)a=∑i,j,k=17v2uiuj∂2f2∂xi∂xj=2v2u1u2∂2f2∂x1∂x2+2v2u3u2∂2f2∂x3∂x2=2v2u2u1∂2f2∂x1∂x2+u3∂2f2∂x3∂x2<2v2u2u1β+ξαu3=−2v2u2βθ1+μβ∗u2+θ1γ1u2,μθ1+μ<0.Thus, the coefficienta is negative.Moreover,(46)b=∑i,k=17vkui∂2fk∂xi∂βECR0=∑i=17v2ui∂2f2∂xi∂β=v2u2∂2f2∂x2∂β=βv2u2>0.Sinceu2>0 and v2>0 then a<0 and b>0, the Racism-Corruption coexistence model (3) exhibits the forward bifurcation, which occurs at RCR=1. That is, if RCR<1, then there is no occurrence of coexistence endemic equilibrium and the coexistence free equilibrium is the only local attractor. But if RCR>1, then the coexistence endemic equilibrium exists. For this reason, there is a forward bifurcation because in the neighborhood of the bifurcation point, the coexistence mind infection prevalence is an increasing function of RCR. Hence, from in Castillo-Chavez and Song [47], the Racism-Corruption coexistence model (3) endemic equilibrium is locally asymptotically stable whenever RCR>1.


## 3.1. Racism Submodel Analysis

We have derived the mathematical model of racism in the absence of corruption from the full racism and corruption coexistence model by makingC=C1=,R1=R3=0, so that we do have the dynamical system
(4)dSdt=Λ+θ2R2−λR+μS,dRdt=λRS−γ2+μR,dR2dt=γ2R−θ2+μR2.Theorem 1 (Positivity).
The solutionsSt, Rt, and R2t of the racism dynamical system (4) are nonnegative for all timet>0.Proof.
Let us defineτ=supt>0:St>0,Rt>0andR2t>0.
Since all,St, Rt, and R2t are continuous so that we can say τ>0. If τ=+∞, then positivity holds. Nevertheless, if 0<τ<+∞, thenSt=0 or t=0R2t=0.
From the first equation of the racism model, we do havedS/dt=Λ+θ2R2−λR+μS. Then, after applying the integrating factor method with some mathematical calculations, we have obtained Sτ=gS0+g∫0τexp∫λR+μdtπλ+θ2R2dt>0, where g=exp−μτ+∫0τλR+μ>0,S0>0,R2t>0. Moreover using the definition of τ, the solution Sτ>0 so that Sτ≠0. Using the same procedure all, the solutions of the dynamical system are nonnegative.


### 3.1.1. Racism-Free Equilibrium Point of the Submodel

Racism-free equilibrium point of the racism model in the absences of corruption is obtained by making the right-hand side of equation is equal to zero providing that the racist class is equal to zero asR=0 which gives following result. Er0=SO,RO,R20=Λ/μ,0,0.


### 3.1.2. Reproduction Number of Racism Model in the Absences of Corruption

The reproduction number is the average number of people that become racist because of the entry of one racial person into a completely susceptible population in the absence of intervention. Moreover, reproduction number utilizes to determining the effect of the control measures and to understand the capability of the corruption to disseminate in the entire community when the control strategies are applied [18].The reproduction number of racism in the absence of corruption model denoted byRr which is manipulated by Van den Driesch, Pauline, and James Warmouth next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiErO/∂xj∂νiErO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another, and Er0 is the disease-free equilibrium point ErO=SO,ROR20=Λ/μ,0,0The general transmission matrixFix and the transition matrix Vix are given by
(5)Fix=λRS00andVix=γ2+μRθ2+μR2−γ2RμS−Λ−θ2R2.Then, we have obtained(6)F=Λαμ000,V=γ1+μ0−γ2θ2+μ,(7)FV−1=Λαγ2+μμ000.Thus the eigenvalues ofFV−1 are 0,αΛ/γ2+μμ.The reproduction number of racism in the absence of corruption model is given by(8)Rr=Λαγ2+μμ.


### 3.1.3. Local Stability of Racism-Free Equilibrium Point

Theorem 2.
The racism-free equilibrium pointEr0=SO,ROR20=Λ/μ,0,0 of the system 4 is locally asymptotically stable if the reproduction number Rr<1, and it is unstable ifRr>1.Proof.
The Jacobean matrix at racism-free equilibrium point isJEr0 of the model 4 given by
(9)JEr0=−μ0θ20αΛμ−γ2+μ00γ2−θ2+μ⟹−μ−λ0θ20αΛμ−γ2+μ−λ00γ2−θ2+μ−λ=0.
By using Software Wolfram Mathematica, we have the eigenvaluesλ1=−μ, λ2=αΛ/μ−γ2+μ, and λ3=−θ2+μ. But λ2=αΛ/μ−γ2+μ can have the form λ2=αΛ/μαΛ/μγ2+μ−1=αΛ/μRr−1. Moreover, λ2=αΛ/μRr−1<0 if and only if Rr<1; hence, all the eigenvalues are negative which implies the racism-free equilibrium point is locally asymptotically stable if and only if Rr<1; otherwise, it is unstable.


### 3.1.4. Global Stability of Racism-Free Equilibrium Point

Theorem 3.
The racism-free equilibrium is globally asymptotically stable ifRr<1.Proof.
To prove the global asymptotic stability (GAS) of the racism-free equilibrium point, we have used the method of Lyapunov functions.
We defined a Lyapunov functionl1 such that; l1=aR, wherea=μ/αΛ.
(10)⇒dl1dt=λRμSαΛ−R.
But we do haveλR=αR and N=S+R+R1=λ/μ.
(11)⇒dl1dt≤Rr−1Rsodl1dt<0,ifRr<1.
Moreover,dl1/dt=0 if R=0 or R=1. From this fact, we do have λ/μ,0,0 is the only singleton set in S,R,R1∈Ω1:dl1/dt=0. Therefore, by the principle of LaSalle (1976), the racism-free equilibrium point is globally asymptotically stable ifRr<1.


### 3.1.5. Existence of Endemic Equilibrium Point of Racism Model in the Absences of Corruption

It is mandatory to be sure about number of endemic equilibrium of the model before investigating the global asymptotic stability of the disease-free equilibrium point (DFE). The endemic equilibrium point of the dynamical system of (4) is solved by making right side of the system equal to zero providing that R≠0. Suppose the endemic equilibrium point of the model is denoted byEr∗=S∗,R∗R2∗.The corresponding force of infection isλRt=αRt, and we have derived the following:
(12)λR=αΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR⇒λRλR+μθ2+μγ2+μ−θ2γ2λR=αΛλRθ2+μ.⇒λR=0 or λR+μθ2+μγ2+μ−θ2γ2λR=αΛθ2+μ after simplification and rearrangement of the terms; we have λR=0 or λR=γ2+μθ2+μ/γ2+θ2+μαΛ/μγ2+μ−1.
(13)⇒λR=0orλR=γ2+μθ2+μγ2+θ2+μRr−1⇒λR>0,ifRr>1.Therefore, there is unique endemic equilibrium point for the racism model in the absence of corruption given byEr∗=S∗,R∗,R2∗ exist when Rr>1 where
(14)S∗=Λθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λR,R∗=ΛλRθ2+μλR+μθ2+μγ2+μS−θ2γ2λR,R2∗=Λγ2γ2+μλR+μθ2+μγ2+μ−θ2γ2λR.Theorem 4.
The racism model in the absence of corruption has a unique endemic equilibrium point wheneverRr>1.Theorem 5.
The endemic equilibrium pointEr∗=S∗,R∗R2∗ is locally asymptotically stable if the Rr>1, otherwise unstable. To deduce the local stability of the endemic equilibrium point, we use the method of Jacobian matrix and Routh Hurwitz stability criteria. The corresponding Jacobian matrix of the dynamical system at the endemic equilibrium point Er∗=S∗,R∗R2∗is
(15)JEr∗=−αR∗+μ−αS∗θ2αR∗αS∗−γ2+μ00γ2−θ2+μ.Then, the characteristic equation of the above Jacobian matrix is given by(16)–αR∗+μ−λ−αS∗θ2αR∗αS∗−γ2+μ−λ00γ2−θ2+μ−λ=0,A−λBθ2CD−λ00EF−λ=0,where A=–αΛλRθ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR+μ, B=−αΛθ2+μγ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR, C=αΛλRθ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR, D=αΛθ2+μγ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR−γ2+μ, E=γ2, and F=−θ2+μ.
(17)⇒A−λD−λF−λ−BCF−λ+θ2CE=0,⇒λ3−A+D+Fλ2+AF+DF−BC+ADλ+ADF+ECθ2−BCF=0,⇒a0λ3+a1λ2+a2λ+a3=0,where a0=1,a1=−A+D+F, a2=AF+DF−BC+AD, and a3=ADF+ECθ2−BCF.Buta0=1>0 and a1=−θ2γ2λR/λR+μθ2+μγ2+μ−θ2γ2λR1−αΛ/γ2+μμ.
(18)a1=−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−Rr>0.Following the same algebraic manipulation, all the coefficients of the characteristic’s polynomial are positives wheneverRr>1. Now, we can determine the local stability of endemic equilibrium point by applying the Routh-Hurwitz criteria ona0λ3+a1λ2+a2λ+a3=0.
(19)λ3λ2λ1λ0a0a2a1a3b10c10,where
(20)b1=−1a1a0a2a1a3=−1a1a0a3−a1a2=b1=−1a1a0a3−a1a2,b1=−1a1a3−a1a2,⇒b1=λR+μθ2+μγ2+μ+θ2γ2λRθ2γ2λR1−RrαΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μαΛθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λRθ2+μ−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−RrαΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μ>0.⇒b1>0 if Rr>1.In the same procedure,(21)c1=−1b1a1a3b1b1=−1b1b1a1−a3b1=a3−a1,c1=αΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μαΛθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λRθ2+μ−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−Rr.⇒c1>0 ifRr>1.We have observed that the first column of the Routh-Hurwitz array has no sign change; thus, the endemic equilibrium point of the dynamical system is locally asymptotically stable forRr>1.


## 3.1.1. Racism-Free Equilibrium Point of the Submodel

Racism-free equilibrium point of the racism model in the absences of corruption is obtained by making the right-hand side of equation is equal to zero providing that the racist class is equal to zero asR=0 which gives following result. Er0=SO,RO,R20=Λ/μ,0,0.


## 3.1.2. Reproduction Number of Racism Model in the Absences of Corruption

The reproduction number is the average number of people that become racist because of the entry of one racial person into a completely susceptible population in the absence of intervention. Moreover, reproduction number utilizes to determining the effect of the control measures and to understand the capability of the corruption to disseminate in the entire community when the control strategies are applied [18].The reproduction number of racism in the absence of corruption model denoted byRr which is manipulated by Van den Driesch, Pauline, and James Warmouth next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiErO/∂xj∂νiErO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another, and Er0 is the disease-free equilibrium point ErO=SO,ROR20=Λ/μ,0,0The general transmission matrixFix and the transition matrix Vix are given by
(5)Fix=λRS00andVix=γ2+μRθ2+μR2−γ2RμS−Λ−θ2R2.Then, we have obtained(6)F=Λαμ000,V=γ1+μ0−γ2θ2+μ,(7)FV−1=Λαγ2+μμ000.Thus the eigenvalues ofFV−1 are 0,αΛ/γ2+μμ.The reproduction number of racism in the absence of corruption model is given by(8)Rr=Λαγ2+μμ.


## 3.1.3. Local Stability of Racism-Free Equilibrium Point

Theorem 2.
The racism-free equilibrium pointEr0=SO,ROR20=Λ/μ,0,0 of the system 4 is locally asymptotically stable if the reproduction number Rr<1, and it is unstable ifRr>1.Proof.
The Jacobean matrix at racism-free equilibrium point isJEr0 of the model 4 given by
(9)JEr0=−μ0θ20αΛμ−γ2+μ00γ2−θ2+μ⟹−μ−λ0θ20αΛμ−γ2+μ−λ00γ2−θ2+μ−λ=0.
By using Software Wolfram Mathematica, we have the eigenvaluesλ1=−μ, λ2=αΛ/μ−γ2+μ, and λ3=−θ2+μ. But λ2=αΛ/μ−γ2+μ can have the form λ2=αΛ/μαΛ/μγ2+μ−1=αΛ/μRr−1. Moreover, λ2=αΛ/μRr−1<0 if and only if Rr<1; hence, all the eigenvalues are negative which implies the racism-free equilibrium point is locally asymptotically stable if and only if Rr<1; otherwise, it is unstable.


## 3.1.4. Global Stability of Racism-Free Equilibrium Point

Theorem 3.
The racism-free equilibrium is globally asymptotically stable ifRr<1.Proof.
To prove the global asymptotic stability (GAS) of the racism-free equilibrium point, we have used the method of Lyapunov functions.
We defined a Lyapunov functionl1 such that; l1=aR, wherea=μ/αΛ.
(10)⇒dl1dt=λRμSαΛ−R.
But we do haveλR=αR and N=S+R+R1=λ/μ.
(11)⇒dl1dt≤Rr−1Rsodl1dt<0,ifRr<1.
Moreover,dl1/dt=0 if R=0 or R=1. From this fact, we do have λ/μ,0,0 is the only singleton set in S,R,R1∈Ω1:dl1/dt=0. Therefore, by the principle of LaSalle (1976), the racism-free equilibrium point is globally asymptotically stable ifRr<1.


## 3.1.5. Existence of Endemic Equilibrium Point of Racism Model in the Absences of Corruption

It is mandatory to be sure about number of endemic equilibrium of the model before investigating the global asymptotic stability of the disease-free equilibrium point (DFE). The endemic equilibrium point of the dynamical system of (4) is solved by making right side of the system equal to zero providing that R≠0. Suppose the endemic equilibrium point of the model is denoted byEr∗=S∗,R∗R2∗.The corresponding force of infection isλRt=αRt, and we have derived the following:
(12)λR=αΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR⇒λRλR+μθ2+μγ2+μ−θ2γ2λR=αΛλRθ2+μ.⇒λR=0 or λR+μθ2+μγ2+μ−θ2γ2λR=αΛθ2+μ after simplification and rearrangement of the terms; we have λR=0 or λR=γ2+μθ2+μ/γ2+θ2+μαΛ/μγ2+μ−1.
(13)⇒λR=0orλR=γ2+μθ2+μγ2+θ2+μRr−1⇒λR>0,ifRr>1.Therefore, there is unique endemic equilibrium point for the racism model in the absence of corruption given byEr∗=S∗,R∗,R2∗ exist when Rr>1 where
(14)S∗=Λθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λR,R∗=ΛλRθ2+μλR+μθ2+μγ2+μS−θ2γ2λR,R2∗=Λγ2γ2+μλR+μθ2+μγ2+μ−θ2γ2λR.Theorem 4.
The racism model in the absence of corruption has a unique endemic equilibrium point wheneverRr>1.Theorem 5.
The endemic equilibrium pointEr∗=S∗,R∗R2∗ is locally asymptotically stable if the Rr>1, otherwise unstable. To deduce the local stability of the endemic equilibrium point, we use the method of Jacobian matrix and Routh Hurwitz stability criteria. The corresponding Jacobian matrix of the dynamical system at the endemic equilibrium point Er∗=S∗,R∗R2∗is
(15)JEr∗=−αR∗+μ−αS∗θ2αR∗αS∗−γ2+μ00γ2−θ2+μ.Then, the characteristic equation of the above Jacobian matrix is given by(16)–αR∗+μ−λ−αS∗θ2αR∗αS∗−γ2+μ−λ00γ2−θ2+μ−λ=0,A−λBθ2CD−λ00EF−λ=0,where A=–αΛλRθ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR+μ, B=−αΛθ2+μγ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR, C=αΛλRθ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR, D=αΛθ2+μγ2+μ/λR+μθ2+μγ2+μ−θ2γ2λR−γ2+μ, E=γ2, and F=−θ2+μ.
(17)⇒A−λD−λF−λ−BCF−λ+θ2CE=0,⇒λ3−A+D+Fλ2+AF+DF−BC+ADλ+ADF+ECθ2−BCF=0,⇒a0λ3+a1λ2+a2λ+a3=0,where a0=1,a1=−A+D+F, a2=AF+DF−BC+AD, and a3=ADF+ECθ2−BCF.Buta0=1>0 and a1=−θ2γ2λR/λR+μθ2+μγ2+μ−θ2γ2λR1−αΛ/γ2+μμ.
(18)a1=−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−Rr>0.Following the same algebraic manipulation, all the coefficients of the characteristic’s polynomial are positives wheneverRr>1. Now, we can determine the local stability of endemic equilibrium point by applying the Routh-Hurwitz criteria ona0λ3+a1λ2+a2λ+a3=0.
(19)λ3λ2λ1λ0a0a2a1a3b10c10,where
(20)b1=−1a1a0a2a1a3=−1a1a0a3−a1a2=b1=−1a1a0a3−a1a2,b1=−1a1a3−a1a2,⇒b1=λR+μθ2+μγ2+μ+θ2γ2λRθ2γ2λR1−RrαΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μαΛθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λRθ2+μ−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−RrαΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μ>0.⇒b1>0 if Rr>1.In the same procedure,(21)c1=−1b1a1a3b1b1=−1b1b1a1−a3b1=a3−a1,c1=αΛλRθ2+μλR+μθ2+μγ2+μ−θ2γ2λR+μαΛθ2+μγ2+μλR+μθ2+μγ2+μ−θ2γ2λRθ2+μ−θ2γ2λRλR+μθ2+μγ2+μ−θ2γ2λR1−Rr.⇒c1>0 ifRr>1.We have observed that the first column of the Routh-Hurwitz array has no sign change; thus, the endemic equilibrium point of the dynamical system is locally asymptotically stable forRr>1.


## 3.2. Mathematical Analysis of the Corruption Model in the Absences of Racism

The mathematical model of corruption in the absence of racism is obtained from the full racism and corruption coexistence model (3) by making R=R2=C1=0 so that we do have the dynamical system.
(22)dSdt=Λ+θ1R1−λC+μS,dCdt=λCS–γ1+μC,dR1dt=γ1C−θ1+μR1.Theorem 6 (Positivity of the submodel solutions).
The solutionsSt, Ct, and R1t of the dynamical system of corruption model (22) are nonnegative for all timet>0.Proof.
Let us defineτ=supt>0:St>0,Ct>0,andR1t>0.
Since all,St, Ct, and R1t are continuous so that we can say τ>0. If τ=+∞, then positivity holds. Nevertheless, if 0<τ<+∞, then St=0 or t=0 andR1t=0.
From the first equation of the racism model, we do havedS/dt=Λ+θ1R1−λC+μS.
Then, applying the integrating factor method with some mathematical calculations, we have obtainedSτ=fS0+f∫0τexp∫λc+μdtΛ+θ1R1dt>0, where f=exp−μτ+∫0τλc+μ>0,S0>0,R1t>0. Moreover, using the definition of τ, the solution Sτ>0 so that Sτ≠0. Using the same procedure, all the solutions of the dynamical system are nonnegative.


### 3.2.1. Corruption-Free Equilibrium Point

Corruption-free equilibrium point of the corruption model in the absences of racism is obtained by making the right-hand side of equation equal to zero providing that the corrupted class is equal to zero asC=0 which gives resultEcO=SO,CO,R1O=Λ/μ,0,0.


### 3.2.2. Reproduction Number of Corruption Model in the Absences of Racism

The reproduction number of corruption in the absence of racism model denoted byRc which is manipulated by next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiEcO/∂xj∂νiEcO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another where 1≤i,j≤m, m is the number of infected compartments, and Ec0 is the disease-free equilibrium point EO=SO,COR1O=Λ/μ,0,0.The general transition matrixFix and the transmission matrix Vix are given by
(23)Fix=λCS00andVix=γ1+μRθ1+μR1−γ1RμS−Λ−θ1R1.Then, we have obtained(24)F=β000,V=γ1+μ0−γ1θ1+μ,(25)FV−1=βγ1+μ000.Thus, the eigenvalues ofFV−1 are 0,β/γ1+μ.Therefore, the reproduction number of the corruption model in the absence of racism is given byRC=β/γ1+μ.


### 3.2.3. Local Stability of Corruption-Free Equilibrium Point

Theorem 7.
The racism-free equilibrium pointECO=SO,COR1O=Λ/μ,0,0 of the system (22) is locally asymptotically stable if the reproduction number Rc<1, and it is unstable ifRc>1.Proof.
The Jacobean matrix at corruption-free equilibrium point isJEcO of the model (3) is given by
(26)JErO=−μ0θ10β−γ1+μ00γ1−θ1+μ,⟹−μ−λ0θ10β−γ1+μ−λ00γ1−θ1+μ−λ=0.
Using Software Wolfram Mathematica, we have obtained the eigenvaluesλ1=−μ, λ2=β−γ1+μ, and λ3=θ1+μ.
Butλ2=β−γ1+μ can have the form λ2=γ2+μβ/γ2+μ−1=γ2+μRc−1. Moreover, λ2=γ2+μRc−1<0 if and only if Rc<1; hence, all the eigenvalues are negative which implies the disease-free equilibrium point is locally asymptotically stable if and only if Rc<1; otherwise, it is unstable.


### 3.2.4. Global Stability of Corruption-Free Equilibrium Point

Theorem 8.
The corruption-free equilibrium is globally asymptotically stable ifRc<1.Proof.
To prove the global asymptotic stability (GAS) of the corruption-free equilibrium point, we have used the method of Lyapunov functions.
We defined a Lyapunov functionl3 such that l3=bC whereb=μ/Λ, ⇒dl3/dt=λCμS/Λ−C.
But we do haveλC=βC/N and N=S+C+R2=λ/μ, ⇒dl2/dt≤Rr−1C, and thus dl2/dt<0 if Rc<1.
Moreover,dl2/dt=0 if C=0 or Rc=1. From this fact, we do have λ/μ,0,0 is the only singleton set in S,C,R2∈Ω:dl2/dt=0. Therefore, by the principle of LaSalle (1976), racism-free equilibrium point is globally asymptotically stable ifRC<1.


### 3.2.5. Existence and Uniqueness of Endemic Equilibrium Point of Corruption Model in the Absences of Racism

It is crucial to be sure about the number of endemic equilibrium of the model before investigating the global asymptotic stability of the DFE. The endemic equilibrium point of the dynamical system of (22) is solved by making right side of the system equal to zero providing that C≠0. Suppose the endemic equilibrium point of the model is denoted by EC∗=S∗,C∗R1∗. The corresponding force of infection is λCt=β/NCt, and we have derived the following.
(27)λC=βCγ1+μ/λCCγ1/θ1+μ+C+Cγ1/θ1+μ,⇒γ1γ1+μθ1+μ+γ1λCθ1+μ+λC=β.After some algebraic simplification and rearrangement of the terms, we have(28)λC=γ1γ1+μθ1+μβγ2+μ−1,⇒λC=γ1γ1+μθ1+μRC−1.Therefore, there is unique endemic equilibrium point for corruption model in the absence of racism given byEc∗=S∗,C∗,R1∗, where
(29)S∗=Λθ1+μγ1+μλC+μθ1+μγ1+μ−θ1γ1λC,C∗=ΛλCθ1+μλC+μθ1+μγ1+μ−θ1γ1λC,R1∗=Λγ1γ1+μλC+μθ1+μγ1+μ−θ1γ1λC.


### 3.2.6. Local Stability of Endemic Equilibrium Point of the Corruption Model in the Absence of Racism

Theorem 8.
The endemic equilibrium pointEc∗=S∗,C∗R1∗ is locally asymptotically stable if the RC>1, otherwise unstable.Proof.
See the Appendix.


## 3.2.1. Corruption-Free Equilibrium Point

Corruption-free equilibrium point of the corruption model in the absences of racism is obtained by making the right-hand side of equation equal to zero providing that the corrupted class is equal to zero asC=0 which gives resultEcO=SO,CO,R1O=Λ/μ,0,0.


## 3.2.2. Reproduction Number of Corruption Model in the Absences of Racism

The reproduction number of corruption in the absence of racism model denoted byRc which is manipulated by next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiEcO/∂xj∂νiEcO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another where 1≤i,j≤m, m is the number of infected compartments, and Ec0 is the disease-free equilibrium point EO=SO,COR1O=Λ/μ,0,0.The general transition matrixFix and the transmission matrix Vix are given by
(23)Fix=λCS00andVix=γ1+μRθ1+μR1−γ1RμS−Λ−θ1R1.Then, we have obtained(24)F=β000,V=γ1+μ0−γ1θ1+μ,(25)FV−1=βγ1+μ000.Thus, the eigenvalues ofFV−1 are 0,β/γ1+μ.Therefore, the reproduction number of the corruption model in the absence of racism is given byRC=β/γ1+μ.


## 3.2.3. Local Stability of Corruption-Free Equilibrium Point

Theorem 7.
The racism-free equilibrium pointECO=SO,COR1O=Λ/μ,0,0 of the system (22) is locally asymptotically stable if the reproduction number Rc<1, and it is unstable ifRc>1.Proof.
The Jacobean matrix at corruption-free equilibrium point isJEcO of the model (3) is given by
(26)JErO=−μ0θ10β−γ1+μ00γ1−θ1+μ,⟹−μ−λ0θ10β−γ1+μ−λ00γ1−θ1+μ−λ=0.
Using Software Wolfram Mathematica, we have obtained the eigenvaluesλ1=−μ, λ2=β−γ1+μ, and λ3=θ1+μ.
Butλ2=β−γ1+μ can have the form λ2=γ2+μβ/γ2+μ−1=γ2+μRc−1. Moreover, λ2=γ2+μRc−1<0 if and only if Rc<1; hence, all the eigenvalues are negative which implies the disease-free equilibrium point is locally asymptotically stable if and only if Rc<1; otherwise, it is unstable.


## 3.2.4. Global Stability of Corruption-Free Equilibrium Point

Theorem 8.
The corruption-free equilibrium is globally asymptotically stable ifRc<1.Proof.
To prove the global asymptotic stability (GAS) of the corruption-free equilibrium point, we have used the method of Lyapunov functions.
We defined a Lyapunov functionl3 such that l3=bC whereb=μ/Λ, ⇒dl3/dt=λCμS/Λ−C.
But we do haveλC=βC/N and N=S+C+R2=λ/μ, ⇒dl2/dt≤Rr−1C, and thus dl2/dt<0 if Rc<1.
Moreover,dl2/dt=0 if C=0 or Rc=1. From this fact, we do have λ/μ,0,0 is the only singleton set in S,C,R2∈Ω:dl2/dt=0. Therefore, by the principle of LaSalle (1976), racism-free equilibrium point is globally asymptotically stable ifRC<1.


## 3.2.5. Existence and Uniqueness of Endemic Equilibrium Point of Corruption Model in the Absences of Racism

It is crucial to be sure about the number of endemic equilibrium of the model before investigating the global asymptotic stability of the DFE. The endemic equilibrium point of the dynamical system of (22) is solved by making right side of the system equal to zero providing that C≠0. Suppose the endemic equilibrium point of the model is denoted by EC∗=S∗,C∗R1∗. The corresponding force of infection is λCt=β/NCt, and we have derived the following.
(27)λC=βCγ1+μ/λCCγ1/θ1+μ+C+Cγ1/θ1+μ,⇒γ1γ1+μθ1+μ+γ1λCθ1+μ+λC=β.After some algebraic simplification and rearrangement of the terms, we have(28)λC=γ1γ1+μθ1+μβγ2+μ−1,⇒λC=γ1γ1+μθ1+μRC−1.Therefore, there is unique endemic equilibrium point for corruption model in the absence of racism given byEc∗=S∗,C∗,R1∗, where
(29)S∗=Λθ1+μγ1+μλC+μθ1+μγ1+μ−θ1γ1λC,C∗=ΛλCθ1+μλC+μθ1+μγ1+μ−θ1γ1λC,R1∗=Λγ1γ1+μλC+μθ1+μγ1+μ−θ1γ1λC.


## 3.2.6. Local Stability of Endemic Equilibrium Point of the Corruption Model in the Absence of Racism

Theorem 8.
The endemic equilibrium pointEc∗=S∗,C∗R1∗ is locally asymptotically stable if the RC>1, otherwise unstable.Proof.
See the Appendix.


## 3.3. The Racism and Corruption Coexistence Model Analysis

### 3.3.1. Basic Properties of the Coexistence Model (3)

The mathematical modeling is the representation of real world phenomena that can be demonstrated by dealing with different quantitative and qualitative attributes. In this newly extended model, we have represented human populations, which cannot be negative. Therefore, we need to show that all the state variables in our model are always nonnegative as well as the solutions of the dynamical system remains positive with positive initial conditions in the bounded region given byΩ=S,C,R,C1,R1,R2,R3∈ℝ+7,N≤Λ/μ.Theorem 9 (Positivity of the model solutions).
The solutionsSt, Ct, Rt, C1t,R1t, R2t,and R3t of the racism and corruption coexistence model (22) are nonnegative for all timet>0.Proof.
By defining,τ=supt>0:St>0,Ct>0,Rt>0,C1t>0,R1t>0,R2t>0,andR3t>0.
AllSt, Ct, Rt, C1t,R1t, R2t,and R3t are continuous so that we can deduce that τ>0. If τ=+∞, then positivity holds. However, if 0<τ<+∞, then St=0 or Ct=0 or Rt=0 or C1t=0 or R1t=0 or R2t=0 orR3t=0.
From the first equation of the racism and corruption coexistence model, we do havedS/dt=Λ+θ2R2+θ3R3+θ1R1−λR+λC+μS, and applying the integrating factor method and after some calculations, we have obtained Sτ=hS0+h∫0τexp∫λR+λC+μdtΛ+θ2R2+θ3R3+θ1R1dt>0, where h=exp−μτ+∫0τλR+λC+μ>0, S0>0, R2t>0, R3t>0 and R1>0. Finally, using the definition of τ, the solution Sτ>0 so that Sτ≠0.Following the same procedure all the solutions of the dynamical system are nonnegative.Theorem 10 (Boundedness of the model solutions).
The regionΩ=S,C,R,C1,R1,R2,R3∈ℝ+7,N≤Λ/μ is bounded inℝ+7.Proof.
Since all the state variables are nonnegative in the absence of infections, we have obtaineddN/dt≤Λ−μN. By incorporating standard comparison theorem, we have obtained ∫dN/Λ−μN≤∫dt and integrating both sides gives −1/μlnΛ−μN≤t+c, where c is some constant, and after some mathematical calculations and simplifications, we have obtained 0≤Nt≤Λ/μ. This result implies all the possible solutions of the given dynamical system with positive initial conditions are bounded.


### 3.3.2. Racism and Corruption-Free Equilibrium Point of the Model (3)

The racism and corruption-free equilibrium point of the model is obtained by making the right-hand side of the system (3) is equal to zero providing that all the infected classes are zero as C=R=C1=0 given byE0=S0,C0,R0,C10,R10,R20,R30=Λ/μ,0,0,0,0,0,0.


### 3.3.3. Basic Reproduction Number of the Coexistence Model (3)

The reproduction number of racism and corruption coexistence model denoted byRcr the Van den Driesch, Pauline, and James Warmouth next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiEO/∂xj∂νiEO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another where 1≤i,j≤m, m is the number of infected compartments, and Ech0 is the disease-free equilibrium point EO=S0,C0,R0,C10,R10R20,R30=Λ/μ,0,0,0,0,0,0.The general transmission matrixFix and the transition matrix Vix are given by
(30)Fix=λCSλRS0000andVix=γ1+ξλR+μC−σ1C1ηλC+γ2+μR−σ2C1σ2+γ3+σ1+μC1−ηλCR−ξλRCθ1+μR1−γ1Cθ2+μR2−γ2Rθ3+μR3−γ3C1.Then, after some calculations, we have obtained(31)F=β0βω10ΛαμΛαω2μ000,V=γ1+μ0−σ10γ2+μ−σ200σ2+γ3+σ1+μ,FV−1=βA0−Bβω1AE0ΛαCμ−DΛαλω2CEμ000,whereA=γ1+μ, B=−σ1, C=γ2+μ, D=−σ2, and E=σ2+γ3+σ1+μ.The eigenvalues ofFV−1 are 0,αΛ/γ2+μμ,β/γ1+μ. Thus, the reproduction number of racism and corruption coexistence denoted by Rrc and ismaxΛα/γ2+μμ,β/γ1+μ.That meansRrc=maxRr,Rc=maxΛα/γ2+μμ,β/γ1+μ.


### 3.3.4. Local Stability of Coexistence-Free Equilibrium Point of the Model (3)

Theorem 11.
The racism and corruption free equilibrium pointEO=SO,CO,RO,C1O,R1O,R2O,R3O=Λ/μ,0,0,0,0,0,0 of the model is locally asymptotically stable if the reproduction number Rrc<1, and it is unstable ifRrc>1.Proof.
The Jacobean matrixJX of the model (3) is given by
(32)JX=−λR+λC+μ−βNS−αS−ω1βNS−ω1αSθ1θ2θ3λCβNS–γ1+ξλR+μ−ξαCσ1000λR−βNηRαS−ηλC+γ2+μσ20000βηRN+ξλRηλC+ξαC−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ,⟹JEO=−μ−β−απλμ−ω1β+ω1απλμθ1θ2θ30β–γ1+μ0σ100000απλμ−γ2+μσ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ.
Using Wolfram Mathematica, we have obtained the eigenvalues ofJEO as λ1=−μ, λ2=β–γ1+μ, λ3=απλ/μ−γ2+μ, λ4=−σ2+γ3+σ1+μ, λ5=−θ1+μ, λ6=−θ2+μ,λ7=−θ3+μ.
Butλ2 and λ3 can be rewritten as follows: λ2=β–γ1+μ=γ1+μβ/γ1+μ−1 and λ3=αΛ/μ−γ2+μ=γ2+μΛα/γ2+μμ−1.
Hence, all the eigenvalues are negative ifRrc<1. Therefore, the racism and corruption coexistence free equilibrium point is locally asymptotically stable if and only if Rrc<1; otherwise, it is unstable.


### 3.3.5. Existence of Racism-Corruption Coexistence Equilibrium Point

The racism-corruption coexistence endemic equilibrium point of the full model (3) is denoted by E∗=S∗,C∗,R∗,C1∗,R1∗,R2∗,R3∗ which occurs when the mind infection persist in the community, and we computed by making the right hand side of the model as zero and obtained as
(33)S∗=Λ+θ1R1∗+θ2R2∗+θ3R3∗λR∗+λC∗+μ,C∗=λC∗S∗+σ1C1∗γ1+ξλR∗+μ,C1∗=ηλC∗R∗+ξλR∗C∗σ2+γ3+σ1+μ,R∗=λR∗S∗+σ2C1∗ηλC∗+γ2+μ,R1∗=γ1C∗θ1+μ,R2∗=γ2R∗θ2+μ,R3∗=γ3C1∗θ3+μ.From the analysis of the corruption only submodel (22) and the racism only submodel (4), we have shown that there is no endemic equilibrium point if Rc<1 and Rr<1, respectively, implying that there is no endemic equilibrium point if Rrc<1 for the coexistence model (3); in other words, the racism-corruption coexistence free equilibrium point is globally stable if Rrc<1.The summary of the racism-corruption mind infection persistence equilibrium points: The explicit computation of the mind infection persistence equilibrium point of the coinfection model (3) in terms of model parameters is difficult analytically since the system is highly nonlinear; however, the model (3) endemic equilibriums corresponds to
(1)
E1∗=S∗,0,R∗,0,R1∗,0,0,if Rr>1 is the corruption-free (racism persistence) equilibrium point. The analysis of the equilibrium E1∗ is similar to the endemic equilibrium Er∗ in the model (2)(2)
E2∗=S∗,C∗,0,0,0,R2∗,0, if Rc>1 is the racism-free (corruption persistence) equilibrium point. The analysis of the equilibrium E2∗ is similar to the endemic equilibrium Ec∗ in equations (3)(3)
Erc∗=S∗,C∗,R∗,C1∗,R1∗,R2∗,R3∗ is the racism-corruption coexistence persistence equilibrium point


### 3.3.6. Bifurcation Analysis of the Racism-Corruption Coexistence Model

In this section, we apply the center manifold theory given by Theorem2 of Castillo-Chavez and Song [47] to ascertain the local asymptotic stability of the endemic equilibrium due to the convolution of the first approach (eigenvalues of the Jacobian). To make use of the center manifold theory, the following change of variables is made by symbolizing: S=x1, C=x2,R=x3, C1=x4, and R1=x5,R2=x6, and R3=x7 such that N=x1+x2+x3+x4+x5+x6+x7. Furthermore, by using vector notation X=x1,x2,x3,x4,x5,,x6,x7T, the Racism-Corruption coexistence model (3) can be written in the form dX/dt=FX with F=f1,f2,f3,f4,f5,f6,f7T, as follows:
(34)dx1dt=f1=Λ+θ1x5+θ2x6+θ3x7−λR+λC+μx1,(35)dx2dt=f2=λCx1+σ1x4–γ1+ξλR+μx2,(36)dx4dt=ηλCx3+ξλRx2−σ2+γ3+σ1+μx4,(37)dx5dt=γ1x2−θ1+μx5,(38)dx6dt=γ2x3−θ2+μx6,withλR=αx3+ω2x4 and λC=β/Nx2+ω1x4.Then, the method entails evaluating the Jacobian of the system (34) at the DFE point ECR0, denoted by JECR0, and this gives us
(39)JECR0=−μ−β−αx10−βω1+αω2x10θ1θ2θ30β−γ1+μ0βω1+σ100000αx10−γ2+μαω2x10+σ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200θ2+μ0000γ300θ3+μ.ConsiderRCR=1, and suppose that β=β∗is chosen as a bifurcation parameter.FromRCR=1 as RCR=β/γ1+μ=1 and solving for β, we have obtainedβ=β∗=γ1+μ.
(40)Jβ∗=−μ−β∗−αx10−β∗ω1+αω2x10θ1θ2θ30β∗−γ1+μ0βω1+σ100000αx10−γ2+μαω2x10+σ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ.After some steps of the calculation, we have obtained the eigenvalues ofJβ∗ as λ1=−μ, λ2=0 or λ3=αx10−γ2+μ=αΛ/μ−γ2+μ=γ2+μαΛ/μγ2+μ−1=γ2+μRR−1<0 if RR<1 or λ4=−σ2+γ3+σ1+μ or λ5=−θ1+μ or λ6=−θ2+μ, and λ7=−θ3+μ.It follows that the JacobianJECR0 of equation (34) at the disease-free equilibrium with β=β∗, denoted by Jβ∗, has a simple zero eigenvalue with all the remaining eigenvalues have negative real part. Hence, the Centre Manifold theory given by Theorem 2 of Castillo-Chavez and Song [47] can be used to analyze the dynamics of the model (3). In particular, it will be used to show that the model (34) undergoes backward bifurcation atRCR=1.In eigenvectors ofJβ∗, for the case RCR=1, it can be shown that the Jacobian of the system (34) at β=β∗ (denoted by Jβ∗) has right eigenvectors associated with the zero eigenvalue given by u=u1,u2,u3,u4,u5,u6,u7T as
(41)−μ−β∗−αx10M1θ1θ2θ30β∗−γ1+μ0β+σ100000M2αω2x10+σ2000000M30000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μu1u2u3u4u5u6u7=0000000,where M1=−β∗ω1−αω2x10, M2=αx10−γ2+μ, and M3=−σ2+γ3+σ1+μ.Then solving equation (41) the right eigenvectors associated with the zero eigenvalue are given by
(42)u1=−θ1+μβ∗u2+θ1γ1u2,μθ1+μ<0,u2=u2>0,u3=0,u4=0,u5=γ1θ1+μu2>0,u6=0,u7=0.Similarly, the left eigenvector associated with the zero eigenvalues atβ=β∗ given by v=v1,v2,v3,v4,v5,v6,v7T as
(43)v1v2v3v4v5v6v7T∗−μβ∗αx10M1θ1θ2θ30β∗−γ1+μ0β+σ100000M2αω2x10+σ2000000M30000γ100−θ1+μ0000γ200θ2+μ0000γ300θ3+μ=00000000.Then, solving equation (43) the left eigenvectors associated with the zero eigenvalue are given by
(44)v1=v3=v4=v5=v6=v7=0andv2=v2>0.After some steps of calculations, the bifurcation coefficientsa and b are obtained as
(45)a=∑i,j,k=17v2uiuj∂2f2∂xi∂xj=2v2u1u2∂2f2∂x1∂x2+2v2u3u2∂2f2∂x3∂x2=2v2u2u1∂2f2∂x1∂x2+u3∂2f2∂x3∂x2<2v2u2u1β+ξαu3=−2v2u2βθ1+μβ∗u2+θ1γ1u2,μθ1+μ<0.Thus, the coefficienta is negative.Moreover,(46)b=∑i,k=17vkui∂2fk∂xi∂βECR0=∑i=17v2ui∂2f2∂xi∂β=v2u2∂2f2∂x2∂β=βv2u2>0.Sinceu2>0 and v2>0 then a<0 and b>0, the Racism-Corruption coexistence model (3) exhibits the forward bifurcation, which occurs at RCR=1. That is, if RCR<1, then there is no occurrence of coexistence endemic equilibrium and the coexistence free equilibrium is the only local attractor. But if RCR>1, then the coexistence endemic equilibrium exists. For this reason, there is a forward bifurcation because in the neighborhood of the bifurcation point, the coexistence mind infection prevalence is an increasing function of RCR. Hence, from in Castillo-Chavez and Song [47], the Racism-Corruption coexistence model (3) endemic equilibrium is locally asymptotically stable whenever RCR>1.


## 3.3.1. Basic Properties of the Coexistence Model (3)

The mathematical modeling is the representation of real world phenomena that can be demonstrated by dealing with different quantitative and qualitative attributes. In this newly extended model, we have represented human populations, which cannot be negative. Therefore, we need to show that all the state variables in our model are always nonnegative as well as the solutions of the dynamical system remains positive with positive initial conditions in the bounded region given byΩ=S,C,R,C1,R1,R2,R3∈ℝ+7,N≤Λ/μ.Theorem 9 (Positivity of the model solutions).
The solutionsSt, Ct, Rt, C1t,R1t, R2t,and R3t of the racism and corruption coexistence model (22) are nonnegative for all timet>0.Proof.
By defining,τ=supt>0:St>0,Ct>0,Rt>0,C1t>0,R1t>0,R2t>0,andR3t>0.
AllSt, Ct, Rt, C1t,R1t, R2t,and R3t are continuous so that we can deduce that τ>0. If τ=+∞, then positivity holds. However, if 0<τ<+∞, then St=0 or Ct=0 or Rt=0 or C1t=0 or R1t=0 or R2t=0 orR3t=0.
From the first equation of the racism and corruption coexistence model, we do havedS/dt=Λ+θ2R2+θ3R3+θ1R1−λR+λC+μS, and applying the integrating factor method and after some calculations, we have obtained Sτ=hS0+h∫0τexp∫λR+λC+μdtΛ+θ2R2+θ3R3+θ1R1dt>0, where h=exp−μτ+∫0τλR+λC+μ>0, S0>0, R2t>0, R3t>0 and R1>0. Finally, using the definition of τ, the solution Sτ>0 so that Sτ≠0.Following the same procedure all the solutions of the dynamical system are nonnegative.Theorem 10 (Boundedness of the model solutions).
The regionΩ=S,C,R,C1,R1,R2,R3∈ℝ+7,N≤Λ/μ is bounded inℝ+7.Proof.
Since all the state variables are nonnegative in the absence of infections, we have obtaineddN/dt≤Λ−μN. By incorporating standard comparison theorem, we have obtained ∫dN/Λ−μN≤∫dt and integrating both sides gives −1/μlnΛ−μN≤t+c, where c is some constant, and after some mathematical calculations and simplifications, we have obtained 0≤Nt≤Λ/μ. This result implies all the possible solutions of the given dynamical system with positive initial conditions are bounded.


## 3.3.2. Racism and Corruption-Free Equilibrium Point of the Model (3)

The racism and corruption-free equilibrium point of the model is obtained by making the right-hand side of the system (3) is equal to zero providing that all the infected classes are zero as C=R=C1=0 given byE0=S0,C0,R0,C10,R10,R20,R30=Λ/μ,0,0,0,0,0,0.


## 3.3.3. Basic Reproduction Number of the Coexistence Model (3)

The reproduction number of racism and corruption coexistence model denoted byRcr the Van den Driesch, Pauline, and James Warmouth next-generation matrix approach [46] is the largest eigenvalue of the next generation matrix FV−1=∂FiEO/∂xj∂νiEO/∂xj−1, where Fi is the rate of appearance of new infection in compartment i, νi is the transfer of infections from one compartment i to another where 1≤i,j≤m, m is the number of infected compartments, and Ech0 is the disease-free equilibrium point EO=S0,C0,R0,C10,R10R20,R30=Λ/μ,0,0,0,0,0,0.The general transmission matrixFix and the transition matrix Vix are given by
(30)Fix=λCSλRS0000andVix=γ1+ξλR+μC−σ1C1ηλC+γ2+μR−σ2C1σ2+γ3+σ1+μC1−ηλCR−ξλRCθ1+μR1−γ1Cθ2+μR2−γ2Rθ3+μR3−γ3C1.Then, after some calculations, we have obtained(31)F=β0βω10ΛαμΛαω2μ000,V=γ1+μ0−σ10γ2+μ−σ200σ2+γ3+σ1+μ,FV−1=βA0−Bβω1AE0ΛαCμ−DΛαλω2CEμ000,whereA=γ1+μ, B=−σ1, C=γ2+μ, D=−σ2, and E=σ2+γ3+σ1+μ.The eigenvalues ofFV−1 are 0,αΛ/γ2+μμ,β/γ1+μ. Thus, the reproduction number of racism and corruption coexistence denoted by Rrc and ismaxΛα/γ2+μμ,β/γ1+μ.That meansRrc=maxRr,Rc=maxΛα/γ2+μμ,β/γ1+μ.


## 3.3.4. Local Stability of Coexistence-Free Equilibrium Point of the Model (3)

Theorem 11.
The racism and corruption free equilibrium pointEO=SO,CO,RO,C1O,R1O,R2O,R3O=Λ/μ,0,0,0,0,0,0 of the model is locally asymptotically stable if the reproduction number Rrc<1, and it is unstable ifRrc>1.Proof.
The Jacobean matrixJX of the model (3) is given by
(32)JX=−λR+λC+μ−βNS−αS−ω1βNS−ω1αSθ1θ2θ3λCβNS–γ1+ξλR+μ−ξαCσ1000λR−βNηRαS−ηλC+γ2+μσ20000βηRN+ξλRηλC+ξαC−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ,⟹JEO=−μ−β−απλμ−ω1β+ω1απλμθ1θ2θ30β–γ1+μ0σ100000απλμ−γ2+μσ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ.
Using Wolfram Mathematica, we have obtained the eigenvalues ofJEO as λ1=−μ, λ2=β–γ1+μ, λ3=απλ/μ−γ2+μ, λ4=−σ2+γ3+σ1+μ, λ5=−θ1+μ, λ6=−θ2+μ,λ7=−θ3+μ.
Butλ2 and λ3 can be rewritten as follows: λ2=β–γ1+μ=γ1+μβ/γ1+μ−1 and λ3=αΛ/μ−γ2+μ=γ2+μΛα/γ2+μμ−1.
Hence, all the eigenvalues are negative ifRrc<1. Therefore, the racism and corruption coexistence free equilibrium point is locally asymptotically stable if and only if Rrc<1; otherwise, it is unstable.


## 3.3.5. Existence of Racism-Corruption Coexistence Equilibrium Point

The racism-corruption coexistence endemic equilibrium point of the full model (3) is denoted by E∗=S∗,C∗,R∗,C1∗,R1∗,R2∗,R3∗ which occurs when the mind infection persist in the community, and we computed by making the right hand side of the model as zero and obtained as
(33)S∗=Λ+θ1R1∗+θ2R2∗+θ3R3∗λR∗+λC∗+μ,C∗=λC∗S∗+σ1C1∗γ1+ξλR∗+μ,C1∗=ηλC∗R∗+ξλR∗C∗σ2+γ3+σ1+μ,R∗=λR∗S∗+σ2C1∗ηλC∗+γ2+μ,R1∗=γ1C∗θ1+μ,R2∗=γ2R∗θ2+μ,R3∗=γ3C1∗θ3+μ.From the analysis of the corruption only submodel (22) and the racism only submodel (4), we have shown that there is no endemic equilibrium point if Rc<1 and Rr<1, respectively, implying that there is no endemic equilibrium point if Rrc<1 for the coexistence model (3); in other words, the racism-corruption coexistence free equilibrium point is globally stable if Rrc<1.The summary of the racism-corruption mind infection persistence equilibrium points: The explicit computation of the mind infection persistence equilibrium point of the coinfection model (3) in terms of model parameters is difficult analytically since the system is highly nonlinear; however, the model (3) endemic equilibriums corresponds to
(1)
E1∗=S∗,0,R∗,0,R1∗,0,0,if Rr>1 is the corruption-free (racism persistence) equilibrium point. The analysis of the equilibrium E1∗ is similar to the endemic equilibrium Er∗ in the model (2)(2)
E2∗=S∗,C∗,0,0,0,R2∗,0, if Rc>1 is the racism-free (corruption persistence) equilibrium point. The analysis of the equilibrium E2∗ is similar to the endemic equilibrium Ec∗ in equations (3)(3)
Erc∗=S∗,C∗,R∗,C1∗,R1∗,R2∗,R3∗ is the racism-corruption coexistence persistence equilibrium point


## 3.3.6. Bifurcation Analysis of the Racism-Corruption Coexistence Model

In this section, we apply the center manifold theory given by Theorem2 of Castillo-Chavez and Song [47] to ascertain the local asymptotic stability of the endemic equilibrium due to the convolution of the first approach (eigenvalues of the Jacobian). To make use of the center manifold theory, the following change of variables is made by symbolizing: S=x1, C=x2,R=x3, C1=x4, and R1=x5,R2=x6, and R3=x7 such that N=x1+x2+x3+x4+x5+x6+x7. Furthermore, by using vector notation X=x1,x2,x3,x4,x5,,x6,x7T, the Racism-Corruption coexistence model (3) can be written in the form dX/dt=FX with F=f1,f2,f3,f4,f5,f6,f7T, as follows:
(34)dx1dt=f1=Λ+θ1x5+θ2x6+θ3x7−λR+λC+μx1,(35)dx2dt=f2=λCx1+σ1x4–γ1+ξλR+μx2,(36)dx4dt=ηλCx3+ξλRx2−σ2+γ3+σ1+μx4,(37)dx5dt=γ1x2−θ1+μx5,(38)dx6dt=γ2x3−θ2+μx6,withλR=αx3+ω2x4 and λC=β/Nx2+ω1x4.Then, the method entails evaluating the Jacobian of the system (34) at the DFE point ECR0, denoted by JECR0, and this gives us
(39)JECR0=−μ−β−αx10−βω1+αω2x10θ1θ2θ30β−γ1+μ0βω1+σ100000αx10−γ2+μαω2x10+σ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200θ2+μ0000γ300θ3+μ.ConsiderRCR=1, and suppose that β=β∗is chosen as a bifurcation parameter.FromRCR=1 as RCR=β/γ1+μ=1 and solving for β, we have obtainedβ=β∗=γ1+μ.
(40)Jβ∗=−μ−β∗−αx10−β∗ω1+αω2x10θ1θ2θ30β∗−γ1+μ0βω1+σ100000αx10−γ2+μαω2x10+σ2000000−σ2+γ3+σ1+μ0000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μ.After some steps of the calculation, we have obtained the eigenvalues ofJβ∗ as λ1=−μ, λ2=0 or λ3=αx10−γ2+μ=αΛ/μ−γ2+μ=γ2+μαΛ/μγ2+μ−1=γ2+μRR−1<0 if RR<1 or λ4=−σ2+γ3+σ1+μ or λ5=−θ1+μ or λ6=−θ2+μ, and λ7=−θ3+μ.It follows that the JacobianJECR0 of equation (34) at the disease-free equilibrium with β=β∗, denoted by Jβ∗, has a simple zero eigenvalue with all the remaining eigenvalues have negative real part. Hence, the Centre Manifold theory given by Theorem 2 of Castillo-Chavez and Song [47] can be used to analyze the dynamics of the model (3). In particular, it will be used to show that the model (34) undergoes backward bifurcation atRCR=1.In eigenvectors ofJβ∗, for the case RCR=1, it can be shown that the Jacobian of the system (34) at β=β∗ (denoted by Jβ∗) has right eigenvectors associated with the zero eigenvalue given by u=u1,u2,u3,u4,u5,u6,u7T as
(41)−μ−β∗−αx10M1θ1θ2θ30β∗−γ1+μ0β+σ100000M2αω2x10+σ2000000M30000γ100−θ1+μ0000γ200−θ2+μ0000γ300−θ3+μu1u2u3u4u5u6u7=0000000,where M1=−β∗ω1−αω2x10, M2=αx10−γ2+μ, and M3=−σ2+γ3+σ1+μ.Then solving equation (41) the right eigenvectors associated with the zero eigenvalue are given by
(42)u1=−θ1+μβ∗u2+θ1γ1u2,μθ1+μ<0,u2=u2>0,u3=0,u4=0,u5=γ1θ1+μu2>0,u6=0,u7=0.Similarly, the left eigenvector associated with the zero eigenvalues atβ=β∗ given by v=v1,v2,v3,v4,v5,v6,v7T as
(43)v1v2v3v4v5v6v7T∗−μβ∗αx10M1θ1θ2θ30β∗−γ1+μ0β+σ100000M2αω2x10+σ2000000M30000γ100−θ1+μ0000γ200θ2+μ0000γ300θ3+μ=00000000.Then, solving equation (43) the left eigenvectors associated with the zero eigenvalue are given by
(44)v1=v3=v4=v5=v6=v7=0andv2=v2>0.After some steps of calculations, the bifurcation coefficientsa and b are obtained as
(45)a=∑i,j,k=17v2uiuj∂2f2∂xi∂xj=2v2u1u2∂2f2∂x1∂x2+2v2u3u2∂2f2∂x3∂x2=2v2u2u1∂2f2∂x1∂x2+u3∂2f2∂x3∂x2<2v2u2u1β+ξαu3=−2v2u2βθ1+μβ∗u2+θ1γ1u2,μθ1+μ<0.Thus, the coefficienta is negative.Moreover,(46)b=∑i,k=17vkui∂2fk∂xi∂βECR0=∑i=17v2ui∂2f2∂xi∂β=v2u2∂2f2∂x2∂β=βv2u2>0.Sinceu2>0 and v2>0 then a<0 and b>0, the Racism-Corruption coexistence model (3) exhibits the forward bifurcation, which occurs at RCR=1. That is, if RCR<1, then there is no occurrence of coexistence endemic equilibrium and the coexistence free equilibrium is the only local attractor. But if RCR>1, then the coexistence endemic equilibrium exists. For this reason, there is a forward bifurcation because in the neighborhood of the bifurcation point, the coexistence mind infection prevalence is an increasing function of RCR. Hence, from in Castillo-Chavez and Song [47], the Racism-Corruption coexistence model (3) endemic equilibrium is locally asymptotically stable whenever RCR>1.


## 4. Sensitivity Analysis and Numerical Simulations

### 4.1. Sensitivity Analysis

The normalized forward sensitivity index of a variable racism and corruption coexistence model (3) basic reproduction number denoted by the symbol Rrc that is differentiable with a parameter ε is defined as SIp=∂Rrc/∂ε∗ε/Rrc as stated in literatures [48–50].The sensitivity indices enable us to examine the relative importance of different parameters in racism and corruption incidence and prevalence. The most sensitive parameter has the magnitude of the sensitivity index larger than that of all other parameters. We can calculate the sensitivity index in terms ofRr and Rc sinceRrc=maxRr,Rc.Using the parameter values given in Table1, the sensitivity indices are listed in Tables 2 and 3 as it is given in Table 2.Table 1
Parameter values for numerical simulation.
ParameterValuesSourceμ0.01[6, 12, 51]Λ50[6, 20, 49]αVariableAssumedσ10.6Assumedσ20.7Assumedη1.3Assumedξ1.2Assumedθ30.2AssumedβVariableAssumedθ10.3[20]θ20.25Assumedγ10.007[39, 50]γ20.006Assumedγ30.008AssumedTable 2
Sensitivity indices ofRr.
Sensitivity indexSensitivity indicesSIΛ+1SIγ2-0.85SIα+1Table 3
Sensitivity indices ofRc.
Sensitivity indexSensitivity indicesSIβ+1SIγ1-0.546SIμ-0.346SIΛ+1Table2 above is the summary sensitivity indices of Rr manipulated with parameters values from Table 1 and provides the result Rr=3.51 at the racism transmission rate α=1.37, which imply that racism spreads throughout the community. Moreover, sensitivity analysis given in Table 2 explains that the human population recruitment rate Λ and racism transmission rate α are highly affecting the racism reproduction numberRr.Here, with the given parameter values in Table3, we have computed Rc=6.03 at the corruption expansion rate β=1.51, which imply that corruption spreads throughout the community. Moreover, sensitivity analysis given in Table 3 explains that the human population recruitment rate Λ and corruption transmission rate β are highly affecting the racism reproduction numberRc.


#### 4.1.1. Sensitivity Analysis Graphical Verification

In this subsection, sensitivity analysis of the racism and corruption coexistence transmission dynamics is performed to identify the most influential parameters for the spread as well as for the control of coexistence mind infection transmission in the community. The results of sensitivity analysis based on the set of parameters values given in Table1 are analyzed in Tables 2 and 3, respectively. Graphical simulation of sensitivity indices given by Figure 2 illustrates the sensitivity index of parameters and has verified the qualitative analysis given in Tables 2 and 3, respectively. It shows that the most sensitive parameters which have direct impacts on the basic reproduction numbers are the racism transmission rate, the corruption transmission rate, and the recruitment rate and the most sensitive parameters which have an indirect impact on the basic reproduction numbers are the racism recovery rate and the corruption recovery rate. Stake holders can minimize the transmission rates and maximize recovery rates to prevent and control the coexistence expansion in the community.Figure 2
Parameter sensitivity analysis.


### 4.2. Numerical Simulations

In this section, numerical simulation has been performed with MATAB ode45 code to analyze the effect of some parameters that causes for conducting this illegal activity. Most specifically, we investigated the stability of the endemic equilibrium point of the main coexistence model (3), parameter effects on the reproduction numbers, and the impact of being honest on racism and corruption co-acting individuals in the community. Moreover, we have used the parameters values stated in Table 1 for numerical simulation.


#### 4.2.1. Behavior of the Coexistence Model Solutions forRrc>1

Figure3 illustrates the result of numerical simulation with ode45 using parameter values that are given in Table 1. From Figure 3, we can observe that after a year the solutions of the racism and corruption coexistence dynamical system (3) will be approaching to the coexistence endemic equilibrium point of the racism-corruption coexistence model depending on the value of Rrc=maxRr,Rc=max3.51,6.03=6.03>1. The numerical simulation justified that the physical phenomenon that can be stated as the expansion and spreading of racism-corruption coexistence activity is consistently occurred throughout the considered community in the study.Figure 3
Behaviors of the coexistence model solutions.


#### 4.2.2. Impact of Racism Recovery Rate onRr

Figure4 illustrates the effect of racism recovery rate γ2 on the racism reproduction number Rr. The plot demonstrates that when the value of γ2 increases, the racism reproduction number is going down, and whenever the value of γ2>0.059 implies that Rr<1. In other words, it mean as the power of treatment rate increases the number of racial individual will decrease. Therefore, the stakeholder shall concentrate on maximizing the value of racism recovery rate γ2 by applying possible interventions strategies to prevent and control the problem of racism.Figure 4
Impact ofγ2 on Rr.


#### 4.2.3. Effect of Corruption Recovery Rateγ1 on Rc

Figure5 illustrates the relation between corruption recovery rate γ1 and corruption reproduction Rc. The plot shows that that when the value of corruption recovery rate γ1 increases, the corruption reproduction number decreases, and whenever the value of γ1>0.508 implies that Rc<1. This means that as the power treatment rate increases, the number of corrupted individual decreases. Moreover, the result tells us the stakeholder shall concentrate on maximizing the values of γ1 to prevent and control the expansion of racism in the community under consideration.Figure 5
Impact ofγ1 onRc.


#### 4.2.4. Impact of Corruption Transmission Rateβ on Rc

Figure6 has been plotted based on the parameters values that given Table 1 which deduces the impact of corruption transmission rate β on the corruption reproduction number Rc. The displayed simulation states that when the value of β increases, the corruption reproduction number is going up and the value of β<0.17 implies that RC<1. Therefore, the stakeholders are expected to minimizing the values of corruption transmission rate β to control corruption expansion in the community.Figure 6
Impact ofβ on Rc.


#### 4.2.5. Impact of Racism Transmission Rateα on Rr

In Figure7, we have simulated on the relation between racism transmission rate and racism reproduction number Rr. The plot demonstrates that whenever the value of α increases, the racism reproduction number increases, and the value of α<0.322 implies that Rr<1. Therefore, the stakeholders are expected to minimizing the values of racism transmission rate α to control the racism expansion in the community.Figure 7
Impact ofα on Rr.


#### 4.2.6. Impact of Coexistence Recovery Rateγ3 onC1

Figure8 investigates the fact that whenever the racism-corruption recovery rate γ3 increases from 0.54 to 0.73 the number of racism-corruption co-occurrence in the community decreases. The figure deduces that when the values of γ3 increases, the number of individual conduction both racism and corruption among population is going down. This means as the power of treatment increases the racial and corrupted coexistence class in the model become decreases. Therefore, the stakeholders shall expect to maximize the values of parameter γ3 to control the expansion of the racism-corruption coexistence in the community.Figure 8
Impact ofγ3 on C1.


## 4.1. Sensitivity Analysis

The normalized forward sensitivity index of a variable racism and corruption coexistence model (3) basic reproduction number denoted by the symbol Rrc that is differentiable with a parameter ε is defined as SIp=∂Rrc/∂ε∗ε/Rrc as stated in literatures [48–50].The sensitivity indices enable us to examine the relative importance of different parameters in racism and corruption incidence and prevalence. The most sensitive parameter has the magnitude of the sensitivity index larger than that of all other parameters. We can calculate the sensitivity index in terms ofRr and Rc sinceRrc=maxRr,Rc.Using the parameter values given in Table1, the sensitivity indices are listed in Tables 2 and 3 as it is given in Table 2.Table 1
Parameter values for numerical simulation.
ParameterValuesSourceμ0.01[6, 12, 51]Λ50[6, 20, 49]αVariableAssumedσ10.6Assumedσ20.7Assumedη1.3Assumedξ1.2Assumedθ30.2AssumedβVariableAssumedθ10.3[20]θ20.25Assumedγ10.007[39, 50]γ20.006Assumedγ30.008AssumedTable 2
Sensitivity indices ofRr.
Sensitivity indexSensitivity indicesSIΛ+1SIγ2-0.85SIα+1Table 3
Sensitivity indices ofRc.
Sensitivity indexSensitivity indicesSIβ+1SIγ1-0.546SIμ-0.346SIΛ+1Table2 above is the summary sensitivity indices of Rr manipulated with parameters values from Table 1 and provides the result Rr=3.51 at the racism transmission rate α=1.37, which imply that racism spreads throughout the community. Moreover, sensitivity analysis given in Table 2 explains that the human population recruitment rate Λ and racism transmission rate α are highly affecting the racism reproduction numberRr.Here, with the given parameter values in Table3, we have computed Rc=6.03 at the corruption expansion rate β=1.51, which imply that corruption spreads throughout the community. Moreover, sensitivity analysis given in Table 3 explains that the human population recruitment rate Λ and corruption transmission rate β are highly affecting the racism reproduction numberRc.


### 4.1.1. Sensitivity Analysis Graphical Verification

In this subsection, sensitivity analysis of the racism and corruption coexistence transmission dynamics is performed to identify the most influential parameters for the spread as well as for the control of coexistence mind infection transmission in the community. The results of sensitivity analysis based on the set of parameters values given in Table1 are analyzed in Tables 2 and 3, respectively. Graphical simulation of sensitivity indices given by Figure 2 illustrates the sensitivity index of parameters and has verified the qualitative analysis given in Tables 2 and 3, respectively. It shows that the most sensitive parameters which have direct impacts on the basic reproduction numbers are the racism transmission rate, the corruption transmission rate, and the recruitment rate and the most sensitive parameters which have an indirect impact on the basic reproduction numbers are the racism recovery rate and the corruption recovery rate. Stake holders can minimize the transmission rates and maximize recovery rates to prevent and control the coexistence expansion in the community.Figure 2
Parameter sensitivity analysis.


## 4.1.1. Sensitivity Analysis Graphical Verification

In this subsection, sensitivity analysis of the racism and corruption coexistence transmission dynamics is performed to identify the most influential parameters for the spread as well as for the control of coexistence mind infection transmission in the community. The results of sensitivity analysis based on the set of parameters values given in Table1 are analyzed in Tables 2 and 3, respectively. Graphical simulation of sensitivity indices given by Figure 2 illustrates the sensitivity index of parameters and has verified the qualitative analysis given in Tables 2 and 3, respectively. It shows that the most sensitive parameters which have direct impacts on the basic reproduction numbers are the racism transmission rate, the corruption transmission rate, and the recruitment rate and the most sensitive parameters which have an indirect impact on the basic reproduction numbers are the racism recovery rate and the corruption recovery rate. Stake holders can minimize the transmission rates and maximize recovery rates to prevent and control the coexistence expansion in the community.Figure 2
Parameter sensitivity analysis.


## 4.2. Numerical Simulations

In this section, numerical simulation has been performed with MATAB ode45 code to analyze the effect of some parameters that causes for conducting this illegal activity. Most specifically, we investigated the stability of the endemic equilibrium point of the main coexistence model (3), parameter effects on the reproduction numbers, and the impact of being honest on racism and corruption co-acting individuals in the community. Moreover, we have used the parameters values stated in Table 1 for numerical simulation.


### 4.2.1. Behavior of the Coexistence Model Solutions forRrc>1

Figure3 illustrates the result of numerical simulation with ode45 using parameter values that are given in Table 1. From Figure 3, we can observe that after a year the solutions of the racism and corruption coexistence dynamical system (3) will be approaching to the coexistence endemic equilibrium point of the racism-corruption coexistence model depending on the value of Rrc=maxRr,Rc=max3.51,6.03=6.03>1. The numerical simulation justified that the physical phenomenon that can be stated as the expansion and spreading of racism-corruption coexistence activity is consistently occurred throughout the considered community in the study.Figure 3
Behaviors of the coexistence model solutions.


### 4.2.2. Impact of Racism Recovery Rate onRr

Figure4 illustrates the effect of racism recovery rate γ2 on the racism reproduction number Rr. The plot demonstrates that when the value of γ2 increases, the racism reproduction number is going down, and whenever the value of γ2>0.059 implies that Rr<1. In other words, it mean as the power of treatment rate increases the number of racial individual will decrease. Therefore, the stakeholder shall concentrate on maximizing the value of racism recovery rate γ2 by applying possible interventions strategies to prevent and control the problem of racism.Figure 4
Impact ofγ2 on Rr.


### 4.2.3. Effect of Corruption Recovery Rateγ1 on Rc

Figure5 illustrates the relation between corruption recovery rate γ1 and corruption reproduction Rc. The plot shows that that when the value of corruption recovery rate γ1 increases, the corruption reproduction number decreases, and whenever the value of γ1>0.508 implies that Rc<1. This means that as the power treatment rate increases, the number of corrupted individual decreases. Moreover, the result tells us the stakeholder shall concentrate on maximizing the values of γ1 to prevent and control the expansion of racism in the community under consideration.Figure 5
Impact ofγ1 onRc.


### 4.2.4. Impact of Corruption Transmission Rateβ on Rc

Figure6 has been plotted based on the parameters values that given Table 1 which deduces the impact of corruption transmission rate β on the corruption reproduction number Rc. The displayed simulation states that when the value of β increases, the corruption reproduction number is going up and the value of β<0.17 implies that RC<1. Therefore, the stakeholders are expected to minimizing the values of corruption transmission rate β to control corruption expansion in the community.Figure 6
Impact ofβ on Rc.


### 4.2.5. Impact of Racism Transmission Rateα on Rr

In Figure7, we have simulated on the relation between racism transmission rate and racism reproduction number Rr. The plot demonstrates that whenever the value of α increases, the racism reproduction number increases, and the value of α<0.322 implies that Rr<1. Therefore, the stakeholders are expected to minimizing the values of racism transmission rate α to control the racism expansion in the community.Figure 7
Impact ofα on Rr.


### 4.2.6. Impact of Coexistence Recovery Rateγ3 onC1

Figure8 investigates the fact that whenever the racism-corruption recovery rate γ3 increases from 0.54 to 0.73 the number of racism-corruption co-occurrence in the community decreases. The figure deduces that when the values of γ3 increases, the number of individual conduction both racism and corruption among population is going down. This means as the power of treatment increases the racial and corrupted coexistence class in the model become decreases. Therefore, the stakeholders shall expect to maximize the values of parameter γ3 to control the expansion of the racism-corruption coexistence in the community.Figure 8
Impact ofγ3 on C1.


## 4.2.1. Behavior of the Coexistence Model Solutions forRrc>1

Figure3 illustrates the result of numerical simulation with ode45 using parameter values that are given in Table 1. From Figure 3, we can observe that after a year the solutions of the racism and corruption coexistence dynamical system (3) will be approaching to the coexistence endemic equilibrium point of the racism-corruption coexistence model depending on the value of Rrc=maxRr,Rc=max3.51,6.03=6.03>1. The numerical simulation justified that the physical phenomenon that can be stated as the expansion and spreading of racism-corruption coexistence activity is consistently occurred throughout the considered community in the study.Figure 3
Behaviors of the coexistence model solutions.


## 4.2.2. Impact of Racism Recovery Rate onRr

Figure4 illustrates the effect of racism recovery rate γ2 on the racism reproduction number Rr. The plot demonstrates that when the value of γ2 increases, the racism reproduction number is going down, and whenever the value of γ2>0.059 implies that Rr<1. In other words, it mean as the power of treatment rate increases the number of racial individual will decrease. Therefore, the stakeholder shall concentrate on maximizing the value of racism recovery rate γ2 by applying possible interventions strategies to prevent and control the problem of racism.Figure 4
Impact ofγ2 on Rr.


## 4.2.3. Effect of Corruption Recovery Rateγ1 on Rc

Figure5 illustrates the relation between corruption recovery rate γ1 and corruption reproduction Rc. The plot shows that that when the value of corruption recovery rate γ1 increases, the corruption reproduction number decreases, and whenever the value of γ1>0.508 implies that Rc<1. This means that as the power treatment rate increases, the number of corrupted individual decreases. Moreover, the result tells us the stakeholder shall concentrate on maximizing the values of γ1 to prevent and control the expansion of racism in the community under consideration.Figure 5
Impact ofγ1 onRc.


## 4.2.4. Impact of Corruption Transmission Rateβ on Rc

Figure6 has been plotted based on the parameters values that given Table 1 which deduces the impact of corruption transmission rate β on the corruption reproduction number Rc. The displayed simulation states that when the value of β increases, the corruption reproduction number is going up and the value of β<0.17 implies that RC<1. Therefore, the stakeholders are expected to minimizing the values of corruption transmission rate β to control corruption expansion in the community.Figure 6
Impact ofβ on Rc.


## 4.2.5. Impact of Racism Transmission Rateα on Rr

In Figure7, we have simulated on the relation between racism transmission rate and racism reproduction number Rr. The plot demonstrates that whenever the value of α increases, the racism reproduction number increases, and the value of α<0.322 implies that Rr<1. Therefore, the stakeholders are expected to minimizing the values of racism transmission rate α to control the racism expansion in the community.Figure 7
Impact ofα on Rr.


## 4.2.6. Impact of Coexistence Recovery Rateγ3 onC1

Figure8 investigates the fact that whenever the racism-corruption recovery rate γ3 increases from 0.54 to 0.73 the number of racism-corruption co-occurrence in the community decreases. The figure deduces that when the values of γ3 increases, the number of individual conduction both racism and corruption among population is going down. This means as the power of treatment increases the racial and corrupted coexistence class in the model become decreases. Therefore, the stakeholders shall expect to maximize the values of parameter γ3 to control the expansion of the racism-corruption coexistence in the community.Figure 8
Impact ofγ3 on C1.


## 5. Discussion

In Section1, we have introduced the backgrounds of the racism, corruption, and racism-corruption coexistence and reviewed some literatures related to the study. In Section 2, we have classified the total human population into seven nonmutual distinct classes and formulated the nonlinear deterministic racism and corruption coexistence dynamical system using a system of ordinary differential equations.In Section3, we have analyzed the qualitative behaviors of the newly developed model such as the positivity of solutions of the model, boundedness of the dynamical system, racism-free equilibrium point, and corruption-free equilibrium point. Additionally, we have analyzed the stability of endemic equilibriums, stability analysis of disease-free equilibrium point, and sensitivity analysis of reproduction numbers. We have also deduced the effect of parameters in the expansion or control of racism and corruption as well as parameter effect on the infected population.


## 6. Conclusion

Nowadays, racism and corruption coexistence is a major problem affecting nations throughout the world, but literatures on prevention and controlling its expansion through the community were rare. In this work, we have developed a first new nonlinear compartmental deterministic mathematical model on the transmission dynamics of racism and corruption coexistence expansion. The developed model has disease-free equilibrium points that are both locally asymptotically and globally asymptotically stable whenever their corresponding basic reproduction number is less than one. All the model mind infection endemic equilibrium points were both locally asymptotically and globally-asymptotically stable whenever their corresponding reproduction number is less than unity. The model did not have the phenomenon of backward bifurcation. The sensitivity analysis of the model showed us the racism transmission and corruption transmission rates are the most sensitive parameters which have a direct effect on the racism and corruption coexistence mind infection transmission in the community. Also the racism recovery rate and corruption recovery rate have high indirect impact on the basic reproduction numbers of the racism model and corruption model, respectively. Using the parameter values given in Table1, we have obtained Rr=3.51 at α=1.37 and Rc=6.03 at β=1.51, i.e., Rcr=maxRr,Rc=max3.51,6.03=6.03. Using numerical simulation, we have verified the qualitative result that the endemic equilibrium point of the racism and corruption coexistence model is locally asymptotically stable when Rcr=maxRr,Rc=max3.51,6.03=6.03>1. Also numerical simulation results showed that whenever the racism transmission rate increases the racism mind infection transmission increases, the corruption transmission rate increases, the corruption mind infection transmission increases, the racism recovery rate increases, the racism mind infection transmission decreases, the corruption recovery rate increases, and the corruption mind infection transmission decreases.Moreover, based on the impact of some changes of parameters on the corresponding reproduction numberRr and Rc, we shall give future directions for the stakeholders in the community. The results we have obtained have a crucial role for stakeholders, as it governs the eradication and/or persistence of racism, corruption, and racism-corruption coexistence which are illegal activities in a community. Stakeholders shall concentrate on decreasing the racism transmission rate, the corruption transmission rate, and increasing or maximizing the values of racism and corruption recovery rates that are used to minimize and possibly to eradicate the problem from the community.Finally, we recommend the governments of nations to introduce, apply and ensure anticorruption and antidiscriminatory laws, and take the bold measures to beak the interconnection of corruption and racism. We want to remark the whole community stay unite to identify common problems and committed to research and advocacy from societies. The international institutions shall be collaborated for better understanding of these two interlinked problems and set up monitoring and investigation bodies. In the limitations of this study, the next potential researchers can incorporate them and extend this study: optimal control approach, stochastic approach, fractional order derivative approach, environmental impacts, age, and spatial structure, whenever possible validating the model by applying appropriate real data.



---

*Source: 9977727-2022-08-11.xml*