This is an introductory example from Olinick's book  An Introduction to Mathematical Models in the Social and Life Sciences . ​  "A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categori
This is an introductory example from Olinick's book An Introduction to Mathematical Models in the Social and Life Sciences. ​

"A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categories: low (below 6 pounds), average (between 6 and 8 pounds), and high (above 8 pounds). Among women whose own birth weights were low, 50 percent of the daughters had low birth weights, 45 percent had average weights, and 5 percent had high weights. Women with average birth weights had daughters with average weights half of the time, while the half was split evenly between low and high categories. Women with high birth weights had female babies with high weights 40 percent of the time, with low and average weights each occuring 30 percent of the time." p. 274-275.

For the Markov chain, you should make sure that you're taking time steps of length 1 in the settings, and Euler. RK-4 effectively looks beyond a single previous step, so it has a sort of memory!

Thanks Mike! Interesting examples, as always....
Andy Long

Next up: an SIR.
This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters .   We're rescaled this SIR model, so that time is given in infection rate-appropriate time units,
This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters.

We're rescaled this SIR model, so that time is given in infection rate-appropriate time units, "rates" are now ratios of rates (with infectivity rate in the denominator), and populations are considered proportions (unfortunately InsightMaker doesn't function properly if I give them all values from 0 to 1, which sum to 1 -- so, at the moment, I give them values that sum to 100, and consider the results percentages).

The new display includes the asymptotics: the three sub-populations will tend to fixed values as time goes to infinity; the infected population goes to zero if the recovery rate is greater than the infectivity rate -- i.e., the disease dies out.

Note the use of a "ghost" stock (for Total Population), which I think is a pretty cool idea. It cuts down on the number of arcs in the model graph.

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-rescaled.nb

This (simplest!) model demonstrates logistic growth.The original differential equation looks like  y'(t) = b y(t) (1 - y(t)/K)   where K is the carrying capacity of the quantity y.       But if we divide each side of the equation by K, we obtain      d(y/K)/dt = b (y/K) (1-y/K)     Defining a new va
This (simplest!) model demonstrates logistic growth.The original differential equation looks like

y'(t) = b y(t) (1 - y(t)/K)

where K is the carrying capacity of the quantity y.

But if we divide each side of the equation by K, we obtain

d(y/K)/dt = b (y/K) (1-y/K)

Defining a new variable w, the population relative to its carrying capacity, we obtain

dw/dt = b w (1 - w)

Finally we divide both sides by b, to write

dw/d(bt) = w (1 - w)

So if we work in dimensionless time units of bt, we have

w' = w (1 - w)

where the derivative is with respect to the variable bt=τ. .
τ=τ
This
       This equation, as simple as possible, contains all the dynamics (all the ways the population can behave), while masking the "trivialities"; but it kind of hides the physical aspects of the problem. So it's easy to study, but harder to interpret: alas, you can't have it all!:) 

τ=1 when t=1b: so if b=.5/year, then τ=1 when t=2.

So the larger b (the greater the birthrate), the shorter the real time t to give τ=1.
τ=τ=

τ=

 MAT375: Non-linear Exam....      This insight implements Newton's method as an InsightMaker model.       It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)      Fun to try a couple of different cases, so I have built four choices into this exa
MAT375: Non-linear Exam....

This insight implements Newton's method as an InsightMaker model.

It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)

Fun to try a couple of different cases, so I have built four choices into this example. You can choose the function ("Function Choice" of 0, 1, 2, or 3) using the slider.

Andy Long
Spring, 2020




This is a simple example of (part of a) simple SIR (Susceptible, Infected, Recovered) model, suggested by De Vries,  et al . in  A Course in Mathematical Biology.    They wanted to illustrate the comparative behavior of differential equations and discrete difference equations. We know that different
This is a simple example of (part of a) simple SIR (Susceptible, Infected, Recovered) model, suggested by De Vries, et al. in A Course in Mathematical Biology.

They wanted to illustrate the comparative behavior of differential equations and discrete difference equations. We know that differential equations are generally solved numerically by discretizing them, so that the comparison is a little bit rigged....

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-w-discrete-version.nb

This is an introductory example from Olinick's book  An Introduction to Mathematical Models in the Social and Life Sciences . ​  "A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categori
This is an introductory example from Olinick's book An Introduction to Mathematical Models in the Social and Life Sciences. ​

"A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categories: low (below 6 pounds), average (between 6 and 8 pounds), and high (above 8 pounds). Among women whose own birth weights were low, 50 percent of the daughters had low birth weights, 45 percent had average weights, and 5 percent had high weights. Women with average birth weights had daughters with average weights half of the time, while the half was split evenly between low and high categories. Women with high birth weights had female babies with high weights 40 percent of the time, with low and average weights each occuring 30 percent of the time." p. 274-275.

For the Markov chain, you should make sure that you're taking time steps of length 1 in the settings, and Euler. RK-4 effectively looks beyond a single previous step, so it has a sort of memory!

Thanks Mike! Interesting examples, as always....
Andy Long

Next up: an SIR.
    Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")  Tags:  Education ,  Chaos ,  Ecology ,  Biology ,  Population   Thanks to Insight Author:  John Petersen       Edits by Andy Long     Everything that follows the dashes was created by John Petersen (or at least came from his Insight model).

Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")
Thanks to Insight Author: John Petersen

Edits by Andy Long

Everything that follows the dashes was created by John Petersen (or at least came from his Insight model). I just wanted to make a few comments.

We are looking at Hare and Lynx, of course. Clone this insight, and change the names.

Then read the text below, to get acquainted with one of the most important and well-known examples of a simple system of differential equations in all of mathematics.

http://www.nku.edu/~longa/classes/mat375/mathematica/Lotka-Volterra.nb
------------------------------------------------------------

Dynamic simulation modelers are particularly interested in understanding and being able to distinguish between the behavior of stocks and flows that result from internal interactions and those that result from external forces acting on a system. 

For some time modelers have been particularly interested in internal interactions that result in stable oscillations in the absence of any external forces acting on a system. 

The model in this last scenario was independently developed by Alfred Lotka (1924) and Vito Volterra (1926).  Lotka was interested in understanding internal dynamics that might explain oscillations in moth and butterfly populations and the parasitoids that attack them.  Volterra was interested in explaining an increase in coastal populations of predatory fish and a decrease in their prey that was observed during World War I when human fishing pressures on the predator species declined. 

Both discovered that a relatively simple model is capable of producing the cyclical behaviors they observed. 

Since that time, several researchers have been able to reproduce the modeling dynamics in simple experimental systems consisting of only predators and prey.  It is now generally recognized that the model world that Lotka and Volterra produced is too simple to explain the complexity of most predator-prey dynamics in nature.  And yet, the model significantly advanced our understanding of the critical role of feedback in predator-prey interactions and in feeding relationships that result in community dynamics.

The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.

As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


    Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")  Tags:  Education ,  Chaos ,  Ecology ,  Biology ,  Population   Thanks to Insight Author:  John Petersen       Edits by Andy Long     Everything that follows the dashes was created by John Petersen (or at least came from his Insight model).

Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")
Thanks to Insight Author: John Petersen

Edits by Andy Long

Everything that follows the dashes was created by John Petersen (or at least came from his Insight model). I just wanted to make a few comments.

We are looking at Hare and Lynx, of course. Clone this insight, and change the names.

Then read the text below, to get acquainted with one of the most important and well-known examples of a simple system of differential equations in all of mathematics.

http://www.nku.edu/~longa/classes/mat375/mathematica/Lotka-Volterra.nb
------------------------------------------------------------

Dynamic simulation modelers are particularly interested in understanding and being able to distinguish between the behavior of stocks and flows that result from internal interactions and those that result from external forces acting on a system. 

For some time modelers have been particularly interested in internal interactions that result in stable oscillations in the absence of any external forces acting on a system. 

The model in this last scenario was independently developed by Alfred Lotka (1924) and Vito Volterra (1926).  Lotka was interested in understanding internal dynamics that might explain oscillations in moth and butterfly populations and the parasitoids that attack them.  Volterra was interested in explaining an increase in coastal populations of predatory fish and a decrease in their prey that was observed during World War I when human fishing pressures on the predator species declined. 

Both discovered that a relatively simple model is capable of producing the cyclical behaviors they observed. 

Since that time, several researchers have been able to reproduce the modeling dynamics in simple experimental systems consisting of only predators and prey.  It is now generally recognized that the model world that Lotka and Volterra produced is too simple to explain the complexity of most predator-prey dynamics in nature.  And yet, the model significantly advanced our understanding of the critical role of feedback in predator-prey interactions and in feeding relationships that result in community dynamics.

The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.

As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


 MAT375: Non-linear Exam....      This insight implements Newton's method as an InsightMaker model.       It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)      Fun to try a couple of different cases, so I have built four choices into this exa
MAT375: Non-linear Exam....

This insight implements Newton's method as an InsightMaker model.

It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)

Fun to try a couple of different cases, so I have built four choices into this example. You can choose the function ("Function Choice" of 0, 1, 2, or 3) using the slider.

Andy Long
Spring, 2020




This is a simple example of (part of a) simple SIR (Susceptible, Infected, Recovered) model, suggested by De Vries,  et al . in  A Course in Mathematical Biology.    They wanted to illustrate the comparative behavior of differential equations and discrete difference equations. We know that different
This is a simple example of (part of a) simple SIR (Susceptible, Infected, Recovered) model, suggested by De Vries, et al. in A Course in Mathematical Biology.

They wanted to illustrate the comparative behavior of differential equations and discrete difference equations. We know that differential equations are generally solved numerically by discretizing them, so that the comparison is a little bit rigged....

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-w-discrete-version.nb

This is an introductory example from Olinick's book  An Introduction to Mathematical Models in the Social and Life Sciences . ​  "A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categori
This is an introductory example from Olinick's book An Introduction to Mathematical Models in the Social and Life Sciences. ​

"A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categories: low (below 6 pounds), average (between 6 and 8 pounds), and high (above 8 pounds). Among women whose own birth weights were low, 50 percent of the daughters had low birth weights, 45 percent had average weights, and 5 percent had high weights. Women with average birth weights had daughters with average weights half of the time, while the half was split evenly between low and high categories. Women with high birth weights had female babies with high weights 40 percent of the time, with low and average weights each occuring 30 percent of the time." p. 274-275.

For the Markov chain, you should make sure that you're taking time steps of length 1 in the settings, and Euler. RK-4 effectively looks beyond a single previous step, so it has a sort of memory!

Thanks Mike! Interesting examples, as always....
Andy Long

Next up: an SIR.
This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters .   We're rescaled this SIR model, so that time is given in infection rate-appropriate time units,
This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters.

We're rescaled this SIR model, so that time is given in infection rate-appropriate time units, "rates" are now ratios of rates (with infectivity rate in the denominator), and populations are considered proportions (unfortunately InsightMaker doesn't function properly if I give them all values from 0 to 1, which sum to 1 -- so, at the moment, I give them values that sum to 100, and consider the results percentages).

The new display includes the asymptotics: the three sub-populations will tend to fixed values as time goes to infinity; the infected population goes to zero if the recovery rate is greater than the infectivity rate -- i.e., the disease dies out.

Note the use of a "ghost" stock (for Total Population), which I think is a pretty cool idea. It cuts down on the number of arcs in the model graph.

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-rescaled.nb

This is an introductory example from Olinick's book  An Introduction to Mathematical Models in the Social and Life Sciences . ​  Next up: and SIR, and his interesting model of female birth weights.
This is an introductory example from Olinick's book An Introduction to Mathematical Models in the Social and Life Sciences. ​

Next up: and SIR, and his interesting model of female birth weights.
 MAT375: Non-linear Exam....      This insight implements Newton's method as an InsightMaker model.       It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)      Fun to try a couple of different cases, so I have built four choices into this exa
MAT375: Non-linear Exam....

This insight implements Newton's method as an InsightMaker model.

It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)

Fun to try a couple of different cases, so I have built four choices into this example. You can choose the function ("Function Choice" of 0, 1, 2, or 3) using the slider.

Andy Long
Spring, 2020




This is an introductory example from Olinick's book  An Introduction to Mathematical Models in the Social and Life Sciences . ​  "A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categori
This is an introductory example from Olinick's book An Introduction to Mathematical Models in the Social and Life Sciences. ​

"A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categories: low (below 6 pounds), average (between 6 and 8 pounds), and high (above 8 pounds). Among women whose own birth weights were low, 50 percent of the daughters had low birth weights, 45 percent had average weights, and 5 percent had high weights. Women with average birth weights had daughters with average weights half of the time, while the half was split evenly between low and high categories. Women with high birth weights had female babies with high weights 40 percent of the time, with low and average weights each occuring 30 percent of the time." p. 274-275.

For the Markov chain, you should make sure that you're taking time steps of length 1 in the settings, and Euler. RK-4 effectively looks beyond a single previous step, so it has a sort of memory!

Thanks Mike! Interesting examples, as always....
Andy Long

Next up: an SIR.
This is a simple example of (part of a) simple SIR (Susceptible, Infected, Recovered) model, suggested by De Vries,  et al . in  A Course in Mathematical Biology.    They wanted to illustrate the comparative behavior of differential equations and discrete difference equations. We know that different
This is a simple example of (part of a) simple SIR (Susceptible, Infected, Recovered) model, suggested by De Vries, et al. in A Course in Mathematical Biology.

They wanted to illustrate the comparative behavior of differential equations and discrete difference equations. We know that differential equations are generally solved numerically by discretizing them, so that the comparison is a little bit rigged....

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-w-discrete-version.nb

This is an introductory example from Olinick's book  An Introduction to Mathematical Models in the Social and Life Sciences . ​  "A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categori
This is an introductory example from Olinick's book An Introduction to Mathematical Models in the Social and Life Sciences. ​

"A recent study focused on the relationship between the birth weights of English women and the birth weights of their daughters. The weights were split into three categories: low (below 6 pounds), average (between 6 and 8 pounds), and high (above 8 pounds). Among women whose own birth weights were low, 50 percent of the daughters had low birth weights, 45 percent had average weights, and 5 percent had high weights. Women with average birth weights had daughters with average weights half of the time, while the half was split evenly between low and high categories. Women with high birth weights had female babies with high weights 40 percent of the time, with low and average weights each occuring 30 percent of the time." p. 274-275.

For the Markov chain, you should make sure that you're taking time steps of length 1 in the settings, and Euler. RK-4 effectively looks beyond a single previous step, so it has a sort of memory!

Thanks Mike! Interesting examples, as always....
Andy Long

Next up: an SIR.
 MAT375: Non-linear Exam....      This insight implements Newton's method as an InsightMaker model.       It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)      Fun to try a couple of different cases, so I have built four choices into this exa
MAT375: Non-linear Exam....

This insight implements Newton's method as an InsightMaker model.

It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)

Fun to try a couple of different cases, so I have built four choices into this example. You can choose the function ("Function Choice" of 0, 1, 2, or 3) using the slider.

Andy Long
Spring, 2020




This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters .   We're rescaled this SIR model, so that time is given in infection rate-appropriate time units,
This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters.

We're rescaled this SIR model, so that time is given in infection rate-appropriate time units, "rates" are now ratios of rates (with infectivity rate in the denominator), and populations are considered proportions (unfortunately InsightMaker doesn't function properly if I give them all values from 0 to 1, which sum to 1 -- so, at the moment, I give them values that sum to 100, and consider the results percentages).

The new display includes the asymptotics: the three sub-populations will tend to fixed values as time goes to infinity; the infected population goes to zero if the recovery rate is greater than the infectivity rate -- i.e., the disease dies out.

Note the use of a "ghost" stock (for Total Population), which I think is a pretty cool idea. It cuts down on the number of arcs in the model graph.

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-rescaled.nb

This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters .   We're rescaled this SIR model, so that time is given in infection rate-appropriate time units,
This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters.

We're rescaled this SIR model, so that time is given in infection rate-appropriate time units, "rates" are now ratios of rates (with infectivity rate in the denominator), and populations are considered proportions (unfortunately InsightMaker doesn't function properly if I give them all values from 0 to 1, which sum to 1 -- so, at the moment, I give them values that sum to 100, and consider the results percentages).

The new display includes the asymptotics: the three sub-populations will tend to fixed values as time goes to infinity; the infected population goes to zero if the recovery rate is greater than the infectivity rate -- i.e., the disease dies out.

Note the use of a "ghost" stock (for Total Population), which I think is a pretty cool idea. It cuts down on the number of arcs in the model graph.

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-rescaled.nb

 MAT375: Non-linear Exam....      This insight implements Newton's method as an InsightMaker model.       It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)      Fun to try a couple of different cases, so I have built four choices into this exa
MAT375: Non-linear Exam....

This insight implements Newton's method as an InsightMaker model.

It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)

Fun to try a couple of different cases, so I have built four choices into this example. You can choose the function ("Function Choice" of 0, 1, 2, or 3) using the slider.

Andy Long
Spring, 2020




 MAT375: Non-linear Exam....      This insight implements Newton's method as an InsightMaker model.       It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)      Fun to try a couple of different cases, so I have built four choices into this exa
MAT375: Non-linear Exam....

This insight implements Newton's method as an InsightMaker model.

It is important to use Euler's method, with step-size of 1. That's what allows us to get away with this!:)

Fun to try a couple of different cases, so I have built four choices into this example. You can choose the function ("Function Choice" of 0, 1, 2, or 3) using the slider.

Andy Long
Spring, 2020




    Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")  Tags:  Education ,  Chaos ,  Ecology ,  Biology ,  Population   Thanks to Insight Author:  John Petersen       Edits by Andy Long     Everything that follows the dashes was created by John Petersen (or at least came from his Insight model).

Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")
Thanks to Insight Author: John Petersen

Edits by Andy Long

Everything that follows the dashes was created by John Petersen (or at least came from his Insight model). I just wanted to make a few comments.

We are looking at Hare and Lynx, of course. Clone this insight, and change the names.

Then read the text below, to get acquainted with one of the most important and well-known examples of a simple system of differential equations in all of mathematics.

http://www.nku.edu/~longa/classes/mat375/mathematica/Lotka-Volterra.nb
------------------------------------------------------------

Dynamic simulation modelers are particularly interested in understanding and being able to distinguish between the behavior of stocks and flows that result from internal interactions and those that result from external forces acting on a system. 

For some time modelers have been particularly interested in internal interactions that result in stable oscillations in the absence of any external forces acting on a system. 

The model in this last scenario was independently developed by Alfred Lotka (1924) and Vito Volterra (1926).  Lotka was interested in understanding internal dynamics that might explain oscillations in moth and butterfly populations and the parasitoids that attack them.  Volterra was interested in explaining an increase in coastal populations of predatory fish and a decrease in their prey that was observed during World War I when human fishing pressures on the predator species declined. 

Both discovered that a relatively simple model is capable of producing the cyclical behaviors they observed. 

Since that time, several researchers have been able to reproduce the modeling dynamics in simple experimental systems consisting of only predators and prey.  It is now generally recognized that the model world that Lotka and Volterra produced is too simple to explain the complexity of most predator-prey dynamics in nature.  And yet, the model significantly advanced our understanding of the critical role of feedback in predator-prey interactions and in feeding relationships that result in community dynamics.

The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.

As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters .   We're rescaled this SIR model, so that time is given in infection rate-appropriate time units,
This is an example of an SIR (Susceptible, Infected, Recovered) model that has been re-parameterized down to the bare minimum, to illustrated the dynamics possible with the fewest number of parameters.

We're rescaled this SIR model, so that time is given in infection rate-appropriate time units, "rates" are now ratios of rates (with infectivity rate in the denominator), and populations are considered proportions (unfortunately InsightMaker doesn't function properly if I give them all values from 0 to 1, which sum to 1 -- so, at the moment, I give them values that sum to 100, and consider the results percentages).

The new display includes the asymptotics: the three sub-populations will tend to fixed values as time goes to infinity; the infected population goes to zero if the recovery rate is greater than the infectivity rate -- i.e., the disease dies out.

Note the use of a "ghost" stock (for Total Population), which I think is a pretty cool idea. It cuts down on the number of arcs in the model graph.

A comparable model in Mathematica is available at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/SIRModel-rescaled.nb