Kasetsart University Workshop. Mathematical modeling using calculus & differential equations concepts
|
|
- Marilyn Atkins
- 6 years ago
- Views:
Transcription
1 Kasetsart University Workshop Mathematical modeling using calculus & differential equations concepts Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA A copy of these slides is available at pardhan/kaset/workshop/
2 Warmup Consider the following models for the rate of change of some population P (t) (a) (b) (c) dp = 0.4P dp = 0.4P dp = 0.4P ( 1 P 200 ( 1 P 200 ) ) ( 1 50 P ) For each model, sketch a rough graph of P vs t for the following 3 initial values: P (0) = 20. P (0) = 100. P (0) = 250. Hint: It is not necessary to solve the problem by integration!
3 Warmup Recall a simple and powerful calculus concept: The sign of the derivative tells us a lot. If f (x) < 0, then f > 0, then f = 0, then f
4 Warmup Solution (a) P = 0.4P Sign of P (as function of P ) P= P The sign graph shows P if P = 0 P if P > 0 P if P < 0 Thus the solutions for P (0) = 20, P (0) = 100, and P (0) = 250 look like: P(t) t
5 (b) P = 0.4P ( 1 P ) 200 Sign of P P= P= The sign graph shows P if P = 0 or 200 P if 0 < P < 200 P if P < 0 or P > 200 P 400 P(t) 200 t
6 (c) P = 0.4P ( 1 P ) ( 1 50 ) 200 P Sign of P P= P= P= The sign graph shows 400 P(t) P if P = 50 or 200 P if 50 < P < 200 P if P < 0, 0 < P < 50 or P > 200 P 200 t
7 Illustration: Predator-prey modeling R(t) = population of rabbits (prey) F (t) = population of fox (predator) Exercise: Fill in reasonable forms for the missing terms (in the boxes) dr = rabbit births rabbit deaths df = fox births fox deaths What assumptions did you make? How did they affect your model?
8 Predator-prey modeling Lotka-Volterra model: Rabbit birth rate = constant Rabbit death rate fox population Fox death rate = constant Fox birth rate rabbit population This leads to models of the form: dr = a R b RF df = c F + d RF where a, b, c, d are positive constant parameters. One way to interpret the parameters: a = per-capita birth rate of rabbits b F = per-capita death rate of rabbits due to F (t) foxes c = per-capita death rate of foxes d R = per-capita birth rate of foxes due to R(t) rabbits Exercise: Find the equilibrium solutions (or, fixed points). Interpret their meaning.
9 Numerical example: dr df Predator-prey modeling = 0.08 R RF = 0.02 F RF Simulate using Maple: Effect of initial conditions (i) R(0) = 500, F (0) = 50 (ii) R(0) = 500, F (0) = 100 Effect of parameter values Carefully vary each parameter and explore the effect. What can you tell about the solution behavior?
10 Brief recap via a short quiz! 1. Suppose the population of some species is given by the model dp = (8 P )2 (12 P ) Find the equilibrium solutions, determine their stability, and sketch a few typical solution curves on a graph of P vs. t. 2. Shown below are different models of two interacting populations. Determine what type of interaction is occurring in each case: (i) predator-prey, or (ii) competing species, or (iii) cooperating species (a) dm = 0.012M MN dn = 0.003N MN (b) dm = 0.003M MN dn = 0.012N MN M, N denote species populations, and t denotes time.
11 Objective: Infectious disease modeling Predict how an infectious disease spreads through a population. Some examples of such diseases include: Swine flu, measles, TB, malaria, HIV, etc. Compartmental modeling strategy: OR
12 Model 1: SIS Assumptions: Infectious disease modeling Total population is constant. No births/deaths. At any time, a person is either infected or susceptible, but not both! Infected persons become susceptible again after they recover. Model: Exercise: Fill in reasonable forms for the missing terms (in the boxes) ds = rate of infection + rate of recovery di = + rate of infection rate of recovery Hint: Rate of infection depends on S(t) and I(t). Rate of recovery depends only on I(t). [why?]
13 Infectious disease modeling A simple form of the terms: Rate of infection = βsi Rate of recovery = γi } β and γ are constants Thus the model becomes ds = βsi + γi di = βsi γi Exercise: What are the units of β and γ? Does the model assume S + I is constant? Let S +I = N. Use this and reduce the model to a single differential equation in I. Does this model look familiar?! What is the long-term behavior of the solution? What do β and γ represent? β how contagious the disease is, and how much contact between S and I γ average time of recovery from disease
14 Infectious disease modeling Basic reproduction number: Let R 0 = βn γ This s a key parameter used in disease modeling. The value R 0 = 1 represents a bifurcation point in this model. The solution behavior changes at this point. Details of how/why: Using S + I = N we get: di = β(n I)I γi = βi(n I γ β ) The two equilibrium solutions are: I = 0 and I = N γ β Exercise: What does each equilibrium solution represent in real life? Explore the effect of R 0 using Maple simulations. Case1: Let N = 50, 000, β = , γ = 1/14 Case2: change β to You can use any initial conditions s.t. S(0) + I(0) = N
15 Model 2: SIR Assumptions: Infectious disease modeling Total population is constant. No births/deaths. At any time, a person is either infected, susceptible or recovered. Recovered persons cannot get the infection again. Model: Exercise: A simple SIR model can be made by rearranging terms in the SIS model. Here it is again for reference ds = βsi + γi di = βsi γi Rearrange terms and construct the SIR model ds = di = dr =
16 Infectious disease modeling Here is one way to do it: ds = βsi di = βsi γi dr = γi Compare with the SIS model: What do the parameters mean? How can we estimate their values? What are the equilibrium solutions? What is the basic reproduction ratio?
17 Infectious disease modeling Simulate an example using Maple: Let N = 50, 000, β = , γ = 1/14, S(0) = 45, 400, I(0) = 2, 100, R(0) = 2, 500 Thus, the equations become S (t) = S I I (t) = S I I 14 Exercise: R (t) = I 14 (1) Why is S + I + R = 0? What does that mean? (2) The sign of which derivative tells you whether the infection is spreading or diminishing? (3) There is an important critical value of S. When S is above this value the infection spreads, and when it is below, the infection decreases. Can you find this value? (4) If a vaccine exists, it is possible (according to this model) to vaccinate enough people to stop the spread of the disease. What is the minimum number of susceptible people that must be vaccinated to stop the spread?
18 A modeling project In 2003 the SARS epidemic broke out in Hong Kong. Data about the number of infected people per week is available for a 10-week period (see Science magazine, vol. 300, 2003) Determine which of the 2 models (SIS or SIR) can be applied to model this situation. Explain why. Construct the model. Simulate it using Maple and check whether it works. [E.g., let S(0) = 1, 000, 000, I(0) = 5, R(0) = 0]
19 A modeling project Some hints: * The given SARS data is a graph of I(t). * Try to estimate the value of γ first. * A rough estimate of I during the initial period can be obtained from the graph. * Can you estimate β if you know γ and I? If yes, how? If not, what else do you think is needed? * Keep γ fixed and try to vary β over some range. Use Maple simulations and compare your graph of I(t) with the given SARS data.
20 Extensions of previous models Predator-prey extensions: Competing species or cooperating species E.g. dx dy = 0.08 X XY = 0.02 Y XY Improved predator-prey models E.g. dx = 0.08 X dy ( 1 X = 0.02 Y XY 3-species hybrid models E.g. dx = 0.08 X dy = 0.05 Y dz ( 1 X ) XY ) XY 0.03 XZ ) XY 0.05 Y Z ( 1 Y = 0.02 Z XZ Y Z Q: Can you explain what interactions are happening here?
21 Extensions of previous models Infectious disease extensions: More complex population compartments E.g., SEIR (Susceptible, Exposed, Infected, Recovered); SEIRS; SIQRS; etc. Modeling the spread of ideas or opinions E.g. Public opinion of government policies S = no opinion I = for the policy Z= against the policy ds = β i(s + Z)I β z (S + I)Z + γ i I + γ z Z di = β i(s + Z)I γ i I dz = β z(s + I)Z γ z Z where β i, β z, γ i, γ z are constant parameters.
22 80 70 Iraq war opinion: no opinion war was wrong decision war was right decision 60 y (percent) t (months)
23 Review what we understand All models are wrong. But some are useful! A famous quote 1. What is mathematical modeling? 2. How do we construct mathematical models? 3. Suppose I want to model a predator-prey system in which the population is a function of position and time. E.g., R(x, y, t) = rabbit population density F (x, y, t) = fox population density How would that change our old model: dr = a R b RF df = c F + d RF
24 The Mathematical Contest in Modeling Annual event for undergraduate students worldwide. It is held every year around the 1st weekend of Feb. The contest consists of open-ended, real-life problems for which there is no known correct solution. Contest participants work in teams. Each team consists of 3 students, plus a faculty advisor. All work is submitted electronically. I strongly encourage Kasetsart students to participate next year! See their website Each year s prize winning entries are published in the UMAP Journal. Examples of prize winning entries from previous years can also be found by doing a Google search with keywords such as mathematical contest in modeling solutions. We will look at examples of some topics and prize winning entries. But, keep in mind, there is no single correct solution to these open-ended problems.
25 Modeling natural phenomena What are natural phenomena? Heating/cooling processes The motion of fluids, electrons, objects, etc. Chemical reactions Climate, weather and geological processes Conservation is a key concept in modeling nature: E.g. Conservation of mass Conservation of momentum Conservation of energy Conservation of (electric) charge When we write the conservation principle in mathematical form, we get a model.
26 Modeling natural phenomena An illustration: Fluid flow Consider an infinitesimal volume V, with surface S, in the flow field (a 2D example is shown in the sketch). Rate of change of mass inside V = d ρ dv V Rate of mass flux through the surface = ρu ds S d By conservation of mass: ρ dv = ρu ds V S ρ t dv = (ρu) dv V V ρ t = (ρu) OR ρ t + (ρu) = 0 This is the famous continuity equation of fluid flow.
27 Modeling natural phenomena Similarly, using conservation of momentum, together with constitutive relations, we get a very complete and accurate mathematical model for fluid flow ρ t + (ρu) = 0 ( ) u ρ t + u u = p + τ + f ρ = density, u = velocity vector, p = pressure, τ = stress tensor, f = body force vector. These are the famous Navier-Stokes equations. In 2D Cartesian coordinates, for incompressible fluids (e.g., water), with no body forces they simplify to u x + v y = 0 ( ) u ρ t + u u x + v u = P ( 2 ) y x + µ u x + 2 u 2 y 2 ρ ( v t + u v x + v v y ) = P y + µ ( 2 v x v y 2 )
28 A simple tsunami model What causes a tsunami? Vertical displacement of large mass of water in deep ocean. Typical characteristics Very long wave-length, small amplitude waves. Wave-length >> ocean depth e.g., wave-length = km, ocean depth = 2-3 km. Simple modeling approach Wave energy and motion is predominantly in 1 direction. So we can use a 1D, depth-averaged approximation of the Navier-Stokes equations. This is called the shallow water model.
29 A simple tsunami model Shallow water wave equations In 1D, they have the form h t + (hu) x = 0 (hu) + (hu gh 2 ) = body forces t x The unknowns are h(x, t) = water level, u(x, t) = velocity. Using these we can show that speed of tsunami wave is v(x) g b(x) where b(x) = depth of ocean at position x. Exercise: Compute the speed of tsunami waves for each case below (1) b(x) = 4300m (average depth of Pacific ocean) (2) b(x) = 3400m (average depth of Indian ocean) (3) b(x) = 11000m (deepest point in the Pacific) Conclusion?
30 A simple tsunami model Tsunamis in deep water travel very fast, but they have very small amplitude. This is captured very well by the shallow water model. Q: As the tsunami nears the shore, what happens according to the above model? This leads to what consequence? The amplitude A of tsunami waves can also be modeled similarly, which leads to the relation Exercise: A 1 b 1/4 Suppose a tsunami in deep ocean has amplitude 1 m when the water depth is 4000 m. Using the above relationship, find the wave amplitude when it is near the shore at water depths of about 2 m. For comparison, the 2011 earthquake-induced tsunami in Japan was about 7 m high when it hit the coast. Let s view a matlab simulation of the shallow water equations for a tsunami runup!
31 Pattern formation in nature Nature is full of fascinating patterns. Some examples include: Fractal patterns in plants, landscapes and crystals Fibonacci patterns in fruit and flowers Animal coat patterns Migration patterns in butterflies, whales and birds Sand dunes, fossils and geological patterns. Mathematics is also full of fascinating patterns! Is there some connection with nature? Let us explore some PDE models and the patterns in their solutions.
32 Pattern formation in nature Consider the following image Suppose the pattern represents the concentration of some chemical. For example, low concentrations correspond to lighter colors, and high concentrations to darker colors. Let u demote the local concentration of this chemical. Exercise: (1) What does the graph of u look like on the square patch? As a first step torwards this, consider a narrow horizontal strip of the square. Make a graph of the pattern (u vs. x) you might see on this horizontal strip. (2) Conjecture a formula for a function whose garph might have similar features. (3) Now return to the 2D case with the square patch. Sketch a qualitatively reasonable graph of u(x, y). (4) Conjecture a formula for a function u(x, y) whose garph might have some similar features.
33 Reaction-diffusion models We consider 2 chemicals in a 2D domain. Let u(t, x, y) = concentration of one chemical v(t, x, y) = concentration of other chemical Example of a typical R-D model [ u u 2 fv ] (u q) (u + q) u t = D 1 2 u + 1 ɛ v t = D 2 2 v + u v Here everything except u, v are constant parameters. This is called the Oregonator model. See typical solution patterns on the website pardhan/kaset/workshop/
34 PDE models Models and research results from various other application areas.
Final 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 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 informationBees and Flowers. Unit 1: Qualitative and Graphical Approaches
Bees and Flowers Often scientists use rate of change equations in their stu of population growth for one or more species. In this problem we stu systems of rate of change equations designed to inform us
More informationMathematical Modeling and Analysis of Infectious Disease Dynamics
Mathematical Modeling and Analysis of Infectious Disease Dynamics V. A. Bokil Department of Mathematics Oregon State University Corvallis, OR MTH 323: Mathematical Modeling May 22, 2017 V. A. Bokil (OSU-Math)
More informationDifferential Equations Handout A
Differential Equations Handout A Math 1b November 11, 2005 These are problems for the differential equations unit of the course, some topics of which are not covered in the main textbook. Some of these
More informationModeling with differential equations
Mathematical Modeling Lia Vas Modeling with differential equations When trying to predict the future value, one follows the following basic idea. Future value = present value + change. From this idea,
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 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 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 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 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 informationMathematical Analysis of Epidemiological Models: Introduction
Mathematical Analysis of Epidemiological Models: Introduction Jan Medlock Clemson University Department of Mathematical Sciences 8 February 2010 1. Introduction. The effectiveness of improved sanitation,
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 informationMathematical modelling and controlling the dynamics of infectious diseases
Mathematical modelling and controlling the dynamics of infectious diseases Musa Mammadov Centre for Informatics and Applied Optimisation Federation University Australia 25 August 2017, School of Science,
More informationFinal Exam Review (and Quiz #12)
Final Exam Review (and Quiz #1) Applied Calculus Math 13.3 Fall 11 Your final exam is scheduled for Saturday, December 17, at 9: a.m. in our usual room. You may bring an 1 11 inch, one-sided, sheet of
More informationDIFFERENCE AND DIFFERENTIAL EQUATIONS IN MATHEMATICAL MODELLING. J. Banasiak
DIFFERENCE AND DIFFERENTIAL EQUATIONS IN MATHEMATICAL MODELLING J. Banasiak ii Contents Basic ideas about mathematical modelling. Introduction: what is mathematical modelling?.........2 Simple difference
More informationModelling of the Hand-Foot-Mouth-Disease with the Carrier Population
Modelling of the Hand-Foot-Mouth-Disease with the Carrier Population Ruzhang Zhao, Lijun Yang Department of Mathematical Science, Tsinghua University, China. Corresponding author. Email: lyang@math.tsinghua.edu.cn,
More informationSection 6: Second Derivative and Concavity Second Derivative and Concavity
Chapter The Derivative Applied Calculus 1 Section 6: Second Derivative and Concavity Second Derivative and Concavity Graphically, a function is concave up if its graph is curved with the opening upward
More informationAE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 1...in the phrase computational fluid dynamics the word computational is simply an adjective to fluid dynamics.... -John D. Anderson 2 1 Equations of Fluid
More informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More information(e) Use Newton s method to find the x coordinate that satisfies this equation, and your graph in part (b) to show that this is an inflection point.
Chapter 6 Review problems 6.1 A strange function Consider the function (x 2 ) x. (a) Show that this function can be expressed as f(x) = e x ln(x2). (b) Use the spreadsheet, and a fine subdivision of the
More informationA Simple Ecological Model
Department of Mathematics Florida State University Florida State University, Tallahassee, FL, Jul 22, 2009 Ecology 1 The branch of biology that deals with the relations of organisms to one another and
More informationStability of SEIR Model of Infectious Diseases with Human Immunity
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 6 (2017), pp. 1811 1819 Research India Publications http://www.ripublication.com/gjpam.htm Stability of SEIR Model of Infectious
More informationMath 142-2, Homework 2
Math 142-2, Homework 2 Your name here April 7, 2014 Problem 35.3 Consider a species in which both no individuals live to three years old and only one-year olds reproduce. (a) Show that b 0 = 0, b 2 = 0,
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationMathematical Analysis of Epidemiological Models III
Intro Computing R Complex models Mathematical Analysis of Epidemiological Models III Jan Medlock Clemson University Department of Mathematical Sciences 27 July 29 Intro Computing R Complex models What
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
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 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 information1.1 GRAPHS AND LINEAR FUNCTIONS
MATHEMATICS EXTENSION 4 UNIT MATHEMATICS TOPIC 1: GRAPHS 1.1 GRAPHS AND LINEAR FUNCTIONS FUNCTIONS The concept of a function is already familiar to you. Since this concept is fundamental to mathematics,
More informationLecture: Wave-induced Momentum Fluxes: Radiation Stresses
Chapter 4 Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Here we derive the wave-induced depth-integrated momentum fluxes, otherwise known as the radiation stress tensor S. These are the 2nd-order
More informationSIR Epidemic Model with total Population size
Advances in Applied Mathematical Biosciences. ISSN 2248-9983 Volume 7, Number 1 (2016), pp. 33-39 International Research Publication House http://www.irphouse.com SIR Epidemic Model with total Population
More informationMath 2930 Worksheet Introduction to Differential Equations. What is a Differential Equation and what are Solutions?
Math 2930 Worksheet Introduction to Differential Equations Week 1 January 25, 2019 What is a Differential Equation and what are Solutions? A differential equation is an equation that relates an unknown
More informationThe SIR Disease Model Trajectories and MatLab
The SIR Disease Model Trajectories and MatLab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 17, 2013 Outline Reviewing the SIR
More informationMultivariable calculus: What is it good for?
: What is it good for? Department of Mathematics and Statistics, UMBC http://www.math.umbc.edu/ rouben/ January 2011 Purpose The purpose of these slides is to show that the material that one learns in
More informationThe death of an epidemic
LECTURE 2 Equilibrium Stability Analysis & Next Generation Method The death of an epidemic In SIR equations, let s divide equation for dx/dt by dz/ dt:!! dx/dz = - (β X Y/N)/(γY)!!! = - R 0 X/N Integrate
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 informationApplied Calculus. Review Problems for the Final Exam
Math135 Study Guide 1 Math 131/135/194, Fall 2004 Applied Calculus Daniel Kaplan, Macalester College Review Problems for the Final Exam Problem 1../DE/102b.tex Problem 3../DE/107.tex Consider the pair
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationModels of Infectious Disease Formal Demography Stanford Spring Workshop in Formal Demography May 2008
Models of Infectious Disease Formal Demography Stanford Spring Workshop in Formal Demography May 2008 James Holland Jones Department of Anthropology Stanford University May 3, 2008 1 Outline 1. Compartmental
More informationAPPM 2360 Lab #3: The Predator Prey Model
APPM 2360 Lab #3: The Predator Prey Model 1 Instructions Labs may be done in groups of 3 or less. One report must be turned in for each group and must be in PDF format. Labs must include each student s:
More informationThree Disguises of 1 x = e λx
Three Disguises of 1 x = e λx Chathuri Karunarathna Mudiyanselage Rabi K.C. Winfried Just Department of Mathematics, Ohio University Mathematical Biology and Dynamical Systems Seminar Ohio University November
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 6
Strauss PDEs 2e: Section 3 - Exercise Page of 6 Exercise Carefully derive the equation of a string in a medium in which the resistance is proportional to the velocity Solution There are two ways (among
More information5.3 Interpretations of the Definite Integral Student Notes
5. Interpretations of the Definite Integral Student Notes The Total Change Theorem: The integral of a rate of change is the total change: a b F This theorem is used in many applications. xdx Fb Fa Example
More informationModule 02 Control Systems Preliminaries, Intro to State Space
Module 02 Control Systems Preliminaries, Intro to State Space Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha August 28, 2017 Ahmad
More information1.2. Introduction to Modeling. P (t) = r P (t) (b) When r > 0 this is the exponential growth equation.
G. NAGY ODE January 9, 2018 1 1.2. Introduction to Modeling Section Objective(s): Review of Exponential Growth. The Logistic Population Model. Competing Species Model. Overview of Mathematical Models.
More 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 informationMathematical Epidemiology Lecture 1. Matylda Jabłońska-Sabuka
Lecture 1 Lappeenranta University of Technology Wrocław, Fall 2013 What is? Basic terminology Epidemiology is the subject that studies the spread of diseases in populations, and primarily the human populations.
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca
More informationCS224W: Analysis of Networks Jure Leskovec, Stanford University
Announcements: Please fill HW Survey Weekend Office Hours starting this weekend (Hangout only) Proposal: Can use 1 late period CS224W: Analysis of Networks Jure Leskovec, Stanford University http://cs224w.stanford.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 informationAARMS Homework Exercises
1 For the gamma distribution, AARMS Homework Exercises (a) Show that the mgf is M(t) = (1 βt) α for t < 1/β (b) Use the mgf to find the mean and variance of the gamma distribution 2 A well-known inequality
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 informationA Note on the Spread of Infectious Diseases. in a Large Susceptible Population
International Mathematical Forum, Vol. 7, 2012, no. 50, 2481-2492 A Note on the Spread of Infectious Diseases in a Large Susceptible Population B. Barnes Department of Mathematics Kwame Nkrumah University
More informationMathematical Models of Biological Systems
Mathematical Models of Biological Systems CH924 November 25, 2012 Course outline Syllabus: Week 6: Stability of first-order autonomous differential equation, genetic switch/clock, introduction to enzyme
More 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 information4.1 Modelling with Differential Equations
Chapter 4 Differential Equations The rate equations with which we began our study of calculus are called differential equations when we identify the rates of change that appear within them as derivatives
More informationSpecies 1 isocline. Species 2 isocline
1 Name BIOLOGY 150 Final Exam Winter Quarter 2002 Before starting please write your name on each page! Last name, then first name. You have tons of time. Take your time and read each question carefully
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 informationApproximation of epidemic models by diffusion processes and their statistical inferencedes
Approximation of epidemic models by diffusion processes and their statistical inferencedes Catherine Larédo 1,2 1 UR 341, MaIAGE, INRA, Jouy-en-Josas 2 UMR 7599, LPMA, Université Paris Diderot December
More informationChap. 3: Kinematics (2D) Recap: Kinematics (1D) 1. Vector Kinematics 2. Projectile Motion 3. Uniform Circular Motion 4.
Chap. 3: Kinematics (2D) Recap: Kinematics (1D) 1. Vector Kinematics 2. Projectile Motion 3. Uniform Circular Motion 4. Relative Velocity 1 Last, This and Next Weeks [Last Week] Chap. 1 and Chap. 2 [This
More informationAlgebra II. Slide 1 / 261. Slide 2 / 261. Slide 3 / 261. Linear, Exponential and Logarithmic Functions. Table of Contents
Slide 1 / 261 Algebra II Slide 2 / 261 Linear, Exponential and 2015-04-21 www.njctl.org Table of Contents click on the topic to go to that section Slide 3 / 261 Linear Functions Exponential Functions Properties
More informationSYMBIOTIC MODELS WITH AN SIR DISEASE
SYMBIOTIC MODELS WITH AN SIR DISEASE A Thesis Presented to the Faculty of California State Polytechnic University, Pomona In Partial Fulfillment Of the Requirements for the Degree Master of Science In
More informationGlobal Analysis of an SEIRS Model with Saturating Contact Rate 1
Applied Mathematical Sciences, Vol. 6, 2012, no. 80, 3991-4003 Global Analysis of an SEIRS Model with Saturating Contact Rate 1 Shulin Sun a, Cuihua Guo b, and Chengmin Li a a School of Mathematics and
More informationMath 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationMATHEMATICAL MODELS Vol. III - Mathematical Models in Epidemiology - M. G. Roberts, J. A. P. Heesterbeek
MATHEMATICAL MODELS I EPIDEMIOLOGY M. G. Roberts Institute of Information and Mathematical Sciences, Massey University, Auckland, ew Zealand J. A. P. Heesterbeek Faculty of Veterinary Medicine, Utrecht
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationFixed Point Analysis of Kermack Mckendrick SIR Model
Kalpa Publications in Computing Volume, 17, Pages 13 19 ICRISET17. International Conference on Research and Innovations in Science, Engineering &Technology. Selected Papers in Computing Fixed Point Analysis
More informationSystems of Differential Equations
WWW Problems and Solutions 5.1 Chapter 5 Sstems of Differential Equations Section 5.1 First-Order Sstems www Problem 1. (From Scalar ODEs to Sstems). Solve each linear ODE; construct and solve an equivalent
More informationA New Mathematical Approach for. Rabies Endemy
Applied Mathematical Sciences, Vol. 8, 2014, no. 2, 59-67 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.39525 A New Mathematical Approach for Rabies Endemy Elif Demirci Ankara University
More informationEpidemics in Complex Networks and Phase Transitions
Master M2 Sciences de la Matière ENS de Lyon 2015-2016 Phase Transitions and Critical Phenomena Epidemics in Complex Networks and Phase Transitions Jordan Cambe January 13, 2016 Abstract Spreading phenomena
More informationIdeas from Vector Calculus Kurt Bryan
Ideas from Vector Calculus Kurt Bryan Most of the facts I state below are for functions of two or three variables, but with noted exceptions all are true for functions of n variables..1 Tangent Line Approximation
More informationU(s) = 0 + 0s 3 + 0s 2 + 0s + 125
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5143 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 2 Andres Rainiero Hernandez Coronado September 6, 2017 The objective of this homework is to test your understanding
More informationSee animations and interactive applets of some of these at. Fall_2009/Math123/Notes
MA123, Chapter 7 Word Problems (pp. 125-153) Chapter s Goal: In this chapter we study the two main types of word problems in Calculus. Optimization Problems. i.e., max - min problems Related Rates See
More informationMath 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
More information1,3. f x x f x x. Lim. Lim. Lim. Lim Lim. y 13x b b 10 b So the equation of the tangent line is y 13x
1.5 Topics: The Derivative lutions 1. Use the limit definition of derivative (the one with x in it) to find f x given f x 4x 5x 6 4 x x 5 x x 6 4x 5x 6 f x x f x f x x0 x x0 x xx x x x x x 4 5 6 4 5 6
More informationMCMC 2: Lecture 3 SIR models - more topics. Phil O Neill Theo Kypraios School of Mathematical Sciences University of Nottingham
MCMC 2: Lecture 3 SIR models - more topics Phil O Neill Theo Kypraios School of Mathematical Sciences University of Nottingham Contents 1. What can be estimated? 2. Reparameterisation 3. Marginalisation
More informationHomework #4 Solutions
MAT 303 Spring 03 Problems Section.: 0,, Section.:, 6,, Section.3:,, 0,, 30 Homework # Solutions..0. Suppose that the fish population P(t) in a lake is attacked by a disease at time t = 0, with the result
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
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 informationFundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering
Fundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering Module No # 05 Lecture No # 25 Mass and Energy Conservation Cartesian Co-ordinates Welcome
More informationMath 310: Applied Differential Equations Homework 2 Prof. Ricciardi October 8, DUE: October 25, 2010
Math 310: Applied Differential Equations Homework 2 Prof. Ricciardi October 8, 2010 DUE: October 25, 2010 1. Complete Laboratory 5, numbers 4 and 7 only. 2. Find a synchronous solution of the form A cos(ωt)+b
More informationMITOCW MITRES18_005S10_DiffEqnsGrowth_300k_512kb-mp4
MITOCW MITRES18_005S10_DiffEqnsGrowth_300k_512kb-mp4 GILBERT STRANG: OK, today is about differential equations. That's where calculus really is applied. And these will be equations that describe growth.
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationAverage rates of change May be used to estimate the derivative at a point
Derivatives Big Ideas Rule of Four: Numerically, Graphically, Analytically, and Verbally Average rate of Change: Difference Quotient: y x f( a+ h) f( a) f( a) f( a h) f( a+ h) f( a h) h h h Average rates
More informationEXERCISES FOR SECTION 1.1
2 CHAPTER 1 FIRST-ORDER DIFFERENTIAL EQUATIONS EXERCISES FOR SECTION 1.1 1. Note that dy/ = 0 if and only if y = 3. Therefore, the constant function y(t) = 3for all t is the only equilibrium solution.
More informationMOMENTS OF INERTIA FOR AREAS, RADIUS OF GYRATION OF AN AREA, & MOMENTS OF INTERTIA BY INTEGRATION
MOMENTS OF INERTIA FOR AREAS, RADIUS OF GYRATION OF AN AREA, & MOMENTS OF INTERTIA BY INTEGRATION Today s Objectives: Students will be able to: a) Define the moments of inertia (MoI) for an area. b) Determine
More informationToday. Logistic equation in many contexts
Today Logistic equation in many contexts Classic example of the power of mathematics - one unifying description for many apparently unrelated phenomena. Rates of change that are proportional to two things
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationMath 2930 Worksheet Introduction to Differential Equations
Math 2930 Worksheet Introduction to Differential Equations Week 1 August 24, 2017 Question 1. Is the function y = 1 + t a solution to the differential equation How about the function y = 1 + 2t? How about
More informationAustralian Journal of Basic and Applied Sciences. Effect of Personal Hygiene Campaign on the Transmission Model of Hepatitis A
Australian Journal of Basic and Applied Sciences, 9(13) Special 15, Pages: 67-73 ISSN:1991-8178 Australian Journal of Basic and Applied Sciences Journal home page: wwwajbaswebcom Effect of Personal Hygiene
More informationC.3 First-Order Linear Differential Equations
A34 APPENDIX C Differential Equations C.3 First-Order Linear Differential Equations Solve first-order linear differential equations. Use first-order linear differential equations to model and solve real-life
More informationEPIDEMIC MODELS I REPRODUCTION NUMBERS AND FINAL SIZE RELATIONS FRED BRAUER
EPIDEMIC MODELS I REPRODUCTION NUMBERS AND FINAL SIZE RELATIONS FRED BRAUER 1 THE SIR MODEL Start with simple SIR epidemic model with initial conditions S = βsi I = (βs α)i, S() = S, I() = I, S + I = N
More informationGrade 7/8 Math Circles Fall October 23/24/25 Qualitative Graph Analysis
Faculty of Mathematics Waterloo, Ontario N2L 3G1 Centre for Education in Mathematics and Computing Grade 7/8 Math Circles Fall 2018 - October 23/24/25 Qualitative Graph Analysis What do we mean when we
More informationArchive of Calculus IV Questions Noel Brady Department of Mathematics University of Oklahoma
Archive of Calculus IV Questions Noel Brady Department of Mathematics University of Oklahoma This is an archive of past Calculus IV exam questions. You should first attempt the questions without looking
More informationMatrix Method of Structural Analysis Prof. Amit Shaw Department of Civil Engineering Indian Institute of Technology, Kharagpur
Matrix Method of Structural Analysis Prof. Amit Shaw Department of Civil Engineering Indian Institute of Technology, Kharagpur Lecture 01 Introduction Hello everyone, welcome to the online course on Matrix
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 informationMATH 280 Multivariate Calculus Fall Integrating a vector field over a surface
MATH 280 Multivariate Calculus Fall 2011 Definition Integrating a vector field over a surface We are given a vector field F in space and an oriented surface in the domain of F as shown in the figure below
More informationFinal Exam December 20, 2011
Final Exam December 20, 2011 Math 420 - Ordinary Differential Equations No credit will be given for answers without mathematical or logical justification. Simplify answers as much as possible. Leave solutions
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More information