0% found this document useful (0 votes)
17 views3 pages

Mathematical Modeling Worksheet

The document is a worksheet for a mathematical modeling course at the University of Zimbabwe, covering various topics such as differential equations, logistic growth models, and systems of equations. It includes exercises on sketching phase lines, explaining mathematical modeling, finding fixed points, and analyzing stability. Additionally, it discusses the SIR model, distinguishes between different types of models, and includes practical applications such as temperature decay in a murder investigation.

Uploaded by

R228351Q
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
17 views3 pages

Mathematical Modeling Worksheet

The document is a worksheet for a mathematical modeling course at the University of Zimbabwe, covering various topics such as differential equations, logistic growth models, and systems of equations. It includes exercises on sketching phase lines, explaining mathematical modeling, finding fixed points, and analyzing stability. Additionally, it discusses the SIR model, distinguishes between different types of models, and includes practical applications such as temperature decay in a murder investigation.

Uploaded by

R228351Q
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

UNIVERSITY OF ZIMBABWE

DEPARTMENT OF MATHEMATICS AND COMPUTATIONAL SCIENCES

HMTHCS324:Mathematical modeling Worksheet

dP
1. Suppose that you wish to model with a differential equation of the form= f (P ) where
dt
P (t) is the population at time t. Experiments have been performed on the population
that give the following information:
(i) The population at P = 0 remains constant,
(ii) A population close to zero will decrease,
(iii) A population of P = 20 will increase,
(v) A population of P > 100 will decrease.

(a) Sketch the simplest phase line that agrees with the experimental information,
(b) Graph a rough sketch of the function f (P ) for the phase line in part (a).

2. (a) Briefly explain what you understand by the term mathematical modeling,
(b) State any three advantages of using mathematics language in modeling,
(c) Outline the steps followed when one is engaging in mathematical modeling.

3. Consider a logistic growth model described by ẋ = ( ab )x(b − x) where a > 0 and b > 0.

(a) Define what a and b represents in the model,


(b) Find the fixed points of the model,
(c) Determine the stability of each fixed point obtained above.

4. A model is given by the following equation ẋ = 2x(x − 1)2 (x − 2)3 . Find the fixed points
of the model

5. A mathematical model is given by the following system of equation


dS
dt = aR − bSI,
dI
dt = bSI − cI,
dR
dt = cI − aR.

(a) Write the state variables of the model,


(b) State the parameters of the model,
(c) Draw a schematic diagram to represent the model.

6. (a) A certain company has 1000 employees. On a particular day a rumor began to spread
among the employees of the company that a senior
supervisor in a department had suddenly resigned. Assuming that the rate of the
spread of the rumor is proportional to the number of possible encounters between
employees who have heard the rumor and those who have not.
i. Find an expression for the possible meetings between employees who have heard
the rumor and those who have not,

1
ii. Write the differential equation that describe the model and solve it.
iii. Given that when t = 0, the number of employees is 50 and when t = 1 the
number of employees who had heard the rumor was 100, Find the time when
800 people have heard the rumor.
(b) When one is engaged in mathematical modeling, they undergo similar steps and pro-
cesses as for conducting a biological experiment. Outline the all the stages involved.
(c) i. A model is described by the following system of equations

x˙1 = −3x1 + 2x2 ,


x˙2 = x1 − 2x2
.
A. Write the Jacobian matrix,
B. Evaluate the eigenvalues and eigenvectors,
C. Draw the phase diagram for the model.
ii. Given the following nonlinear system

ẋ = x(1 − x) − xz,
ẏ = ry(1 − y) − yz,
ż = z(2x + 2y − 1).
.
A. Find the Jacobian matrix,
B. Show that the equilibrium in which all populations are zero is always
unstable,
C. Show that the equilibrium points (1, 0, 0) and (0, 1, 0) are unstable.
(d) Consider a model represented by the following system

x˙1 = x1 (3 − x1 − x2 ),
x˙2 = x2 (x1 − 1).
i. Find the fixed points of the system,
ii. Evaluate the Jacobian matrix for the system,
iii. At each equilibrium point find the eigenvalue and the associated eigenvetors and
classify the point.
i. Consider the following SIR model:
dS
dt = −ρβSI,
dI
dt = ρβSI − κI,
dR
dt = κI.

where β > 0, κ > 0 and 0 < ρ < 1


A. State the biological meaning of all the state variables and the parameters in
the model,
B. Explain why the change in S proportional to both S and I,
C. Explain the reason why the terms ρβSI and κI appear in two different
equations and why are they positive in one instance and negative in the
other

2
ii. Consider the following nonlinear system:

ẋ = x(3 − x − 2y),
ẏ = y(4 − 3x − y).

A. Find the stationary points of the system,


B. Compute the linearization of the system around the nontrivial stationary
point,
C. Solve the linear system from (b) and sketch its phase portrait.
(e) Distinguish between:
(a) empirical vs mechanistic models,
(b) deterministic vs Stochastic models,
(c) discrete vs continuous models.
2
(f) The equation of motion of a linear spring with no damping is m ddt2x + kx = 0. What
are the units of the spring constant.
(g) Consider a logistic growth model described by ẋ = ( ab )x(b − x) where a > 0 and
b > 0.
i. Define what a and b represents in the model,
ii. Find the fixed points of the model,
iii. Determine the stability of each fixed point obtained above.
(h) i. A detective discovers a murder victim in a hotel room at 9 : 00am one morning.
The temperature of the body is 800 F . One hour later at 10 : 00am the body
has cooled to 750 F . The room is kept at a temperature of 700 F . Assume that
the victim had a normal temperature of 98.60 F at the time of death. Find the
time when the murder took place. [10]
ii. A model is described by the following system of equations

x˙1 = 4x1 − x2 ,
x˙2 = −x1 + 4x2
.
A. Write the Jacobian matrix,
B. Evaluate the eigenvalues and eigenvectors,
C. Draw the phase diagram for the model.
iii. A model is given by the following equation ẋ = 2x(x − 1)2 (x − 2)3 .
A. Find the fixed points of the model,
B. Determine the stability of each of the fixed points obtained.
(i) Consider a model represented by the following system

x˙1 = x1 (3 − x1 − x2 ),
x˙2 = x2 (x1 − 1).
i. Find the fixed points of the system,
ii. Evaluate the Jacobian matrix for the system,
iii. At each equilibrium point find the eigenvalue and the associated eigenvetors and
classify the point.

Common questions

Powered by AI

Mathematical modeling involves the following steps: defining the problem, formulating a model, analyzing the model, testing the model, and validating it against real-world data before using it for predictions or further applications. These steps closely resemble the processes in conducting a biological experiment, where one starts with a hypothesis, conducts experiments, analyzes results, compares them with predictions, and refines the hypothesis based on the discrepancies observed .

The Jacobian matrix plays a crucial role in analyzing a system of differential equations by providing a linear approximation of the system near equilibrium points. It contains partial derivatives of the system's functions with respect to its variables, offering a way to evaluate how small changes in state variables affect the system behavior. Eigenvalues of the Jacobian reveal stability characteristics: if all eigenvalues have negative real parts, the equilibrium is stable; if any have positive real parts, it is unstable .

In the SIR model, the change in the susceptible population S is proportional to the product of S and I because the rate of new infections depends on the frequency of contacts between susceptible and infected individuals. Each contact has a probability of transmission, hence the rate of change in susceptible individuals is directly proportional to both the number of susceptible individuals who can potentially be infected and the number of infective individuals who can spread the disease .

A logistic growth model often offers more realistic population dynamics modeling compared to exponential growth because it incorporates a carrying capacity, reflecting resource limitations. While exponential growth assumes unlimited resources leading to perpetual increase, logistic growth approaches a saturation point where growth slows as the population nears the environmental carrying capacity. This makes the logistic model more suitable for populations in constrained environments and highlights how intrinsic biological or environmental factors curb growth .

Continuous models, which use differential equations, are suited for representing processes that evolve smoothly over time, such as population growth or chemical reactions, and can model changes at any infinitesimal interval. Discrete models, however, are appropriate for systems where changes occur at distinct intervals, such as digital systems or population models with clear generations. The choice between them affects computational approaches, with continuous models requiring calculus-based solutions and discrete models often analyzed via difference equations, highlighting different trade-offs between precision and ease of computation in different scenarios .

The stability of fixed points in a logistic growth model can be determined by evaluating the derivative of the function at the fixed points. In the model ˙x = (a/b)x(b − x), the fixed points occur when x = 0 and x = b. The stability is determined by the eigenvalue test: if the derivative at a fixed point is negative, the fixed point is stable (attracting), and if positive, it is unstable (repelling). For x = 0, the derivative is positive, indicating instability. For x = b, the derivative is negative, indicating stability .

To sketch the phase diagram of a 2-variable linear system of differential equations such as ˙x1 = 4x1 − x2 and ˙x2 = −x1 + 4x2, plot the direction of the vector field representing system dynamics in a plane. Find eigenvalues and eigenvectors of the system's Jacobian matrix to identify eigen-directions. The trajectory's orientation and type (nodes, saddle points) depend on eigenvalue signs and directions; positives indicate expansion, negatives indicate contraction. Drawing several trajectories that highlight these behaviors will visualize the system's dynamic behaviour .

In the SIR (Susceptible-Infective-Recovered) model, the term ρβSI represents the rate of change due to infections. In the equation for susceptible individuals, it appears as negative because it represents the decrease in susceptible individuals as they become infected. Conversely, in the equation for infective individuals, it is positive since it corresponds to the increase in the number of infective individuals due to new infections. This capture of interacting population dynamics explains why the same term can have opposing signs in different equations .

To determine the fixed points of the equation ˙x = 2x(x −1)²(x −2)³, set the right-hand side to zero: 2x(x −1)²(x −2)³ = 0. The solutions are x = 0, x = 1, and x = 2. By analyzing the derivative of the equation at these points or by examining the power of terms involved (even powers indicate stability due to the function's local structure), the stability can be deduced. Generally, the fixed point at x = 0 is unstable, x = 1 is stable, and x = 2 is unstable .

Deterministic models produce a specific outcome from a given set of initial conditions without randomness, offering precise predictions. In contrast, stochastic models incorporate random variables to account for uncertainties and variations in predictions. This difference significantly impacts their usage; deterministic models are suitable for systems with predictable interactions, while stochastic models are better for systems where uncertainty plays a critical role .

You might also like