BioMed Research International
Volume 2015 (2015), Article ID 584234, 10 pages
http://dx.doi.org/10.1155/2015/584234

Research Article

Modeling the Effects of Multiple Intervention Strategies on Controlling Foot-and-Mouth Disease

Steady Mushayabasa and Gift Tapedzesa

Department of Mathematics, University of Zimbabwe, P.O. Box MP 167, Harare, Zimbabwe

Received 26 March 2015; Revised 4 September 2015; Accepted 13 September 2015

Academic Editor: Klaus Wimmers

Copyright © 2015 Steady Mushayabasa and Gift Tapedzesa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Foot-and-mouth disease (FMD) is a threat to economic security and infrastructure as well as animal health, in both developed and developing countries. We propose and analyze an optimal control problem where the control system is a mathematical model for FMD that incorporates vaccination and culling of infectious animals. The control functions represent the fraction of animals that are vaccinated during an outbreak, infectious symptomatic animals that are detected and culled, and infectious nonsymptomatic animals that are detected and culled. Our aim was to study how these control measures should be implemented for a certain time period, in order to reduce or eliminate FMD in the community, while minimizing the interventions implementation costs. A cost-effectiveness analysis is carried out, to compare the application of each one of the control measures, separately or in combination.

1. Introduction

Foot-and-mouth disease (FMD) is a highly contagious viral disease of cloven-hoofed animals and is one of the most economically important diseases of livestock [1]. The causative agent of the disease is a small icosahedral nonenveloped RNA virus classified within the Aphthovirus genus, as a member of the Picornaviridae [2]. The virus is airborne and can also be transmitted through physical contact with infected animals’ expired air, saliva, milk, urine, semen, animal feed and bedding, and so forth [3]. Direct or indirect contact with FMD-infected animals can result in susceptible animals becoming diseased or subclinically infected [2]. The incubation period for FMD can vary with the species of animal, the dose of virus, the viral strain, and the route of inoculation [4]. It is reported to be 1 to 12 days in sheep, with most infections appearing in 2–8 days, and 2 to 14 days in cattle [5]. After the incubation, a fraction of infected animals progress to symptomatic stages and the remainder become FMD carriers. FMD carriers are defined as animals in which either viral nucleic acids or live virus can be found for more than 28 days after infection [4]. How long an animal can remain a carrier varies with the species [6]. Most cattle carry foot-and-mouth disease virus (FMDV) for six months or less, but some animals can remain persistently infected for up to 3.5 years [4, 5]. The virus or its nucleic acids have been found for up to 12 months in sheep, up to 4 months in goats, for a year in water buffalo, and up to 8 months in yaks [7]. Carrier animals can only be identified by collecting esophageal-pharyngeal fluids for virus isolation and/or the detection of nucleic acids. Repeated sampling may be necessary to identify a carrier, as the amount of virus is often low and fluctuates [4].

Conventional control measures of the disease include movement restriction—for example, through construction of “veterinary boundaries”—that is, cordon fences erected to divide a country into multiple subregions to restrict the movement of animals across the borders; educational awareness; quarantine; vaccination and culling of detected infected animals [3]. Practising import regulations can be essential to prevent FMDV from being introduced from endemic regions in infected animals or contaminated foodstuffs fed to animals [4].

Since the 2001 FMD outbreak in the United Kingdom mathematical modeling of FMD has been an interesting topic for a number of researchers; see, for example, [3, 8–11]. Tildesley et al. [8] proposed a probabilistic FMD transmission model to explore an optimal deployment strategy of limited reactive ring vaccination of cattle in a single epidemic outbreak. Their work suggested that optimal ring size is highly dependent upon logistic constraints. In 2011, Hansen and Day [10] constructed an SIR, Susceptible-Infectious-Removed, model to assess the impact of limited isolation resources and limited vaccination resources. Results from their study highlighted a number of areas that warrant further study and also emphasized the impact of time-optimal control on controlling the generation of new infections in resource-limited settings. More recently, Ringa and Bauch [3] developed an SEIRVC, Susceptible-Exposed-Infectious-Removed-Vaccinated-Culled, pair-approximation model of FMD transmission in a near-endemic population. Their work suggested that the optimal long-term control of FMD by vaccination in near-endemic settings can be achieved by rolling out a prophylactic vaccine as much as possible, especially if resources are limited.

The primary goal of this paper is to formulate a model for FMD dynamics that includes relevant biological detail and accounts for multiple intervention strategies. The FMD intervention strategies to be considered are pre- and postexposure vaccines and culling of symptomatic and nonsymptomatic infectious animals. During an FMD outbreak, a vaccine may be protective or suppressive. Suppressive vaccination reduces the potential of FMDV production in herds and flocks that may already have been exposed to infection [12]. Vaccination of exposed animals often arose due to the fact that it is difficult to detect an animal in this stage of disease as the acute phase of virus replication may be transient [12, 13]. By vaccinating all the exposed animals, it is hoped that those not already infected will develop sufficient immunity to provide at least partial protection against clinical disease [12]. Protective vaccination is used on herds and flocks that are in the vicinity of an outbreak but are thought not to have been exposed [12]. Prior studies suggest that a vaccinated exposed animal is highly likely to progress to be an FMD carrier [14]. Apart from suppressive and protective vaccines, our model explores the role of culling infectious symptomatic and nonsymptomatic animals.

2. Methods and Results

2.1. Model Framework

The total population of livestock is subdivided into proportions of susceptible animals , vaccinated animals , latently infected animals not vaccinated , exposed animals that have been administered a suppressive vaccine , infectious animals displaying clinical signs of the disease , and infectious animals not displaying clinical signs of the disease , also known as FMD carriers. Thus the total population is

Susceptible animals acquire FMD infection at rate , where denotes the FMD transmission parameter, which is considered to be a product of the between-animal contact rate —that is, the average number of contacts between animals per unit of time—and the FMDV transmission probability per each contact ; that is, . Factor , accounts for the reduced probability of FMD transmission by animals in class If it implies that animals in class do not transmit FMDV, while implies that animals in class have less chance of infecting susceptible animals. Animals in class acquire FMD infection at rate , where captures the impact of vaccination on reducing the susceptibility of animals in this class to FMD infection, that is, the vaccine efficacy. The model takes the following form:

where the upper dot represents the derivative of the component with respect to time. The constant parameter denotes birth or natural permanent exit of animals from the community. Susceptible animals acquire protective vaccination at per capita rate , denotes the incubation period of animals in class and this is usually in the range of 2–14 days [13], and is the proportion of exposed animals in class that receive suppressive vaccination and the remainder that do not receive suppressive vaccination progress to infectious carrier or symptomatic infectious animals . It is assumed that a fraction of nonvaccinated exposed animals progress to the infectious carrier population and the complementary proportion becomes symptomatic and infectious. Further, we assume that vaccinated exposed animals progress only to the infectious carrier population at rate ; accounts for the permanent exit rate of animals in class due to FMD infection, with a fraction progressing to infectious carrier population and the remainder succumbing to disease-related mortality; is the average life-span of an infectious carrier animal. Prior studies suggest that in cattle population an FMD carrier may exist for a period of 3.5 years [4]; denotes the detection and culling rate of infectious FMD carrier animals.

Our assumptions on the transfer of the population are demonstrated in Figure 1.

Figure 1: Model structure.

For biological reasons, we study the solutions of system (2) in the closed set:

is positively invariant. Thus all solutions of system (2) with nonnegative initial data will remain nonnegative for all time.

2.2. Equilibrium Points and Stability Analysis

Infection-Free Equilibrium Point. In the absence of FMD in the community, system (2) admits an equilibrium point known as infection-free, and it is given by

Following the next-generation method [16], the reproductive number of system (2) is

is a threshold for disease invasion or eradication, under suitable conditions, such as the absence of a backward bifurcation. See [17] for more discussion. Theorems 1 and 2 are based on computations in Appendix A.

Theorem 1. The FMD-free equilibrium point is globally asymptotically stable whenever .

Theorem 2. The unique endemic equilibrium of system (2) is locally asymptotically stable for but close to one.

Sensitivity Analysis of the Reproductive Number. Sensitivity analysis of model parameters is very important to design and control strategies as well as being a direction to future research. Local sensitivity indices allow us to measure the relative change in a state variable when a parameter changes. In computing the sensitivity analysis, we adopt the approach described by Arriola [18]. The normalized forward sensitivity index of a variable to a parameter is the ratio of the relative change in the variable to the relative change in the parameter. When the variable is a differentiable function of the parameter, the sensitivity index may be alternatively defined using partial derivatives.

Definition 3. The normalized forward sensitivity index of a variable, , that depends differentiably on a parameter, , is defined as

Model parameters whose sensitivity index values are near −1 or +1 suggest that a change in their magnitude has a significant impact on either increasing or decreasing the size of the reproductive number . From Table 2, it is clear that is most sensitive to , , , and . An increase in by 10% would increase by 10%. An increase in by 10% would decrease 9.9%. Similarly, if increases by 10%, then decreases by 9.9%. In summary, the numerical estimation of the indices suggests that reduction in disease-transmission rate coupled by an increase in preexposure vaccine efficacy and culling of infectious FMD carriers can lead to a significant reduction in new FMD cases.

2.3. Optimal Control

In this section, we formulate an optimal control problem for the transmission dynamics of FMD by extending the autonomous system (2). The goal here is to study the best strategies to curtail the epidemic. Three intervention methods, called controls, are included in our initial model. Controls are represented as functions of time and assigned reasonable upper and lower bounds. The first control attempts to strengthen the impact of vaccination and the second control attempts to strengthen the impact of detection and culling of infectious symptomatic infectious animals while the third control attempts to strengthen the impact of detection and culling of infectious nonsymptomatic animals. Using the same parameter and class names as in system (2) and Table 1, the system of differential equations describing our model with controls is

The objective functional, , is used to formulate the relevant optimization problem: finding the most effective strategy that reduces or eliminates the levels of FMD at minimal cost. This minimization goal will be achieved through the implementation of controls , , and over the preselected time interval . Mathematically, this corresponds to the minimization of the functional over a set of feasible, , strategies on . The functional is defined as follows:

where the constants , , are a measure of the relative cost of the interventions associated with the controls , , and , respectively.

Table 1: Model parameters and their interpretations.

Description
Symbol
Units
Baseline value
Ref.

Vaccine efficacy

—
0.5
[2]

Modification factor

—
0.4
Estimate

FMD transmission rate

Day⁻¹
0.6
[3]

Non-FMD related exit rate

Day⁻¹
0.001
[15]

Latency period for animals in class

Day⁻¹
0.26
[11]

Latency period for animals in class

Day⁻¹
0.07
[11]

Fraction of exposed animals vaccinated

—
0.65
Estimate

FMD-related death for animals in class

Day⁻¹
0.0008
[5]

FMD-related permanent exit from class

Day⁻¹
0.143
[3]

Detection and culling of infectious FMD carriers

Day⁻¹
0.25 (0–0.25)
[3]

Rate of vaccination for the susceptible population

Day⁻¹
0.25
[3]

Proportion of infectious animals detected and culled

—
0.85
Estimate

Proportion of infected animals that become carriers

—
0.45 (0.15–0.5)
[5]

Table 2: Sensitivity indices of to parameters for the FMD model (2), evaluated at the baseline parameter values given in Table 1.

Parameter
Sensitivity index

We consider state system (7) of ordinary differential equations in with the set of admissible control functions given by

We consider the optimal control problem of determining , , , , , and associated with an admissible control treble on the time interval , satisfying (7), given the initial conditions , , , , , and and minimizing the cost functional (8); that is,

The existence of optimal controls follows from standard results in optimal control theory [19]. Theorem 4 follows from Appendix B.

Theorem 4. Problem (7)–(10) with given initial conditions , , , , , and and fixed final time admits a unique optimal solution associated to an optimal control treble on

The optimal control treble predicted by Theorem 4 represents the optimal intervention strategy, given cost constraints, and can be found by the application of the Pontryagin maximum principle [19].

Numerical Results. In this section, we explore the role of optimal control on minimizing cumulative FMD infections in the community. To achieve this, we solve system (7) with a guess for the controls over the time interval using a forward fourth-order Runge-Kutta scheme and the transversality conditions , Then, system (B.2) is solved by a backward fourth-order Runge-Kutta scheme using the current iteration solution of system (7). The controls are updated by using a convex combination of the previous controls and the values from (B.3). The iteration is stopped when the values of the unknowns at the previous iteration are very close to the ones at the present iteration.

In all the simulations performed in this section, all parameters are fixed according to Table 1, together with the following estimated initial population levels and controls: , , , , , , , , and . For the weights, we assumed that vaccination is the most expensive intervention due to the fact that it involves a large number of animals; hence the cost associated with control is assumed to be higher than the cost tied to culling infectious animals, whether symptomatic or nonsymptomatic; thus . Further, we assumed that detection and culling of infectious nonsymptomatic animals involves a number of procedures such as taking blood samples for examination in laboratory since these animals do not display clinical signs of the disease, compared to detection and culling of symptomatic infectious animals; thus . Specifically, we have assigned our weights , , and .

Figure 2 illustrates the impact of time dependent intervention strategies on controlling FMD prevalence in the community over a period of 60 days. Results from these simulations demonstrate that the presence of optimal intervention strategies can lead to a significant decrease in cumulative FMD cases; the epidemiological classes and die out after 30 days and 40 days, respectively, of implementing the methods.

(a)

(b)

(c)

(d)

Figure 2: Time-series plots showing the impact of optimal intervention strategies on controlling FMD prevalence in the community, over a period of 60 days.

The solution for the optimal control problem (7)–(10) is illustrated in Figure 3. From these simulation results, we note that vaccination control, , and culling of infectious nonsymptomatic animals control, , must be maximum for the entire period of 60 days during which the interventions last, while the detection and culling of infectious symptomatic animals control, , should stay at its maximum intensity, for approximately 50 days. The total relative cost for all the strategies is estimated by computing the approximate area under control functions with daily costs computed by dividing the total costs by the number of days where nonzero controls are applied [20]. In this case, the total relative cost for the 60-day strategy—, , and —is $0.6125 with a daily cost of $0.0001.

(a)

(b)

(c)

Figure 3: Control profiles for control function: (a) , (b) , and (c) , over a period of 60 days.

Figure 4 shows the long-term dynamics of FMD in epidemiological classes and We note that, in the presence of controls, the epidemiological class dies out after 350 days of implementing these methods, while dies out after 200 days. The overall conclusion of the simulations is that all FMD cases can be eliminated from the community after 350 days of implementing the aforementioned control methods.

(a)

(b)

Figure 4: Time-series plots showing the impact of optimal intervention strategies on controlling FMD prevalence in the community, over a period of 500 days.

Figure 5 displays the time-dependent controls, , , and , over a period of 500 days. The simulations suggest that vaccination control, , and culling of infectious nonsymptomatic animals control, , must be maximized for the entire period of 500 days during which the interventions lasts, while detection and culling of infectious symptomatic animals control, , should stay at its maximum intensity, for approximately 65 days, and then can be progressively reduced. The total relative cost for the three control strategies over a long time horizon, 500 days, is computed as follows: , , and to get $0.045 with a daily cost of $0.00009. The costs for the short and long time horizons seem “roughly” comparable; in fact the daily costs are approximately the same.

(a)

(b)

(c)

Figure 5: Control profiles for control function: (a) , (b) , and (c) , over a period of 500 days.

Efficacy of Optimal Intervention Strategy. The efficacy of an intervention strategy on curbing new infections reflects the strength of the strategy to effectively control the epidemic. In this section, we explore the effectiveness of the aforementioned optimal intervention methods on reducing cumulative infectious FMD cases. We define the efficacy function by

where and denote the optimal solutions associated with the optimal control of the corresponding variable and and denote the corresponding initial condition. The function (11) measures the proportional decrease in the number of FMD infectious animals imposed by the intervention with controls, , , , by comparing the number of FMD infectious animals at time with the initial conditions for which there are no controls implemented; that is, . By construction, for all time . Thus the upper bound of is one.

Figure 6 illustrates the effectiveness of optimal intervention methods aimed at curtailing FMD in the community over a period of 500 days. We note that, after 50 days of implementing the strategies, the efficiency level will be above 90% and will reach the 100% mark after 300 days. This demonstrates that optimal intervention strategies can be effective in reducing or eliminating new FMD infections in the community.

Figure 6: Time-series plot demonstrating the efficacy of optimal intervention strategies over a period of 500 days.

3. Concluding Remarks

We have developed a dynamic model for foot-and-mouth disease (FMD). In this work, our research is focused on proposing the “optimal prevention and control strategy of FMD” from mathematical modeling. We have introduced three control mechanisms representing pre- and postexposure vaccination, culling of symptomatic animals, and culling of infectious nonsymptomatic animals into our model. The various strategies associated with these three controls have been investigated. Numerical simulations demonstrate that vaccination and case finding of infectious nonsymptomatic animals are the most effective controls. Hence, for effective control of FMD during an outbreak, these two controls should be maximized for the entire period.

If all three controls are used, then the number of latently infected (not vaccinated) and symptomatic will be almost zero in a period of 40 days. However, the number of animals in epidemiological classes and will be negligible in a period of 350 and 250 days, respectively. This shows that the optimal control strategy of FMD elimination can be effective in controlling the disease during an outbreak. Further, our analysis demonstrated that control, , (culling of infectious symptomatic animals) can be sustainable for a 100-day strategy; thereafter the control can be dropped.

Our model has several limitations, which should be acknowledged. We have assumed that infections can be transmitted through contact between an infectious and a susceptible animal, although airborne foot-and-mouth disease virus transmission has been documented [4]. Incorporating this aspect may bring a new dimension to our results. In practice, the movement of animals can be influenced by a number of factors such as seasonal variations. The current paper did not include such factors, though these might be as well worthwhile to model and analyze mathematically. Further, we assumed mass action incidence, although some researchers believe that there is little evidence that in ecology any form of contact among animals or individuals abides closely to this law [21].

Appendices

A. Stability Analysis of Model Steady States

A.1. Global Stability of the Disease-Free Equilibrium

Following Castillo-Chavez et al. (2002) [22], system (2) is presented in the form

where and . Here denotes the uninfected population and denotes the infected population. The infection-free equilibrium point is now denoted by where . We have to prove that following the two conditions

(H1)For , is globally asymptotically stable,(H2) for ,

are satisfied, where is a positively invariant attracting domain. Consider

Thus,

Since and at , it follows that . Hence is globally asymptotically stable whenever .

A.2. Endemic Equilibrium and Its Local Stability Analysis

Solving system (2) in terms of yields

Substituting and into , one gets

where

such that there exists an endemic equilibrium for model (2) if and only if there is a positive solution to . Because

there exists an endemic equilibrium if . Further, global stability of implies that is a unique endemic equilibrium point.

We now use the Centre Manifold Theory [23] to investigate the local stability of whenever but close to one.

The Jacobian matrix of system (2) evaluated about is

Let be the bifurcation parameter. From (A.8), we can deduce that the right eigenvectors are

Further, the left eigenvector of associated with the zero eigenvalue at is given by , where

It can easily be deduced that the bifurcation coefficients and are given by

Since and , then Theorem 2 holds.

B. Proof of Theorem 1

Proof. The Hamiltonian associated with the problem (7)–(10) is given by

where , , denotes the adjoint functions associated with the states , , , , , respectively. Note that, in , each adjoint function multiplies the right-hand side of the differential equation of its corresponding state function. The first term in comes from the integrand of the objective functional.
Given an optimal control triple , , and corresponding states , there exist adjoint functions satisfying

with transversality conditions , for . Furthermore, the optimal controls are characterised by

Conflict of Interests

The authors declare that there is no conflict of interests.

Acknowledgments

The authors would like to thank the anonymous reviewer and the editor for their suggestions that improved this paper.

References

G. J. Hughes, V. Mioulet, R. P. Kitching, M. E. J. Woolhouse, S. Alexandersen, and A. I. Donaldson, “Foot-and-mouth disease virus infection of sheep: implications for diagnosis and control,” The Veterinary Record, vol. 150, no. 23, pp. 724–727, 2002.

S. J. Cox and P. V. Barnett, “Experimental evaluation of foot-and-mouth disease vaccines for emergency use in ruminants and pigs: a review,” Veterinary Research, vol. 40, no. 3, article 13, 2009.

N. Ringa and C. T. Bauch, “Impacts of constrained culling and vaccination on control of foot and mouth disease in near-endemic settings: a pair approximation model,” Epidemics, vol. 9, pp. 18–30, 2014.

Aftosa Fever, “Foot-and-mouth disease,” May 2014, http://www.cfsph.iastate.edu

S. Alexandersen, Z. Zhang, and A. I. Donaldson, “Aspects of the persistence of foot-and-mouth disease virus in animals-the carrier problem,” Microbes and Infection, vol. 4, no. 10, pp. 1099–1110, 2002.

R. P. Kitching, “Identification of foot and mouth disease virus carrier and subclinically infected animals and differentiation from vaccinated animals,” OIE Revue Scientifique et Technique, vol. 21, no. 3, pp. 531–538, 2002.

H. Chang, Y. Ma, T. Lin, G. Cong, J. Du, and J. Ma, “Foot-and-mouth disease virus carrier status in Bos grunniens yaks,” Virology Journal, vol. 10, article 81, 2013.

M. J. Tildesley, N. J. Savill, D. J. Shaw et al., “Optimal reactive vaccination strategies for a foot-and-mouth outbreak in the UK,” Nature, vol. 440, no. 7080, pp. 83–86, 2006.

M. J. Keeling, M. E. J. Woolhouse, R. M. May, G. Davies, and B. T. Grenfell, “Modelling vaccination strategies against foot-and-mouth disease,” Nature, vol. 421, no. 6919, pp. 136–142, 2003.

E. Hansen and T. Day, “Optimal control of epidemics with limited resources,” Journal of Mathematical Biology, vol. 62, no. 3, pp. 423–451, 2011.

S. Mushayabasa, C. P. Bhunu, and M. Dhlamini, “Impact of vaccination and culling on controlling foot and mouth disease: a mathematical modelling approach,” World Journal of Vaccines, vol. 1, no. 4, pp. 156–161, 2011.

S. Kahn, D. W. Geale, P. R. Kitching, A. Bouffard, D. G. Allard, and R. J. Duncan, “Vaccination against foot-and-mouth disease: the implications for Canada,” Canadian Veterinary Journal, vol. 43, no. 5, pp. 349–354, 2002.

S. Parida, “Vaccination against foot-and-mouth disease virus: strategies and effectiveness,” Expert Review of Vaccines, vol. 8, no. 3, pp. 347–365, 2009.

J. Beale, “Foot and mouth disease: why not vaccinate?” Journal of the Royal Society of Medicine, vol. 94, no. 6, pp. 263–264, 2001.

G. E. Lahodny Jr., R. Gautam, and R. Ivanek, “Estimating the probability of an extinction or major outbreak for an environmentally transmitted infectious disease,” Journal of Biological Dynamics, vol. 9, supplement 1, pp. 128–155, 2014.

P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.

J. Li, D. Blakeley, and R. J. Smith?, “The failure of ,” Computational and Mathematical Methods in Medicine, vol. 2011, Article ID 527610, 17 pages, 2011.

L. Arriola, “Sensitivity analysis for uncertainty quantification in mathematical models,” in Mathematical and Statistical Estimation Approaches in Epidemiology, G. Chowell, J. M. Hyman, L. M. A. Bettencourt, and C. Castillo-Chavez, Eds., pp. 195–248, Springer, New York, NY, USA, 2009.

L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, and E. F. Mishchenko, The Mathematical Theory of Optimal Processes, Wiley, New Jersey, NJ, USA, 1962.

S. Lee, E. Jung, and C. Castillo-Chavez, “Optimal control intervention strategies in low- and high-risk problem drinking populations,” Socio-Economic Planning Sciences, vol. 44, no. 4, pp. 258–265, 2010.

J. Zhou and H. W. Hethcote, “Population size dependent incidence in models for diseases without immunity,” Journal of Mathematical Biology, vol. 32, no. 8, pp. 809–834, 1994.

C. Castillo-Chavez, Z. Feng, and W. Huang, “On the computation of and its role on global stability,” in Mathematical Approaches for Emerging and Re-emerging Infectious Diseases: An Introduction, C. Castillo- Chavez, S. Blower, P. Driessche, D. van Kirschner, and A.-A. Yakubu, Eds., vol. 125, pp. 229–250, IMA, 2002.

C. Castillo-Chavez and B. Song, “Dynamical models of tuberculosis and their applications,” Mathematical Biosciences and Engineering, vol. 1, no. 2, pp. 361–404, 2004.