Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt
|
|
- Damian Dorsey
- 5 years ago
- Views:
Transcription
1 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 of the solution. Periodicity in biology: life proceeds in a rhythmic and periodic style... the periodic patterns are not easily disrupted or changed by a noisy and random environment. 1
2 Examples: Predator-prey models (lynx-hare population, sharkfish population,), nerve conduction (Hodgkin-Huxley equations, Fitzhugh-Nagumo model), chemical reactions Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt Critiques of the equation: there are infinitely many periodic solutions, and all of them are neutrally stable. Ideally, the system should have a limit cycle which all solutions tend to. 2
3 Designing better models of predator-prey: More general predation equation: (prey P, predator Q) dp dt = kp ( 1 P ) N h(p )Q h(p ) is the functional response of predators to prey population Assumptions: 1. h(0) = 0: when the population is zero, there is nothing to harvest 2. h (P ) 0: when there is more can be harvested, then the harvest rate is higher In Lotka-Volterra model: h(p ) = kp (linear predation) 3
4 Holling s type II functional response h(p ) = dp dt = kp ( 1 P ) N AP 1 + BP Q AP 1 + BP Assumptions: The number of predator is assumed to be constant, and they cannot consume more preys when P is large. It takes the predator a certain amount of time to kill and eat each prey. So suppose that in one hour, the predator (a wolf) can catch AP number of prey (rabbits) (it is proportional to P since when P is larger, the wolf has better chance to meet rabbits,) but it needs T hour to handle and eat each rabbit caught. So for all AP rabbits, it takes AT P hours, and in fact the wolf spends 1 + AT P hours on these AP rabbits. So in 1 hour, the wolf AP actually only eats rabbits. We use B = AT as a new 1 + AT P parameter in the equation. 4
5 Holling s type I model dp dt = kp ( 1 P ) N h(p ) where h(p ) = ap when 0 P P 0 and h(p ) = ap 0 when P > P 0. Assumptions: The number of predator is assumed to be constant, and they cannot consume more preys when P is large. Type I is unusual, except for the case of filter-feeding crustacea feeding on algar cells. Holling s type III model dp dt = kp spruce budworm model (lecture 8) ( 1 P ) AP n N 1 + BP n 5
6 Predator-prey model with Holling type II predation (saturating interaction) (non-dimensionalized) du ds = u(1 mu) uv u + 1, dv ds = cv + puv u + 1. u-nullcline: u = 0, v = (1 mu)(u + 1) v-nullcline: v = 0, u = c p c Equilibrium points: (0, 0), (m 1 c p(p c mc), 0), (, p c (p c) 2 ). Case 1: (p < (m + 1)c) Two equilibrium points at (0, 0) and (m 1, 0), and (m 1, 0) is a sink which attracts any initial values. Case 2: (p > (m + 1)c) Three equilibrium points at (0, 0), (m 1 c p(p c mc), 0), (u 0, v 0 ) = (, p c (p c) 2 ). Stability of (u 0, v 0 )? 6
7 A brief theory of periodic orbits Theorem 1 Inside a periodic orbit, there is at least one equilibrium point. Corollary 2 Let (x 0, y 0 ) be an equilibrium point. If there is a stable or unstable orbit of this equilibrium going to infinity as t ±, then there is no periodic orbit around this equilibrium point. Corollary 3 Let (x 0, y 0 ) be an equilibrium point. If the direction of vector field on the nullclines near (x 0, y 0 ) is not clockwise or counterclockwise, then there is no periodic orbit around this equilibrium point. 7
8 Theorem 4 (Poincaré-Bendixon) If there is a donut region (a region without equilibrium point) that the vector field is pointing inside at any point on the boundary, then there is a periodic orbit in the donut region. Corollary 5 If in a region, there is only one equilibrium point which is a source or spiral source, and the vector field is pointing inside at any point on the boundary, then there is a periodic orbit around this equilibrium point. A typical case: the solutions from outside are spiraling in, and the equilibrium point is a spiral source. 8
9 Theorem 6 (Poincaré-Andronov-Hopf bifurcation) If there is a bifurcation occurring that an equilibrium point changes from a spiral sink to spiral source, and the solutions from outside are always spiraling in, then there is a periodic orbit around the spiral source. Jules Henri Poincaré ( ): One of greatest mathemati- 9
10 cians in 20th century history/mathematicians/poincare.ht
11 Hopf Bifurcation 10
12 du ds = u(1 mu) uv u + 1, dv ds = cv + puv u + 1. Stability of (u 0, v 0 ): c(p c mp mc) J = p(p c) c p p (m + 1)c 0 Eigenvalue equation: λ 2 c(p c mp mc) λ+ c p(p c) p (p c mc) = 0 Spiral sink or sink: p c mp mc < 0 (or m > p c p + c ) Spiral source or source: p c mp mc > 0 (or 0 < m < p c p + c ) 11
13 Conclusion: m = p c p + c m > p c p + c is where Hopf bifurcation occurs. When but near p c p + c, (u 0, v 0 ) is a spiral sink, and all solutions tend to (u 0, v 0 ) in an oscillating fashing; when m passes p c p + c but near p c p + c, (u 0, v 0 ) becomes a spiral source, but all solutions away from (u 0, v 0 ) still tends toward (u 0, v 0 ) (spiral inward), thus a periodic solution (with small oscillation) emerges around (u 0, v 0 ), and it is a limit cycle. Examples: p = 1, c = 0.5, the Hopf bifurcation point is m = 1/
14 Hopf bifurcation destablizes the equilibrium (from stable spiral sink to unstable spiral source) du ds = u(1 mu) uv u + 1, dv ds = cv + puv u + 1. What can cause Hopf bifurcation: m is small! m 1 is carrying capacity, so that m is smaller implies carrying capacity is larger. Larger carrying capacity for the prey is called enrichment, which seems lead to stabilizing. But instead it leads to Hopf bifurcation and destablizing of the equilibrium. This is called paradox of enrichment 13
15 Hodgkin-Huxley model: v (t) = 1 C [g Nam 3 h(v v Na ) + g K n 4 (v v K ) + g L (v v L )] n (t) = α n (v)(1 n) β n (v)n, m (t) = α m (v)(1 m) β m (v)m, h (t) = α h (v)(1 h) β h (v)h, where α n (v) = 0.01(v + 10)(e (v+10)/10 1) 1, β n (v) = 0.125e v/80 α m (v) = 0.01(v + 25)(e (v+25)/10 1) 1, β m (v) = 4e v/18 α h (v) = 0.07e v/20, β h (v) = (e (v+30)/10 + 1) 1 and g Na = 120, g K = 36, g L = 0.3, v Na = 115, v K = 12, and v L = Numerical simulation of HH-model was a great success! 14
16 Hodgkin, A., and Huxley, A. (1952): A quantitative description of membrane current and its application to conduction and excitation in nerve. Journal of Physiology 117: Alan Lloyd Hodgkin( ), and Andrew Fielding Huxley (1917-) won the 1963 Nobel Prize in Physiology or Medicine for their work on the basis of nerve action potentials, the electrical impulses that enable the activity of an organism to be coordinated by a central nervous system. Hodgkin and Huxley shared the prize that year with John Carew Eccles, who was cited for research on synapses. Hodgkin and Huxley s findings led the pair to hypothesize the existence of ion channels, which were isolated only decades later. Confirmation of ion channels came with the development of the patch clamp, which led to a Nobel prize in 1991 to Erwin Neher and Bert Sakmann. 15
17 FitzHugh-Nagumo equation: a simpler system but retain many of qualitative features of Hodgkin-Huxley system (and the simplified 2-d system of FitzHugh) FitzHugh (1961), Nagumo et. al. (1962) ɛ dv dw = v(v a)(1 v) w, = v bw c. dt dt v(t) is the excitability of the system (voltage), w(t) is recovery variable representing the force that tends to return the resting state Excitability: small perturbation will return to resting state quickly, and a larger perturbation will cause a full circle going around before returning to the resting state. 16
18 Case 1: ɛ = 0.01, a = 0.1, b = 0.5 and c = 0 Equilibrium (0, 0) stable Case 2: ɛ = 0.01, a = 0.1, b = 0.5 and c = 0.1 Equilibrium (0, 0) unstable, and a periodic solution 17
Math 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction
Math 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction Junping Shi College of William and Mary November 8, 2018 Neuron Neurons Neurons are cells in the brain and other subsystems
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
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 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 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 informationModeling Action Potentials in Cell Processes
Modeling Action Potentials in Cell Processes Chelsi Pinkett, Jackie Chism, Kenneth Anderson, Paul Klockenkemper, Christopher Smith, Quarail Hale Tennessee State University Action Potential Models Chelsi
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 informationNonlinear Dynamics. Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna.
Nonlinear Dynamics Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna http://www.moreno.marzolla.name/ 2 Introduction: Dynamics of Simple Maps 3 Dynamical systems A dynamical
More informationTHE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL
THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd
More information1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model
1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model Alan Hodgkin and Andrew Huxley developed the first quantitative model of the propagation of an electrical signal (the action potential)
More informationFrom neuronal oscillations to complexity
1/39 The Fourth International Workshop on Advanced Computation for Engineering Applications (ACEA 2008) MACIS 2 Al-Balqa Applied University, Salt, Jordan Corson Nathalie, Aziz Alaoui M.A. University of
More informationLast lecture: Recurrence relations and differential equations. The solution to the differential equation dx
Last lecture: Recurrence relations and differential equations The solution to the differential equation dx = ax is x(t) = ce ax, where c = x() is determined by the initial conditions x(t) Let X(t) = and
More informationDynamical Systems in Neuroscience: Elementary Bifurcations
Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................
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 informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
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 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 informationStability and bifurcation in a two species predator-prey model with quintic interactions
Chaotic Modeling and Simulation (CMSIM) 4: 631 635, 2013 Stability and bifurcation in a two species predator-prey model with quintic interactions I. Kusbeyzi Aybar 1 and I. acinliyan 2 1 Department of
More informationOn predator-prey models
Predator-prey On models ddd Department of Mathematics College of William and Mary Math 41/CSUMS Talk February 3, 1 Collaborators Sze-Bi Hsu (Tsinghua University, Hsinchu, Taiwan) Junjie Wei (Harbin Institute
More informationFinal Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations
Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger Project I: Predator-Prey Equations The Lotka-Volterra Predator-Prey Model is given by: du dv = αu βuv = ρβuv
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 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 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 informationLIMIT CYCLE OSCILLATORS
MCB 137 EXCITABLE & OSCILLATORY SYSTEMS WINTER 2008 LIMIT CYCLE OSCILLATORS The Fitzhugh-Nagumo Equations The best example of an excitable phenomenon is the firing of a nerve: according to the Hodgkin
More informationMathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued
Mathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland
More informationΝευροφυσιολογία και Αισθήσεις
Biomedical Imaging & Applied Optics University of Cyprus Νευροφυσιολογία και Αισθήσεις Διάλεξη 5 Μοντέλο Hodgkin-Huxley (Hodgkin-Huxley Model) Response to Current Injection 2 Hodgin & Huxley Sir Alan Lloyd
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 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 informationSignal processing in nervous system - Hodgkin-Huxley model
Signal processing in nervous system - Hodgkin-Huxley model Ulrike Haase 19.06.2007 Seminar "Gute Ideen in der theoretischen Biologie / Systembiologie" Signal processing in nervous system Nerve cell and
More informationMotivation and Goals. Modelling with ODEs. Continuous Processes. Ordinary Differential Equations. dy = dt
Motivation and Goals Modelling with ODEs 24.10.01 Motivation: Ordinary Differential Equations (ODEs) are very important in all branches of Science and Engineering ODEs form the basis for the simulation
More informationPresented by Sarah Hedayat. Supervised by Pr.Cappy and Dr.Hoel
1 Presented by Sarah Hedayat Supervised by Pr.Cappy and Dr.Hoel Outline 2 Project objectives Key elements Membrane models As simple as possible Phase plane analysis Review of important Concepts Conclusion
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 informationComputational Neuroscience. Session 4-2
Computational Neuroscience. Session 4-2 Dr. Marco A Roque Sol 06/21/2018 Two-Dimensional Two-Dimensional System In this section we will introduce methods of phase plane analysis of two-dimensional systems.
More informationCPS 5310: Parameter Estimation. Natasha Sharma, Ph.D.
Example Suppose our task is to determine the net income for year 2019 based on the net incomes given below Year Net Income 2016 48.3 million 2017 90.4 million 2018 249.9 million Last lecture we tried to
More informationGetting Started With The Predator - Prey Model: Nullclines
Getting Started With The Predator - Prey Model: Nullclines James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline The Predator
More informationPredator-Prey Population Models
21 Predator-Prey Population Models Tools Used in Lab 21 Hudson Bay Data (Hare- Lynx) Lotka-Volterra Lotka-Volterra with Harvest How can we model the interaction between a species of predators and their
More informationBifurcation and Stability Analysis of a Prey-predator System with a Reserved Area
ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
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 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 informationSlow Manifold of a Neuronal Bursting Model
Slow Manifold of a Neuronal Bursting Model Jean-Marc Ginoux 1 and Bruno Rossetto 2 1 PROTEE Laboratory, Université du Sud, B.P. 2132, 83957, La Garde Cedex, France, ginoux@univ-tln.fr 2 PROTEE Laboratory,
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 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 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 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 informationWave Pinning, Actin Waves, and LPA
Wave Pinning, Actin Waves, and LPA MCB 2012 William R. Holmes Intercellular Waves Weiner et al., 2007, PLoS Biology Dynamic Hem1 waves in neutrophils Questions? How can such waves / pulses form? What molecular
More informationSystems Biology: Theoretical Biology. Kirsten ten Tusscher, Theoretical Biology, UU
Systems Biology: Theoretical Biology Kirsten ten Tusscher, Theoretical Biology, UU 1 Chapter 9 Hodgkin-Huxley model 2 Introduction Hodgkin-Huxley model : detailed model for nerve action potential Perfect
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 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 informationDynamical systems in neuroscience. Pacific Northwest Computational Neuroscience Connection October 1-2, 2010
Dynamical systems in neuroscience Pacific Northwest Computational Neuroscience Connection October 1-2, 2010 What do I mean by a dynamical system? Set of state variables Law that governs evolution of state
More informationFish Population Modeling
Fish Population Modeling Thomas Wood 5/11/9 Logistic Model of a Fishing Season A fish population can be modelled using the Logistic Model, P = kp(1 P/a) P is the population at time t k is a growth constant
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 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 informationNonlinear Dynamics of Neural Firing
Nonlinear Dynamics of Neural Firing BENG/BGGN 260 Neurodynamics University of California, San Diego Week 3 BENG/BGGN 260 Neurodynamics (UCSD) Nonlinear Dynamics of Neural Firing Week 3 1 / 16 Reading Materials
More informationJournal of the Vol. 36, pp , 2017 Nigerian Mathematical Society c Nigerian Mathematical Society
Journal of the Vol. 36, pp. 47-54, 2017 Nigerian Mathematical Society c Nigerian Mathematical Society A CLASS OF GENERALIZATIONS OF THE LOTKA-VOLTERRA PREDATOR-PREY EQUATIONS HAVING EXACTLY SOLUBLE SOLUTIONS
More informationdv dt Predator-Prey Models
Predator-Prey Models This is a diverse area that includes general models of consumption: Granivores eating seeds Parasitoids Parasite-host interactions Lotka-Voterra model prey and predator: V = victim
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 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 informationIntegral Equation and PDE Neuronal Models:
Integral Equation and PDE Neuronal Models: 1952-2015 W. C. Troy Integral Equation and PDE Neuronal Models: 1952-2015 p.1/51 Part 1. 1D Single Neuron Models. (i) Hodgkin-Huxley (1952), Fitzhugh-Nagumo (1961)
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 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 informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
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 informationStability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations
Stability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations Sreelakshmi Manjunath Department of Electrical Engineering Indian Institute of Technology Madras (IITM), India JTG Summer
More informationTopics in Neurophysics
Topics in Neurophysics Alex Loebel, Martin Stemmler and Anderas Herz Exercise 2 Solution (1) The Hodgkin Huxley Model The goal of this exercise is to simulate the action potential according to the model
More informationDynamical Systems for Biology - Excitability
Dynamical Systems for Biology - Excitability J. P. Keener Mathematics Department Dynamical Systems for Biology p.1/25 Examples of Excitable Media B-Z reagent CICR (Calcium Induced Calcium Release) Nerve
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 information3 Action Potentials - Brutal Approximations
Physics 172/278 - David Kleinfeld - Fall 2004; Revised Winter 2015 3 Action Potentials - Brutal Approximations The Hodgkin-Huxley equations for the behavior of the action potential in squid, and similar
More informationOutline. Calculus for the Life Sciences. What is a Differential Equation? Introduction. Lecture Notes Introduction to Differential Equa
Outline Calculus for the Life Sciences Lecture Notes to Differential Equations Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu 1 Department of Mathematics and Statistics Dynamical Systems Group Computational
More informationCellular Electrophysiology. Cardiac Electrophysiology
Part 1: Resting and Action Potentials Cardiac Electrophysiology Theory Simulation Experiment Scale The membrane: structure, channels and gates The cell: resting potential, whole cell currents, cardiac
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 informationChapter 7. Nonlinear Systems. 7.1 Introduction
Nonlinear Systems Chapter 7 The scientist does not study nature because it is useful; he studies it because he delights in it, and he delights in it because it is beautiful. - Jules Henri Poincaré (1854-1912)
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 informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
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 informationBME 5742 Biosystems Modeling and Control
BME 5742 Biosystems Modeling and Control Hodgkin-Huxley Model for Nerve Cell Action Potential Part 1 Dr. Zvi Roth (FAU) 1 References Hoppensteadt-Peskin Ch. 3 for all the mathematics. Cooper s The Cell
More informationEnrichment in a Producer-Consumer Model with varying rates of Stoichiometric Elimination
Enrichment in a Producer-Consumer Model with varying rates of Stoichiometric Elimination Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth
More information1 Introduction and neurophysiology
Dynamics of Continuous, Discrete and Impulsive Systems Series B: Algorithms and Applications 16 (2009) 535-549 Copyright c 2009 Watam Press http://www.watam.org ASYMPTOTIC DYNAMICS OF THE SLOW-FAST HINDMARSH-ROSE
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 informationHopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models
Vol.37 (SUComS 6), pp.57-66 http://dx.doi.org/.457/astl.6.37. Hopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models Shuang Guo, Xiuli Li, Jian Yu, Xiu Ren Department of Mathematics, Teacher
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 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 information2013 NSF-CMACS Workshop on Atrial Fibrillation
2013 NSF-CMACS Workshop on A Atrial Fibrillation Flavio H. Fenton School of Physics Georgia Institute of Technology, Atlanta, GA and Max Planck Institute for Dynamics and Self-organization, Goettingen,
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 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 informationLecture Notes 8C120 Inleiding Meten en Modelleren. Cellular electrophysiology: modeling and simulation. Nico Kuijpers
Lecture Notes 8C2 Inleiding Meten en Modelleren Cellular electrophysiology: modeling and simulation Nico Kuijpers nico.kuijpers@bf.unimaas.nl February 9, 2 2 8C2 Inleiding Meten en Modelleren Extracellular
More informationMA 138 Calculus 2 with Life Science Applications Handout
.. MA 138 Calculus 2 with Life Science Applications Handout Alberto Corso alberto.corso@uky.edu Department of Mathematics University of Kentucky February 17, 2017 . Example 4 (Lotka-Volterra Predator-Prey
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 informationGlobal Stability Analysis on a Predator-Prey Model with Omnivores
Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani
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 informationInteractions between predators and prey
Interactions between predators and prey What is a predator? Predator An organism that consumes other organisms and inevitably kills them. Predators attack and kill many different prey individuals over
More informationLecture 10 : Neuronal Dynamics. Eileen Nugent
Lecture 10 : Neuronal Dynamics Eileen Nugent Origin of the Cells Resting Membrane Potential: Nernst Equation, Donnan Equilbrium Action Potentials in the Nervous System Equivalent Electrical Circuits and
More informationElectrophysiology of the neuron
School of Mathematical Sciences G4TNS Theoretical Neuroscience Electrophysiology of the neuron Electrophysiology is the study of ionic currents and electrical activity in cells and tissues. The work of
More informationLocal Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)
Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with
More informationGerardo Zavala. Math 388. Predator-Prey Models
Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked
More informationModeling Prey-Predator Populations
Modeling Prey-Predator Populations Alison Pool and Lydia Silva December 13, 2006 Alison Pool and Lydia Silva () Modeling Prey-Predator Populations December 13, 2006 1 / 25 1 Introduction 1 Our Populations
More informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA November 5, 2011 1d systems To understand more complex
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More information