Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations
|
|
- Kelley Webb
- 5 years ago
- Views:
Transcription
1 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 γv where αu models the birth of the prey, βuv models death of the prey by being captured by predators, ρβuv models the rate of birth of the predators (which depends on how many prey are captured), and γv models the death of the predators. The basic model exhibits solutions which have periodic oscillating predator and prey populations. There is one non-trivial steady-state which is a center which is stable, but is not asymptotically stable. Consequently, the initial conditions for the predator and prey population size play an important role in determining the overall dynamics. To obtain the basic Lotka-Volterra Predator-Prey Model a number of assumptions were made. In particular, the predators are assumed to die at a constant per capita rate and the prey are assumed to reproduce at a constant per capita rate. The predator birth rate is assumed to depend only linearly on the rate at which predators catch prey. While the basic model may be a good starting point in modeling predator-prey populations, it neglects many features of populations which may be potentially important in practice. For example, the effect of the different seasons of the year on the birth rate of the prey and the death rate of the predators was not taken into account. For instance, in the winter we might expect the birth rate of the prey population to be greatly reduced and the death rate of the predators to increase. An interesting question is how these periodic variations effect the population dynamics, which is already prone to its own intrinsic periodic oscillations through the predator-prey interactions. In this project you will explore modifications of the basic predator-prey model using both analytical and computational techniques. One way to incorporate seasonality into the model is to modify the constant rates above, say β, to be β(1 + ǫ sin(2πωt)) where ǫ controls the influence of the seasonality and ω controls the frequency of the seasons. Use matlab codes to simulate the model and explore the influence of seasonality separately for each rate constant, then try combinations. In your report carefully address the following questions: (i) How does the dynamics of the system change qualitatively when seasonality is incorporated into the various rate constants? (ii) Which rate constants have the most dramatic effect on the dynamics when modified to incorporate seasonality? (Birth of prey, death of predators, rate of predation, or are they all comparable?)
2 (iii) What influence does ω have on the effects of the seasonality? If ω is very large is this any different than when ω is small? What about the effect of ǫ small verse very large? Be careful to keep ǫ small enough so that all rates remain positive. (iv) Are there steady-states or limit cycles for the system? If so, try to determine their stability either by linearization or by arguing in terms of phase-plane portrait. Give a biological interpretation of why the steady-states or limit cycles occur and how they are stabilised. On the class website you can find numerical codes for computing trajectories of the system and for the computation of Poincare maps.
3 Project II: Ecological Modeling Species of plankton grow on the surface of the oceans and play a fundamental role in the oceanic food chain. They can be classified into two broad categories, phytoplankton and zooplankton. Phytoplankton are plant species that use photosynthesis to sustain themselves and are often limited in their growth by available nutrients in the environment. Zooplankton are animal species and feed on the phytoplankton. In this project you will explore the effects of the nutrients in the environment and the interaction of the two species with one another. A basic model for the plankton populations and nutrients is the following: dp dz dn = cnp dpz ap = dpz bz = ap + bz cnp where P, Z model the phytoplanton and zooplankton population sizes and N models the amount of the available nutrient nitrogen. The ap, bz terms model the rate of nitrogen production which occurs when plankton die and cnp models the rate of consumption of the nitrogen by phytoplanton. The terms appearing in the equations for P and Z follow the usual predator-prey interpretation (see Project I). Some species of phytoplankton form massive colonies of cocolithiphore clusters on the surface of the ocean and produce the gas dimethyl sulfide which is thought to play an important role in controlling the buoyancy of the clusters. When dimethyl sulfide is released from the ocean surface into the atmosphere it becomes a good base for condensation of water and the formation of clouds. One hypotheses is that the phytoplankton, which depend on sunlight for photosynthesis, could in principle form a regulating loop with the climate through the role dimethyl sulfide plays in cloud formation. In this project the hypothesized feedback loop and its effects on the dynamics of the plankton populations will be explored. In the model the reproduction of the phytoplankton was modeled as depending on the size of the phytoplanton population and the amount of the nutrient nitrogen with a constant coefficient c. Since the rate of reproduction of the phytoplanton population is expected to depend on the amount of sunlight available for photosynthesis we can incorporate effects of cloud cover by using a non-constant coefficient c. One simple model is to modify the photosynthetic rate constant c to c(1 Q), where Q is the average density of clouds in the sky having a range from 0 for no clouds to 1 for an overcast sky saturated with dark clouds. To model how the clouds depend on the phytoplankton population we could use a logistic equation for 1 Q with an effective carrying capacity of 1. This yields the model equations P for Q: dq = r(1 Q)(1 P(1 Q)). In the project use the analytic methods discussed in class and the matlab codes to address the following questions:
4 (i) Does the dynamics change significantly when the cloud-phytoplankton feedback loop is incorporated into the model? (ii) If we introduce a delay τ into the equation for the cloud formation so that P(t) is replaced by P(t τ) how does this change the dynamics? (iii) What are the steady-states or limit-cycles for the model with and without delays? Which are stable? Argue using linearization theory and plots of the dynamical trajectories. (iv) What happens in the model if the cloud formation is influenced by the seasons? In particular, what happens if the seasons amplify or diminish the effects of the phytoplankton on cloud formation in the sense that P is modified to (1+ǫ sin(2πωt))p. For example, what happens when ω is small? When it is large? On the class website can be found numerical codes to solve the ODEs.
5 Project III: Epidemiology of Diseases The spread of many diseases can be modeled by classifying individuals in a population into one of three classes. An individual can either be disease-free but susceptible (S), infected with the disease and contagious (I), or immune to the disease and recovered (R). Many features of how the disease behaves in the population can be described in terms of how the size of these classes changes over time and by the steady-states of the system. We showed in class how to predict whether a disease will lead to an epidemic and how to estimate the total number of individuals infected over the course of the epidemic. We also discussed how to predict whether a disease will be endemic and how to estimate the fraction of the population that will be chronically infected. The effectiveness of various vaccination strategies was also analyzed. In many such disease models the law of mass action is used to model the interactions between susceptible and infected members of the population. This is based on the underlying assumption that the classes are well mixed so that the members of each of the classes are uniformly spread-out over the regions occupied by the population. For many diseases this assumption is unrealistic and in fact does not give good results, for example, some diseases depend importantly on social interactions and lifestyle habits which play a central role in transmission. In this project you will explore the role that non-uniform interactions play in the dynamics of diseases. One approach to modeling non-uniform interactions of the population is to use a graph (social network) to model the interactions between individuals. The individuals are modeled as the nodes of a graph where an edge is formed between any pair of individuals who are in contact. The population members then interact using the specific structure of the graph. To model the network of interactions we shall use a random graph where an edge is formed between a pair of individuals with probability 0 p 1. Using the random graph one infected individual is introduced into the population. Over the time interval δτ the following cycle of events occur to update the population. (i) Each individual who is connected to an infected member has probability βδτ of becoming infected. (ii) Infected individuals have a probability of γδτ of recovering from the disease over the time step. We exclude the newly infected individuals from consideration in this latter step. Using matlab codes and the analytic techniques discuss in class explore features of this model. In your report the following questions should be addressed: (i) What are the sizes of the S, I, R classes as the disease spreads throughout the population? How does the parameter p influence the occurrence of and time course of the epidemic? Is there a critical value of p? Is there an endemic state? What is its size? (ii) How does the random network model compare to the S, I, R models discussed in class when the same parameters are chosen for β,γ, etc..? Does the size N of the population play a role? What happens when N is small? What appears to be the range of values for p that gives good agreement with the S,I,R models in class?
6 (iii) How does the model change qualitatively when disease-unrelated death and disease-related death are added to the model? Does the model agree well with the results derived in class? What range of p leads to good agreement? For what values of p does the model behave differently? On the class website can be found codes to generate the random graphs and to perform the updates for the dynamics of the model.
7 Project IV: Neuroscience The Hodgkin-Huxley-Katz equations originally derived for experiments on giant squid axons provides a mathematical model for how electrical signals are propagated in many neural tissues. The equations model how individual neurons respond to a stimulating depolarising pulse of current and how the subsequent action potential is propagated along the axon. The propagation of the action potential occurs in four phases: (i) a positive feedback phase where the depolarisation is increased by the opening of sodium ion channels, (ii) an excited phase of reversed polarisation where potassium ion channels are opened to counterbalance the reverse polarisation, (iii) a hyperpolarised phase where the sodium ion channels are closed and the neuron begins its recovery to the rest-state, (iv) a refractory period where the rest-polarisation has been substantially recovered but the sodium channels remain inactive. In this project you will explore how a neuron responds to various types of stimuli using a variant of the Hodgkin-Huxley-Katz model. To make the problem more tractable we shall work with a reduced system of equations derived by FitzHugh and Nagumo which model only two variables, an effective conductance and voltage associated with an action potential. The FitzHugh-Nagumo equations are: ǫ dv dw = f(v,w) + I(t) = g(v, w) where f(v,w) = v(v a)(1 v) w, g(v,w) = v c bw, and the value of ǫ > 0 is typically small. The function I(t) represents the input stimuli to the neuron. Using matlab codes and the analytic techniques discussed in class you will explore the response of the neuron for various types of stimulating currents I(t). In your report address the following questions: (i) What does the time series for the voltage look like when the input stimuli is a single pulse of current? What magnitude and duration is required to get a response from the neuron? Analyze this in terms of the phase-plane of the equations. (ii) What does the time series for the voltage look like when the input stimuli is a constant current? How does this depend on the magnitude of the current? Is there a critical magnitude? How does the constant current change the phase-plane of the system? (iii) How does the time series change if we make the input current a sequence of pulses of fixed duration? What types of behavior does the model exhibit as the frequency and duration of the pulses are varied? How might this behavior arise? Try to explain it in terms of the phase-plane. (iv) What happens if we make the input current white noise? How does the behavior depend on the magnitude of the noise? What if we make the stimulating currents pulsed again but with an exponential waiting time between pulses? Give an explanation of what happens in terms of trajectories in the phase-plane.
8 On the class website can be found numerical codes to integrate the FitzHugh-Nugamo equations and to generate the currents I(t).
Predation is.. The eating of live organisms, regardless of their identity
Predation Predation Predation is.. The eating of live organisms, regardless of their identity Predation 1)Moves energy and nutrients through the ecosystem 2)Regulates populations 3)Weeds the unfit from
More information6.3.4 Action potential
I ion C m C m dφ dt Figure 6.8: Electrical circuit model of the cell membrane. Normally, cells are net negative inside the cell which results in a non-zero resting membrane potential. The membrane potential
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 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 informationKasetsart University Workshop. Mathematical modeling using calculus & differential equations concepts
Kasetsart University Workshop Mathematical modeling using calculus & differential equations concepts Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu
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 informationNeural Modeling and Computational Neuroscience. Claudio Gallicchio
Neural Modeling and Computational Neuroscience Claudio Gallicchio 1 Neuroscience modeling 2 Introduction to basic aspects of brain computation Introduction to neurophysiology Neural modeling: Elements
More informationEcology Regulation, Fluctuations and Metapopulations
Ecology Regulation, Fluctuations and Metapopulations The Influence of Density on Population Growth and Consideration of Geographic Structure in Populations Predictions of Logistic Growth The reality of
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 informationDiscrete and Indiscrete Models of Biological Networks
Discrete and Indiscrete Models of Biological Networks Winfried Just Ohio University November 17, 2010 Who are we? What are we doing here? Who are we? What are we doing here? A population of interacting
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 informationProject 1 Modeling of Epidemics
532 Chapter 7 Nonlinear Differential Equations and tability ection 7.5 Nonlinear systems, unlike linear systems, sometimes have periodic solutions, or limit cycles, that attract other nearby solutions.
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 informationPopulations Study Guide (KEY) All the members of a species living in the same place at the same time.
Populations Study Guide (KEY) 1. Define Population. All the members of a species living in the same place at the same time. 2. List and explain the three terms that describe population. a. Size. How large
More informationBo Deng University of Nebraska-Lincoln UNL Math Biology Seminar
Mathematical Model of Neuron Bo Deng University of Nebraska-Lincoln UNL Math Biology Seminar 09-10-2015 Review -- One Basic Circuit By Kirchhoff's Current Law 0 = I C + I R + I L I ext By Kirchhoff s Voltage
More informationQualitative Analysis of a Discrete SIR Epidemic Model
ISSN (e): 2250 3005 Volume, 05 Issue, 03 March 2015 International Journal of Computational Engineering Research (IJCER) Qualitative Analysis of a Discrete SIR Epidemic Model A. George Maria Selvam 1, D.
More informationPopulation Ecology and the Distribution of Organisms. Essential Knowledge Objectives 2.D.1 (a-c), 4.A.5 (c), 4.A.6 (e)
Population Ecology and the Distribution of Organisms Essential Knowledge Objectives 2.D.1 (a-c), 4.A.5 (c), 4.A.6 (e) Ecology The scientific study of the interactions between organisms and the environment
More informationNonlinear Systems in Applications
Chapter 10 Nonlinear Systems in Applications 375 This chapter is a survey of a few nonlinear systems that arise in biology, chemistry, and mechanics. Our intent is to describe to some extent both the modeling
More informationNon-Linear Models Cont d: Infectious Diseases. Non-Linear Models Cont d: Infectious Diseases
Cont d: Infectious Diseases Infectious Diseases Can be classified into 2 broad categories: 1 those caused by viruses & bacteria (microparasitic diseases e.g. smallpox, measles), 2 those due to vectors
More informationYakın Doğu Üniversitesi Mimarlık Fakültesi Peyzaj Mimarlığı Bölümü. PM 317 Human and Environment Assoc. Prof. Dr. Salih GÜCEL
Yakın Doğu Üniversitesi Mimarlık Fakültesi Peyzaj Mimarlığı Bölümü PM 317 Human and Environment Assoc. Prof. Dr. Salih GÜCEL Ecology & Ecosystems Principles of Ecology Ecology is the study of the interactions
More informationOrganizing Diversity Taxonomy is the discipline of biology that identifies, names, and classifies organisms according to certain rules.
1 2 3 4 5 6 7 8 9 10 Outline 1.1 Introduction to AP Biology 1.2 Big Idea 1: Evolution 1.3 Big Idea 2: Energy and Molecular Building Blocks 1.4 Big Idea 3: Information Storage, Transmission, and Response
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 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 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 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 informationthey give no information about the rate at which repolarization restores the
497 J. Physiol. (1952) ii6, 497-506 THE DUAL EFFECT OF MEMBRANE POTENTIAL ON SODIUM CONDUCTANCE IN THE GIANT AXON OF LOLIGO BY A. L. HODGKIN AND A. F. HUXLEY From the Laboratory of the Marine Biological
More informationIntroduction. Chapter What is this book about?
Chapter 1 Introduction 1.1 What is this book about? This book is about how to construct and use computational models of specific parts of the nervous system, such as a neuron, a part of a neuron or a network
More informationIntroduction to SEIR Models
Department of Epidemiology and Public Health Health Systems Research and Dynamical Modelling Unit Introduction to SEIR Models Nakul Chitnis Workshop on Mathematical Models of Climate Variability, Environmental
More informationBIO S380T Page 1 Summer 2005: Exam 2
BIO S380T Page 1 Part I: Definitions. [5 points for each term] For each term, provide a brief definition that also indicates why the term is important in ecology or evolutionary biology. Where I ve provided
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 informationSocial Influence in Online Social Networks. Epidemiological Models. Epidemic Process
Social Influence in Online Social Networks Toward Understanding Spatial Dependence on Epidemic Thresholds in Networks Dr. Zesheng Chen Viral marketing ( word-of-mouth ) Blog information cascading Rumor
More informationSpotlight on Modeling: The Possum Plague
70 Spotlight on Modeling: The Possum Plague Reference: Sections 2.6, 7.2 and 7.3. The ecological balance in New Zealand has been disturbed by the introduction of the Australian possum, a marsupial the
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best) Koch, Chapters 6, 8, 9 See also Clay, J. Neurophysiol. 80:903-913 (1998) (for a recent version of the HH squid axon model) Rothman
More informationThe dynamics of disease transmission in a Prey Predator System with harvesting of prey
ISSN: 78 Volume, Issue, April The dynamics of disease transmission in a Prey Predator System with harvesting of prey, Kul Bhushan Agnihotri* Department of Applied Sciences and Humanties Shaheed Bhagat
More informationINTERPRETING POPULATION DYNAMICS GRAPH
INTERPRETING POPULATION DYNAMIS GRAPH OJETIVES TASKS Name: To learn about three types of population dynamics graphs To determine which type of graph you constructed from the Pike and Perch Game To interpret
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 informationAPES Fall Final REVIEW
Class: Date: APES Fall Final REVIEW Short Answer 1. The difference between chemical and physical weathering of rock is that 2. The difference between weathering and erosion is that 3. Select the correct
More informationBIOS 3010: ECOLOGY. Dr Stephen Malcolm. Laboratory 6: Lotka-Volterra, the logistic. equation & Isle Royale
BIOS 3010: ECOLOGY Dr Stephen Malcolm Laboratory 6: Lotka-Volterra, the logistic equation & Isle Royale This is a computer-based activity using Populus software (P), followed by EcoBeaker analyses of moose
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 informationMath 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 informationA NUMERICAL STUDY ON PREDATOR PREY MODEL
International Conference Mathematical and Computational Biology 2011 International Journal of Modern Physics: Conference Series Vol. 9 (2012) 347 353 World Scientific Publishing Company DOI: 10.1142/S2010194512005417
More informationEcosystems. 2. Ecosystem
1. Studying our living Planet The biosphere consist of all life on Earth and all parts of the Earth in which life exists, including land, water, and the atmosphere. Ecology is the scientific study of interactions
More informationCHAPTER 52 Study Questions (An Introduction to Ecology and the Biosphere)
WLHS / AP Bio / Monson Name CHAPTER 52 Study Questions (An Introduction to Ecology and the Biosphere) 52.1: Earth s climate varies by latitude and season and is changing rapidly (p. 1144-1150) 1) Distinguish
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 informationMATH3203 Lecture 1 Mathematical Modelling and ODEs
MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1
More informationMarine Resources Development Foundation/MarineLab Grades: 9, 10, 11, 12 States: AP Biology Course Description Subjects: Science
Marine Resources Development Foundation/MarineLab Grades: 9, 10, 11, 12 States: AP Biology Course Description Subjects: Science Highlighted components are included in Tallahassee Museum s 2016 program
More information3.3 Simulating action potentials
6 THE HODGKIN HUXLEY MODEL OF THE ACTION POTENTIAL Fig. 3.1 Voltage dependence of rate coefficients and limiting values and time constants for the Hodgkin Huxley gating variables. (a) Graphs of forward
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 informationThursday. Threshold and Sensitivity Analysis
Thursday Threshold and Sensitivity Analysis SIR Model without Demography ds dt di dt dr dt = βsi (2.1) = βsi γi (2.2) = γi (2.3) With initial conditions S(0) > 0, I(0) > 0, and R(0) = 0. This model can
More information1. Population dynamics of rabbits and foxes
1. Population dynamics of rabbits and foxes (a) A simple Lotka Volterra Model We have discussed in detail the Lotka Volterra model for predator-prey relationships dn prey dt = +R prey,o N prey (t) γn prey
More informationEssential Ideas of Mathematical Modeling in Population Dynamics
Essential Ideas of Mathematical Modeling in Population Dynamics Toward the application for the disease transmission dynamics Hiromi SENO Research Center for Pure and Applied Mathematics Department of Computer
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 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 informationUnit 8: Ecology: Ecosystems and Communities
Unit 8: Ecology: Ecosystems and Communities An ecosystem consists of all the plants and animals that interact with the nonliving things in an area. Biosphere = area on Earth where living things are found
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 informationDelay SIR Model with Nonlinear Incident Rate and Varying Total Population
Delay SIR Model with Nonlinear Incident Rate Varying Total Population Rujira Ouncharoen, Salinthip Daengkongkho, Thongchai Dumrongpokaphan, Yongwimon Lenbury Abstract Recently, models describing the behavior
More informationName Student ID. Good luck and impress us with your toolkit of ecological knowledge and concepts!
Page 1 BIOLOGY 150 Final Exam Winter Quarter 2000 Before starting be sure to put your name and student number on the top of each page. MINUS 3 POINTS IF YOU DO NOT WRITE YOUR NAME ON EACH PAGE! You have
More informationUnit 1 Lesson 3 Population Dynamics. Copyright Houghton Mifflin Harcourt Publishing Company
Movin Out How can a population grow or get smaller? If new individuals are added to the population, it grows. If individuals are removed from a population, it gets smaller. The population stays at about
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
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 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 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 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 informationA.P. Biology CH Population Ecology. Name
1 A.P. Biology CH. 53 - Population Ecology Name How many ants (shown below - 6 ants / cm 2 ) would there be in an ant colony that is flat and one meter long on each side? Dispersion Patterns Matching A
More informationSingle-Compartment Neural Models
Single-Compartment Neural Models BENG/BGGN 260 Neurodynamics University of California, San Diego Week 2 BENG/BGGN 260 Neurodynamics (UCSD) Single-Compartment Neural Models Week 2 1 / 18 Reading Materials
More informationENVE203 Environmental Engineering Ecology (Nov 05, 2012)
ENVE203 Environmental Engineering Ecology (Nov 05, 2012) Elif Soyer Ecosystems and Living Organisms Population Density How Do Populations Change in Size? Maximum Population Growth Environmental Resistance
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best Koch, Chapters 6, 8, 9 See also Hodgkin and Huxley, J. Physiol. 117:500-544 (1952. (the source Clay, J. Neurophysiol. 80:903-913
More informationAPPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics
APPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics Due: Thursday, December 6, 2018 by 4:59 p.m. Submit as a PDF to Assignments on Canvas 1 Introduction Cardiac Arrhythmia, or irregular
More informationPopulation Ecology. Text Readings. Questions to Answer in the Chapter. Chapter Reading:
Population Ecology Text Readings Chapter Reading: Chapter # 26 in Audesirk, Audesirk and Byers: Population Growth and Regulation Pg. # 513-534. Questions to Answer in the Chapter How Does Population Size
More informationDynamical Analysis of a Harvested Predator-prey. Model with Ratio-dependent Response Function. and Prey Refuge
Applied Mathematical Sciences, Vol. 8, 214, no. 11, 527-537 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/12988/ams.214.4275 Dynamical Analysis of a Harvested Predator-prey Model with Ratio-dependent
More informationConverting energy into nerve impulses, resting potentials and action potentials Sensory receptors
D 1.3 s Converting energy into nerve impulses, resting potentials and action potentials Sensory receptors A receptor converts an external or internal stimulus into an electrical signal. Sensory receptors
More information14.1. KEY CONCEPT Every organism has a habitat and a niche. 38 Reinforcement Unit 5 Resource Book
14.1 HABITAT AND NICHE KEY CONCEPT Every organism has a habitat and a niche. A habitat is all of the living and nonliving factors in the area where an organism lives. For example, the habitat of a frog
More informationEnergy, Producers, and Consumers. Lesson Overview. Lesson Overview. 4.1 Energy, Producers, and Consumers
4.1 Energy, Producers, and Consumers THINK ABOUT IT At the core of every organism s interaction with the environment is its need for energy to power life s processes. Where does energy in living systems
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 informationApplications in Biology
11 Applications in Biology In this chapter we make use of the techniques developed in the previous few chapters to examine some nonlinear systems that have been used as mathematical models for a variety
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 informationPopulation Ecology Density dependence, regulation and the Allee effect
2/22/15 Population Ecology Density dependence, regulation and the Allee effect ESRM 450 Wildlife Ecology and Conservation Wildlife Populations Groups of animals, all of the same species, that live together
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 informationC2 Differential Equations : Computational Modeling and Simulation Instructor: Linwei Wang
C2 Differential Equations 4040-849-03: Computational Modeling and Simulation Instructor: Linwei Wang Part IV Dynamic Systems Equilibrium: Stable or Unstable? Equilibrium is a state of a system which does
More informationChapter 9 Population Dynamics, Carrying Capacity, and Conservation Biology
Chapter 9 Population Dynamics, Carrying Capacity, and Conservation Biology 9-1 Population Dynamics & Carrying Capacity Populations change in response to enviromental stress or changes in evironmental conditions
More informationBUNDLE 9: ENERGY AND ECOLOGY Review
BUNDLE 9: ENERGY AND ECOLOGY Review 1. Describe Cellular Respiration, what happens, where does it happen, what type of organism does it take place in? What is the equation for respiration? Happens in the
More informationEcology Test Biology Honors
Do Not Write On Test Ecology Test Biology Honors Multiple Choice Identify the choice that best completes the statement or answers the question. 1. The study of the interaction of living organisms with
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 informationCh 5. Evolution, Biodiversity, and Population Ecology. Part 1: Foundations of Environmental Science
Ch 5 Evolution, Biodiversity, and Population Ecology Part 1: Foundations of Environmental Science PowerPoint Slides prepared by Jay Withgott and Heidi Marcum Copyright 2006 Pearson Education, Inc., publishing
More information1 R.V k V k 1 / I.k/ here; we ll stimulate the action potential another way.) Note that this further simplifies to. m 3 k h k.
1. The goal of this problem is to simulate a propagating action potential for the Hodgkin-Huxley model and to determine the propagation speed. From the class notes, the discrete version (i.e., after breaking
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 information11/10/13. How do populations and communities interact and change? Populations. What do you think? Do you agree or disagree? Do you agree or disagree?
Chapter Introduction Lesson 1 Populations Lesson 2 Changing Populations Lesson 3 Communities Chapter Wrap-Up How do populations and communities interact and change? What do you think? Before you begin,
More informationPopulation Ecology. Chapter 44
Population Ecology Chapter 44 Stages of Biology O Ecology is the interactions of organisms with other organisms and with their environments O These interactions occur in different hierarchies O The simplest
More informationPULSE-SEASONAL HARVESTING VIA NONLINEAR DELAY DIFFERENTIAL EQUATIONS AND APPLICATIONS IN FISHERY MANAGEMENT. Lev V. Idels
PULSE-SEASONAL HARVESTING VIA NONLINEAR DELAY DIFFERENTIAL EQUATIONS AND APPLICATIONS IN FISHERY MANAGEMENT Lev V. Idels University-College Professor Mathematics Department Malaspina University-College
More information2. What is a phytoplankton bloom and when does it generally occur in the North Atlantic?
Name: Date: Guiding Question: Seasonal Cycles: the North Atlantic Phytoplankton Bloom What are the factors that control the patterns/cycles of phytoplankton growth in the North Atlantic Ocean? Introduction
More informationMATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev
MATH 723 Mathematical Neuroscience Spring 28 Instructor: Georgi Medvedev 2 Lecture 2. Approximate systems. 2. Reduction of the HH model to a 2D system. The original HH system consists of 4 differential
More informationExam 2. Principles of Ecology. March 25, Name
Exam 2. Principles of Ecology. March 25, 2009. Name N t = N o λ t N t = N o e rt N t+1 = N t + r o N t (1-N t /K) N t = K/(1 + [(K N o )/N o ] * e rt ) dn/dt = rn(1-n/k) N captured and marked initially
More informationUnit 2 Ecology Study Guide. Niche Autotrophs Heterotrophs Decomposers Demography Dispersion
Vocabulary to know: Ecology Ecosystem Abiotic Factors Biotic Factors Communities Population Unit 2 Ecology Study Guide Niche Autotrophs Heterotrophs Decomposers Demography Dispersion Growth Rate Carrying
More informationGlobal Stability of a Computer Virus Model with Cure and Vertical Transmission
International Journal of Research Studies in Computer Science and Engineering (IJRSCSE) Volume 3, Issue 1, January 016, PP 16-4 ISSN 349-4840 (Print) & ISSN 349-4859 (Online) www.arcjournals.org Global
More informationTransmission in finite populations
Transmission in finite populations Juliet Pulliam, PhD Department of Biology and Emerging Pathogens Institute University of Florida and RAPIDD Program, DIEPS Fogarty International Center US National Institutes
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 informationLecture 11 : Simple Neuron Models. Dr Eileen Nugent
Lecture 11 : Simple Neuron Models Dr Eileen Nugent Reading List Nelson, Biological Physics, Chapter 12 Phillips, PBoC, Chapter 17 Gerstner, Neuronal Dynamics: from single neurons to networks and models
More informationInterspecific Competition
Interspecific Competition Intraspecific competition Classic logistic model Interspecific extension of densitydependence Individuals of other species may also have an effect on per capita birth & death
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationDYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS. Wei Feng
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 7 pp. 36 37 DYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS Wei Feng Mathematics and Statistics Department
More information