- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Abstract and Applied Analysis

Volume 2014 (2014), Article ID 169609, 12 pages

http://dx.doi.org/10.1155/2014/169609

## Nonlinear Dynamic in an Ecological System with Impulsive Effect and Optimal Foraging

^{1}School of Life and Environmental Science, Wenzhou University, Wenzhou, Zhejiang 325035, China^{2}Zhejiang Provincial Key Laboratory for Water Environment and Marine Biological Resources Protection, Wenzhou University, Wenzhou, Zhejiang 325035, China

Received 13 January 2014; Accepted 15 May 2014; Published 16 June 2014

Academic Editor: Weiming Wang

Copyright © 2014 Min Zhao and Chuanjun Dai. 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

The population dynamics of a three-species ecological system with impulsive effect are investigated. Using the theories of impulsive equations and small-amplitude perturbation scales, the conditions for the system to be permanent when the number of predators released is less than some critical value can be obtained. Furthermore, because the predator in the system follows the predictions of optimal foraging theory, it follows that optimal foraging promotes species coexistence. In particular, the less beneficial prey can support the predator alone when the more beneficial prey goes extinct. Moreover, the influences of the impulsive effect and optimal foraging on inherent oscillations are studied using simulation, which reveals rich dynamic behaviors such as period-halving bifurcations, a chaotic band, a periodic window, and chaotic crises. In addition, the largest Lyapunov exponent and the power spectra of the strange attractor, which can help analyze the chaotic dynamic behavior of the model, are investigated. This information will be useful for studying the dynamic complexity of ecosystems.

#### 1. Introduction

In recent years, interest in studying nonlinear dynamic systems has exploded. In the 1970s, since the pioneering work of May on the relationship between food-web complexity and stability and the chaotic phenomenon [1–3], more and more researchers have become interested in dynamic behavior involving ecological mechanisms that promote species diversity [4–20]. More recently, dynamic systems' studies have benefited from an infusion of interest and new techniques in ecology.

It is known that when a predator is shared by two noncompeting species, predator-mediated apparent competition often leads to competitive exclusion of one prey population [21]. This phenomenon is related to optimal foraging and adaptive foraging. A two-prey-one-predator population model with optimal predator foraging behavior has been studied in a fine-grained environment [22–24]. On this basis, Křivan and Eisner considered a system composed of two prey types and an optimally foraging predator [25] in a system described by the following model:

This paper considers an impulsive differential-equation model based on model (1), which assumes that predators forage according to optimal foraging theory [23, 24]. This system can be expressed by the following equations: where , , are, respectively, the densities of two prey types and one predator at time , is the per capita prey intrinsic growth rate, , are the respective carrying capacities of the prey, and , are the corresponding values of available resources or in the ideal case (i.e., where no resources are wasted) the carrying capacity. However, the ideal case is impossible in reality. The ratios , express the relative efficiency of nutrient utilization in species , . At any time, if , , the efficiency is high as long as , are close to one; when the values are lower, this indicates that resource limitations are restricting the population increase [23]. are the rate of intraspecific competition of the prey, is the cropping rate of a predator feeding on the th prey type, is the conversion factor relating predator reproduction to prey consumption, and is the per capita mortality rate for the forager. In this paper, it is assumed that prey type is more beneficial than the other and hence [26, 27]. To study optimal foraging, a control parameter is introduced [25], which represents the probability that the alternative second prey type is included in the predator’s diet. is the period of the impulsive effect, , and is the number of predators released at . To achieve a set of conditions which can guarantee that the system will be permanent and that the numbers of the two prey types are not so large that they go extinct because of exceeding the carrying capacity of the environment, the model will release a certain number of predators only at because the predator is assumed to be a versatile and advanced predator.

The rest of this paper is organized as follows. Section 2 will review the effect of impulsive perturbations, establish conditions for extinction, and obtain the conditions for permanence of System (2) using the Floquet theory of impulsive equations at small-amplitude perturbation scales. In Section 3, the results of computer-based numerical analysis are shown and discussed briefly. In addition, the largest Lyapunov exponent, which also indicates the chaotic dynamic behavior of the model, is computed, and the Fourier spectra, which illustrate the qualitative nature of strange attractors, are plotted. Finally, conclusions and remarks are stated.

#### 2. Analysis of the System

Let , , , and let be the set of all nonnegative integers. The map is defined by the right-hand side of the first three equations of System (2).

Let ; then is said to belong to class if(1) is continuous in , and for each , , exists;(2) is locally Lipschitzian in .

*Definition 1. *Let ; then, for , the upper right derivative of with respect to the impulsive differential System (2) can be defined as

The solution of System (2) is a piecewise continuous function , where is continuous on , , and exists. Obviously, the smoothness property of guarantees the global existence and uniqueness of a solution of System (2) (for details, see [28–30]).

*Definition 2. *System (2) is said to be permanent if there exists a compact region such that every solution of System (2) will eventually enter and remain in the region .

The following lemma will now be presented.

Lemma 3. *Let be a solution of System (2) with ; then for all , and furthermore , if .*

An important comparison theorem will now be used on the impulsive differential equation.

Lemma 4 (see [28–30]). *Let and assume that
**
where is continuous in and for , , exists and is nondecreasing. Let be the maximal solution to the scalar impulsive differential equation
**
existing on . Then implies that , , where is any solution to System (2).*

Theorem 5. *There exists a constant such that , , and for each solution of System (2) for all large enough.*

*Proof. *Define such that
where . Since and , then , , and the upper right derivative of can be calculated along a solution of System (2), yielding the following impulsive differential equation:
Obviously,

Let ; then is bounded. Select , such that
where , are two positive constants.

According to Lemma 4,
where . Hence,

Therefore, is ultimately bounded, and it follows that each positive solution of System (2) is uniformly ultimately bounded. This completes the proof.

Next, some basic properties of the following subsystem of System (2), in which the two prey types are absent, will be defined:

Clearly, , , is a positive periodic solution of System (12). Hence, is a solution of System (12) with initial value , where , .

Lemma 6. *For a positive periodic solution of System (12) and every solution of System (12) with , .*

Hence, when only the predator is present, it is possible to obtain the complete expression for the periodic solution of System (2).

Based on these discussions, the following theorems can be proved.

Theorem 7. *Let be any solution of System (2). Then*(1)* is said to be locally asymptotically stable if ;*(2)* is said to be globally asymptotically stable if and
*

*Proof. *The local stability of the periodic solution may be determined by considering the behavior of small-amplitude perturbations of the solution. Define
Then substitute (15) into System (2). The linearization of the system becomes
Therefore,
where satisfies
and , the identity matrix, and
The stability of the periodic solution is determined by the eigenvalues of
which are

According to Floquet theory, is locally asymptotically stable if and ; that is, .

If is locally asymptotically stable and a global attractor, then is globally asymptotically stable. In the following, global attractiveness will be demonstrated.

Let ; then

By Theorem 5, there exists a constant such that , , and for each solution of System (2) with all large enough. Therefore,

By Lemmas 4 and 6, it is known that there exists and it is possible to select small enough so that . Therefore, for all ,
Define

Then and . Thus, for ,

So , , as . Note that the limiting case of System (2) is exactly System (12) together with Lemma 6. It follows that the periodic solution is a global attractor. This completes the proof.

Theorem 8. *System (2) is permanent if
*

*Proof. *Let be any solution of System (2) with . From Theorem 5, assume that .

From System (2), .

Consider the following equation:
It is possible to obtain and as . Hence, for any , for all sufficiently large. For simplification, it may be assumed that holds for all . The same arguments can be made for any , for all . Let . Note that , and consider the following equation:

From Lemmas 4 and 6, , and hence for all sufficiently large. Therefore, it is necessary to find and such that , for all large enough. This can be done as shown in the following two steps.*First*, choose , such that . Then there exist and such that , . Otherwise,(1)there exists a such that , but for all ;(2)there exists a such that , but for all ;(3)there are and for all .

For case (1), according to Theorem 8, select small enough so that

From case (1),
where . Therefore, and as , where is the solution of the following equation:
and . Therefore, there exists a , when ,

Let and . Integrating (34) on , , the following result can be obtained:
Then as , which is a contradiction.

For case (2), the same arguments can be used.

Now consider case (3). Choose small enough so that

From case (3),
where . Therefore, and as , where is the solution of the following equation:
and . Then there exists a , when , such that

Let and . Integrating (40) on , , the following result can be obtained:
Then as , which is a contradiction.

In conclusion, there exist and such that .*Second*, if , for all , then the objective has been attained. Otherwise, there exists such that , for . Let . Then , for and , , and , because is continuous. Choose such that
where . Let ; then there must be a such that ; otherwise , . Considering (32) with ,
for and . Then, for ,
It can be concluded that (34) holds for . As in the first step above, it is possible to obtain . There are two possible cases for .*Case (**)*. If for , then for . From System (2),
Integrating (45) on , .

Thus
which is a contradiction.

Let , so that and (45) holds on . Integrating (45) on ,
For , the same arguments can be used because .*Case (**)*. There exists a such that ; let ; then for and . For , (45) holds true. Integrating (45) on ,
For , the same arguments can be used because .

In summary, can be obtained for all . In the same way, it can be proved that for all . This completes the proof.

#### 3. Numerical Analysis

##### 3.1. The Impulsive Effect and Optimal Foraging

To study the population dynamics of a three-species ecological model with impulsive effect, the solution of System (2) with initial conditions in the first quadrant is obtained numerically for a biologically feasible range of parametric values, and the bifurcation diagram provides a summary of the basic population dynamic behavior of the system.

Now two different control parameters will be discussed, the number of predators released, , and the probability . Other parameters are set to

From Theorem 7, it is known that the prey-eradication periodic solution is locally asymptotically stable provided that . Figure 1 shows a typical prey-eradication periodic solution of System (2), in which it can be observed that the variable oscillates in a stable cycle. At the same time, the prey types and rapidly diminish and go to zero beyond . If the number of predators released, , is less than , the prey-eradication solution becomes unstable. It is, however, possible that the two prey types and the predator can coexist in a stable positive periodic solution. In other words, the system can be permanent when the number of predators released, , is less than .

Next, the bifurcation diagrams for the control parameter will be examined. Figure 2 is plotted as a function of the bifurcation parameter and shows that the system has rich population dynamic behavior consistent with the theoretical analysis, such as period-halving bifurcation (see Figure 3), a chaotic band, a periodic window, and chaotic crises. Furthermore, Theorem 8 indicates that the system is permanent when the value of is less than some critical value. When the value of is in the interval [0.009, 3.257895], the two prey types and one predator can coexist. When the value of is in the interval [3.257895, 4.098648], the prey will become extinct rapidly, but the prey and the predator can coexist. These results may show that prey is inferior to prey in its ability to reproduce or prey is a favorite food of predator . When the number of predators released is greater than some critical value, all species in the system will become extinct. All these results demonstrate the effectiveness of mathematical analysis for understanding such systems.

The next question is how impacts the complex population dynamics. In Figure 4, when prey and predator populations are plotted as a function of the probability , the value of is 1.45. In the former case, it is assumed that the foraging behavior of predator follows optimal foraging theory [16–18] and prey is more beneficial for predator than prey . In other words, the more beneficial prey is always included in the predator’s diet, but if the density of prey falls below a critical threshold or goes to zero, prey is included with probability one. From Figure 4, it can be clearly observed that the two prey types and the predator can coexist in the intervals [0, 0.1136] and [0.7215, 0.8792], where the system dynamics can be chaotic, periodic, or nonperiodic. In the interval (0.1136, 0.7215), prey goes extinct, while prey and predator can coexist stably. This means that if the more beneficial prey disappears, prey alone can support the population of predator . As increases, prey goes extinct, but prey and predator can coexist stably.

##### 3.2. The Largest Lyapunov Exponent

Deterministic chaos is an important problem that is solved by measuring the largest Lyapunov exponent [31–36]. Based on research by various investigators, these results have confirmed the importance of detecting and exploring chaos. In this paper, the largest Lyapunov exponents for chaotic system (2) are examined. The largest Lyapunov exponents take into account the average exponential rates of divergence or convergence of nearby orbits in phase space [31, 32]. If the attractor is chaotic, the largest Lyapunov exponent must be positive, which implies a stable or a periodic state. In Figure 2, the corresponding largest Lyapunov exponent ([]) can be calculated for System (2) (see Figure 5).

##### 3.3. The Strange Attractor and Power Spectra

To study the properties of strange attractors, commonly used methods such as power spectra can be used [35, 36]. A power spectrum was calculated using 4096 points corresponding to the time series of the variable with time increment [35, 36]. For strange attractors (a) and (b), it is known that the value of the largest Lyapunov exponent for the strange attractor (a) is 0.25603, while for (b) the computed largest Lyapunov exponent is 0.30567. Therefore, strange attractors (a) and (b) are chaotic attractors. Moreover, the strange attractor (b) displays more chaotic dynamics than (a) because its positive exponent is larger than that of (a). In addition, the spectra of strange attractors (a) and (b) consist of strong broadband components and sharp peaks (Figures 6(c) and 6(d)) These results conform to the observation that strange attractors (a) and (b) arise from some weak limit cycles which can lose stability due to noise.

#### 4. Conclusions and Remarks

Complex population dynamics of a three-species ecological model with optimal foraging and impulsive control strategy have been investigated both numerically and analytically. The periodic solution has been shown to be globally asymptotically stable by use of the Floquet theorem and small-amplitude perturbations, if and . At the same time, using the method of comparison involving multiple Lyapunov functions, the permanence of the system can be proved. Bifurcation diagrams of the impulsive perturbation and the probability parameter have also been obtained. The bifurcation diagrams of have shown that dynamic complexity exists in System (2), including chaotic behavior, periodic windows, chaotic bands, chaotic crises, and period-halving bifurcations. The bifurcation diagrams of indicate that optimal foraging promotes species coexistence and that if the more beneficial prey goes extinct, the less beneficial prey can support the predator so that it will not die out. In addition, the presence of chaotic dynamics was confirmed, and the qualitative nature of strange attractors was investigated using computer simulations of the largest Lyapunov exponents and Fourier spectra. All these results may be useful in the study of the dynamic complexity of ecosystems.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 31370381 and 31170338), by the Key Program of Zhejiang Provincial Natural Science Foundation of China (Grant no. LZ12C03001), and by the National Key Basic Research Program of China (973 Program, Grant no. 2012CB426510).

#### References

- R. M. May, “Will a large complex system be stable?”
*Nature*, vol. 238, no. 5364, pp. 413–414, 1972. View at Publisher · View at Google Scholar · View at Scopus - R. M. May, “Simple mathematical models with very complicated dynamics,”
*Nature*, vol. 261, no. 5560, pp. 459–467, 1976. View at Google Scholar · View at Scopus - R. M. May, “Biological populations with nonoverlapping generations: stable points, stable cycles, and chaos,”
*Science*, vol. 186, no. 4164, pp. 645–647, 1974. View at Google Scholar · View at Scopus - H. Zhu, S. A. Campbell, and G. S. K. Wolkowicz, “Bifurcation analysis of a predator-prey system with nonmonotonic functional response,”
*SIAM Journal on Applied Mathematics*, vol. 63, no. 2, pp. 636–682, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S.-B. Hsu, T.-W. Hwang, and Y. Kuang, “Global analysis of the Michaelis-Menten-type ratio-dependent predator-prey system,”
*Journal of Mathematical Biology*, vol. 42, no. 6, pp. 489–506, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Do, H. Baek, Y. Lim, and D. Lim, “A three-species food chain system with two types of functional responses,”
*Abstract and Applied Analysis*, vol. 2011, Article ID 934569, 16 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Lv and M. Zhao, “The dynamic complexity of a three species food chain model,”
*Chaos, Solitons and Fractals*, vol. 37, no. 5, pp. 1469–1480, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. Wang, J. Shen, and J. J. Nieto, “Permanence and periodic solution of predator-prey system with Holling type functional response and impulses,”
*Discrete Dynamics in Nature and Society*, vol. 2007, Article ID 81756, 15 pages, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Zhang and M. Zhao, “Dynamic complexities in a hyperparasitic system with prolonged diapause for host,”
*Chaos, Solitons and Fractals*, vol. 42, no. 2, pp. 1136–1142, 2009. View at Publisher · View at Google Scholar · View at Scopus - R. K. Upadhyay and R. K. Naji, “Dynamics of a three species food chain model with Crowley-Martin type functional response,”
*Chaos, Solitons and Fractals*, vol. 42, no. 3, pp. 1337–1346, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Zhao, L. Zhang, and J. Zhu, “Dynamics of a host-parasitoid model with prolonged diapause for parasitoid,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 1, pp. 455–462, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. Wei and L. Chen, “A delayed epidemic model with pulse vaccination,”
*Discrete Dynamics in Nature and Society*, vol. 2008, Article ID 746951, 12 pages, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Jiao, S. Cai, and L. Chen, “Analysis of a stage-structured predatory-prey system with birth pulse and impulsive harvesting at different moments,”
*Nonlinear Analysis: Real World Applications*, vol. 12, no. 4, pp. 2232–2244, 2011. View at Publisher · View at Google Scholar · View at MathSciNet - R. Shi and L. Chen, “Stage-structured impulsive $SI$ model for pest management,”
*Discrete Dynamics in Nature and Society*, vol. 2007, Article ID 97608, 11 pages, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H.-F. Huo, Z.-P. Ma, and C.-Y. Liu, “Persistence and stability for a generalized Leslie-Gower model with stage structure and dispersal,”
*Abstract and Applied Analysis*, vol. 2009, Article ID 135843, 17 pages, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Yu, S. Zhong, and R. P. Agarwal, “Mathematics analysis and chaos in an ecological model with an impulsive control strategy,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 2, pp. 776–786, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Wang, L. Chen, and J. J. Nieto, “The dynamics of an epidemic model for pest control with impulsive effect,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 3, pp. 1374–1386, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Yu, S. Zhong, R. P. Agarwal, and S. K. Sen, “Effect of seasonality on the dynamical behavior of an ecological system with impulsive control strategy,”
*Journal of the Franklin Institute: Engineering and Applied Mathematics*, vol. 348, no. 4, pp. 652–670, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Luo, “Permanence and extinction of a generalized Gause-type predator-prey system with periodic coefficients,”
*Abstract and Applied Analysis*, vol. 2010, Article ID 845606, 24 pages, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Wang, W. Wang, and X. Lin, “Dynamics of a two-prey one-predator system with Watt-type functional response and impulsive control strategy,”
*Chaos, Solitons and Fractals*, vol. 40, no. 5, pp. 2392–2404, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. D. Holt, “Predation, apparent competition and the structure of prey communities,”
*Theoretical Population Biology*, vol. 12, no. 2, pp. 197–229, 1977. View at Publisher · View at Google Scholar · View at MathSciNet - H. Werner, “Optimal foraging and the size selection of prey by the bluegill sunfish (Lepomis macrochirus),”
*Ecology*, vol. 55, pp. 1042–1052, 1974. View at Publisher · View at Google Scholar - E. L. Charnov, “Optimal foraging: attack strategy of a mantid,”
*The American Naturalist*, vol. 110, no. 971, pp. 141–151, 1976. View at Publisher · View at Google Scholar - D. W. Stephens and J. R. Krebs,
*Foraging Theory*, Princeton University Press, Princeton, NJ, USA, 1986. - V. Křivan and J. Eisner, “Optimal foraging and predator-prey dynamics III,”
*Theoretical Population Biology*, vol. 63, no. 4, pp. 269–279, 2003. View at Publisher · View at Google Scholar · View at Scopus - V. Křivan, “Optimal foraging and predator-prey dynamics,”
*Theoretical Population Biology*, vol. 49, no. 3, pp. 265–290, 1996. View at Publisher · View at Google Scholar · View at Scopus - V. Křivan and A. Sikder, “Optimal foraging and predator-prey dynamics, II,”
*Theoretical Population Biology*, vol. 55, no. 2, pp. 111–126, 1999. View at Publisher · View at Google Scholar · View at Scopus - V. Lakshmikantham, D. D. Baĭnov, and P. S. Simeonov,
*Theory of Impulsive Differential Equations*, vol. 6, World Scientific, Singapore, 1989. View at Publisher · View at Google Scholar · View at MathSciNet - D. D. Baĭnov and P. S. Simeonov,
*Impulsive Differential Equations: Asymptotic Properties of the Solutions*, World Scientific, Singapore, 1993. View at Publisher · View at Google Scholar · View at MathSciNet - D. D. Bainov and V. C. Covachev,
*Impulsive Differential Equations with a Small Parameter*, vol. 24, World Scientific, Singapore, 1994. View at Publisher · View at Google Scholar · View at MathSciNet - S. Lv and M. Zhao, “The dynamic complexity of a host-parasitoid model with a lower bound for the host,”
*Chaos, Solitons and Fractals*, vol. 36, no. 4, pp. 911–919, 2008. View at Publisher · View at Google Scholar · View at MathSciNet - M. Zhao and S. Lv, “Chaos in a three-species food chain model with a Beddington-DeAngelis functional response,”
*Chaos, Solitons and Fractals*, vol. 40, no. 5, pp. 2305–2316, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Zhao and L. Zhang, “Permanence and chaos in a host-parasitoid model with prolonged diapause for the host,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 12, pp. 4197–4203, 2009. View at Publisher · View at Google Scholar · View at Scopus - L. Zhu and M. Zhao, “Dynamic complexity of a host-parasitoid ecological model with the Hassell growth function for the host,”
*Chaos, Solitons and Fractals*, vol. 39, no. 3, pp. 1259–1269, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Yu, S. Zhong, R. P. Agarwal, and S. K. Sen, “Three-species food web model with impulsive control strategy and chaos,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 2, pp. 1002–1013, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Yu, S. Zhong, and M. Ye, “Dynamic analysis of an ecological model with impulsive control strategy and distributed time delay,”
*Mathematics and Computers in Simulation*, vol. 80, no. 3, pp. 619–632, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet