Ordinary Differential Equations
|
|
- Hugo O’Connor’
- 5 years ago
- Views:
Transcription
1 Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005
2 Overview What are Ordinary Differential Equations (ODEs)? Equilibria: existence and properties. What techniques exist to solve ODEs? How do we draw up ODEs: A case study on Predator-Prey systems. A variant: Delay Differential Equations (DDEs)
3 Ordinary Differential Equations ODEs can be used to model the behaviour of a system through time The systems studied range from spring-mass systems to ecology and economics. ODEs have the general form: d x dt = f( x, t). (1) Quite frequently, we deal with time independent ODEs, which have the form The simplest form they can take are linear ODEs: d x dt = f( x). (2) d x dt = A x. (3) with A a matrix.
4 Equilibria Equilibria occur where Three types of equilibria occur: d x dt = f( x, t) = 0. (4) Stable: the system returns to the equilibrium after small perturbation. Unstable: the system diverges from the equilibrium state after small perturbation. Neutrally Stable: the system neither diverges from, nor converges towards the equilibrium after small perturbation. In the case of time independent ODEs we can determine the stability by local stability analysis quite easily. In the first step, we linearize the ODE around the equilibrium.
5 Equilibria Let the equilibrium point be x 0. The ODE given by d x dt = f( x) (5) can be approximated by Taylor expansion. We first introduce x = x x 0. We then approximate (5) by in which J is the Jacobian matrix: J = f 1 d x dt = J x (6) f 1 x n x f n x 1 f n x n x= x0 (7)
6 Equilibria The solution to such a set of linear differential equations is just x (t) = a i Ei e λ it, (8) with E i the eigenvectors of J, λ i the corresponding eigenvalues, and a i the amplitude of that particular eigenmode of the system. Clearly, if all the eigenvalues have a negative real part, the system returns to the equilibrium position and the system is stable. If one or more eigenvalues have positive real parts then it is unstable. If none have positive and one or more zero real parts it is neutrally stable.
7 Equilibria To determine eigenvalues we need to solve the equation f 1 f x 1 λ 1 f x 2 1 x n f 2 f 2 f det x 1 x 2 λ 2 x n = 0 (9) f n f x 2 n x n λ x= x0 f n x 1 This yields an n th order polynomial c 0 + c 1 λ + c 2 λ c n λ n = 0 (10) the roots of which are the eigenvalues. Note that if all c i are real and have the same sign, the system is stable.
8 Solving ODEs Usually, we are concerned with systems which are too large for analytical treatment, but can be handled numerically. Most often we are concerned with initial-value problems: given the state of the system at time t 0, compute the state at a number of points in time t n > t 0. The simplest method is the Euler method x(t n+1 ) = x(t n ) + h f( x(t n ), t n ) + O(h 2 ), (11) in which h is the time step and O(h 2 ) indicates an error which is proportional to h 2. An improvement is the midpoint method: k1 = h f( x(t n ), t) (12) x(t n+1 ) = x(t n ) + h f( x(t n ) + k 1 /2, t n + h/2) + O(h 3 ) (13)
9 Solving ODEs Better still is the 4 th order Runge-Kutta method: k1 = h f( x(t n ), t) (14) k2 = h f( x(t n ) + k 1 /2, t n + h/2) (15) k3 = h f( x(t n ) + k 2 /2, t n + h/2) (16) k4 = h f( x(t n ) + k 3, t n + h) (17) x(t n+1 ) = x(t n ) + k k k k O(h5 ) (18) Various strategies for computing an optimal time step h exist. Problems occur in so-called stiff ODEs. In many cases (chemistry, biology) we need to check for zero crossings!
10 Predator-Prey Systems One of the best studied systems is the so-called predator-prey (or host-parasite) system. It is fundamental to understanding food webs in ecology. Its general form is in which dx 1 dt = F (x 1) G(x 1, x 2 ) (19) dx 2 dt = ηg(x 1, x 2 ) H(x 2 ) (20) F (x 1 ) is the growth rate of x 1, G(x 1, x 2 ) is the rate of predation of x 1 by x 2, η is the efficiency with which prey biomass is turned predator biomass, and H(x 2 ) is the starvation rate of x 2.
11 Predator-Prey Systems In its simplest form, we have the Lotka-Volterra system in which ( dx 1 dt = r 1 1 x ) 1 x 1 fx 1 x 2 (21) K 1 dx 2 dt = ηfx 1x 2 d 2 x 2 (22) f determines the rate at which predator and prey encounter eachother, r 1 is the maximum relative growth rate of x 1, K 1 is the carrying capacity of the ecosystem for x 1, and d 2 is the starvation rate of x 2. This system can be stable, neutrally stable, or unstable, depending on the parameters.
12 Predator-Prey Systems Other models for F include the Monod-model for bacteria S F (x 1, S) = µ max K S + S x 1 (23) in which S is the substrate (food) concentration, µ max the maximum growth rate, and K S a saturation constant. Other models for G include: Holling Type II Holling Type III Jost f f x 1 K 2 +x 1 x 2 x 2 1 x K 2 +x x 2 1 x K 2 +k 2 x 1 +x Ivlev f(1 e kx 1)x 2 f same as Monod for vertebrates similar to above All these models try to limit maximum growth in some way.
13 A Bacterial Predator-Prey System Suppose we have a bacterium X 1 using resource X 0, and which is preyed upon by predator Y. We assume that these species are present in a chemostat-like environment, in which all substances are well mixed, and which has a dilution rate D. Food (X 0 ) enters at a concentration S. Uptake of food by X 1 is modelled using the Monod-form. Uptake of X 1 by Y is modelled using the Holling-type-II form. Starvation of Y takes the usual form, with d y the starvation rate.
14 A Bacterial Predator-Prey System The differential equations for this model are dx 0 dt = D(S X 0 ) V 1 X 0 K 1 + X 0 X 1 (24) dx 1 X 0 X 1 = µ 1 X 1 V y Y DX 1 dt K 1 + X 0 K X + X 1 (25) dy dt = µ X 1 y Y (D + d y )Y K X + X 1 (26) with V x and µ x maximum uptake and growth rates respectively, and K x saturation constants.
15 A Bacterial Predator-Prey System This system can be in four phases Phase 0: Insufficient food is available to allow survival of X 1 ; stable equilibrium: X 0 = S, X 1 = 0, and Y = 0. Phase I: Sufficient food for survival of X 1, but insufficient number survive to allow survival of Y. Phase II: Both X 1 and Y coexist stably. Phase III: unstable coexistence of X 1 and Y. Important question: which parameters determine outcome?
16 A Bacterial Predator-Prey System The equilibrium with all species present is given by X 0 = 1 2 ( S K 1 V 1 K X µ y D d y ± ( S K1 V 1 K X µ y D d y ) 2 + 4K 1 S) (27) X 1 = Y = D + d y µ y D d y D D + d y µ y V y ( µ1 V 1 (S X 0 ) X 1 ) (28) (29) Only if all three solutions are positive do we have either a phase II or phase III system. The boundary between these phases II and III is determined by local stability analysis.
17 A Bacterial Predator-Prey System The dynamical behaviour of the predator prey system after introduction of predator Prey Predator Prey Predator Biomass (mg/l) Biomass (mg/l) time (h) time (h)
18 A Bacterial Predator-Prey System The boundaries of the phases I 0.18 I D (h 1 ) II III D (h 1 ) II III S (mg/l) S (mg/l)
19 Decoys in Predator-Prey Systems In a simple complication, we can add a third species X 2, which does not compete with X 1 but can collide with the predator. Assume that the predator can be in three states: free, bound to X 1 and bound to X 2. We can then draw up the following set of differential equations: d[x 1 Y ] dt d[x 2 Y ] dt dy free dt = k 1 [X 1 Y ] + rx 1 Y free (30) = k 2 [X 2 Y ] + rx 2 Y free (31) = (y x + 1)k 1 [X 1 Y ] + k 2 [X 2 Y ] r(x 1 + X 2 )Y free (32)
20 Decoys in Predator-Prey Systems We will assume that these reactions take place at a much faster rate than the predator prey dynamics We can then in quasi-steady state analysis say that k 1 [X 1 Y ] = rx 1 Y free and k 2 [X 2 Y ] = rx 2 Y free (33) Adding all versions of Y together yields dy dt = y x k 1 X 1 Y k 1 /r + X 1 + k 1 X 2 /k 2 = µ y X 1 Y K X + X 1 + K inh X 2 (34) Thus the presence of X 2 leads to an increase in the apparent saturation constant K X = K X + K inh X 2.
21 Decoys in Predator-Prey Systems The dynamical behaviour of the predator prey system after introduction of predator with decoys: K X = 2K X Prey Predator Prey Predator Biomass (mg/l) time (h) Biomass (mg/l) time (h)
22 Decoys in Predator-Prey Systems The boundaries of the phases with decoys: K X = 2K X I 0.18 I D (h 1 ) II D (h 1 ) II 0.06 III 0.06 III S (mg/l) S (mg/l)
23 Delay Differential Equations (DDEs) In ODEs the time evolution depends only on the current state of the system In DDEs the time evolution of the system depends on the current state, plus the state of the system at one or more points some time in the past. This type of equation is useful to model systems in which there are time delays between cause and effect. Examples are: incubation times of infections in models of epidemics, delay times between penetration of host by parasites and emergence of the new parasite generation delay between mating and birth, etc.
24 Delay Differential Equations (DDEs) DDEs have the general form: d x dt = f( x(t), x(t τ 1 ), x(t τ 2 ),..., x(t τ N )). (35) with τ 1, τ 2,..., τ N the set of delays in the system. The stability condition becomes f( x(t), x(t), x(t),..., x(t)) = 0. (36) Finding equilibria is easy, analysis of their properties generally done numerically. Numerical treatment requires different techniques than for ODEs. Initial value problems require information about the past for initialization!
25 Example DDE we use the bacterial predator-prey model as before, but use a delay equation to model the predator behaviour. The prey is converted to infected prey [X 1 Y ] as in the decoy model After a time delay τ the new generation of predators is released. The differential equations for this model are dx 0 dt = D(S X 0 ) V 1 X 0 K 1 + X 0 X 1 (37) dx 1 dt = µ 1 X 0 K 1 + X 0 X 1 rx 1 Y DX 1 (38) d[x 1 Y ] = e Dτ rx dt 1Y + rx 1 Y D[X 1 Y ] (39) dy dt = e Dτ (y x + 1)rX 1Y (D + d y )Y (40) in which X 1 and Y denote the concentrations of X 1 and Y at t τ.
26 Example DDE The delay effect can sometimes be approaximated by an average rate of formation of new prey, leading to a Holling type II version as in the original bacterial predator-prey model density (cm 3 ) density (cm 3 ) time (h) (a) time (h) (b) Part (a) shows the DDE appraoch, part (b) the ODE approximation.
Physics: 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 informationWorkshop on Theoretical Ecology and Global Change March 2009
2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More informationInteractions. Yuan Gao. Spring Applied Mathematics University of Washington
Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,
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 informationContinuous time population models
Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
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 informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
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 informationMATH 1014 N 3.0 W2015 APPLIED CALCULUS II - SECTION P. Perhaps the most important of all the applications of calculus is to differential equations.
MATH 1014 N 3.0 W2015 APPLIED CALCULUS II - SECTION P Stewart Chapter 9 Differential Equations Perhaps the most important of all the applications of calculus is to differential equations. 9.1 Modeling
More informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
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 informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp DOI: /m2an:
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp. 339 344 DOI: 10.1051/m2an:2003029 PERSISTENCE AND BIFURCATION ANALYSIS
More informationModeling Microbial Populations in the Chemostat
Modeling Microbial Populations in the Chemostat Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Modeling Microbial Populations in the Chemostat MBI, June 3, 204 / 34 Outline Why
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More informationFIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material
More information1 2 predators competing for 1 prey
1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2
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 information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More informationAN EXTENDED ROSENZWEIG-MACARTHUR MODEL OF A TRITROPHIC FOOD CHAIN. Nicole Rocco
AN EXTENDED ROSENZWEIG-MACARTHUR MODEL OF A TRITROPHIC FOOD CHAIN Nicole Rocco A Thesis Submitted to the University of North Carolina Wilmington in Partial Fulfillment of the Requirements for the Degree
More informationMath 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
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 informationNon-Linear Models. Non-Linear Models Cont d
Focus on more sophistiated interaction models between systems. These lead to non-linear, rather than linear, DEs; often not soluble exactly in analytical form so use Phase-Plane Analysis. This is a method
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 informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1
ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate
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 informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Transient Response of a Chemical Reactor Concentration of a substance in a chemical reactor
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
More informationThe Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor
, pp. 35-46 http://dx.doi.org/10.14257/ijbsbt.2017.9.3.04 The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor Alemu Geleta Wedajo
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationEECS 700: Exam # 1 Tuesday, October 21, 2014
EECS 700: Exam # 1 Tuesday, October 21, 2014 Print Name and Signature The rules for this exam are as follows: Write your name on the front page of the exam booklet. Initial each of the remaining pages
More informationROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 6, Number 1, Winter 1997 ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM J. CHATTOPADHYAY, E. BERETTA AND F. SOLIMANO ABSTRACT. The present paper deals with
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
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 informationNumerical bifurcation analysis of a tri-trophic food web with omnivory
Mathematical Biosciences 177&178 (2002) 201 228 www.elsevier.com/locate/mbs Numerical bifurcation analysis of a tri-trophic food web with omnivory B.W. Kooi *, L.D.J. Kuijper, M.P. Boer 1, S.A.L.M. Kooijman
More informationLab 5: Nonlinear Systems
Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting
More informationSolution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)
Math 380 Practice Final Solutions This is longer than the actual exam, which will be 8 to 0 questions (some might be multiple choice). You are allowed up to two sheets of notes (both sides) and a calculator,
More informationSystems Biology in Photosynthesis INTRODUCTION
1 / 26 Systems Biology in Photosynthesis INTRODUCTION Rainer Machné Institute for Theoretical Chemistry, University of Vienna, Austria PSI - Photon System Instruments, Czech Republic Brno, April, 2011
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationLesson 9: Predator-Prey and ode45
Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and
More informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
More informationDynamics Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses
Journal of Physics: Conference Series PAPER OPEN ACCESS Dynamics Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses To cite this article: D Savitri 2018
More informationLokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt
Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period
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 informationCalculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y
Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point
More informationIntroduction to standard and non-standard Numerical Methods
Introduction to standard and non-standard Numerical Methods Dr. Mountaga LAM AMS : African Mathematic School 2018 May 23, 2018 One-step methods Runge-Kutta Methods Nonstandard Finite Difference Scheme
More informationPreservation of local dynamics when applying central difference methods: application to SIR model
Journal of Difference Equations and Applications, Vol., No. 4, April 2007, 40 Preservation of local dynamics when applying central difference methods application to SIR model LIH-ING W. ROEGER* and ROGER
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 informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationMath 341 Fall 2006 Final Exam December 12th, Name:
Math 341 Fall 2006 Final Exam December 12th, 2006 Name: You may use a calculator, your note card and something to write with. You must attach your notecard to the exam when you turn it in. You cannot use
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationFunctional Response to Predators Holling type II, as a Function Refuge for Preys in Lotka-Volterra Model
Applied Mathematical Sciences, Vol. 9, 2015, no. 136, 6773-6781 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.53266 Functional Response to Predators Holling type II, as a Function Refuge
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 informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationAssignment 1 Logistic equation and Lotka Volterra two species competition model
Assignment 1 Logistic equation and Lotka Volterra two species competition model Edvin Listo Zec edvinli@student.chalmers.se 920625-2976 August 28, 2014 Co-workers: Emma Ekberg, Sofia Toivonen 1 Question
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 informationPREDATOR-PREY MODELS WITH DISTRIBUTED TIME DELAY
PREDATOR-PREY MODELS WITH DISTRIBUTED TIME DELAY PREDATOR-PREY MODELS WITH TIME DISTRIBUTED DELAY By ALEXANDRA TESLYA, M.SC. A Thesis Submitted to the School oif Graduate Studies in Partial Fulfillment
More informationWorksheet 8 Sample Solutions
Technische Universität München WS 2016/17 Lehrstuhl für Informatik V Scientific Computing Univ.-Prof. Dr. M. Bader 19.12.2016/21.12.2016 M.Sc. S. Seckler, M.Sc. D. Jarema Worksheet 8 Sample Solutions Ordinary
More informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationTHETA-LOGISTIC PREDATOR PREY
THETA-LOGISTIC PREDATOR PREY What are the assumptions of this model? 1.) Functional responses are non-linear. Functional response refers to a change in the rate of exploitation of prey by an individual
More informationNonstandard Finite Difference Methods For Predator-Prey Models With General Functional Response
Nonstandard Finite Difference Methods For Predator-Prey Models With General Functional Response Dobromir T. Dimitrov Hristo V. Kojouharov Technical Report 2007-0 http://www.uta.edu/math/preprint/ NONSTANDARD
More informationHarvesting Model for Fishery Resource with Reserve Area and Modified Effort Function
Malaya J. Mat. 4(2)(2016) 255 262 Harvesting Model for Fishery Resource with Reserve Area and Modified Effort Function Bhanu Gupta and Amit Sharma P.G. Department of Mathematics, JC DAV College, Dasuya
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationInteracting Populations.
Chapter 2 Interacting Populations. 2.1 Predator/ Prey models Suppose we have an island where some rabbits and foxes live. Left alone the rabbits have a growth rate of 10 per 100 per month. Unfortunately
More informationA Stability Analysis on Models of Cooperative and Competitive Species
Research Journal of Mathematical and Statistical Sciences ISSN 2320 6047 A Stability Analysis on Models of Cooperative and Competitive Species Abstract Gideon Kwadzo Gogovi 1, Justice Kwame Appati 1 and
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
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 informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory
More informationThe major interactions between two organisms in a mixed culture are : Competition Neutralism Mutualism Commensalism Amensalism Prey-predator
1 Introduction Major classes of interaction in mixed cultures Simple models describing mixed-cultures interactions Mixed culture in nature Industrial utilization of mixed cultures 2 1 The dynamic of mixed
More informationThird In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix
Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.
More information1.2. Introduction to Modeling. P (t) = r P (t) (b) When r > 0 this is the exponential growth equation.
G. NAGY ODE January 9, 2018 1 1.2. Introduction to Modeling Section Objective(s): Review of Exponential Growth. The Logistic Population Model. Competing Species Model. Overview of Mathematical Models.
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationNONSTANDARD NUMERICAL METHODS FOR A CLASS OF PREDATOR-PREY MODELS WITH PREDATOR INTERFERENCE
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 67 75. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More informationA Discrete Numerical Scheme of Modified Leslie-Gower With Harvesting Model
CAUCHY Jurnal Matematika Murni dan Aplikasi Volume 5(2) (2018), Pages 42-47 p-issn: 2086-0382; e-issn: 2477-3344 A Discrete Numerical Scheme of Modified Leslie-Gower With Harvesting Model Riski Nur Istiqomah
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 informationDiscrete time dynamical systems (Review of rst part of Math 361, Winter 2001)
Discrete time dynamical systems (Review of rst part of Math 36, Winter 2) Basic problem: x (t);; dynamic variables (e.g. population size of age class i at time t); dynamics given by a set of n equations
More informationCOMPETITION OF FAST AND SLOW MOVERS FOR RENEWABLE AND DIFFUSIVE RESOURCE
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number, Spring 22 COMPETITION OF FAST AND SLOW MOVERS FOR RENEWABLE AND DIFFUSIVE RESOURCE SILOGINI THANARAJAH AND HAO WANG ABSTRACT. In many studies of
More information2015 Holl ISU MSM Ames, Iowa. A Few Good ODEs: An Introduction to Modeling and Computation
2015 Holl Mini-Conference @ ISU MSM Ames, Iowa A Few Good ODEs: An Introduction to Modeling and Computation James A. Rossmanith Department of Mathematics Iowa State University June 20 th, 2015 J.A. Rossmanith
More informationWe have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),
Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationExam in TMA4195 Mathematical Modeling Solutions
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Exam in TMA495 Mathematical Modeling 6..07 Solutions Problem a Here x, y are two populations varying with time
More informationCoupled differential equations
Coupled differential equations Example: dy 1 dy 11 1 1 1 1 1 a y a y b x a y a y b x Consider the case with b1 b d y1 a11 a1 y1 y a1 ay dy y y e y 4/1/18 One way to address this sort of problem, is to
More informationAnalysis of bacterial population growth using extended logistic Growth model with distributed delay. Abstract INTRODUCTION
Analysis of bacterial population growth using extended logistic Growth model with distributed delay Tahani Ali Omer Department of Mathematics and Statistics University of Missouri-ansas City ansas City,
More informationx 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 20, Number 4, Fall 1990 ON THE STABILITY OF ONE-PREDATOR TWO-PREY SYSTEMS M. FARKAS 1. Introduction. The MacArthur-Rosenzweig graphical criterion" of stability
More informationStable Coexistence of a Predator-Prey model with a Scavenger
ANTON DE KOM UNIVERSITEIT VAN SURINAME INSTITUTE FOR GRADUATE STUDIES AND RESEARCH Stable Coexistence of a Predator-Prey model with a Scavenger Thesis in partial fulfillment of the requirements for the
More informationField experiments on competition. Field experiments on competition. Field experiments on competition
INTERACTIONS BETWEEN SPECIES Type of interaction species 1 species 2 competition consumer-resource (pred, herb, para) mutualism detritivore-detritus (food is dead) Field experiments on competition Example
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationPredator-Prey Model with Ratio-dependent Food
University of Minnesota Duluth Department of Mathematics and Statistics Predator-Prey Model with Ratio-dependent Food Processing Response Advisor: Harlan Stech Jana Hurkova June 2013 Table of Contents
More informationMathematical analysis of three species model and introduction of the canonical model. Isao Kawaguchi NIRS, Japan
Mathematical analysis of three species model and introduction of the canonical model. Isao Kawaguchi NIRS, Japan Suggestions from Tatiana a) could you transform your experimental model into a generic model
More informationDynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs
Dynamics of a opulation Model Controlling the Spread of lague in rairie Dogs Catalin Georgescu The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationA NEW SOLUTION OF SIR MODEL BY USING THE DIFFERENTIAL FRACTIONAL TRANSFORMATION METHOD
April, 4. Vol. 4, No. - 4 EAAS & ARF. All rights reserved ISSN35-869 A NEW SOLUTION OF SIR MODEL BY USING THE DIFFERENTIAL FRACTIONAL TRANSFORMATION METHOD Ahmed A. M. Hassan, S. H. Hoda Ibrahim, Amr M.
More information