30
Dynamic Modeling with Difference Equations
replace P + r P(1 P) by these approximations in Pt+1 = Pt +
r Pt (1 Pt ). Use this to determine the stability of the equilibria. Your
answer should agree
64
Linear Models of Structured Populations
in this matrix is saying about the population. Be careful in trying to
explain the meaning of the .11 in the upper left corner.
2.2.9. a. Show that Ax = Ay d
60
Linear Models of Structured Populations
Although we will not prove it here, it is possible to show that, for square
matrices, if Q P = I , then P Q = I . So, if Q is the inverse of P, then P is
the
4
Dynamic Modeling with Difference Equations
It may seem odd to call Pt+1 = (1 + f d)Pt a difference equation, when
the difference !P does not appear. However, the equations
Pt+1 = (1 + f d)Pt
and
!P
viii
Preface
Our writing style is intentionally informal. We have not tried to offer definitive coverage of any topic, but rather draw students into an interesting field.
In particular, we often only
32
Dynamic Modeling with Difference Equations
get started, try the MATLAB commands:
t=[0:50]
K=5+sin(2*pi/12)*t)
p=.1; pops=p
for i=1:50 p=p+.2*p*(1-p/K(i); pops=[pops p]; end
plot(t,K,t,pops)
Suggest
66
Linear Models of Structured Populations
Table 2.2. Linear
Model Simulation
with Eigenvector as
Initial Values
t
xt
0
1
2
3
.
.
v
Av = v
Av = 2 v
A2 v = 3 v
.
.
The practical consequence of this is
28
Dynamic Modeling with Difference Equations
the flour beetle tribolium, that exhibits chaotic dynamics (see (Cushing et al.,
2001).
Problems
1.3.1. The equilibrium points of a model are located wher
74
Linear Models of Structured Populations
the eigenvalues i are complex, if 1 is strictly dominant
so |1 | > |i | for
! !
! !
i = 2, 3, . . . , n, then by part (c) of the theorem, ! 1i ! < 1 as befor
12
Dynamic Modeling with Difference Equations
P/P
r
P
K
Figure 1.1. Per-capita growth rate as a function of population size.
Of course we cannot say exactly what a graph of !P/P should look
like with
46
Linear Models of Structured Populations
Matrices (the plural of matrix) are usually denoted by capital letters,
such as A, M, or P. For instance, we might say
!
"
.9925 .0125
P=
.0075 .9875
is the
80
Linear Models of Structured Populations
While our argument above also applies to larger matrices (provided you
learn to compute determinants of larger matrices in some other course), solving the ch
70
Linear Models of Structured Populations
Definition. An eigenvalue of A that is largest in absolute value is called a
dominant eigenvalue of A. An eigenvector corresponding to it is called a
dominan
68
Linear Models of Structured Populations
Technical remark: Not every matrix has eigenvectors that can be used as
columns to form an invertible matrix. However, it is possible to prove that if a
matr
54
Linear Models of Structured Populations
To capture the effects on population growth, we might begin modeling a
human population by creating five age classes with:
x1 (t) = no. of individuals age 0
16
Dynamic Modeling with Difference Equations
next_p = p+.7*p*(1-p/10)
20
Population at time t+1
18
16
14
12
10
8
6
4
2
0
0
2
4
6
8
10
12
14
16
18
20
Population at time t
Figure 1.4. Cobweb plot of a
6
Dynamic Modeling with Difference Equations
1.1.2. In the early stages of the development of a frog embryo, cell division
occurs at a fairly regular rate. Suppose you observe that all cells
divide, a
24
Dynamic Modeling with Difference Equations
When linearizing to determine stability, it is vital that you are focusing
on an equilibrium. Do not attempt to decide if a point is a stable or unstable
48
Linear Models of Structured Populations
W (Dx0 ), which we could compute relatively easily by matrix multiplication:
!
"!
" !
"
.9925 .0125
106.45
116.82
x2 =
.
.0075 .9875
893.55
883.18
A more int
14
Dynamic Modeling with Difference Equations
next_p = p+.7*p*(1p/10)
Population P
12
10
8
6
4
2
0
0
5
10
15
Time
Figure 1.2. Population values from a nonlinear model.
If we measure population size in
20
Dynamic Modeling with Difference Equations
intuitively? (Note that r will be very small, because we are using a
small time interval.) The logistic growth model is sometimes also
referred to as the
86
Nonlinear Models of Interactions
3.1. A Simple PredatorPrey Model
Imagine two species, one of which, the predator, preys on the other, the prey.
To keep things simple, we imagine that the predatorp
40
Dynamic Modeling with Difference Equations
Problems
1.5.1. (Calculus) The logistic differential equation is
dN
= r N (1 N /K ).
dt
a. Show that
N (t) =
K
1 + Cer t
where C =
K N0
N0
is a solution w
52
Linear Models of Structured Populations
find the following without a computer. Then check your answers with
MATLAB.
a. A + B
b. AB
c. B A
d. A2 = A A
e. 2A
f. Show C(A + B) = C A + C B.
! "
! "
r s
To J., R., and K.,
may reality live up to the model
MATHEMATICAL MODELS IN BIOLOGY
AN INTRODUCTION
ELIZABETH S. ALLMAN
Department of Mathematics and Statistics,
University of Southern Maine
JOHN A. RH