4 Insect outbreak model
|
|
- Horatio Wiggins
- 5 years ago
- Views:
Transcription
1 4 Insect outbreak model In this lecture I will put to a good use all the mathematical machinery we discussed so far. Consider an insect population, which is subject to predation by birds. It is a very well known and very well documented fact that often it is possible to observe so-called outbreaks of the insect population, when the number of insects increase from almost undetectable to extremely high. Here is the question: Is it possible to suggest a minimalist mathematical model in the form of one parameterdependent ODE, which would exhibit this phenomenon? Assume that the dynamics of the insect population is governed by the ODE of the form Ṅ = F N) P N), where Nt) is the population size at time t, F is the growth rate of the population in the absence of predators, and P describes the effect of the birds on the rate of change of the population. I am free to choose any reasonable functional form for F and P. Since I am looking for a model as simple as possible, let me take the logistic equation for F : ) F N) = r B N 1 NKB, where r B > 0 is the per capita growth rate when N is close to zero, and K B is the population carrying capacity. I know that if there are no predators, then Nt) K B as t for any initial condition N 0 > 0. To choose P, I note that it is reasonable to expect that if N then P N) B = const since the birds can eat only some amount of insect. It is usually said that function P in this case should exhibit a saturation effect. Another qualitative feature of P is that it should decrease to zero faster than the linear function. The rational for this is that if the insect population is low, then birds prefer looking for food somewhere else. Putting these two assumptions together, I obtain that the graph of P should look similar to the one shown in the figure. PN) Figure 1: Qualitative form of P N) in the insect outbreak model Math 484/684: Mathematical modeling of biological processes by Artem Novozhilov artem.novozhilov@ndsu.edu. Fall 2015 N 1
2 It is possible to pick various analytical expressions for P, but one of the simplest with the desired properties is possibly P N) = BN 2 A 2 + N 2, where B > 0 is the saturation constant and A > 0 determines the inflection point of the graph of P N). Putting everything together, ) Ṅ = r B N 1 NKB BN 2 A 2 + N 2. 1) Model 1) has four parameters. It is usually quite a difficult task to fully describe the behavior of solutions depending on more than two parameters. Luckily for our case, it is frequently possible to reduce the number of model parameters by properly rescaling the variables. Let me make the change of variables and introduce new parameters x = N A, τ = Bt A, r = Ar B B, q = K B A. Note that new variables and parameters are dimensionless I will explain how to actually figure out a possible change of the variables at the end of the lecture). In the new variables I have ẋ = rx 1 x ) x2 =: fx, r, q). 2) q 1 + x2 To find equilibria of 2) I need to solve fx, r, q) = 0. Obviously, ˆx = 0 is always an equilibrium, which is unstable for any parameter values r, q actually, f 0, r, q) = r > 0). Other possible equilibria, if they exist, should satisfy r 1 x ) = x q 1 + x 2. Finding the zeros of the last expression boils down to solving a cubic equation, which, albeit formally possible, usually leads to messy expressions, whose use is dubious at best. It is better to look for the equilibria graphically. The equilibria are the points of intersections of g 1 x) = r1 x/q) and g 2 x) = u/1 + x 2 ). Since g 2 x) does not have any parameters, I can plot it once and for all. g 1 x) defines a straight line, which intersect x-axis at the point q and y-axis at the point r. By changing r for a fixed q I can see that it is possible to have additionally from one to three equilibria see the figure), which I will denote ˆx 1, ˆx 2, ˆx 3. How to determine stability of these equilibria? I can use the standard technique by analyzing the sign of the derivative f evaluated at the equilibria, but an easier method is to note that if all the equilibria are hyperbolic, then any two equilibria without other equilibria between them must be such that one is unstable and another one is asymptotically stable. I know that ˆx = 0 is unstable. Therefore, if there is only one additional equilibrium, say ˆx 1 and it is hyperbolic i.e., f ˆx 1 ) 0), then it must be asymptotically stable. If there are all three equilibria then ˆx 1 and ˆx 3 are asymptotically stable and ˆx 2, which is between them, is unstable see the figure). In the same figure I can also see that the bifurcations in the system happen when either two equilibria collide and disappear, or when two equilibria appear as the parameters change. This is 2
3 g1x), g2x) r 1 r 2 g 1 x) = r 1 x ) q r 3 g 2 x) = x 1+x 2 x q Figure 2: Nontrivial equilibria of 2) given as the points of intersection of g 1 x) and g 2 x). It is possible to have one r 1 and r 3 cases) or three r 2 case) equilibria. For some parameter values the graph of g 1 x) touches the graph of g 2 x) and there are two equilibria, one of which is non-hyperbolic exactly the fold bifurcation. It is interesting to find the bifurcation values of r, q. They can be found as the solutions to the following system fx, r, q) = rx 1 x ) x2 q 1 + x 2 = 0, df x, r, q) = r 1 2x ) 2x dx q 1 + x 2 ) 2 = 0. The first equation here is the condition for x to be an equilibrium, and the second equation is the condition for this equilibrium to be non-hyperbolic. By solving this system with respect to r and q I get 2x 3 r = 1 + x 2 ) 2, q = 2x3 x 2 1. fx) 0 ˆx 1 ˆx 2 ˆx 3 x Figure 3: The case of bistability in model 2). asymptotically stable Equilibria 0, ˆx 2 are unstable, and ˆx 1 and ˆx 3 are 3
4 The last expressions is a parametrically defined curve on the plane q, r). The shape of this curve can be analyzed by the usual means I used Mathematica to plot it, see the figure). Note that there is a singularity on this curve, which is called the cusp singularity. In the area inside the curve the model 2) has three nontrivial equilibria, and outside the curve there is only one equilibrium different from zero. Everywhere on the curve, except for the cusp point, crossing of this curve corresponds to the fold bifurcation, when either two equilibria appear or disappear. D 1 2 C r B q A Figure 4: The bifurcation curve in the space of parameters r, q) of model 2) Model 2) and, correspondingly, 1)) exhibits the phenomenon of hysteresis. Hysteresis can be defined as the dependence of the system state not only on the current environment but also on the previous history how the system found itself in this particular state. To illustrate this phenomenon with the help of the model 2), consider the surface in the space r, q, x), which is defined by the equilibrium condition fx, r, q) = 0. This surface is shown in the figure below, here for each fixed r, q) the x-coordinate on the surface gives the equilibria of 2). The previous figure can be obtained by the projection of this surface onto r, q) plane. Now I assume that parameter q is fixed, but somehow I can change the value of parameter r in the model. I start at the point A see the figures), where there is only one nontrivial equilibrium ˆx 1 this corresponds to the case in Figure 3 for r 1 ). This point is asymptotically stable, and the insect population settles at this point. Now I assume that r starts growing such that I reach point B in Figure 4. For an observer not much changed, since the population is still at the equilibrium ˆx 1, which is slightly greater than before, but quite close. On the other hand, for the whole system actually I have that the fold bifurcation happened, and aside of ˆx 1 two new equilibria ˆx 2 and ˆx 3 appeared, the former is unstable, and the latter one is asymptotically stable. If I increase r to the point C in Figure 4, I have observed continuous growth of ˆx 1 -equilibrium, in which the population resided. At this point actually another bifurcation happens, points ˆx 1 and ˆx 2 collide and disappear. What is going to happen with the population? It must find another stable state, which is actually present in the system and given by ˆx 3. Note however, that the coordinate of ˆx 3 is much larger than that of ˆx 1, and at this point the system experiences the outbreak! Now the insect population at ˆx 3 -equilibrium for any larger values 4
5 q C D B C ˆx A B r Figure 5: This surface is defined by the condition fx, r, q) = 0 in model 2) of r point D). The path for the population state in Figure 5 is given by A, lower B, lower C, upper C, D. Now I assume that somehow I am able to reduce the value of r, the insect population will start declining, but staying still at the point ˆx 3, now I will need to reduce the parameter value to the point B to actually make sure that the insect population is again at ˆx 1. Hence, the path on Figure 5 when we reduce the values of r is given by D, upper C, upped B, lower B, A. 4.1 Non-dimensional variables In the mathematical model of the insect outbreak a change of variables was presented, which allowed me to reduce the number of parameters. Here is the explanation how to actually find such a change of variables. In physics this process is frequently referred to as non-dimensionalization. As an example consider the logistic equation in the form Ṅ = rn 1 N ), K where r > 0 and K > 0 are parameters of the problem. The former is the rate of growth and the latter is the carrying capacity. Nt) is the size of population at time moment t. Let me assume that here I am talking about elephants, hence Nt) is the number of elephants in the population at the time moment t. What are the dimensions of the parameters in this problem? First note that on the left hand side I have the derivative Ṅ, i.e., the rate of change, which means that the units on the left 5
6 [ ] have to be [Ṅ] = elephants time the square parenthesis mean dimensions ), which implies that the right hand side has to have the same units, and I have no choice other than consider that [K] = [elephants] and [r] = [ 1 time]. Now I make the substitutions Nt) = Anτ), t = T τ, where τ is a new independent variable, nτ) is a new dependent variable, and A and T are constants to be determined. Using the chain rule hence A T dnt) dt = A dnτ) dt dn dτ = ran 1 An ) K Since I am free to choose A and T, I can set = A dnτ) dτ dτ dt = A T dnτ) dτ = dn dτ = rt N 1 An ). K T = 1, [T ] = [time], A = K, [A] = [elephants], r and the change of variables nτ) = Nt) K, τ = rt, reduces the logistic equation to the form without any parameters: ṅ = n1 n), where now ṅ means the derivative with respect to τ. Note that the new variables are dimensionless: nτ) is the proportion of the population at time τ, and time is measured in some abstract nondimensional units. It is usually possible to suggest several alternative changes of the variables that lead to dimensionless form. Which one to choose is up to the researcher. For example, in the insect outbreak model the choice was made so that all the remaining parameters are in the linear function, whose behavior is easy to figure out. Here is another example for the system of two autonomous equations: Ṅ = an bnp, P = dp + cnp. Here N, P are dependent variables, that depend on the independent variable t time), and a, b, c, d are parameters, all of which are supposed to be nonnegative. Let Nt) = Axτ), P t) = Byτ), t = T τ be my change of variables such that A, B, T are some constants I would like to determine. Using again the chain rule, I find or, after canceling, Aẋ = aat x babt xy, Bẏ = dt By + cabt xy, ẋ = at x bbt xy, ẏ = dt y + cat xy. 6,
7 Putting at = 1, bbt = 1, cat = 1, I find that the change of variables and parameters leads to xτ) = cnt) a, yτ) = bp t) a, τ = at, α = d a, ẋ = x xy, ẏ = αy + xy, which has only one dimensionless parameter. You should convince yourself that the change yields xτ) = cnt) d, yτ) = bp t) a, τ = at, γ = d a, ẋ = x xy, ẏ = γyx 1), which is in some cases more convenient for subsequent analysis. 7
then the substitution z = ax + by + c reduces this equation to the separable one.
7 Substitutions II Some of the topics in this lecture are optional and will not be tested at the exams. However, for a curious student it should be useful to learn a few extra things about ordinary differential
More informationHomework 2. Due Friday, July We studied the logistic equation in class as a model of population growth. It is given by dn dt = rn 1 N
Problem 1 (10 points) Homework Due Friday, July 7 017 We studied the logistic equation in class as a model of population growth. It is given by dn dt = rn 1 N, (1) K with N(0) = N 0. (a) Make the substitutions
More informationONE DIMENSIONAL FLOWS. Lecture 3: Bifurcations
ONE DIMENSIONAL FLOWS Lecture 3: Bifurcations 3. Bifurcations Here we show that, although the dynamics of one-dimensional systems is very limited [all solutions either settle down to a steady equilibrium
More information8 Autonomous first order ODE
8 Autonomous first order ODE There are different ways to approach differential equations. Prior to this lecture we mostly dealt with analytical methods, i.e., with methods that require a formula as a final
More information14 Periodic phenomena in nature and limit cycles
14 Periodic phenomena in nature and limit cycles 14.1 Periodic phenomena in nature As it was discussed while I talked about the Lotka Volterra model, a great deal of natural phenomena show periodic behavior.
More informationLecture 3. Dynamical Systems in Continuous Time
Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver Yue-Xian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a
More information1 Mathematical models of the population growth
1 Mathematical models of the population growth 1.1 How many people can the Earth support? (Or a quick refresher on the least square method) I would like to start with a very simple and yet interesting
More informationHysteresis. Lab 6. Recall that any ordinary differential equation can be written as a first order system of ODEs,
Lab 6 Hysteresis Recall that any ordinary differential equation can be written as a first order system of ODEs, ẋ = F (x), ẋ := d x(t). (6.1) dt Many interesting applications and physical phenomena can
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationContinuous-Time Dynamical Systems: Sketching solution curves and bifurcations
Math 360 Winter 07 Section 0 HW#4 Continuous-Time Dynamical Systems: Sketching solution curves and bifurcations (Due in one week! Thursday Nov 9, 07) Materials contained in this assignment are of great
More information1.2. Introduction to Modeling
G. NAGY ODE August 30, 2018 1 Section Objective(s): Population Models Unlimited Resources Limited Resources Interacting Species 1.2. Introduction to Modeling 1.2.1. Population Model with Unlimited Resources.
More informationLogistic Model with Harvesting
Logistic Model with Harvesting Model Here, we will talk about a model for the growth and harvesting of a fish population. Suppose that, absent any fishermen, the population grows logistically. That is,
More information1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy
.. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution
More informationMA 137 Calculus 1 with Life Science Application A First Look at Differential Equations (Section 4.1.2)
MA 137 Calculus 1 with Life Science Application A First Look at Differential Equations (Section 4.1.2) Alberto Corso alberto.corso@uky.edu Department of Mathematics University of Kentucky October 12, 2015
More information3 Single species models. Reading: Otto & Day (2007) section 4.2.1
3 Single species models 3.1 Exponential growth Reading: Otto & Day (2007) section 4.2.1 We can solve equation 17 to find the population at time t given a starting population N(0) = N 0 as follows. N(t)
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More information2 One-dimensional models in discrete time
2 One-dimensional models in discrete time So far, we have assumed that demographic events happen continuously over time and can thus be written as rates. For many biological species with overlapping generations
More informationCHAPTER 2 POLYNOMIALS KEY POINTS
CHAPTER POLYNOMIALS KEY POINTS 1. Polynomials of degrees 1, and 3 are called linear, quadratic and cubic polynomials respectively.. A quadratic polynomial in x with real coefficient is of the form a x
More informationSpring /30/2013
MA 138 - Calculus 2 for the Life Sciences FINAL EXAM Spring 2013 4/30/2013 Name: Sect. #: Answer all of the following questions. Use the backs of the question papers for scratch paper. No books or notes
More informationContinuous Threshold Policy Harvesting in Predator-Prey Models
Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July
More informationKey ideas about dynamic models
GSTDMB 202: DYNAMICAL MODELLING FOR BIOLOGY AND MEDICINE Lecture 2 Introduction to modelling with differential equations Markus Owen ey ideas about dynamic models Dynamic models describe the change in
More informationObjective. Single population growth models
Objective Single population growth models We are given a table with the population census at different time intervals between a date a and a date b, and want to get an expression for the population. This
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More information1 Simple one-dimensional dynamical systems birth/death and migration processes, logistic
NOTES ON A Short Course and Introduction to Dynamical Systems in Biomathematics Urszula Foryś Institute of Applied Math. & Mech. Dept. of Math., Inf. & Mech. Warsaw University 1 Simple one-dimensional
More informationThe functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z).
Review of Multi-variable calculus: The functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z). The spatial variable represents the environment where the species
More informationUCI MATH CIRCLE: POPULATION DYNAMICS
UCI MATH CIRCLE: POPULATION DYNAMICS ADAPTED FROM A SESSION OF THE UNIVERSITY OF UTAH MATH CIRCLE When Europeans first arrived in the Americas, Passenger Pigeons were perhaps the most populous bird in
More informationUse separation of variables to solve the following differential equations with given initial conditions. y 1 1 y ). y(y 1) = 1
Chapter 11 Differential Equations 11.1 Use separation of variables to solve the following differential equations with given initial conditions. (a) = 2ty, y(0) = 10 (b) = y(1 y), y(0) = 0.5, (Hint: 1 y(y
More informationMA 138: Calculus II for the Life Sciences
MA 138: Calculus II for the Life Sciences David Murrugarra Department of Mathematics, University of Kentucky. Spring 2016 David Murrugarra (University of Kentucky) MA 138: Section 11.4.2 Spring 2016 1
More informationQualitative analysis of differential equations: Part I
Qualitative analysis of differential equations: Part I Math 12 Section 16 November 7, 216 Hi, I m Kelly. Cole is away. Office hours are cancelled. Cole is available by email: zmurchok@math.ubc.ca. Today...
More informationMA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, Exam Scores. Question Score Total
MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, 2016 Exam Scores Question Score Total 1 10 Name: Section: Last 4 digits of student ID #: No books or notes may be used. Turn off all
More informationAutonomous Equations and Stability Sections
A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Autonomous Equations and Stability Sections 2.4-2.5 Dr. John Ehrke Department of Mathematics Fall 2012 Autonomous Differential
More informationIntroduction to Biomathematics. Problem sheet
Introction to Biomathematics Problem sheet 1. A model for population growth is given in non-dimensional units in the form = u1 u ) u0) > 0. a) Sketch the graph of the function fu) = u1 u ) against u for
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationMath 3B: Lecture 14. Noah White. February 13, 2016
Math 3B: Lecture 14 Noah White February 13, 2016 Last time Accumulated change problems Last time Accumulated change problems Adding up a value if it is changing over time Last time Accumulated change problems
More informationMathematical models in life sciences. Lecture notes of a course held by Marco Di Francesco based on the book by J. D. Murray
Mathematical models in life sciences Lecture notes of a course held by Marco Di Francesco based on the book by J. D. Murray December 31, 2010 Contents I Continuous time dependent models 3 1 Ordinary differential
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More information5 Alfred Lotka, Vito Volterra, and Population Cycles
5 Alfred Lotka, Vito Volterra, and opulation Cycles Dr. Umberto D Ancona entertained me several times with statistics that he was compiling about fishing during the period of the war and in periods previous
More informationMAT 311 Midterm #1 Show your work! 1. The existence and uniqueness theorem says that, given a point (x 0, y 0 ) the ODE. y = (1 x 2 y 2 ) 1/3
MAT 3 Midterm # Show your work!. The existence and uniqueness theorem says that, given a point (x 0, y 0 ) the ODE y = ( x 2 y 2 ) /3 has a unique (local) solution with initial condition y(x 0 ) = y 0
More informationLecture 9. Systems of Two First Order Linear ODEs
Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form
More information1. The growth of a cancerous tumor can be modeled by the Gompertz Law: dn. = an ln ( )
1. The growth of a cancerous tumor can be modeled by the Gompertz Law: ( ) dn b = an ln, (1) dt N where N measures the size of the tumor. (a) Interpret the parameters a and b (both non-negative) biologically.
More information8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds
c Dr Igor Zelenko, Spring 2017 1 8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds to section 2.5) 1.
More informationMA 777: Topics in Mathematical Biology
MA 777: Topics in Mathematical Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.math.uky.edu/~dmu228/ma777/ Spring 2018 David Murrugarra (University of Kentucky) Lecture
More informationIntroduction to Dynamical Systems
Introduction to Dynamical Systems Autonomous Planar Systems Vector form of a Dynamical System Trajectories Trajectories Don t Cross Equilibria Population Biology Rabbit-Fox System Trout System Trout System
More informationThe Big, Big Picture (Bifurcations II)
The Big, Big Picture (Bifurcations II) Reading for this lecture: NDAC, Chapter 8 and Sec. 10.0-10.4. 1 Beyond fixed points: Bifurcation: Qualitative change in behavior as a control parameter is (slowly)
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationDynamical Systems: Lecture 1 Naima Hammoud
Dynamical Systems: Lecture 1 Naima Hammoud Feb 21, 2017 What is dynamics? Dynamics is the study of systems that evolve in time What is dynamics? Dynamics is the study of systems that evolve in time a system
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationMath 2930 Worksheet Introduction to Differential Equations. What is a Differential Equation and what are Solutions?
Math 2930 Worksheet Introduction to Differential Equations Week 1 January 25, 2019 What is a Differential Equation and what are Solutions? A differential equation is an equation that relates an unknown
More informationChaos and Liapunov exponents
PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic
More informationAgenda Sections 2.4, 2.5
Agenda Sections 2.4, 2.5 Reminders Read 3.1, 3.2 Do problems for 2.4, 2.5 Homework 1 due Friday Midterm Exam I on 1/23 Lab on Friday (Shapiro 2054) Office hours Tues, Thurs 3-4:30 pm (5852 East Hall) Theorem:
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationAndronov Hopf and Bautin bifurcation in a tritrophic food chain model with Holling functional response types IV and II
Electronic Journal of Qualitative Theory of Differential Equations 018 No 78 1 7; https://doiorg/10143/ejqtde018178 wwwmathu-szegedhu/ejqtde/ Andronov Hopf and Bautin bifurcation in a tritrophic food chain
More informationA Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More information1 (t + 4)(t 1) dt. Solution: The denominator of the integrand is already factored with the factors being distinct, so 1 (t + 4)(t 1) = A
Calculus Topic: Integration of Rational Functions Section 8. # 0: Evaluate the integral (t + )(t ) Solution: The denominator of the integrand is already factored with the factors being distinct, so (t
More informationf (x) = x 2 Chapter 2 Polynomial Functions Section 4 Polynomial and Rational Functions Shapes of Polynomials Graphs of Polynomials the form n
Chapter 2 Functions and Graphs Section 4 Polynomial and Rational Functions Polynomial Functions A polynomial function is a function that can be written in the form a n n 1 n x + an 1x + + a1x + a0 for
More informationAn Application of Perturbation Methods in Evolutionary Ecology
Dynamics at the Horsetooth Volume 2A, 2010. Focused Issue: Asymptotics and Perturbations An Application of Perturbation Methods in Evolutionary Ecology Department of Mathematics Colorado State University
More informationSample Questions, Exam 1 Math 244 Spring 2007
Sample Questions, Exam Math 244 Spring 2007 Remember, on the exam you may use a calculator, but NOT one that can perform symbolic manipulation (remembering derivative and integral formulas are a part of
More informationHomework 2 Modeling complex systems, Stability analysis, Discrete-time dynamical systems, Deterministic chaos
Homework 2 Modeling complex systems, Stability analysis, Discrete-time dynamical systems, Deterministic chaos (Max useful score: 100 - Available points: 125) 15-382: Collective Intelligence (Spring 2018)
More informationA Comparison of Two Predator-Prey Models with Holling s Type I Functional Response
A Comparison of Two Predator-Prey Models with Holling s Type I Functional Response ** Joint work with Mark Kot at the University of Washington ** Mathematical Biosciences 1 (8) 161-179 Presented by Gunog
More informationNon-Dimensionalization of Differential Equations. September 9, 2014
Non-Dimensionalization of Differential Equations September 9, 2014 1 Non-dimensionalization The most elementary idea of dimensional analysis is this: in any equation, the dimensions of the right hand side
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationMath 266: Autonomous equation and population dynamics
Math 266: Autonomous equation and population namics Long Jin Purdue, Spring 2018 Autonomous equation An autonomous equation is a differential equation which only involves the unknown function y and its
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationDifferential Equations ( DEs) MA 137 Calculus 1 with Life Science Application A First Look at Differential Equations (Section 4.1.
/7 Differential Equations DEs) MA 37 Calculus with Life Science Application A First Look at Differential Equations Section 4..2) Alberto Corso alberto.corso@uky.edu Department of Mathematics University
More informationProblem Set Number 1, j/2.036j MIT (Fall 2014)
Problem Set Number 1, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept., Cambridge, MA 02139) September 9, 2014 Due Fri., September 19, 2014. 1 Get equation from phase line portrait problem
More information2 Discrete Dynamical Systems (DDS)
2 Discrete Dynamical Systems (DDS) 2.1 Basics A Discrete Dynamical System (DDS) models the change (or dynamics) of single or multiple populations or quantities in which the change occurs deterministically
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationIntroduction to bifurcations
Introduction to bifurcations Marc R. Roussel September 6, Introduction Most dynamical systems contain parameters in addition to variables. A general system of ordinary differential equations (ODEs) could
More information7.2. Differential Equations. Phase Plots. Example 1 The logistic equation. Phase Plots. Phase Plots. Phase Plots, Equilibria, and Stability
Differential Equations 7 7.2, Equilibria, and Stability Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. Phase plots provide a way to visualize the dynamics
More information2.2 Average vs. Instantaneous Description
2 KINEMATICS 2.2 Average vs. Instantaneous Description Name: 2.2 Average vs. Instantaneous Description 2.2.1 Average vs. Instantaneous Velocity In the previous activity, you figured out that you can calculate
More informationNonlinear Autonomous Dynamical systems of two dimensions. Part A
Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous
More informationAsymptotic and Perturbation Methods: Lecture Notes.
Asymptotic and Perturbation Methods: Lecture Notes. Blerta Shtylla University of Utah Mathematics Department Fall 2008 1 Lecture 1 and Lecture 2 Week of August 26, 2008 Introduction: For many problems
More informationLecture 6, September 1, 2017
Engineering Mathematics Fall 07 Lecture 6, September, 07 Escape Velocity Suppose we have a planet (or any large near to spherical heavenly body) of radius R and acceleration of gravity at the surface of
More informationFigure 5: Bifurcation diagram for equation 4 as a function of K. n(t) << 1 then substituting into f(n) we get (using Taylor s theorem)
Figure 5: Bifurcation diagram for equation 4 as a function of K n(t)
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationMath 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations
Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function
More informationMathematical Models of Biological Systems
Mathematical Models of Biological Systems CH924 November 25, 2012 Course outline Syllabus: Week 6: Stability of first-order autonomous differential equation, genetic switch/clock, introduction to enzyme
More informationMath 2930 Worksheet Equilibria and Stability
Math 2930 Worksheet Equilibria and Stabilit Week 3 September 7, 2017 Question 1. (a) Let C be the temperature (in Fahrenheit) of a cup of coffee that is cooling off to room temperature. Which of the following
More information22: Applications of Differential Calculus
22: Applications of Differential Calculus A: Time Rate of Change The most common use of calculus (the one that motivated our discussions of the previous chapter) are those that involve change in some quantity
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationDynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems
. Ordinary Differential Equations. An ordinary differential equation (ODE, for short) is an equation involves a dependent variable and its derivatives with respect to an independent variable. When we speak
More informationv n+1 = v T + (v 0 - v T )exp(-[n +1]/ N )
Notes on Dynamical Systems (continued) 2. Maps The surprisingly complicated behavior of the physical pendulum, and many other physical systems as well, can be more readily understood by examining their
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More informationIntroduction to Dynamical Systems Basic Concepts of Dynamics
Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More information1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos.
Dynamical behavior of a prey predator model with seasonally varying parameters Sunita Gakkhar, BrhamPal Singh, R K Naji Department of Mathematics I I T Roorkee,47667 INDIA Abstract : A dynamic model based
More informationtutorial ii: One-parameter bifurcation analysis of equilibria with matcont
tutorial ii: One-parameter bifurcation analysis of equilibria with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht February 13, 2018 1 This session
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More information(e) Use Newton s method to find the x coordinate that satisfies this equation, and your graph in part (b) to show that this is an inflection point.
Chapter 6 Review problems 6.1 A strange function Consider the function (x 2 ) x. (a) Show that this function can be expressed as f(x) = e x ln(x2). (b) Use the spreadsheet, and a fine subdivision of the
More informationCalculus. Weijiu Liu. Department of Mathematics University of Central Arkansas 201 Donaghey Avenue, Conway, AR 72035, USA
Calculus Weijiu Liu Department of Mathematics University of Central Arkansas 201 Donaghey Avenue, Conway, AR 72035, USA 1 Opening Welcome to your Calculus I class! My name is Weijiu Liu. I will guide you
More informationChapter 4: Growth & Decay. Introduction
Linear Difference Equations Non-Linear Models with Applications Chapter 4: Growth & Decay 70 / 117 Introduction Here, model systems with population so big (or growth so fine-grained) that continuity can
More informationSolutions to Math 53 First Exam April 20, 2010
Solutions to Math 53 First Exam April 0, 00. (5 points) Match the direction fields below with their differential equations. Also indicate which two equations do not have matches. No justification is necessary.
More informationAnalysis of Plankton Populations
Analysis of Plankton Populations Sean Li September 21, 25 Abstract I use numerical and analytical techniques to study the dynamics of various models of plankton food webs with and without resource fluctuations.
More informationPopulation Dynamics. Max Flöttmann and Jannis Uhlendorf. June 12, Max Flöttmann and Jannis Uhlendorf () Population Dynamics June 12, / 54
Population Dynamics Max Flöttmann and Jannis Uhlendorf June 12, 2007 Max Flöttmann and Jannis Uhlendorf () Population Dynamics June 12, 2007 1 / 54 1 Discrete Population Models Introduction Example: Fibonacci
More informationNon-Linear Dynamics Homework Solutions Week 1
Non-Linear Dynamics Homework Solutions Week 1 Chris Small March, 007 Please email me at smachr09@evergreen.edu with any corrections or concerns reguarding these solutions..3. We are given that ẋ = k 1
More information