MATH 221 Biocalculus II Project 5 Nonlinear Systems of Difference Equations. BIO: To explore some modifications of the Nicholson-Bailey Model
|
|
- Melinda Stevens
- 5 years ago
- Views:
Transcription
1 Goals MATH 221 Biocalculus II Project 5 Nonlinear Systems of Difference Equations MATH: To analyze nonlinear systems of difference equations BIO: To explore the Nicholson-Bailey Model BIO: To explore some modifications of the Nicholson-Bailey Model COMP: Implement nonlinear difference equations in EXCEL dynamically Computational Tools: EXCEL (Make sure the Data Analysis Add-In loaded.) and Maple as needed Create the Nicholson Bailey Spread Sheet The Nicholson Bailey Model is given by N t+1 = bn t e apt = f(n t,p t ) P t+1 = cn t (1 e apt )=g(n t,p t ), where N t is the size of the host population in the tth generation, P t is the size of the parasitoid population in the tth generation, b and c are growth parameters, and a is parameter that characterizes the fraction of the the host population that escapes parasitism. (Specifically, the factor e apt is the fraction of the N t that escapes parasitism.) We denote the initial populations of the host by N 0 and the parasitoid by P Label the current worksheet, Sheet1, Home. 2. Create Sliders for the parameters a, b, andc, and the initial population values N 0 and P 0 in the model. 3 Initial Number of Hosts N0 Slider from 0 to 20, incremental change 1, linked to D3 4 Initial Number of Parasitoids P0 Slider from 0 to 20, incremental change 1, linked to D4 5 Host Escape Rate a =E5/100 Slider from 0 to 100, incremental change 1, linked to E5 6 Host Growth Parameter b =E6/10 Slider from 0 to 100, incremental change 1, linked to E6 7 Parasitoid Growth Parameter c =E7/10 Slider from 0 to 100, incremental change 1, linked to E7 Note: Cells D3 and D4 will be filled with slider values. 3. Click on Sheet2, and change its name to Data. 4. Label cell A4 T, cell B4 N, cell C4 P. 5. Fill the A column from A5 to A205 with the integers 0, 1,2,..., Enter =Home!$D$3 in B5 and =Home!D4 in C5. 7. In B6, enter =Home!$D$6*B5*EXP(-Home!$D$5*C5) and in C6, enter =Home!$D$7*B5*(1 - EXP(-Home!$D$5*C5)) 8. Copy the formulas in B6 and C6 down their respective columns through row 205. At this point, the times series data has been generated. 9. In the Home worksheet, create a phase plane diagram of P versus N. Label the axes and give the graph a title. Use an XY scatter plot without connecting the points. (Does this figure seem helpful?) 10. By copying the phase plane diagram, create a new phase plane diagram in which the points are connected by smoothed lines. (Does this figure seem helpful?) 1
2 11. In the Home worksheet, create a chart containing plots of P versus T and N versus T. Make the horizontal axis range from 0 to 50 and the vertical axis range from 0 to 100. Again, use XY scatter plots and connect the points with smoothed lines. (Does this plot seem helpful?) 12. Play with your parameter sliders and look at this third plot to gain some insight into the behavior of this model. Will the populations approach equilibrium values? 13. Find the equilibria of the Nicholson Bailey Model. Let ( ˆN, ˆP ) denote the nontrivial equilibrium of the system. Enter the expression for ˆN in terms of the cells containing the corresponding parameter values into cell L4 and the expression for ˆP into cell L5. Label these cells by entering appropriate titles in cells H4 and H5, respectively. 14. We need to compute the Jacobian of this system evaluated at and its eigenvalues to determine whether the nontrivial equilibrium is locally stable or unstable. Calculuate the matrix: J(N,P) = [ F N G N F P G P ] Now evaluate this matrix at ( ˆN, ˆP ). What conditions are on the parameter b so that ˆN is a positive quantity? What is the biological significance of this nontrivial equilibrium. Explain. Add the nontrivial equilibrium point ot the phase plane plots. Enter the Jacobian matrix in terms of the cells corresponding to appropriate parameter values and into cells B12, C12, B13, and C13. (You will either need to calculuate the Jacobian by hand or using a computer algebra system. 15. For an equilibrium point of system of two nonlinear difference equations, if the magnitudes of both eigenvalues of the Jacobian matrix are less than one, the equilibrium is locally stable. If at least one eigenvalue has magnitude greater than 1, the equilibrium is unstable. The spreadsheet will determine stability by directly calculating the eigevalues and by using the Jury Conditions below. We begin this process by placing the determinant of the Jacobian in cell F12 and the trace of the Jacobian in cell F13. Name cell F12 Det and cell F13 Trace. 16. Excel does not compute in terms of complex numbers, but we need to be able to work with complex eigenvalues. To do this, we perform the following calculations that return the real and imaginary parts to solutions of quadratic equations separately. Enter the following in K21: =SQRT((Trace^2-4*Det+ABS(Trace^2-4*Det))/2) Enter the following in K22: =SQRT((-(Trace^2-4*Det)+ABS(Trace^2-4*Det))/2) Name K21 alpha and K22 beta. 17. Enter Lambda1 and Lambda2 into H12 and H13, respectively. Make sure the Data Analysis Add-In has been installed so that the following commands will work. Enter the value of the eigenvalue Lambda1 in cell I12: =COMPLEX(ROUND((Trace+alpha)/2,4),ROUND(beta/2,4)) Enter the value of eigenvalue Lambda2 in cell I13: =COMPLEX(ROUND((Trace-alpha)/2,4),ROUND(-beta/2,4)) 2
3 Enter the magnitude of these eigenvalues in cells J12 and J13. The formula for the norm of Lamdba1 is given by: =IMABS(COMPLEX((Trace+alpha)/2,beta/2)) By further playing with several different parameter values, what do you predict about the stability of the nontrivial equilibrium? 18. (Jury Conditions.) For a system of two nonlinear difference equations, there are necessary and sufficient conditions that indicate when both eigenvalues of the Jacobian are less than one in magnitude. The conditions are the following: If the characteristic equation for J = J( ˆN, ˆP )is λ 2 Trace(J)λ +Det(J) =0, then both roots of the equation (the eigenvalues) will have magnitude less than one if and only if 2 > 1+Det(J) > Trace(J). Implement this test in Excel as follows: Enter 1 + Det(J) in N6, Trace(J), in N7, =1 + Det in O6, and =ABS(Trace) in O7. Enter Jury Test in N9, Stability in N10, =IF(AND(2 > $O$6,$O$6> $O$7), SATISFIED, FAILED ) in O9, and =IF($O$9= SATISFIED, LOCALLY STABLE, UNSTABLE ) in O10. What you observe here should be consistent with the observed magnitudes of the eigenvalues of the Jacobian. 19. Calculate Trace(J) anddet(j) in general and show that Det(J) is greater than 1 for the values of b satisfying the condition above. (Note that you can accomplish this by showing that the function S(b) =b 1 b ln(b) < 0, on the appropriate domain for b. WhatisS(1) and what is S (b)?) Draw a conclusion about the stability of ( ˆN, ˆP ). 20. Determine the stability of the trivial equilibrium as well. 21. Use the spreadsheet to observe how the Nicholson Bailey Model behaves with the following parameter choices. Describe the dynamical behavior in biological terms in each case. Print each time series graph. (a) N 0 = 15, P 0 =8,a =0.02, b =1.5, c =1.5 (b) N 0 = 15, P 0 =0,a =0.02, b =1.5, c =1.5 (c) N 0 = 15, P 0 =8,a =0.02, b =0.5, c =1.5 (d) N 0 = 15, P 0 =8,a =0.068, b =2,c =1 Negative Binomial Model With a few minor changes, we can implement the spreadsheet for the Negative Binomial model: ( N t+1 = bn t 1+ ap ) k = f(n t,p t ) k ( P t+1 = cn t (1 1+ ap ) ) k = g(n t,p t ), k where all the parameters are the same as in the previous model, with the addition of the binomial power k. In this model, ( ) 1+ ap k k serves as the fraction of the host populaton that escapes parasitism. You will be able to observe how this model may improve upon the Nicholson Bailey model. 3
4 1. Save a copy of your original spreadsheet as Negbinomial. 2. We first add the new parameter k to the Home worksheet in the spreadsheet: 8 binomial power k = E8/10 Slider from 0 to 100, incremental change 1, linked to E8 3. In the Data Worksheet, you must update the cells from rows 6 to 205 with the new formula. 4. After the data values have been updated with the new formula, experiment with different parameter values to begin to observethe behavior of this model. Modify the axes on your three graphs as necessary to obtain meaningful figures. 5. Determine the equilibria for the system and update cells L4 and L5 accordingly. 6. Compute the Jacobian matrix J = J(N,P) and evaluate J at the nontrivial equilibrium ( ˆN, ˆP ). Update the cells B12, C12, B13, C13 accordingly. 7. Using the new spreadsheet investigate the model for each of the parameter value combinations provided above for the Nicholson Bailey model. Take k =0.7 andthenk =0.5. Describe the behavior in each case and compare to the Nicholson Bailey Model. Print each graph. Make a prediction about the stability behavior of the nontrivial equilbrium in general. Can you prove your assertion? A Density Dependent Variation We will now look at another variation of the Nicholson-Bailey model in which we assume that the host population grows to some limiting density. The model becomes: N t+1 = N t exp(r(1 N t /K) ap t ) (1) P t+1 = N t (1 e apt ), (2) where r>0 is a growth constant and K is the carrying capacity. The quantity q = ˆN/K is indicates to what extent the equilibrium population ˆN is depressed by the presence of parasitoids. 1. Save another copy of the Nicholson Bailey file as HostParDep. 2. Update the parameter sliders as follows: 3 Initial Number of Hosts N0 Slider from 0 to 20, incremental change 1, linked to D3 4 Initial Number of Parasitoids P0 Slider from 0 to 20, incremental change 1, linked to D4 5 Host Escape Rate a =E5/100 Slider from 0 to 100, incremental change 1, linked to E5 6 Host Growth Parameter b =E6/10 Slider from 0 to 100, incremental change 1, linked to E6 7 Parasitoid Growth Parameter c =E7/10 Slider from 0 to 100, incremental change 1, linked to E7 8 Growth Parameter r =E8/10 Slider from 0 to 100, incremental change 1, linked to E8 9 Population Depression Factor q =E9/10 Slider from 0 to 100, incremental change 1, linked to E9 Note: we keep b = c = 1 throughout this spreadsheet. Simply leave existing sliders for these parameters set to 1. The nontrivial equilibrium can be expressed as: ˆP = r (1 q) (3) a Update cells L4 and L5 accordingly. ˆN = ˆP 1 e a ˆP (4) 4
5 3. The carrying capacity, K, will be calculuated using the formula K = q ˆN. Input this quantity into cell L3 as: =$L$4/$D$9 Label this cell by entering an appropriate title in cell H3. 4. In the Data Worksheet, you must update the cells from rows 6 to 205 with the new formula. 5. Compute the Jacobian matrix J = J(N,P) and evaluate J at the nontrivial equilibrium ( ˆN, ˆP ). Update the cells B12, C12, B13, C13 accordingly. We will investigate this model in the NP phase plane. Use the new spreadsheet with the following parameter choices. Determine the nontrivial equilibrium ( ˆN, ˆP )andk in each case. Report the moduli of the two eigenvalues of the Jacobian matrix at the equilibrium point. Print each graph of the NP phase plane and the plots of the two populations. Describe the behavior of the dynamics and the stability in each case. Assume that a =0.2 andq =0.4 in each case. (a) r =0.5 (b) r =2 (c) r =2.2 (d) r =2.65 References [1] J.R. Beddington, C.A. Free, and J.H. Lawton, Dynamic complexity in predator-prey models framed in difference equations, Nature 255, 58-60, [2] L. Edelstein-Keshet, Mathematical Models in Biology, Classics in Applied Mathematics 46, SIAM, [3] K.J. Griffiths, The importance of conicidence in the functional and numerical responses of two parasites of the European pine sawfly, Neodiprion sertifer, Canadian Entomologist , [4] R.M. May, Host-parasitoid systems in patchy environments: a phenomenological model, Journal of Animal Ecology , [5] A.J. Nicholson and V.A. Bailey, The balance of animal populations, Proceedings of the Zoological Society of London ,
EXCELLING WITH BIOLOGICAL MODELS FROM THE CLASSROOM T0 RESEARCH
EXCELLING WITH BIOLOGICAL MODELS FROM THE CLASSROOM T0 RESEARCH Timothy D. Comar Benedictine University Department of Mathematics 5700 College Road Lisle, IL 60532 tcomar@ben.edu Introduction Computer
More informationComputer simulation of radioactive decay
Computer simulation of radioactive decay y now you should have worked your way through the introduction to Maple, as well as the introduction to data analysis using Excel Now we will explore radioactive
More informationModeling Prey-Predator Populations
Modeling Prey-Predator Populations Alison Pool and Lydia Silva December 13, 2006 Alison Pool and Lydia Silva () Modeling Prey-Predator Populations December 13, 2006 1 / 25 1 Introduction 1 Our Populations
More informationRichter Scale and Logarithms
activity 7.1 Richter Scale and Logarithms In this activity, you will investigate earthquake data and explore the Richter scale as a measure of the intensity of an earthquake. You will consider how numbers
More informationAn area chart emphasizes the trend of each value over time. An area chart also shows the relationship of parts to a whole.
Excel 2003 Creating a Chart Introduction Page 1 By the end of this lesson, learners should be able to: Identify the parts of a chart Identify different types of charts Create an Embedded Chart Create a
More informationGeology Geomath Computer Lab Quadratics and Settling Velocities
Geology 351 - Geomath Computer Lab Quadratics and Settling Velocities In Chapter 3 of Mathematics: A simple tool for geologists, Waltham takes us through a brief review of quadratic equations and their
More informationMATH 3MB3 final sample questions
MATH 3MB3 final sample questions Special Instructions: Casio FX-991 MS or MS Plus calculator allowed, no other external aids Stability condition for 2D discrete-time systems: stable if T < 1 + < 2 Stability
More information1. The Basic X-Y Scatter Plot
1. The Basic X-Y Scatter Plot EXCEL offers a wide range of plots; however, this discussion will be restricted to generating XY scatter plots in various formats. The easiest way to begin is to highlight
More informationCanimals. borrowed, with thanks, from Malaspina University College/Kwantlen University College
Canimals borrowed, with thanks, from Malaspina University College/Kwantlen University College http://commons.wikimedia.org/wiki/file:ursus_maritimus_steve_amstrup.jpg Purpose Investigate the rate heat
More informationConformational Analysis of n-butane
Conformational Analysis of n-butane In this exercise you will calculate the Molecular Mechanics (MM) single point energy of butane in various conformations with respect to internal rotation around the
More informationLinear Motion with Constant Acceleration
Linear Motion 1 Linear Motion with Constant Acceleration Overview: First you will attempt to walk backward with a constant acceleration, monitoring your motion with the ultrasonic motion detector. Then
More informationModule 2A Turning Multivariable Models into Interactive Animated Simulations
Module 2A Turning Multivariable Models into Interactive Animated Simulations Using tools available in Excel, we will turn a multivariable model into an interactive animated simulation. Projectile motion,
More informationReview of Dynamic complexity in predator-prey models framed in difference equations
Review of Dynamic complexity in predator-prey models framed in difference eqations J. Robert Bchanan November 10, 005 Millersville University of Pennsylvania email: Robert.Bchanan@millersville.ed Review
More informationModeling Prey and Predator Populations
Modeling Prey and Predator Populations Alison Pool and Lydia Silva December 15, 2006 Abstract In this document, we will explore the modeling of two populations based on their relation to one another. Specifically
More informationIntroduction to Spectroscopy: Analysis of Copper Ore
Introduction to Spectroscopy: Analysis of Copper Ore Using a Buret and Volumetric Flask: 2.06 ml of solution delivered 2.47 ml of solution delivered 50.00 ml Volumetric Flask Reading a buret: Burets are
More informationIntroduction to Spectroscopy: Analysis of Copper Ore
Introduction to Spectroscopy: Analysis of Copper Ore Using a Buret and Volumetric Flask: 2.06 ml of solution 2.47 ml of solution 50.00 ml delivered delivered Volumetric Flask Reading a buret: Burets are
More informationCalculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y
Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point
More informationα m ! m or v T v T v T α m mass
FALLING OBJECTS (WHAT TO TURN IN AND HOW TO DO SO) In the real world, because of air resistance, objects do not fall indefinitely with constant acceleration. One way to see this is by comparing the fall
More informationChapter 9 Ingredients of Multivariable Change: Models, Graphs, Rates
Chapter 9 Ingredients of Multivariable Change: Models, Graphs, Rates 9.1 Multivariable Functions and Contour Graphs Although Excel can easily draw 3-dimensional surfaces, they are often difficult to mathematically
More informationI. Pre-Lab Introduction
I. Pre-Lab Introduction Please complete the following pages before the lab by filling in the requested items. A. Atomic notation: Atoms are composed of a nucleus containing neutrons and protons surrounded
More informationTitrator 3.0 Tutorial: Calcite precipitation
Titrator 3.0 Tutorial: Calcite precipitation November 2008 Steve Cabaniss A. Introduction This brief tutorial is intended to acquaint you with some of the features of the program Titrator. It assumes that
More informationDynamics of a plant-herbivore model
Journal of Biological Dynamics Vol., No., Month-Month x, 3 Dynamics of a plant-herbivore model Yun Kang, Dieter Armbruster and Yang Kuang (Received Month x; revised Month x; in final form Month x) We formulate
More informationLab #10 Atomic Radius Rubric o Missing 1 out of 4 o Missing 2 out of 4 o Missing 3 out of 4
Name: Date: Chemistry ~ Ms. Hart Class: Anions or Cations 4.7 Relationships Among Elements Lab #10 Background Information The periodic table is a wonderful source of information about all of the elements
More informationExperiment: Oscillations of a Mass on a Spring
Physics NYC F17 Objective: Theory: Experiment: Oscillations of a Mass on a Spring A: to verify Hooke s law for a spring and measure its elasticity constant. B: to check the relationship between the period
More informationLab 4: The Simple Pendulum
Page 1 Technical Math II Lab 4: Simple Pendulum Lab 4: The Simple Pendulum Purpose: To investigate the relationship between the length of a simple pendulum and the time it takes to complete a full swing.
More informationLocal Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)
Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two
More informationDISCRETE RANDOM VARIABLES EXCEL LAB #3
DISCRETE RANDOM VARIABLES EXCEL LAB #3 ECON/BUSN 180: Quantitative Methods for Economics and Business Department of Economics and Business Lake Forest College Lake Forest, IL 60045 Copyright, 2011 Overview
More informationa. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules.
Lab #1 - Free Vibration Name: Date: Section / Group: Procedure Steps (from lab manual): a. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules. b. Locate the various springs and
More informationAdvanced Forecast. For MAX TM. Users Manual
Advanced Forecast For MAX TM Users Manual www.maxtoolkit.com Revised: June 24, 2014 Contents Purpose:... 3 Installation... 3 Requirements:... 3 Installer:... 3 Setup: spreadsheet... 4 Setup: External Forecast
More informationSME 864 Mark Urban-Lurain
SME 864 Mark Urban-Lurain 1 Import data from non-excel sources Probe Software Web sites Other sources Organize data Structure file for analysis Clean values Analyze Summarize Statistics Graph 2 Get files
More informationKey words and phrases. Bifurcation, Difference Equations, Fixed Points, Predator - Prey System, Stability.
ISO 9001:008 Certified Volume, Issue, March 013 Dynamical Behavior in a Discrete Prey- Predator Interactions M.ReniSagaya Raj 1, A.George Maria Selvam, R.Janagaraj 3.and D.Pushparajan 4 1,,3 Sacred Heart
More informationBoyle s Law: A Multivariable Model and Interactive Animated Simulation
Boyle s Law: A Multivariable Model and Interactive Animated Simulation Using tools available in Excel, we will turn a multivariable model into an interactive animated simulation. Projectile motion, Boyle's
More informationMath 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations
Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function
More information+ i. cos(t) + 2 sin(t) + c 2.
MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,
More information18-Dec-12 PHYS Simple Pendulum. To investigate the fundamental physical properties of a simple pendulum.
Objective Simple Pendulum To investigate the fundamental physical properties of a simple pendulum. Equipment Needed Simple Pendulum Apparatus with Meter Scale and Protractor Bobs 4 (Aluminum, Brass, Lead,
More informationChapter 3 Engineering Solutions. 3.4 and 3.5 Problem Presentation
Chapter 3 Engineering Solutions 3.4 and 3.5 Problem Presentation Organize your work as follows (see book): Problem Statement Theory and Assumptions Solution Verification Tools: Pencil and Paper See Fig.
More informationEXAM. Exam 1. Math 5316, Fall December 2, 2012
EXAM Exam Math 536, Fall 22 December 2, 22 Write all of your answers on separate sheets of paper. You can keep the exam questions. This is a takehome exam, to be worked individually. You can use your notes.
More informationDetermination of Density 1
Introduction Determination of Density 1 Authors: B. D. Lamp, D. L. McCurdy, V. M. Pultz and J. M. McCormick* Last Update: February 1, 2013 Not so long ago a statistical data analysis of any data set larger
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More informationH-Pre-Calculus Targets Chapter I can write quadratic functions in standard form and use the results to sketch graphs of the function.
H-Pre-Calculus Targets Chapter Section. Sketch and analyze graphs of quadratic functions.. I can write quadratic functions in standard form and use the results to sketch graphs of the function. Identify
More informationarxiv: v1 [math.ds] 11 Feb 2011
Journal of Biological Dynamics Vol. 00, No. 00, October 2011, 1 20 RESEARCH ARTICLE Global Dynamics of a Discrete Two-species Lottery-Ricker Competition Model arxiv:1102.2286v1 [math.ds] 11 Feb 2011 Yun
More informationMA 138: Calculus II for the Life Sciences
MA 138: Calculus II for the Life Sciences David Murrugarra Department of Mathematics, University of Kentucky. Spring 2016 David Murrugarra (University of Kentucky) MA 138: Section 11.4.2 Spring 2016 1
More 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 informationError Analysis, Statistics and Graphing Workshop
Error Analysis, Statistics and Graphing Workshop Percent error: The error of a measurement is defined as the difference between the experimental and the true value. This is often expressed as percent (%)
More informationYour work from these three exercises will be due Thursday, March 2 at class time.
GEO231_week5_2012 GEO231, February 23, 2012 Today s class will consist of three separate parts: 1) Introduction to working with a compass 2) Continued work with spreadsheets 3) Introduction to surfer software
More informationIntroduction to Computer Tools and Uncertainties
Experiment 1 Introduction to Computer Tools and Uncertainties 1.1 Objectives To become familiar with the computer programs and utilities that will be used throughout the semester. To become familiar with
More informationStability Analysis of Predator- Prey Models via the Liapunov Method
Stability Analysis of Predator- Prey Models via the Liapunov Method Gatto, M. and Rinaldi, S. IIASA Research Memorandum October 1975 Gatto, M. and Rinaldi, S. (1975) Stability Analysis of Predator-Prey
More informationPhotometry of Supernovae with Makali i
Photometry of Supernovae with Makali i How to perform photometry specifically on supernovae targets using the free image processing software, Makali i This worksheet describes how to use photometry to
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 informationMath 132. Population Growth: Raleigh and Wake County
Math 132 Population Growth: Raleigh and Wake County S. R. Lubkin Application Ask anyone who s been living in Raleigh more than a couple of years what the biggest issue is here, and if the answer has nothing
More informationA Discrete Model of Three Species Prey- Predator System
ISSN(Online): 39-8753 ISSN (Print): 347-670 (An ISO 397: 007 Certified Organization) Vol. 4, Issue, January 05 A Discrete Model of Three Species Prey- Predator System A.George Maria Selvam, R.Janagaraj
More informationEXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY. S. H. Saker
Nonlinear Funct. Anal. & Appl. Vol. 10 No. 005 pp. 311 34 EXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY S. H. Saker Abstract. In this paper we derive
More information371 Lab Rybolt Data Analysis Assignment Name
Data Analysis Assignment 1 371 Lab Rybolt Data Analysis Assignment Name You wake up one morning and feel you may have a fever. You have an oral thermometer marked in Celsius degrees and find your temperature
More informationEach copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.
Host-Parasitoid Systems in Patchy Environments: A Phenomenological Model Author(s): Robert M. May Source: Journal of Animal Ecology, Vol. 47, No. 3 (Oct., 1978), pp. 833-844 Published by: British Ecological
More informationDeterministic Changes - Chapter 5 of Heinz
Deterministic Changes - Chapter 5 of Heinz Mathematical Modeling, Spring 2019 Dr Doreen De Leon 1 Introduction - Section 51 of Heinz Our plan now is to use mathematics to explain changes in variables There
More informationDynamics of a plant-herbivore model
Journal of Biological Dynamics Vol., No., Month-Month 2x, 4 Dynamics of a plant-herbivore model Yun Kang, Dieter Armbruster and Yang Kuang (Received Month 2x; revised Month 2x; in final form Month 2x)
More informationSpring /30/2013
MA 138 - Calculus 2 for the Life Sciences FINAL EXAM Spring 2013 4/30/2013 Name: Sect. #: Answer all of the following questions. Use the backs of the question papers for scratch paper. No books or notes
More informationSome hints for the Radioactive Decay lab
Some hints for the Radioactive Decay lab Edward Stokan, March 7, 2011 Plotting a histogram using Microsoft Excel The way I make histograms in Excel is to put the bounds of the bin on the top row beside
More informationEXPERIMENT: THE SPRING II Subtitle: Time Dependence of Oscillations
EXPERIMENT: THE SPRING II Subtitle: Time Dependence of Oscillations OBJECTIVES : A stretched spring is a good example of an oscillator, which oscillates around an equilibrium position ( sum of all forces
More informationThe Coefficient of Friction
The Coefficient of Friction OBJECTIVE To determine the coefficient of static friction between two pieces of wood. To determine the coefficient of kinetic friction between two pieces of wood. To investigate
More informationKEELE UNIVERSITY SCHOOL OF CHEMICAL AND PHYSICAL SCIENCES Year 1 ASTROPHYSICS LAB. WEEK 1. Introduction
KEELE UNIVERSITY SCHOOL OF CHEMICAL AND PHYSICAL SCIENCES Year 1 ASTROPHYSICS LAB WEEK 1. Introduction D. E. McLaughlin January 2011 The purpose of this lab is to introduce you to some astronomical terms
More informationSUMMER MATH PACKET ALGEBRA TWO COURSE 229
SUMMER MATH PACKET ALGEBRA TWO COURSE 9 MATH SUMMER PACKET INSTRUCTIONS MATH SUMMER PACKET INSTRUCTIONS Attached you will find a packet of exciting math problems for your enjoyment over the summer. The
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationPhysics 4C Simple Harmonic Motion PhET Lab
Physics 4C Simple Harmonic Motion PhET Lab Scott Hildreth Chabot College Goal: Explore principles of Simple Harmonic Motion through both hanging masses and pendula. Then, verify your understanding of how
More informationAP PHYSICS C Mechanics - SUMMER ASSIGNMENT FOR
AP PHYSICS C Mechanics - SUMMER ASSIGNMENT FOR 2018-2019 Dear Student: The AP physics course you have signed up for is designed to prepare you for a superior performance on the AP test. To complete material
More informationTests for Two Coefficient Alphas
Chapter 80 Tests for Two Coefficient Alphas Introduction Coefficient alpha, or Cronbach s alpha, is a popular measure of the reliability of a scale consisting of k parts. The k parts often represent k
More informationDOC / GRAPH QUADRATIC FUNCTIONS CALCULATOR
04 April, 2018 DOC / GRAPH QUADRATIC FUNCTIONS CALCULATOR Document Filetype: PDF 232.11 KB 0 DOC / GRAPH QUADRATIC FUNCTIONS CALCULATOR Let's explore how to solve quadratic equations by looking. It can
More informationLooking hard at algebraic identities.
Looking hard at algebraic identities. Written by Alastair Lupton and Anthony Harradine. Seeing Double Version 1.00 April 007. Written by Anthony Harradine and Alastair Lupton. Copyright Harradine and Lupton
More informationPurpose: Materials: WARNING! Section: Partner 2: Partner 1:
Partner 1: Partner 2: Section: PLEASE NOTE: You will need this particular lab report later in the semester again for the homework of the Rolling Motion Experiment. When you get back this graded report,
More informationInterdisciplinary Lively Application Project Spread of Disease Activity
Interdisciplinary Lively Application Project Spread of Disease Activity Title: Spread of Disease Activity Authors: Bruce MacMillan Lynn Bennethum University of Colorado at Denver University of Colorado
More informationUsing Microsoft Excel
Using Microsoft Excel Objective: Students will gain familiarity with using Excel to record data, display data properly, use built-in formulae to do calculations, and plot and fit data with linear functions.
More information'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. EGR 103L Fall Test 2. Michael R. Gustafson II
'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ EGR 103L Fall 2017 Test 2 Michael R. Gustafson II Name (please print) NET ID (please print): In keeping with the Community Standard, I have neither provided
More informationAugust 7, 2007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION
August 7, 007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION PURPOSE: This experiment illustrates the numerical solution of Laplace's Equation using a relaxation method. The results of the relaxation method
More informationRemember that C is a constant and ë and n are variables. This equation now fits the template of a straight line:
CONVERTING NON-LINEAR GRAPHS INTO LINEAR GRAPHS Linear graphs have several important attributes. First, it is easy to recognize a graph that is linear. It is much more difficult to identify if a curved
More informationLet's look at some higher order equations (cubic and quartic) that can also be solved by factoring.
GSE Advanced Algebra Polynomial Functions Polynomial Functions Zeros of Polynomial Function Let's look at some higher order equations (cubic and quartic) that can also be solved by factoring. In the video,
More informationGetting to the Roots of Quadratics
NAME BACKGROUND Graphically: The real roots of a function are the x-coordinates of the points at which the graph of the function intercepts/crosses the x-axis. For a quadratic function, whose graph is
More informationM61 1 M61.1 PC COMPUTER ASSISTED DETERMINATION OF ANGULAR ACCELERATION USING TORQUE AND MOMENT OF INERTIA
M61 1 M61.1 PC COMPUTER ASSISTED DETERMINATION OF ANGULAR ACCELERATION USING TORQUE AND MOMENT OF INERTIA PRELAB: Before coming to the lab, you must write the Object and Theory sections of your lab report
More informationSpeed of waves. Apparatus: Long spring, meter stick, spring scale, stopwatch (or cell phone stopwatch)
Name: Speed of waves Group Members: Date: TA s Name: Apparatus: Long spring, meter stick, spring scale, stopwatch (or cell phone stopwatch) Objectives 1. To directly calculate the speed of waves in a stretched
More informationChemical Kinetics: Integrated Rate Laws. ** updated Procedure for Spec 200 use **
Chemical Kinetics: Integrated Rate Laws ** updated Procedure for Spec 200 use ** *DISCLAIMER: It is highly recommended that students bring in their own computers to lab this week to use excel. There may
More informationConnecticut Common Core Algebra 1 Curriculum. Professional Development Materials. Unit 8 Quadratic Functions
Connecticut Common Core Algebra 1 Curriculum Professional Development Materials Unit 8 Quadratic Functions Contents Activity 8.1.3 Rolling Ball CBR Activity 8.1.7 Galileo in Dubai Activity 8.2.3 Exploring
More informationA Scientific Model for Free Fall.
A Scientific Model for Free Fall. I. Overview. This lab explores the framework of the scientific method. The phenomenon studied is the free fall of an object released from rest at a height H from the ground.
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 information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationMath 2930 Worksheet Equilibria and Stability
Math 2930 Worksheet Equilibria and Stabilit Week 3 September 7, 2017 Question 1. (a) Let C be the temperature (in Fahrenheit) of a cup of coffee that is cooling off to room temperature. Which of the following
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More information1 Overview of Simulink. 2 State-space equations
Modelling and simulation of engineering systems Simulink Exercise 1 - translational mechanical systems Dr. M. Turner (mct6@sun.engg.le.ac.uk 1 Overview of Simulink Simulink is a package which runs in the
More information7. STREAMBED TEXTURE ANALYSIS
Geomorphology 7. Streambed Texture Analysis 7. STREAMBED TEXTURE ANALYSIS 50 Points The purpose of this exercise is to analyze the size characteristics of the sediment samples collected in the field. We
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Experiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.01 Fall Term 2010 Experiment 03: Work and Energy Purpose of the Experiment: In this experiment you allow a cart to roll down an inclined
More informationTaguchi Method and Robust Design: Tutorial and Guideline
Taguchi Method and Robust Design: Tutorial and Guideline CONTENT 1. Introduction 2. Microsoft Excel: graphing 3. Microsoft Excel: Regression 4. Microsoft Excel: Variance analysis 5. Robust Design: An Example
More informationDynamics of a plant-herbivore model Yun Kang a ; Dieter Armbruster a ; Yang Kuang a a
This article was downloaded by: On: 30 September 2010 Access details: Access Details: Free Access Publisher Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More informationPHY222 Lab 2 - Electric Fields Mapping the Potential Curves and Field Lines of an Electric Dipole
Print Your Name PHY222 Lab 2 - Electric Fields Mapping the Potential Curves and Field Lines of an Electric Dipole Print Your Partners' Names Instructions January 23, 2015 Before lab, read the Introduction,
More informationFIBONACCI S NUMBERS, A POPULATION MODEL, AND POWERS OF MATRICES
FIBONACCI S NUMBERS, A POPULATION MODEL, AND POWERS OF MATRICES The goal of these notes is to illustrate an application of large powers of matrices Our primary tools are the eigenvalues and eigenvectors
More informationUnit 2, Ongoing Activity, Little Black Book of Algebra II Properties
Unit 2, Ongoing Activity, Little Black Book of Algebra II Properties Little Black Book of Algebra II Properties Unit 2 - Polynomial Equations & Inequalities 2.1 Laws of Exponents - record the rules for
More information5-Sep-15 PHYS101-2 GRAPHING
GRAPHING Objectives 1- To plot and analyze a graph manually and using Microsoft Excel. 2- To find constants from a nonlinear relation. Exercise 1 - Using Excel to plot a graph Suppose you have measured
More informationChemical Kinetics I: The Dry Lab. Up until this point in our study of physical chemistry we have been interested in
Chemical Kinetics I: The Dry Lab Up until this point in our study of physical chemistry we have been interested in equilibrium properties; now we will begin to investigate non-equilibrium properties and
More informationPurpose: This lab is an experiment to verify Malus Law for polarized light in both a two and three polarizer system.
Purpose: This lab is an experiment to verify Malus Law for polarized light in both a two and three polarizer system. The basic description of Malus law is given as I = I 0 (cos 2 θ ) Where I is the transmitted
More informationSystems of Algebraic Equations and Systems of Differential Equations
Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices
More informationLAB 5 INSTRUCTIONS LINEAR REGRESSION AND CORRELATION
LAB 5 INSTRUCTIONS LINEAR REGRESSION AND CORRELATION In this lab you will learn how to use Excel to display the relationship between two quantitative variables, measure the strength and direction of the
More informationDescribing the Relationship between Two Variables
1 Describing the Relationship between Two Variables Key Definitions Scatter : A graph made to show the relationship between two different variables (each pair of x s and y s) measured from the same equation.
More information