Epidemic Model of Leptospirosis Containing Fractional Order

and Applied Analysis 3 The Riemann Liouville derivative has some disadvantages when applying to real worldmodel with fractional differential equations. Therefore, here we use a modified fractional differential operatorD a which is used in his work by Caputo on the theory of viscoelasticity. Definition 3. The Caputo fractional derivative of g(x) of order α > 0 with a ≥ 0 is defined as (D α a g) (x) = (J m−α a g (m) ) (x)


Introduction
Leptospirosis disease is an important infectious disease.This kind of infection occurs in urban areas of industrialized and developed countries and also in the rural areas.The people of the city who walk in dirty water are mostly infected.Workers planting rice, sewer cleaners, cleaning canals workers, and agriculture labor get the disease easily.The disease flourishes due to delay in diagnosis and unavailability of clinical infrastructure.The cause of the disease is bacteria.It is potentially fatal infection of brain, kidney, liver, heart, and lung.The people who can get infection are those who have contact with infected animals, soil, or water in which the bacteria is present.The outdoor people, who work with animals, face the risk of leptospirosis infection, similarly workers in farms, sewer, mine, slaughter houses, dairy farmers, and animal caretakers and those who work with fishes and military personnel.Those people who work outdoors like swimming, rafting, and kayaking also face the risk of infection [1].
The mathematical modeling dealing with the complex biological structures is a great challenge for the researchers.For simple biological models, the study of integer order differential equations is enough for the description of their dynamics while for those systems contains complexities are characterized due to its variability structured properties, like nonlinearity of the systems, multiscale behavior, and the mathematical relation between the parameters [2].To study such nonlinear biological models, with complex mathematical structures, the fractional derivatives provide the tool to handle the dynamical behavior of such complex systems.The fundamental property related to these models is the nonexistence of the differential operator of integer order.These characteristics show information about the present but also show its historical states.
The fractional derivative and fractional integrals were first used by Magin [3] in his model which is related to the relation between stress-strain and biomaterials.A new fractional analytical approach via a modified Riemann Liouville derivative has been used by [4].For solution of the 2D region with obstacles, the author [5] used the fractional differential equations.For the giving up smoking model, the fractional differential equation was used by [6].The author used the differential transforms method and Pade approximants for a fractional population growth model [7].For more information about fractional calculus, see [2,[8][9][10].The author had made many efforts for modeling of leptospirosis epidemic disease since 2011; see [11][12][13][14][15][16].The author proposed 2 Abstract and Applied Analysis a mathematical model that describes the epidemic leptospirosis disease [11].The model proposed by [11] is given by and the related initial conditions: Here  ℎ (),  ℎ (),  ℎ (),  V (), and  V (), respectively, represent the population of susceptible human, infected human, recovered human, susceptible vector, and infected vector at time .The rate at which the population of human increases is shown by  1 .The natural mortality rate for the human population is  ℎ ;  1 ,  2 , and  3 represent the transmission coefficients.The parameter  ℎ shows the individuals who become susceptible again.The death from the disease that occurs to humans is shown by  ℎ .The rate of recovery from infection for the human is denoted by  ℎ .The growth rate of the vector population is represented by  2 ;  V is the natural death rate for vector and disease related death rate for the vector is  V .The model presented by [11] is given by system (1) which represents a system of nonlinear ODE; their fractional order differential equation is given as follows: where    represents the fractional derivative in Caputo derivative sense and the order of the fractional derivative is shown by the parameter  and 0 <  < 1, with the related initial conditions (2).The parameter which describes the general response expression in representing the order of fractional derivative gives different results for different values.Obviously, the integer-order system can be viewed as a special case of the fractional-order system by putting the time-fractional order of the derivative equal to one.To put it simple, for the higher in order, the behavior of the fractional order system is the same in the case of integer order.According to research, this is the first work, which is available on the epidemic model of leptospirosis in fractional order.Due to this reason, we express the approximate solutions to the problems (3) by the technique of (MSGDTM).The solution to such nonlinear fractional equation can also be available by using an efficient approach [17].
The MSGDTM is a modified form of the GDTM [6,8,[18][19][20], in which it is treated as an algorithm in a sequence of small steps, to obtain the accurate approximate solution to the desire models.By using GDTM, the obtained solution is valid for short interval of time while the solution obtained from multistep generalized differential transform method (MSGDTM) is valid for a long time.To obtain the solution, using MSGDTM is more accurate and valid for long interval of time and agrees well with the classical Runge-Kutta numerical solution method, with the unity order derivative [6,21].
The paper has been organized/structured as follows.Section 1 is introduction; Section 2 contains basic definitions of fractional calculus and their notations; Section 3 gives nonnegative solution of the model; Section 4 contains application of MSGDTM; Section 5 gives numerical results; and Section 6 is the conclusion.

Basic Definitions and Notations
This section contains some fundamental definitions and notations of fractional calculus that will help us in the proceeding sections [3,[22][23][24].Definition 1.A function () with the positive values of  is known in the space   ( ∈ ) if it is written in the form () =    1 () and for some  > , where  1 () is continuous in [0, ∞), and it is known to be in the space    if  () ∈   ,  ∈ N. Definition 2. The Riemann Liouville integral operator of order  > 0 with  ≥ 0 is defined as For the properties of the operator [22], we need only the following.For  ∈   ,  > 0,  > 0,  ∈ , and  > −1, one has where   (,  + 1) represents the incomplete beta function defined as Abstract and Applied Analysis 3 The Riemann Liouville derivative has some disadvantages when applying to real world model with fractional differential equations.Therefore, here we use a modified fractional differential operator    which is used in his work by Caputo on the theory of viscoelasticity.Definition 3. The Caputo fractional derivative of () of order  > 0 with  ≥ 0 is defined as for  − 1 <  ≤ ,  ∈ ,  ≥ , () ∈   −1 .Many authors investigated the Caputo fractional order derivative; for  − 1 <  ≤ , () ∈    , and  ≥ −1, one gets

Numerical Results
In this section, we present the numerical solution of the proposed model (3) numerically by the classical Runge-Kutta order four method for the integer case derivative and analytically we use the MSGDTM.We apply the algorithm on the interval [0, 30] for the approximate solution of nonlinear fractional differential equations (3).For  = 10 and  = 3000, the results are obtained.The results are obtained by using computer algebra package Mathematica.The initial conditions  1 = 40,  2 = 20,  3 = 10,  4 = 50, and  5 = 30 are used.The values for the parameters we assumed are presented in Table 2.In Figures 1, 2, 3, 4, and 5, the approximate solution is obtained by MSGDTM and classical Runge-Kutta method order four scheme for  = 1.The approximate solution of MSGDTM has good agreement with  the Runge-Kutta order four method.Figures 6, 7, 8, 9, and 10 represent the approximate solution obtained by MSGDTM and classical Runge-Kutta order four scheme considering different values of .We conclude from our graphical results that the MSGDTM results behave as closely as Runge-Kutta method.

Conclusion
In this paper, we have presented an approximate solution to an epidemic model of leptospirosis in fractional order.First, we obtained the analytical solution to the fractional order of epidemic model.Then, we solved the fractional order        system and the results were compared with classical Runge-Kutta method.The results obtained from MGDTM have a good agreement with Runge-Kutta order four method.The MGDTM results are valid for a larger .

Figure 1 :
Figure 1: The dotted line and the solid line connectively show the Runge-Kutta method and MSGDTM.

Table 2 :
Parameter values used in the numerical simulations of the optimal control.