Modeling of infectious diseases with stochastic differential equation (SDE) has increased foundation lately due to its extensive variety of applications and its aptitude to reflect actuality in epidemiology  . The diseases outbreaks in a population of susceptibles rationally go behind stochastic processes  . Stochastic process occurs naturally in lots of physical applications where randomness is to be incorporated in the mathematical model   . In recent years, main studies on stochastic model that have been published by researchers have recognized the growing significance of study the stability of stochastic positive equilibrium, as well as the global stability of the endemic equilibrium    . In this paper we approached by using deterministic and stochastic model. Briefly deterministic models are model processes which are often described by differential equations, with a unique input leading to a unique output for well-defined linear models and with multiple outputs possible for non-linear models. Throughout this paper, let be a complete probability space with a filtration satisfying the usual conditions (i.e., it is right continuous and increasing while contains all -null sets).
Considering the general n dimensional stochastic differential equation
on with initial value , the solution is denoted by . Assume that and for all , so (1.0) has the solution . This solution is called the trivial solution.
Definition 1.1. The trivial solution of (1.0) is said to be as follows:
1) stable in probability if for all ,
2) asymptotically stable if it is stable in probability and, moreover,
3) asymptotically stable in the large if it is stable in probability and, moreover, for all
In this paper, we consider the epidemic model in Meta-population setting. From the proposed schematics of the compartment model shown in  , we will extract a metapopulation model for HIV dynamics among the youth coupled with awareness/education i.e., we extended the single patch disease model to include multiple patches (see Figure 1).
Figure 1. A schematic of the metapopulation model. A schematic of the metapopulation model for HIV transmission in the youths coupled with awareness/education in each patch i, i = 1, ・・・, n.
Differential Equation of the Model
States variables: , susceptible children, , susceptible uneducated female youth, , susceptible educated female youth, , susceptible uneducated male youth, , susceptible educated male youth, , infected children(children who were infected either during pregnancy or at childbirth), , infected uneducated female youth, , infected educated female youth, , infected uneducated male youth, , infected educated male youth, and , antiretroviral therapy treatment among the youths. Parameters: , death rate of the youth and children, , disease induced rate in the youth and children before ART, , disease induced deaths in the youth and children after ART, , birth rate, , rate of vertical transmission, , probability a susceptible female youth gets infected by infected male youth, , probability a susceptible male youth gets infected by infected female youth, , rate at which children grow to become uneducated female youth, , rate at which children grow to become educated female youths, , rate at which children grow to become uneducated male youths, , rate at which children grow to become educated male youth, , awareness/education rate, , rate at which infected youth take the ART. The solid lines represent movement between classes and the dashed lines represent rate at which a susceptible individual moves in to infected class.
The model answers one important underlying research subjects; the determination of the existence of the threshold parameter which hints on the spreading or dying out of an invading epidemic into a population of susceptible. In this research article, we first study the positivity and boundedness of the system (1.0). The basic reproduction ratio is determined. Applying the hypothetical theorem of the Lyapunov functional, we find out the global stability of the two equilibria for system (1.0). We extend our stability analysis to the stochastic system (5.0), which is obtained by random perturbation of the deterministic system (1.0) and find the stability of its positive equilibrium. Finally, numerical examples which shows the dynamics of systems (1.0) and (5.0) are given, which gives the explicit difference in the dynamics of the models.
2. Basic Properties of the Model
In this section, the basic properties of model system (1.1) which are useful in the proofs of stability are studied. These are the invariant region and positivity of solutions. The former describes the region in which the solutions of system (1.1) makes biological sense while the latter describes non-negativity of solutions of system (1.1). The model under consideration monitors a human population and as such, we need to have that all the parameters and the variables of the model are positive for all .
2.1. Positivity of Solutions
The theory of ordinary differential equations requires that, for every set of initial conditions
the state variables of the solution must remain non-negative.
Proposition 2.1 Let the initial data be
Then, the solution set
of system (1.1) is positive for all
Proof. Let From the first equation of model
Integrating (2.0) by separation of variables gives
This proves that for all . Similarly, it can be shown that the remaining variables of system (1.1) are also positive .
Remark 2.1. for all .
2.2. Invariant Region
We now apply Birkhoff and Rotas theorem on differential inequality (2.1). By separation of variables of differential inequality (2.1), we get
Integrating (2.2) on both sides gives,
where A is a constant. Now, applying the initial condition in (2.3), we get
Substituting (2.4) into (2.3) gives
Making N the subject in (2.5) we have,
As in (2.6) above, the population size N, approaches
Therefore, the feasible solutions set of system (2.7) enters the region
In this case, whenever , then which means that . On the other hand, whenever , every solution with initial
condition in remains in that region for . Thus, the region is positively-invariant.
3. Basic Reproduction Number
The basic reproduction ratio ( ) is defined as an infections originating from an infected individual that invades a population originally of susceptible individuals.
The above system can be represented in matrix form as where f is the matrix of the infection rates and v is the matrix of the transition rates.
The spectral radius of the Metzler Matrix, , is defined as the largest eigenvalue of the Metzler Matrix. Thus:
If for j = 1, 2, 3, 4, then each infectious individual in Sub-Population j infects on average less than one other person and the disease is likely to die out Otherwise, If for j = 1, 2, 3, 4, then each infectious individual in Sub-Population j infects on average more than one other person; the infection could therefore establish itself in the population and become endemic. An SIR epidemic model, where the presence or absence of an epidemic wave is characterized by the value of.
4. The Global Stability of the Endemic Equilibrium
In this part, we analyse the global stability of the endemic equilibrium point by construction a appropriate Lyapunov function. For simplicity, we consider the reduced model system (6) to prove for global stability. We use the come up to of  as it is used for several complicated epidemiological models. We consider the Lyapunov function of the form
where (for i = 1, 2, ・・・, 6 is a properly chosen positive constant in the given region. is a population of compartment i and is the equilibrium level. So we define the Lyapunov function as
The time derivative of L is
At an endemic equilibrium point we have
F is non-positive by following the approach of  . Thus, for
. Hence and is zero when,
, , , ,. Therefore, the largest
invariant set in such that is the singleton which is our
endemic equilibrium point. By LaSalles invariant principle  we conclude that is globally asymptotically stable (g.a.s). Thus, we establish the following theory.
Theorem 3.1 When, the endemic equilibrium point is globally asymptotically stable in.
5. The Stochastic Model
Stochastic perturbations were bring in some of the major parameters involved in the model equations.
Here, we bring in stochastic perturbations in the major parameters of the deterministic model (1.1). Thus we permit stochastic perturbations of the variable
around their values at positive equilibrium. Hence, we assume that the white noise of the stochastic perturbations of the variable around values of are proportional to the distances of
Hence the stochastic version of model (1.1) is
With where are real constants, and
are independent wiener processes. We examine the asymptotic stability behavior of the equilibrium of the stochastic equation (5.0) and contrast results with the deterministic model (1.1).
Stochastic Stability of the Positive Equilibrium
It can be shown clearly that, the deterministic model (1.1) has one disease-free equilibrium
which is globally asymptotically stable when. However, when, the disease-free equilibrium is unstable. Obviously, there is also a unique positive endemic equilibrium
This equilibrium is globally asymptotically stable. The stochastic system (5.0) has the similar equilibria as the deterministic system (1.1). Assuming that, we examine the stability of the endemic equilibrium of (5.0). The stochastic differential equation (5.0) can be centered at its positive equilibrium by the change of variables
The linearized system of the stochastic model (5.0) around takes the form
Clearly, the endemic equilibrium corresponds to the trivial solution in (5.2). We denote L to be the differential operator associated with (5.2), defined for the family of nonnegative functions such that it is continuously differentiable with respect to t and twice with respect to x.
According to Afanas ev et al.  , the differential operator L for a function is given by
where, “T” and “Tr” are the transposition and trace respectively. With reference to Afanas?? ev et al.  , the following results hold.
Theorem 5.1: Suppose a function exist, satisfying the following inequalities
where and. Then the trivial solution of (5.2) is pth moment exponentially stable. Again, given that p = 2 the trivial solution is supposed to be exponentially stable in mean square and the equilibrium x = 0 is globally asymptotically stable. From theorem 5.1, the conditions for stochastic asymptotic stability of trivial solution of (5.0) are given theorem 5.2.
Theorem 5.2: Suppose
and hold, then the zero solution of (5.0) is asymptotically mean square stable.
Proof: We consider the Lyapunov
where non-negative constants that will be chosen in the course of the proof. It can be easily ascertained that inequality (5.4) hold true when Applying the operator L on gives
Now remark that
Now, from Equation (5.5), if we choose
and then from Equation (5.5), it is easy to verify that,
From the assumptions of the theorem, we deduce that and. Hence, D is a symmetric positive definite matrix. Let denote the minimum of its eleven positive eigenvalues, , , , , , , , , and; then,we can easily get
According to Theorem 5.1, we conclude that the trivial solution of stochastic system is globally asymptotically stable.
Hence, according to theorem 5.1, the proof is completed.
In this article, the dynamics of deterministic epidemic model and its stochastic variant are presented. The stability analyses of the deterministic model were investigated. Suitable Lyapunov functions were constructed for the global stability of the two equilibria. Mathematical analysis was done and it was established that in the absence of the disease a disease free equilibrium will always exist if for j = 1, 2, 3, 4,. We also established that the endemic equilibrium exists in the presence of the disease that is when for j = 1, 2, 3, 4, with the infectious population greater than zero. Reducing the infection in the vector population reduces for j = 1, 2, 3, 4, greatly. Thus the best methods of controlling HIV transmission is to target the Infected uneducated female youth, Infected educated female youth, Infected uneducated male youth, Infected educated male youth. is a threshold that completely determines the global dynamics of disease transmission. Our major purpose of the study was to examine the asymptotic stability behavior of the endemic equilibrium of the stochastic version of the deterministic epidemic model in Metapopulation setting.
Conflict of Interest
The author(s) declare(s) that there is no conflict of interest regarding the publication of this paper.
 Liu, X.Q., Zhong, S.M., Tian, B.D. and Zheng, F.X. (2013) Asymptotic Properties of a Stochastic Predator-Prey Model with Crowley-Martin Functional Response. Journal of Applied Mathematics and Computing, 43, 479-490.