Grey, G, and red, R, squirrel populations undergo the following interactions:
Grey squirrels are born at a rate proportional to their population, where the proportionality constant is b1.
Red squirrels are born at a rate proportional to their population, where the proportionality constant is b2.
Both squirrel populations compete internally with members of their own population (i.e. grey compete with grey
and red compete with red, respectively). Namely, whenever two members of the same population interact one of the
members is suppressed. The constant of proportionality is d1 for the grey squirrels and d2 for the red squirrels.
The squirrels also compete across species. Namely, a grey and a red squirrel interaction can lead to a reduction in
red squirrels, or a reduction in grey squirrels. The rate of proportionality for the reduction of grey squirrels is c1,
whilst the rate of proportionality for the reduction of red squirrels is c2.
1.
Write down the interaction equations specified by the above description.
Hint: You should create six equations over all and they should be symmetric in G and R (i.e. if you swap G and R
the dynamics should be the same).
2.
Show that the interaction equations lead to the following ODEs
Ġ
= b1G - d1G2- c1RG,
(1)
Ṙ
= b2R - d2R2- c2RG.
(2)
3.
What are the dimensions of b1, b2, d1, d2, c1 and c2 in terms of density and time?
4.
Non-dimensionalise the system to produce the following equations
=
u(1 - u - α12v),
(3)
=
ρv(1 - v - α21u).
(4)
Specify the population scales and the time scales and show that the scales have the appropriate dimensions.
5.
What are the parameter groupings ρ, α12 and α21 in terms of b1, b2, d1, d2, c1 and c2? Show that ρ, α12 and α21 are
dimensionless.
6.
Find all steady states of the system, noting any parameter dependencies that need to be satisfied for the steady state to
exist.
7.
Calculate the stability of the trivial steady states. Namely, those in which one of the populations is zero, noting parameter
dependencies of the stability criterion.
Calculating the stability of the fourth steady state,
(5)
leads to a complicated eigenvalue calculation and we will use curve sketching in the next question to consider the stability of
this steady state.
Hint: you should find there are four cases:
(a)
α12< 1 and α21< 1;
(b)
α12< 1 and α21> 1;
(c)
α12> 1 and α21< 1;
(d)
α12> 1 and α21> 1.
8.
Sketch the (u,v) phase plane in each of the cases (a)-(d). You should include: the nullclines; annotations regarding the signs
of du∕ dτ and dv∕ dτ in each region delineated by the nullclines; accompanying directional arrows in each region;
directional arrows on the nullclines; and, finally, trajectories illustrating the expected evolution with initial conditions taken
from each region.
Use these sketches to determine the stability of (us,vs), when it exists.
9.
Describe what happens to the grey and red squirrels in each of the cases (a)-(d).
2 Fitzhugh-Nagumo
In 1952 Alan Hodgkin and Andrew Huxley developed the first quantitative model of the propagation of an electrical signal along a
squid giant axon. They received the 1963 Nobel Prize in Physiology or Medicine for this work. Fitzhugh and Nagumo simplified
this work and showed how the essentials of the excitable process could be distilled into a simpler model, which is analytically
tractable.
Let v be membrane voltage potential (note this can go negative) and let n be the cellular recovery variable. The
non-dimensional equation controlling voltage spikes in a nerve are given by
ϵ
= v(0.5 - v)(v - 1) - n,
(6)
ṅ
= v - γn,
(7)
where γ and ϵ are positive.
1.
Derive the steady states noting any parameter restrictions.
2.
Assume γ < 16, characterise the stability of any steady states that exist.
3.
Draw, with as much detail as possible, the (v,n) phase diagrams in the γ < 16 case.
3 Computer simulation
Simulate equations (6) and (7) with parameters γ = 2 and ϵ = 0.001. Specifically, plot the trajectories in the (v,n) space along
with the nullclines. Further, initialise the trajectories near the (0,0) steady state but with negative v initially. Hint: for this
simulation you will have to use a ‘stiff’ solver. Matlab’s standard solver is ode23s.
Suppose we have a trajectory that starts at the origin and we increase γ past 16 and then lower it again. Could this
model exhibit hysteresis?
Figure 2: Simulation of equations (6) and (7) with parameters γ = 2 and ϵ = 0.001.
Exam revision
4 Damped pendulum
The equation of a damped pendulum is
(8)
where k and η are positive constants.
1.
Write the second order differential equation (8) as two first order ODEs.
2.
Identify the steady states of your two variable system, noting any parameter dependencies.
3.
Characterise the stability of the steady states, noting any parameter dependencies.
4.
Draw the phase diagram of the system in the (θ,) plane. Include nullclines, trajectory directions, stability
information of the steady states and draw a number of trajectories to illustrate potential dynamics.
5 Similar systems
Consider the two following systems:
1
= u2,
(9)
2
= -u1,
(10)
and
1
= v1v2,
(11)
2
= -v12.
(12)
1.
Show that the quantities of u12+ u22 and v12+ v22 are conserved across their respective systems. What does this
mean?
2.
Illustrate the differences between the two systems by calculating the steady states, characterising their stability and
drawing the phase planes with as much information as possible.
=
=
) plane. Include nullclines, trajectory directions, stability
information of the steady states and draw a number of trajectories to illustrate potential dynamics.
1
2
1
2