Global Stability of Seiqrs Computer Virus Propagation Model with Non-linear Incidence Function

In this paper, we present an SEIQRS epidemic model with non-linear incidence function. The proposed model exhibits two equilibrium points, the virus free equilibrium and viral equilibrium. The model stability is connected with the basic reproduction number R0. If R0 < 1 then the virus free equilibrium point is stable locally and globally. In the opposite case R0 > 1, then the model is locally and globally stable at viral equilibrium point. Numerical methods are used for supporting the analytical work.


Introduction
The Malicious objects are harmful codes that reproduce and spread by way of internet [1].A number of malicious objects have been noticed to appear on the internet for the last decade.They are powerful enough to confront and disfunction the infected computer system.As the internet has been used for a wide range of function so, the malicious objects have become a serious threat to man's work.They have hampered the economic and financial growth of man.Now it has become an issue of great concern and therefore it is vital to suggest such plans that are effective in combating the malicious objects.The correspondence between the propagation of biological virus and the malicious objects compel the researchers to advocate epidemic models that could explain and overcome the propagation of malicious objects.Re F. Wang et al. in [2] present the SEIQRS model with graded infection rates for internet worms.The infection ratio of exposed class have less then infected ratio and Q. Badshah show the worm free and viral equilibrium stability locally as well as globally, which is connected to threshold quantity.In the previous decade or something like that, various epidemic models (see [3]- [6]) were essentially obtained to depict the spread of computer virus.B. K. Mishrand G. M. Ansari present an E-SIRS epidemic model for virus and worms in a computer network.In which they consider the immune as well as the latent period and self-replication time.They stated the global stability of the proposed model in term of reproduction number [7].M. Kumar et.al.studied the quarantine and vaccination role in the control of virus and worms in the computer network [8].S. T. Ge et al. presented the SEIQR epidemic model with the discrete time in computer networks.The stability of the proposed model was shown by Lyapunove method in [9].M. Kummar and A. Kumar analyzes the SI j RS E-Epidemic model for various groups of infection in computer networking [10].
In this paper, we present a Propagation Model with non-linear incidence function, (susceptible, exposed, infected, quarantined, and recovered) which exhibits two equilibria the virus free and viral equilibrium point.The stability of both equilibrium points is connected with the threshold quantity.On the off chance that its value is less then unity then the virus free equilibrium point is locally and additionally globally asymptotically stable and in the opposite case of the threshold quantity, the same condition of stability is satisfied for viral equilibrium point.
The rest of the paper is set as follows, in Section 2 we shall formulate the new model, in Section 3 the basic propertied is discussed like reproduction number, virus free and viral equilibrium.In Section 4, we shall study the local and global stability of the virus free equilibrium point.In Section 5, we shall examine the local and global stability at viral equilibrium point.Finally, we support our analytical work with numerical simulations.

Model Formulation
As usual the computer it either internal or external (i.e.connected with the internet or not).In this model we divide the internal computers into five states.Susceptible, exposed, infected, quarantined and recovered.Where the variable S(t), E(t), I(t), Q(t) and R(t) denote the sizes of nodes at time t in the states susceptible, exposed, infectious, quarantined and recovered respectively.N(t) the total number of computer at time t, where where ∂ is the recruitment rate of the computers.The ratio of crashing nodes without attack of malicious objectsis ϑ .The transmission rate is denoted by β .While  is the ratio at which the recovered computers loss their immunity.The contact ratio between exposed and infected computers is 1 η .The rate of recovery is 2 η and 3 η is the contact ratio between infected and recovered computers.ξ is the contact ratio between infected and recovered computers.ϕ is the fraction of computers to be immune from virus.
We suppose that the like transmission rate in the form of ( ) , where f is a positive function with, ( ) ψ ′ ≥ , used by [11].This is the generalized form of mass action law, that is, ( )

Basic Reproduction Number and Equilibrium Points
In this subsection, we shall investigate the feasible region and the threshold quantity for the proposed model (1).Where the size of total population is represented by N satisfy the following equation.

(
) and then ( ) Therefore, for the system (1) the positively invariant feasible region is given below.
For viral free equilibrium point we take 0 E I Q = = = , thus from the system (1) we get ( To investigate the reproduction Number for system (1) we use the next generation method [12].Assume that ( ) , z E I = then from system (1) we can write.
The spectral radius of the matrix Now we examine the virus equilibrium point for the given system (1) and denoted by ( ) )( ) Next to take the derivative of ( ) Therefore, F has a unique positive zero if and only if ( ) It can be stated as below.Proposition 3.1.Suppose that ( ) F I holds the conditions imposed on it.Then, the system (1) has virus free equilibrium point which is true for all parameter values.When, the system (1) admits also a unique viral equilibrium point ( )

The Stability Analysis of the Virus Free Equilibrium Point
In this section, we shall study the local and global stability of the given system (1) at virus free equilibrium point 0 E .Theorem 4.1.The given system (1) is locally asymptotically stable if 0 1 R < , otherwise unstable.Proof: To examine the local stability of the given system (1) we construct the jacobian matrix at 0 E , which is given below.
( ) With row operation we can get the characteristic equation of the above jacobian matrix is The first three roots of the above equation is and the other roots will be find on the following way.

(
)( ) After a little algebraic calculation we can get the following equation T > if 0 1 R < then according to Routh-Hurtwiz criteria the virus free equilibrium point 0 E is locally asymptotically stable.If 0 1 R > then 2 0 T < and hence one of the eigenvalue have positive real part.So the virus free equilibrium point is unstable.

Global Stability of Viral Free Equilibrium Point
Theorem 4.2.The system (1) at 0 E is globally stable when 0 1 R > otherwise unstable.Proof: To examine the global stability of the proposed model consider a Lyapunov function as below: ( ) ( ) Taking the derivative i.e. d d

V t
, we get the following and after a little simplification we get the following results

( ) ( )( )( )
then according to Lasalle's invariance principle the virus free equilibrium point 0 E is globally stable and in case of 0 1 R > the mention result is failed and hence the system is unstable [13].

The Local Stability Analysis of Viral Equilibrium Point
In this section we will discuss the local as well as the global stability of the system (1) for the viral equilibrium Q. Badshah point.
Theorem 5.1.The given system (1) at 1 E is locally asymptotically stable if 0 1 R > .Proof: To investigate the local stability of the given system (1) we find the jacobian matrix at 1 E , which is given below.
And after some row operation we get the following matrix [ ] Thus, the system (1) at 1 E has eigenvalues, that contains negative real part.So, we conclude that the model (1) is locally asymptotically stable.

Global Stability of Viral Equilibrium Point
In this segment we should look at the global stability of the given system (1) at viral equilibrium point.We use the method presented in [14] established by Li and Mouldowney.According to which we find sufficient condition for the global stability of the model (1) at viral equilibrium point.First we shall discuss the method briefly.Consider ( ), where : , We know that equilibrium point * y is globally stable in K when it locally stable as well as all the orbits in K converges to * y .Bendixson criterion for 2 c ≥ we mean that g satisfied a condition which prohibit the existence of non-constant periodic solutions of (4) The classical Bendixson's condition ( ) div 0 g y < for 2 c = , is robust under 1  C robust properties are discussed in [14] [15].If a neighborhood V exist for 0 y K ∈ for 0 I > such that i.e.
( ) Then the point 0 y is wandering for the Equation ( 4).The following principle for global stability is established in [14] for any finite dimension autonomous systems.
Lemma 5.1.Assume that conditions ( ) i L and ( ) ii L holds.Such that (4) satisfies a Bendixson criterion that is robust under 1  C local perturbations of g at all non-equilibrium non-wandering points for (4).Then * y is globally stable in K provided it is stable.
The accompanying Bandixson criterion is displayed in [14] and indicates to have the heartiness needed by be a nonsingular matrix and ( ) Suppose that for a compact absorbing set 1 , O P − not only exist but also continues for h O ∈ .A q quantity is defined as.
( ) ( ) ( ) And [ ] ** J represents the second additive compound matrix of the Jacobian matrix J, i.e., ( ) ( ) and ( ) B ℵ be the Lozinskii measure of B with respect to a vector norm . in N  , ( ) It is shown in [14] that, if K is simply connected, the condition 0 q < rules out the presence of any trajec- tory that fives rise to a simple closed rectifiable curve, like periodic orbits, and heteroclinic cycles it is invariant for system (4).The accompanying result for the global stability is displayed in [14] by Li and Muldowney.Lemma 5.2.Let the simply connected set K satisfy the conditions ( ) i L and ( ) ii L .Then the system (4) is globally stable in K at a unique equilibrium * h if 0 q < .Now for the analysis of global stability at viral equilibrium we follow the method presented by Li and Muldowney in [14].The unstable viral equilibrium point 1  E means that the viruses will persists [16], i.e., for any solution , , Θ Θ Θ in the trajectory of the model (1) there exist a constant 0 s > which satisfy the following Proof: For the global stability of the system (1) we find the second additive compound matrix [ ] . 2 Suppose the norm in 3  as: where ( ) , , v v v be the vector in 3   and ℵ signify the Lozinskii measure as regard to this norm follows [17].
where, ( ) ( ) using second equation of system (1) ( ) where, ( ) Thus the viral equilibrium 1 E is globally asymptotically stable according to [14].Next we talk about the given subsystem of system (1) Based on (6), we get * 0 e 0 e d [18] the system (1) at 1  E is globally asymptotically stable.

Discussion
The aim of this work is to study and analyze the dynamic behavior of an epidemic model SEIQRS with a nonlinear incidence function.We consider a mathematical model of the type SEIQRS and obtained the basic reproduction number, to determine its dynamical behavior.In epidemiology, the reproduction ratio is very important, because the stability of the proposed model is associated with reproduction ratio.For virus free equilibrium point, the model is stable locally and globally if 0 1 R < .In the case of when 0 R exceeds then unity then the system (1) is asymptotically stable at viral equilibrium point.
The purpose of this section is to support the analytic results mentioned in above work are supported through numerical results.Numerical results we choose different values of parameters, which we discuss below with the help of graph.For 0 ϕ = , the reproduction ratio for immune free system is ( )( ) R ≤ , then clearly the virus is vanished.But if 1 1 R > (see Figure 1), then the vaccination ( ) v ϕ is needed so that 0 1 R ≤ or equivalently to Nothing that when 0 1 R > , since v ϕ is an increasing function of  , and the value of the vaccination cover    Q. Badshah simply connected, open set and ( )1 g C K ∈.The solution of Equation (4) is denoted by ( ) y t y i e g y = Let us suppose that the assumptions given below are true:

Figure 3 . 1 v
Figure 3.In the opposite case ( ) 1 v ϕ > , the virus will persists, even if we can vaccinate all recruitment Figure 4.Nothing that when 0 1 R > , since v ϕ is an increasing function of  , and the value of the vaccination cover