The hare/lynx interaction model is a typical two-species predator/prey model. This example has been borrowed from SystemDynamics 1.0. The original version of this example can be found in the STELLA tutorial.

---

References:

1. Richmond, B., S. Peterson, and P. Vescuso (1987), An Academic User's Guide to STELLA, High Performance Systems, Inc., Lyme, N.H.

---

Simulate the model across 60 time units and plot the two animal populations against simulation time:

You need to reduce the accuracy to 1e-10 in order for the simulation to run correctly.

Also plot the two population against each other:

You shall notice that the limit cycle is mildly unstable.

The simulation results obtained are not exactly the same as with SystemDynamics 1.0. The reason is that the previous version used tabular interpolation functions with linear extrapolation, whereas the new version uses tabular interpolation functions without extrapolation capability (a more prudent approach). In line with STELLA philosophy, the support values of the two tabular functions were manually extended to reflect a constant rather than a linear extrapolation.

hare	Type: Real Description: Hare population
lynx	Type: Real Description: Lynx population