Mathematical Models for Erythropoiesis
|
|
- Lilian Pitts
- 6 years ago
- Views:
Transcription
1 Mathematical Models for Erythropoiesis Joseph M. Mahaffy Department of Mathematical and Statistical Sciences San Diego State University January 2003 p.1/37
2 Collaborators Jacques Bélair - Université de Montréal Mike Mackey - McGill University NSF-REU Students Sam Polk Roland Roeder p.2/37
3 Outline Discuss the Physiology Develop the Age-Structured Model Reduce the Model to Delay Equations Bifurcation Analysis Compare to Examples Rabbit with Induced Auto-Immune Hemolytic Anemia Human Subject following a Phlebotomy Summary p.3/37
4 Hematopoiesis p.4/37
5 Erythropoiesis p.5/37
6 Features of Erythropoiesis BFU-E and CFU-E differentiate and proliferate in response to EPO Maturation requires about 6 days EPO accelerates maturation Lack of EPO causes apoptosis Cell divisions every 8 hours for about 4 days Reticulocytes do not divide - increase hemoglobin Erythrocytes lose nucleus - live 120 days Macrophages degrade RBCs EPO released near kidneys with half-life of 6 hours p.6/37
7 Age-Structured Model p.7/37
8 Detailed Model p.8/37
9 Active Degradation of RBCs RBCs age - Cell membrane breaks down Membrane marked with antibodies Macrophages destroy least pliable cells Model assumes constant supply macrophages Saturated consumption of Erythrocytes - Satiated predator Constant flux of RBCs being destroyed p.9/37
10 Constant Flux Boundary Condition Let Q be rate of removal of erythrocytes Erythrocytes lost are Q t Mean Value Theorem - average number RBCs m(ξ, ν F (ξ)) for ξ (t, t + t) Balance law Q t = W t m(ξ, ν F (ξ)) [ν F (t + t) ν F (t)]m(ξ, ν F (ξ)) As t 0, Q = [W ν F (t)]m(t, ν F (t)) p.10/37
11 Constant Flux - Diagram p.11/37
12 Simplifying Assumptions Several simplifying assumptions allow reduction of the age-structured model to delay differential equations Assume that V (E) = W = 1. Assume the birth rate β satisfies: β(µ, E) = { β, µ < µ1, 0, µ µ 1, Assume that γ is constant. p.12/37
13 Reduced PDEs The model satisfies the partial differential equations: p t + p µ = β(µ)p m t + m ν = γm p.13/37
14 Reduced PDEs The model satisfies the partial differential equations: p t + p µ = β(µ)p m t + m ν = γm with the boundary conditions: p(t, 0) = S 0 (E) and p(t, µ F ) = m(t, 0) (1 ν F (t))m(t, ν F (t)) = Q p.13/37
15 Negative Control The negative feedback by EPO satisfies the equation: E = a 1 + KM r ke where the total mature erythrocyte population is M(t) = νf (t) 0 m(t, ν)dν p.14/37
16 Method of Characteristics Applied to a Simplified Model The method of characteristics can be used to simplify the partial differential equations given above. Let P (s) = p(t(s), µ(s)), dp ds = p t dt ds + p µ dµ ds = β(µ(s), E(t(s)))P (s), which has the solution, P (s) = p(t, µ) = P (0) exp provided [ s 0 ] β(µ(r), E(t(r)))dr, dt ds = 1 and dµ ds = V (E(t(s))) = 1 p.15/37
17 Characteristics Diagram p.16/37
18 Evaluating p(t, µ) and m(t, ν) To find p at µ F, p(t, µ F ) = p(t 0, 0)exp [ µf 0 ] β(r)dr = p(t µ F, 0)e βµ 1 = e βµ 1 S 0 (E(t µ F )) Similar use of the characteristics gives m(t, ν) = m(t ν, 0)e γν, p.17/37
19 Finding M(t) M(t) = = = νf (t) 0 νf (t) 0 νf (t) 0 m(t ν, 0)e γν dν p(t µ F ν, 0)e γν dν e βµ 1 S 0 (E(t µ F ν))e γν dν, = e γ(t µ F ) e βµ 1 t µf t µ F ν F (t) S 0 (E(w))e γw dw, p.18/37
20 Leibnitz s Rule We apply Leibnitz s rule for differentiating an integral: M(t) = γe γ(t µ F ) e βµ 1 t µf t µ F ν F (t) S 0 (E(w))e γw dw +e βµ 1 [S 0 (E(t µ F )) ] S 0 (E(t µ F ν F (t)))e γν F (t) (1 ν F (t)) = γm(t) + e βµ 1 S 0 (E(t µ F )) Q, p.19/37
21 Model with Delays After reduction of PDEs, the state variables become total mature erythrocytes, M, EPO, E, and age of RBCs, ν F. dm(t) dt de(t) dt dν F (t) dt = e βµ 1 S 0 (E(t T 1 )) γm(t) Q = f(m(t)) ke(t) = 1 Qe βµ 1 eγν F (t) S 0 (E(t T 1 ν F (t)) where T 1 = µ F. This is a state-dependent delay differential equation. p.20/37
22 Properties of the Model State-dependent delay model has a unique positive equilibrium p.21/37
23 Properties of the Model State-dependent delay model has a unique positive equilibrium Delay T 1 accounts for maturing time p.21/37
24 Properties of the Model State-dependent delay model has a unique positive equilibrium Delay T 1 accounts for maturing time State-dependent delay in equation for ν F, the varying age of maturation p.21/37
25 Properties of the Model State-dependent delay model has a unique positive equilibrium Delay T 1 accounts for maturing time State-dependent delay in equation for ν F, the varying age of maturation The ν F differential equation is uncoupled from the differential equations for M and E p.21/37
26 Properties of the Model State-dependent delay model has a unique positive equilibrium Delay T 1 accounts for maturing time State-dependent delay in equation for ν F, the varying age of maturation The ν F differential equation is uncoupled from the differential equations for M and E Stability determined by equations for M and E p.21/37
27 Linear Analysis of the Model Linearizing about the unique equilibrium ( M, Ē, ν F ), M(t) = e βµ 1 S 0(Ē)E(t T 1) γm(t) E(t) = f ( M)M(t) ke(t) ν F (t) = 1 Ē E(t T 1 ν F ) γν F (t) p.22/37
28 Linear Analysis of the Model Linearizing about the unique equilibrium ( M, Ē, ν F ), M(t) = e βµ 1 S 0(Ē)E(t T 1) γm(t) E(t) = f ( M)M(t) ke(t) ν F (t) = 1 Ē E(t T 1 ν F ) γν F (t) The characteristic equation is given by (λ + γ) [ (λ + γ)(λ + k) + Āe λt 1] = 0, where A e βµ 1 S 0(Ē)f ( M) > 0. p.22/37
29 Linear Analysis of the Model Linearizing about the unique equilibrium ( M, Ē, ν F ), M(t) = e βµ 1 S 0(Ē)E(t T 1) γm(t) E(t) = f ( M)M(t) ke(t) ν F (t) = 1 Ē E(t T 1 ν F ) γν F (t) The characteristic equation is given by (λ + γ) [ (λ + γ)(λ + k) + 1] Āe λt = 0, where A e βµ 1 S 0(Ē)f ( M) > 0. One solution is λ = γ, which shows the stability of the ν F equation. p.22/37
30 Analysis of Characteristic Equation Remains to analyze (λ + γ)(λ + k) = Ae λt 1 p.23/37
31 Analysis of Characteristic Equation Remains to analyze (λ + γ)(λ + k) = Ae λt 1 A Hopf bifurcation occurs when λ = iω solves the characteristic equation p.23/37
32 Analysis of Characteristic Equation Remains to analyze (λ + γ)(λ + k) = Ae λt 1 A Hopf bifurcation occurs when λ = iω solves the characteristic equation From complex variables, we match the magnitudes and arguments: (iω + γ)(iω + k) = A ( ) ω ( ω ) Θ(ω) arctan + arctan γ k = π ωt 1 p.23/37
33 Analysis of Characteristic Equation Remains to analyze (λ + γ)(λ + k) = Ae λt 1 A Hopf bifurcation occurs when λ = iω solves the characteristic equation From complex variables, we match the magnitudes and arguments: (iω + γ)(iω + k) = A ( ) ω ( ω ) Θ(ω) arctan + arctan γ k = π ωt 1 Solve for ω by varying parameters such as γ p.23/37
34 Hopf - Argument Principle Link to Var. Vel. Anal. p.24/37
35 Variable Velocity of Maturing Physiologically, erythrocytes mature more rapidly with increasing EPO p.25/37
36 Variable Velocity of Maturing Physiologically, erythrocytes mature more rapidly with increasing EPO Velocity of maturation, V (E), is a non-decreasing function, e.g. V (E) = κ 1E κ 2 + E p.25/37
37 Variable Velocity of Maturing Physiologically, erythrocytes mature more rapidly with increasing EPO Velocity of maturation, V (E), is a non-decreasing function, e.g. V (E) = κ 1E κ 2 + E Method of characteristics leaves a threshold-type functional equation rather than simpler delay equation p.25/37
38 Variable Velocity of Maturing Physiologically, erythrocytes mature more rapidly with increasing EPO Velocity of maturation, V (E), is a non-decreasing function, e.g. V (E) = κ 1E κ 2 + E Method of characteristics leaves a threshold-type functional equation rather than simpler delay equation Linear analysis of the age-structured model relatively simple p.25/37
39 Linear Analysis with Variable Velocity The modified characteristic equation for the threshold-type functional equation becomes (λ + γ)(λ + k) e βµ 1 f ( M)(V 1 + V 2 e λµ F ) = 0 where V 2 V 1 V (Ē)S 0(Ē) V (Ē) V (Ē)S 0(Ē) V (Ē)S 0(Ē) V (Ē) p.26/37
40 Linear Analysis (continued) This has the form (λ + γ)(λ + k) + α 1 + (α 2 α 1 )e λµ F = 0 with α 1 = V 1 and α 2 positive. p.27/37
41 Linear Analysis (continued) This has the form (λ + γ)(λ + k) + α 1 + (α 2 α 1 )e λµ F = 0 with α 1 = V 1 and α 2 positive. Since α 2 = A from previous characteristic equation, the α 1 shifts our geometric diagram above to the right with a smaller radius circle. p.27/37
42 Linear Analysis (continued) This has the form (λ + γ)(λ + k) + α 1 + (α 2 α 1 )e λµ F = 0 with α 1 = V 1 and α 2 positive. Since α 2 = A from previous characteristic equation, the α 1 shifts our geometric diagram above to the right with a smaller radius circle. This can be readily seen to stabilize the model. p.27/37
43 Hopf - Variable Velocity Link to Anal. Link to Hopf p.28/37
44 Auto-Immune Induced Anemia Rabbits were injected with antibodies to their Red Blood Cells p.29/37
45 Identify Parameters Physiological parameters found for rabbit Extensive literature search Some gaps using other species studied Difficult, but essential process Increasing the destruction of Erythrocytes, γ, causes a Hopf Bifurcation Model most sensitive to parameters γ and µ F Variable velocity of maturation stabilizes the model p.30/37
46 Simulation p.31/37
47 Phlebotomy (Blood Donation) Normal blood donation is about 8% of blood O 2 sensors near kidneys probably sense concentration, not M(t) Blood donation loses erythrocytes and plasma, no concentration change Plasma recovers quickly p.32/37
48 Model for Phlebotomy Define the hemoglobin concentration h(t) = H M(t) M(t) + ρ(t) The plasma function chosen to fit data ρ(t) = α [ 1 + (β 1 t 0.08)e β 2t ] The age-structured model remains the same except E(t) = f(h(t)) ke(t) Fit model to data of Maeda et al and Wadsworth p.33/37
49 Model a Phlebotomy p.34/37
50 Model Compared to Data p.35/37
51 Summary The mathematical model provides a good example of how age-structured models are related to models with delays. The modeling of satiation for destruction of erythrocytes could prove valuable in other population models. The model for erythropoiesis can be fit to existing data and can hopefully provide insight into the study of some hematopoietic diseases. It is unlikely to aid in the study of normal individuals for improved blood donation schemes. p.36/37
52 Summary (cont) Our numerical simulations and analytical study show how a variable velocity of maturation stabilizes the model. This implies that plasticity in the precursor compartment may be an important evolutionary adaptation. Current studies have identified the most significant parameters in the model, which could give insight to the likely causes of the disease states and possible therapeutic approaches. New studies examine thrombopoietic systems using a multi-compartment model to account for size structures of megakaryocytes. p.37/37
A Structured Population Model for Erythropoiesis
A Structured Population Model for Erythropoiesis Doris H. Fuertinger 1, Franz Kappel 1, Peter Kotanko 2, Nathan W. Levine 2, Stephan Thijssen 2 1 Institute for Mathematics and Scientific Computation, University
More informationAdding Self-Renewal in Committed Erythroid Progenitors Improves the Biological Relevance of a Mathematical Model of Erythropoiesis
Adding Self-Renewal in Committed Erythroid Progenitors Improves the Biological Relevance of a Mathematical Model of Erythropoiesis Fabien Crauste a, Laurent Pujo-Menjouet a, Stéphane Génieys a, Clément
More informationAge-Structured Modeling of Hematopoiesis
Age-Structured Modeling of Hematopoiesis Joseph M. Mahaffy CRM-269 May 1999 Department of Mathematical Sciences, San Diego State University, San Diego, CA 92182. This work was supported in part by NSF
More information8. MATHEMATICAL MODELS OF HEMATOPOIETIC CELL REPLICATION AND CONTROL
8. MATHEMATICAL MODELS OF HEMATOPOIETIC CELL REPLICATION AND CONTROL Michael C. Mackey 8.1 INTRODUCTION The biological sciences offer an abundance of examples in which the dynamics are dependent on the
More informationAdding self-renewal in committed erythroid progenitors improves the biological relevance of a mathematical model of erythropoiesis
Journal of Theoretical Biology 25 (28) 322 338 www.elsevier.com/locate/yjtbi Adding self-renewal in committed erythroid progenitors improves the biological relevance of a mathematical model of erythropoiesis
More informationAsynchronous oscillations due to antigenic variation in Malaria Pf
Asynchronous oscillations due to antigenic variation in Malaria Pf Jonathan L. Mitchell and Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX SIAM LS, Pittsburgh, 21 Outline
More informationModeling Red Blood Cell and Iron Dynamics in Patients Undergoing Periodic EPO and Iron Treatments
Modeling Red Blood Cell and Iron Dynamics in Patients Undergoing Periodic EPO and Iron Treatments H T Banks 1, Karen M Bliss 2, D F Kapraun 1, and Hien Tran 1 1 Center for Research in Scientific Computation
More informationHEMATOLOGY/ HEMATOPOIESIS. Introduction
HEMATOLOGY/ HEMATOPOIESIS Introduction HEMATOLOGY Introduction Study of blood & its components Window of rest of body BLOOD Raison d etre Delivery of nutrients Oxygen Food Vitamins Removal of wastes Carbon
More informationc Copyright 2011 by Karen M. Bliss All Rights Reserved
ABSTRACT BLISS, KAREN M. Modeling of Red Blood Cell Dynamics in Patients with Chronic Kidney Disease. (Under the direction of H.T. Banks and H.T. Tran.) Kidneys are the main site of production of the hormone
More informationStability and Hopf bifurcation in a mathematical model of pluripotent stem cell dynamics
Nonlinear Analysis: Real World Applications 6 (2005) 651 670 www.elsevier.com/locate/na Stability and Hopf bifurcation in a mathematical model of pluripotent stem cell dynamics Mostafa Adimy a, Fabien
More informationBLOOD. YONG-MEI CHEN ( 陈咏梅 ) Dept. of Anatomy, Histology & Embryology Peking Union Medical College Tel: address:
BLOOD YONG-MEI CHEN ( 陈咏梅 ) Dept. of Anatomy, Histology & Embryology Peking Union Medical College Tel:69156461 E-mail address: pumc_he@yahoo.com.cn Content Composition & functions of plasma Structure &
More informationDelay Differential Equations and Autonomous Oscillations in Hematopoietic Stem Cell Dynamics Modeling
Math. Model. Nat. Phenom. Vol. 7, No. 6, 212, pp. 1 22 DOI: 1.151/mmnp/212761 Delay Differential Equations and Autonomous Oscillations in Hematopoietic Stem Cell Dynamics Modeling M. Adimy 1,2, F. Crauste
More informationGene Regulatory Network and Functional Analysis: Hopf Bifurcation causes by molecular movement
Gene Regulatory Network and Functional Analysis: Hopf Bifurcation causes by molecular movement Mariya Ptashnyk Division of Mathematics, University of Dundee joined work with Mark Chaplain and Marc Sturrock
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationDELAY DIFFERENTIAL EQUATION MODELING FOR THE ANALYSIS OF HEMATOLOGICAL DISEASES BY ALTERNATIVE G-CSF TREATMENT
DELAY DIFFERENTIAL EQUATION MODELING FOR THE ANALYSIS OF HEMATOLOGICAL DISEASES BY ALTERNATIVE G-CSF TREATMENT A THESIS REPORT Submitted by BALAMURALITHARAN S. Under the guidance of Dr. S.RAJASEKARAN in
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 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 informationBLOOD. PUMC Dept. of Anat. Hist. & Embry. 陈咏梅 YONG-MEI CHEN Dept. of Anatomy, Histology & Embryology Peking Union Medical College
BLOOD PUMC Dept. of Anat. Hist. & Embry. 陈咏梅 YONG-MEI CHEN Dept. of Anatomy, Histology & Embryology Peking Union Medical College Content Composition & functions of plasma PUMC Dept. of Anat. Hist. & Embry.
More informationSWINGING UP A PENDULUM BY ENERGY CONTROL
Paper presented at IFAC 13th World Congress, San Francisco, California, 1996 SWINGING UP A PENDULUM BY ENERGY CONTROL K. J. Åström and K. Furuta Department of Automatic Control Lund Institute of Technology,
More informationSabrina Rogg Doris H. Fuertinger Stefan Volkwein Franz Kappel Peter Kotanko
Universität Konstanz Optimal EPO dosing in hemodialysis patients using a non-linear model predictive control approach Sabrina Rogg Doris H. Fuertinger Stefan Volkwein Franz Kappel Peter Kotanko Konstanzer
More informationarxiv: v2 [q-bio.to] 27 Jan 2017
Normal and pathological dynamics of platelets in humans arxiv:168.286v2 [q-bio.to] 27 Jan 217 Gabriel P. Langlois Morgan Craig Antony R. Humphries Michael C. Mackey Joseph M. Mahaffy Jacques Bélair Thibault
More informationStoichiometry Of Tumor Dynamics: Models And Analysis
Stoichiometry Of Tumor Dynamics: Models And Analysis Yang Kuang Department of Mathematics and Statistics Arizona State University Supported by NSF grant DMS-0077790 Y. Kuang, J. Nagy and J. Elser: Disc.
More informationDynamic Feedback and the Design of Closed-loop Drug Delivery Systems
Dynamic Feedback and the Design of Closed-loop Drug Delivery Systems John Milton 1,2, Sue Ann Campbell 2,3,4 and Jacques Bélair 2,4,5 1) Department of Neurology, The University of Chicago Hospitals, MC
More informationSupplement to TB in Canadian First Nations at the turn-of-the twentieth century
Supplement to TB in Canadian First Nations at the turn-of-the twentieth century S. F. Ackley, Fengchen Liu, Travis C. Porco, Caitlin S. Pepperell Equations Definitions S, L, T I, T N, and R give the numbers
More informationOutline. Calculus for the Life Sciences. Crow Predation on Whelks. Introduction. Lecture Notes Optimization. Joseph M. Mahaffy,
Calculus for the Life Sciences Lecture Notes Optimization Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationIntroduction to Mathematical Physiology I - Biochemical Reactions
Introduction to Mathematical Physiology I - Biochemical Reactions J. P. Keener Mathematics Department Math Physiology p.1/28 Introduction The Dilemma of Modern Biology The amount of data being collected
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationVariable mass oscillator
Variable mass oscillator José Flores Departamento de Física J. J. Giambiagi, Universidad de Buenos Aires, Argentina and Facultad de Ingeniería y Ciencias de la Universidad Favaloro, Buenos Aires, Argentina
More informationReaction-Diffusion Models and Bifurcation Theory Lecture 11: Bifurcation in delay differential equations
Reaction-Diffusion Models and Bifurcation Theory Lecture 11: Bifurcation in delay differential equations Junping Shi College of William and Mary, USA Collaborators/Support Shanshan Chen (PhD, 2012; Harbin
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
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 informationPeriodic Solutions in a Mathematical Model for the Treatment of Chronic Myelogenous Leukemia
Math. Model. Nat. Phenom. Vol. 7, No.,, pp. 35-44 DOI:.5/mmnp/7 Periodic Solutions in a Mathematical Model for the Treatment of Chronic Myelogenous Leukemia A. Halanay Department of Mathematics I, Politehnica
More informationClassical Description of NMR Parameters: The Bloch Equations
Classical Description of NMR Parameters: The Bloch Equations Pascale Legault Département de Biochimie Université de Montréal 1 Outline 1) Classical Behavior of Magnetic Nuclei: The Bloch Equation 2) Precession
More informationMathematical Modeling in the Kidney. Characterizing a Long-Looped Nephron
: Characterizing a Long-Looped Nephron Quinton Neville Department of MSCS, St. Olaf College October 13, 2016 Northfield Undergraduate Math Symposium Kidney and Nephron Tubuloglomerular Feedback (TGF) System
More informationSpace-time kinetics. Alain Hébert. Institut de génie nucléaire École Polytechnique de Montréal.
Space-time kinetics Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE613: Week 1 Space-time kinetics 1/24 Content (week 1) 1 Space-time kinetics equations
More informationMathematical physiology. 1.1 Derive a suitably scaled form of the Michaelis-Menten model for the reaction E + P,, λ = k 2
Problem sheet 1. 1.1 Derive a suitably scaled form of the Michaelis-Menten model for the reaction S + E and show that it depends on the parameters k 1 k 1 C k 2 E + P K = k 1 + k 2 k 1 S 0 λ = k 2 k 1
More informationModels of Infectious Disease Formal Demography Stanford Summer Short Course James Holland Jones, Instructor. August 15, 2005
Models of Infectious Disease Formal Demography Stanford Summer Short Course James Holland Jones, Instructor August 15, 2005 1 Outline 1. Compartmental Thinking 2. Simple Epidemic (a) Epidemic Curve 1:
More informationPARAMETER ESTIMATION IN EPIDEMIC MODELS: SIMPLIFIED FORMULAS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 19, Number 4, Winter 211 PARAMETER ESTIMATION IN EPIDEMIC MODELS: SIMPLIFIED FORMULAS Dedicated to Herb Freedman on the occasion of his seventieth birthday
More informationName: Date: Answer: Answer:
Name: Date: 5 6 7 8 9 0 Scoring Guide: Scoring Guide: 5 6 7 8 9 0 5 6 7 8 9 0 Scoring Guide: Scoring Guide: 5 Scoring Guide: 6 7 8 9 0 5 6 7 8 9 50 Scoring Guide: 5 Scoring Guide: Standard(s):..0.F,...F,..0.D,...D,..0.C,...C,..0.E,...E,.5.0.F,.5..F
More informationarxiv: v1 [math.ds] 6 Aug 2013
Regions of Stability for a Linear Differential Equation with Two Rationally Dependent Delays Joseph M. Mahaffy Nonlinear Dynamical Systems Group Department of Mathematics San Diego State University San
More informationClassical Description of NMR Parameters: The Bloch Equations
Classical Description of NMR Parameters: The Bloch Equations Pascale Legault Département de Biochimie Université de Montréal 1 Outline 1) Classical Behavior of Magnetic Nuclei: The Bloch Equation 2) Precession
More informationCloud modeling and Lotka-Voltera Spence Lunderman December 9, 2017 University of Arizona Department of Mathematics
1 Low dimensional cloud modeling Cloud modeling and Lotka-Voltera Spence Lunderman December 9, 17 University of Arizona Department of Mathematics Cloud modeling is a very difficult task: one must understand
More informationS-Duality for D3-Brane in NS-NS and R-R Backgrounds
S-Duality for D3-Brane in NS-NS and R-R Backgrounds Chen-Te Ma Collaborator: Pei-Ming Ho National Taiwan University arxiv:1311.3393 [hep-th] January 17, 2014 Reference P. -M. Ho and Y. Matsuo, M5 from
More informationTutorial on Control and State Constrained Optimal Control Pro. Control Problems and Applications Part 3 : Pure State Constraints
Tutorial on Control and State Constrained Optimal Control Problems and Applications Part 3 : Pure State Constraints University of Münster Institute of Computational and Applied Mathematics SADCO Summer
More informationSpring Force and Power
Lecture 14 Chapter 9 Physics I Spring Force and Power Course website: http://faculty.uml.edu/andriy_danylov/teaching/physicsi IN THIS CHAPTER, you will learn how to solve problems using two new concepts:
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 informationModelling of cell choice between differentiation and apoptosis on the basis of intracellular and extracellular regulations and stochasticity
Modelling of cell choice between differentiation and apoptosis on the basis of intracellular and extracellular regulations and stochasticity M. Banerjee 1, V. Volpert 2 arxiv:163.8712v1 [q-bio.cb] 29 Mar
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 informationarxiv: v1 [math.ap] 16 Apr 2009
A Mathematical Study of the Hematopoiesis Process with Applications to Chronic Myelogenous Leukemia arxiv:94.2493v1 [math.ap] 16 Apr 29 Mostafa Adimy, Fabien Crauste Year 24 and Shigui Ruan Laboratoire
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationAlgorithmische Parameterschätzung und optimale Versuchsplanung am Beispiel der Erythropoese. Masterarbeit
Algorithmische Parameterschätzung und optimale Versuchsplanung am Beispiel der Erythropoese An der Fakultät für Mathematik der Otto-von-Guericke-Universität Magdeburg zur Erlangung des akademischen Grades
More informationCellular lifespan based pharmacodynamic analysis of erythropoiesis
University of Iowa Iowa Research Online Theses and Dissertations Summer 2009 Cellular lifespan based pharmacodynamic analysis of erythropoiesis Kevin Jay Freise University of Iowa Copyright 2009 Kevin
More informationConnecting Caputo-Fabrizio Fractional Derivatives Without Singular Kernel and Proportional Derivatives
Connecting Caputo-Fabrizio Fractional Derivatives Without Singular Kernel and Proportional Derivatives Concordia College, Moorhead, Minnesota USA 7 July 2018 Intro Caputo and Fabrizio (2015, 2017) introduced
More informationOutline. Calculus for the Life Sciences. Malthusian Growth Model. Introduction. Lecture Notes Separable Differential Equations
Calculus for the Life Sciences Lecture Notes Outline 1 Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationA Parameter Identification Technique for Structured Population Equations Modeling Erythropoiesis in Dialysis Patients
, 23-25 October, 2013, San Francisco, USA A Parameter Identification Technique for Structured Population Equations Modeling Erythropoiesis in Dialysis Patients Doris H. Fuertinger, Franz Kappel Abstract
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 informationNumerical qualitative analysis of a large-scale model for measles spread
Numerical qualitative analysis of a large-scale model for measles spread Hossein Zivari-Piran Department of Mathematics and Statistics York University (joint work with Jane Heffernan) p./9 Outline Periodic
More informationDouble Hopf Bifurcation of a Hematopoietic System with State Feedback Control
Journal of Mathematics Research; Vol. 1, No. 4; August 218 ISSN 1916-9795 E-ISSN 1916-989 Published by Canadian Center of Science and Education Double Hopf Bifurcation of a Hematopoietic System with State
More informationPolynomial Chaos Approach for Maxwell s Equations in Dispersive Media
Polynomial Chaos Approach for Maxwell s Equations in Dispersive Media Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 15, 2013 Prof. Gibson (OSU) PC-FDTD
More informationDynamics of Erythroid Progenitors and Erythroleukemia
Math. Model. Nat. Phenom. Vol. 4, No. 3, 2009, pp. 210-232 DOI: 10.1051/mmnp/20094309 Dynamics of Erythroid Progenitors and Erythroleukemia N. Bessonov a,b, F. Crauste b1, I. Demin b, V. Volpert b a Institute
More informationI. INTRODUCTION AND HISTORICAL PERSPECTIVE
I. INTRODUCTION AND HISTORICAL PERSPECTIVE A. Failures of Classical Physics At the end of the 19th century, physics was described via two main approaches. Matter was described by Newton s laws while radiation
More informationThe interaction of light and matter
Outline The interaction of light and matter Denise Krol (Atom Optics) Photon physics 014 Lecture February 14, 014 1 / 3 Elementary processes Elementary processes 1 Elementary processes Einstein relations
More informationStability and Hopf Bifurcation Analysis of the Delay Logistic Equation
1 Stability and Hopf Bifurcation Analysis of the Delay Logistic Equation Milind M Rao # and Preetish K L * # Department of Electrical Engineering, IIT Madras, Chennai-600036, India * Department of Mechanical
More informationOSCILLATIONS IN A MATURATION MODEL OF BLOOD CELL PRODUCTION
SIAM J. APPL. MATH. Vol. 66, No. 6, pp. 227 248 c 26 Society for Industrial and Applied Mathematics OSCILLATIONS IN A MATURATION MODEL OF BLOOD CELL PRODUCTION IVANA DROBNJAK, A. C. FOWLER, AND MICHAEL
More information7. Non-LTE basic concepts
7. Non-LTE basic concepts LTE vs NLTE occupation numbers rate equation transition probabilities: collisional and radiative examples: hot stars, A supergiants 1 Equilibrium: LTE vs NLTE LTE each volume
More informationLINEAR DISPERSIVE WAVES
LINEAR DISPERSIVE WAVES Nathaniel Egwu December 16, 2009 Outline 1 INTRODUCTION Dispersion Relations Definition of Dispersive Waves 2 SOLUTION AND ASYMPTOTIC ANALYSIS The Beam Equation General Solution
More informationsections June 11, 2009
sections 3.2-3.5 June 11, 2009 Population growth/decay When we model population growth, the simplest model is the exponential (or Malthusian) model. Basic ideas: P = P(t) = population size as a function
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 informationLink Models for Packet Switching
Link Models for Packet Switching To begin our study of the performance of communications networks, we will study a model of a single link in a message switched network. The important feature of this model
More informationExercise 1 (Formula for connection 1-forms) Using the first structure equation, show that
1 Stokes s Theorem Let D R 2 be a connected compact smooth domain, so that D is a smooth embedded circle. Given a smooth function f : D R, define fdx dy fdxdy, D where the left-hand side is the integral
More informationBrownian Motion: Fokker-Planck Equation
Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential
More informationStability and bifurcation of a simple neural network with multiple time delays.
Fields Institute Communications Volume, 999 Stability and bifurcation of a simple neural network with multiple time delays. Sue Ann Campbell Department of Applied Mathematics University of Waterloo Waterloo
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationExam Thermodynamics 12 April 2018
1 Exam Thermodynamics 12 April 2018 Please, hand in your answers to problems 1, 2, 3 and 4 on separate sheets. Put your name and student number on each sheet. The examination time is 12:30 until 15:30.
More informationAP Physics 2012 Practice Quiz 4, Conduction & Electric Fields
Name: Class: Date: ID: A AP Physics 01 Practice Quiz 4, Conduction & Electric Fields Multiple Choice Identify the choice that best completes the statement or answers the question. 1. ( points) A repelling
More information1 Isotropic Covariance Functions
1 Isotropic Covariance Functions Let {Z(s)} be a Gaussian process on, ie, a collection of jointly normal random variables Z(s) associated with n-dimensional locations s The joint distribution of {Z(s)}
More informationGA A22443 STUDY OF H MODE THRESHOLD CONDITIONS IN DIII D
GA A443 STUDY OF H MODE THRESHOLD CONDITIONS IN DIII D by R.J. GROEBNER, T.N. CARLSTROM, K.H. BURRELL, S. CODA, E.J. DOYLE, P. GOHIL, K.W. KIM, Q. PENG, R. MAINGI, R.A. MOYER, C.L. RETTIG, T.L. RHODES,
More informationAtherosclerosis Initiation Modeled as an Inflammatory Process
Math. Model. Nat. Phenom. Vol. 2, No. 2, 2007, pp. 126-141 Atherosclerosis Initiation Modeled as an Inflammatory Process N. El Khatib 1, S. Génieys and V. Volpert Université Lyon 1, Institut Camille Jordan,
More informationA comparison of delayed SIR and SEIR epidemic models
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 2, 181 190 181 A comparison of delayed SIR and SEIR epidemic models Abdelilah Kaddar a, Abdelhadi Abta b, Hamad Talibi Alaoui b a Université
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 informationA Lévy-Fokker-Planck equation: entropies and convergence to equilibrium
1/ 22 A Lévy-Fokker-Planck equation: entropies and convergence to equilibrium I. Gentil CEREMADE, Université Paris-Dauphine International Conference on stochastic Analysis and Applications Hammamet, Tunisia,
More informationFast neural network simulations with population density methods
Fast neural network simulations with population density methods Duane Q. Nykamp a,1 Daniel Tranchina b,a,c,2 a Courant Institute of Mathematical Science b Department of Biology c Center for Neural Science
More informationOpen Quantum Systems and Markov Processes II
Open Quantum Systems and Markov Processes II Theory of Quantum Optics (QIC 895) Sascha Agne sascha.agne@uwaterloo.ca July 20, 2015 Outline 1 1. Introduction to open quantum systems and master equations
More informationOutline. Formic Acid
Outline Calculus for the Life Sciences Lecture Notes and Functions Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu 1 2 Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences
More informationx R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ).
Continuation of equilibria Problem Parameter-dependent ODE ẋ = f(x, λ), x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x,
More informationCalculus for the Life Sciences
Calculus for the Life Sciences Lecture Notes and Functions Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationDelay systems Batz-sur-Mer October
Delay systems Batz-sur-Mer October 17-21 2011 Thomas Erneux Université Libre de Bruxelles Optique Nonlinéaire Théorique T. Erneux, Applied Delay Differential Equations, Springer (2009) PLAN 1. Familiar
More informationUniversity of Bristol - Explore Bristol Research. Early version, also known as pre-print
Erzgraber, H, Krauskopf, B, & Lenstra, D (2004) Compound laser modes of mutually delay-coupled lasers : bifurcation analysis of the locking region Early version, also known as pre-print Link to publication
More informationCell Division and Reproduction
Cell Division and Reproduction What do you think this picture shows? If you guessed that it s a picture of two cells, you are right. In fact, the picture shows human cancer cells, and they are nearing
More informationPost Publication Errata April 27, 2011
1 Applied Stochastic Processes and Control for Jump-Diffusions: Modeling, Analysis, and Computation Floyd B. Hanson, University of Illinois at Chicago, Chicago, Illinois Post Publication Errata April 27,
More informationOutline. Math Numerical Analysis. Intermediate Value Theorem. Lecture Notes Zeros and Roots. Joseph M. Mahaffy,
Outline Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationOscillation Onset In. Neural Delayed Feedback
Oscillation Onset In Neural Delayed Feedback Andre Longtin Complex Systems Group and Center for Nonlinear Studies Theoretical Division B213, Los Alamos National Laboratory Los Alamos, NM 87545 Abstract
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 informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation
High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y
More informationNumerical bifurcation analysis of delay differential equations
Numerical bifurcation analysis of delay differential equations Dirk Roose Dept. of Computer Science K.U.Leuven Dirk.Roose@cs.kuleuven.be Acknowledgements Many thanks to Koen Engelborghs Tatyana Luzyanina
More informationQuantum Light-Matter Interactions
Quantum Light-Matter Interactions QIC 895: Theory of Quantum Optics David Layden June 8, 2015 Outline Background Review Jaynes-Cummings Model Vacuum Rabi Oscillations, Collapse & Revival Spontaneous Emission
More informationThe Fractional-order SIR and SIRS Epidemic Models with Variable Population Size
Math. Sci. Lett. 2, No. 3, 195-200 (2013) 195 Mathematical Sciences Letters An International Journal http://dx.doi.org/10.12785/msl/020308 The Fractional-order SIR and SIRS Epidemic Models with Variable
More informationBIBO STABILITY AND ASYMPTOTIC STABILITY
BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to
More information9231 FURTHER MATHEMATICS
UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Advanced Level MARK SCHEME for the October/November 2011 question paper for the guidance of teachers 9231 FURTHER MATHEMATICS 9231/21 Paper 2, maximum
More information