Optimal Control Applied to the Spread of Influenza A(H1N1)
|
|
- Caren Richardson
- 6 years ago
- Views:
Transcription
1 Applied Matematical Sciences, Vol. 6, 2012, no. 82, Optimal Control Applied to te Spread of Influenza AH11 M. El ia 1, O. Balatif 2, J. Bouyagroumni, E. Labriji, M. Racik Laboratoire d Analyse Modélisation et Simulation, Département de Matématiques et d Informatique, Faculté des Sciences Ben M Sik, Université Hassan II Moammedia, B.P 7955, Sidi Otman, Casablanca, Maroc Abstract Te aim of tis work is te application of optimal control to te study of te impact of vaccination on te spread of influenza A H11. Seeking to reduce te infected group, we use a control simulating a vaccination program. Te Pontryagin s maximum principle is used to caracterize te optimal control. Te optimality system is derived and solved numerically. Keywords: Optimal control, influenza A H11, Pontryagin s maximum principle 1 Introduction Influenza A H is an acute respiratory disease of umans appeared in Te contamination is mainly by air, tat is to say, couging and sneezing. Te virus can survive 8 to 48 ours in te open air, depending on te nature of te surface on wic it rests. It caused an influenza epidemic in te monts tat followed its appearance. Given te scale of te epidemic, te WHO described te pandemic in June In August 2010, te world came into a post-pandemic stage according to WHO. However, it is only a warning, te virus is still circulating widely around te globe. To deal wit tis epidemic, many governments ave decided to implement toug measures; for example, in several countries, most governments ave 1 Corresponding autor. elia moamed@yaoo.fr 2 Corresponding autor. balatif.mats@gmail.com
2 4058 M. El ia et al opted for mass vaccination plans as a precaution of te epidemic. In tis context and given tat te optimal control teory as been used successfully to make decisions involving biological or medical models, we propose to study, in tis paper, te impact of a vaccination campaign on te spread of te influenza AH11. Te paper is organized as follows. In section 2, we present a matematical model describing te spread of te influenza AH11 wit a control term. Te analysis of optimization problem is presented in section 3. In section 4, we give a numerical appropriate metod and te simulation corresponding results. Finally, te conclusion are summarized in section 5. 2 Matematical model In tis article we consider te simple model SIR used in [1]. Te ost population is divided into tree epidemiological classes: namely susceptible S, infective I and removed R. denotes te total population. We assume tat an individual can be infected only troug contacts wit infectious individuals. Tis model is governed by te following system of ordinary differential equations ds dt = Λ µs βs I di dt = βs I µ + d + r I dr = ri µr dt were all parameters are non negative and defined as follows Parameter Definition β Effective contact rate Λ Recruitment rate µ atural mortality rate d H11 induced mortality rate r Recovery rate Table 1: Parameter definitions We introduce into te above mentioned model a control ut representing te vaccination rate at time t. Te control ut is te fraction of susceptibles individuals being vaccinated per unit time. We assume tat all susceptibles vaccinees are transferred directly to te removed class. Te matematical
3 Control of AH system wit control is given by te nonlinear differential equations ds dt = Λ µs βs I us di dt = βs I µ + d + r I dr = ri µr + us dt 1 wit S0 0, I0 0 and R0 0 are given and t = St + It + Rt for all t. 3 Te optimal control problem In tis section we use te optimal control teory to analyze te beavior of te model 1. Our goal is to reduce te infected individuals and te cost of vaccination. Matematically, te problem is to minimize te objective functional Ju = t 0 fit + A 2 u2 tdt 2 were te parameter A 0 denotes te weigt on cost. And t f represents te duration of te vaccination program. In oter words, we seek te optimal control u suc tat J u = min {Ju : u U} 3 were U is te set of admissible controls defined by U = {ut : 0 u 0.9, t [0, t f ], u is Lebesgue mesurable} Te Pontryagin s maximum principle [2] converted 1, 2, 3 into problem of minimizing an Hamiltonian, H, defined by H = It + A 2 u2 t + 3 λ i f i 4 i=1 were f i is te rigt side of te differential equation of te i t state variable. By applying te Pontryagin s maximum principle [2] and te existence result of optimal control from [3], we obtain te following teorem:
4 4060 M. El ia et al Teorem 1 Tere exists an optimal control u, and corresponding solution S, I and R, tat minimizes Ju over U. Moreover, tere exists adjoint functions, λ 1, λ 2 and λ 3 verifying λ 1 = λ 1 µ + λ 1 λ 3 u + λ 1 λ 2 β I λ 2 = 1 + λ 1 λ 2 β S + λ 2µ + r + d λ 3 r λ 3 = λ 3 µ wit te transversality conditions λ 1 t f = λ 2 t f = λ 3 t f = 0 Futermore, te optimal control u is given by Proof. u = min0.9, max0, λ 1 λ 3 S 5 A Te existence of optimal control can be proved by using te results from [3] see Teorem 2.1. Te adjoint equations and transversality conditions can be obtained by using Pontryagin s Maximum Principle suc tat λ 1 = H, λ S 1t f = 0 λ 2 = H I, λ 2t f = 0 λ 3 = H, λ R 3t f = 0 Te optimal control u can be solve from te optimality condition, H u = 0 tat is H u = Au λ 1S + λ 3 S = 0 By te bounds in U of te controls, it is easy to obtain u in te form of 5 4 umerical simulations In tis section we first present an iterative metod for te numerical solution of te optimality system. ext, we present numerical results obtained using Matlab.
5 Control of AH Algoritm Te numerical algoritm presented below is a semi-implicit finite difference metod. We discretize te interval [t 0, t f ] at te points t i = t 0 + i i = 0, 1,..., n, were is te time step suc tat t n = t f, [4]. ext, we define te state and adjoint variables St, It, Rt, λ 1, λ 2, λ 3 and te control u in terms of nodal points S i, I i, R i, λ i 1, λ i 2, λ i 3 and u i. ow a combination of forward and backward difference approximation is used as follows: Te Metod, developed by [5] and presented in [6] and [7], is ten read as: S i+1 S i I i = Λ µs i+1 βs i+1 ui S i+1 I i+1 I i I i+1 = βs i+1 µ + d + r I i+1 R i+1 R i = ri i+1 µr i+1 + u i S i+1 By using a similar tecnique, we approximate te time derivative of te adjoint variables by teir first-order backward-difference and we use te appropriate sceme as follows = 1 1 µ u β I = 1 + λ1 n i 1 λ2 n i 1 β S i+1 + λn i 1 2 µ + r + d 3 r = 1 3 µ Te algoritm describing te approximation metod to obtain te optimal control is te following Algoritm 2 Step1: S0 = S 0, I0 = I 0, R0 = R 0, λ i t f = 0i = 1,..., 3, u0 = 0. Step 2: for i = 0,..., n 1, do: S i+1 = S i + Λ 1 + µ + β I i + ui I i+1 = 1 + µ + r + d β S i+1 R i+1 = R i + ri i+1 + u i S i µ I i
6 4062 M. El ia et al 1 1 = 1 2 = β I i µ + β I i = λn i µ λ n i 1 1 λ3 n i 1 + λn i + ui 1 λ1 n i 1 β S i+1 3 u i 3 r + λn i β S i+1 µ + r + d Si+1 T i+1 = A u i+1 = min0.9, max0, T i+1 end for Step 3: for i = 0,..., n, write: S t i = S i, I t i = I i, R t i = R i, u t i = u i. end for 4.2 umerical results Te numerical simulations are carried out using Matlab and using te following parameter values and initial conditions: Initial conditions S I0 30 R0 28 Parameter Values β Λ µ d 0.63% r 0.2 ote tat te initial conditions and te parameter values Λ, µ, d and r are taken from [1] and te parameter β is calculated from β = R 0 µ + d + r wit R 0 = 1.9 Te graps below, allow us to compare canges in te number of infected, susceptible and removed individuals before and after te introduction of control. In figure 1 we notice tat in presence of control te number I decreases greatly. Tus te maximum number of I would be in te presence of control against oterwise. It s nearly 87% of te effectiveness of te vaccination campaign.
7 Control of AH Te figure 2 indicates tat te number S decreases more rapidly during te vaccination campaign. It reaces at te end of tis campaign against in te absence of control, i.e a reduction of cases. Te figure 3 sows tat te number of people removed begins to grow notably from 22 nd day instead of 41 st days in te absence of control. Moreover te number R at te end of te vaccination period is instead of , wic represents an increase of cases. Finally, figure 4 gives a representation of te optimal control u. Figure 1: te function I wit and witout control
8 4064 M. El ia et al Figure 2: te function S wit and witout control Figure 3: te function R wit and witout control Figure 4: te optimal control u
9 Control of AH Conclusion In tis work we discussed an efficient numerical metod based on optimal control to study te impact of a vaccincation strategy on te spread of te influenza AH11. Te numerical simulation of bot te systems i.e wit control and witout control, sows tat tis strategy elps to reduce te number of infected and susceptible individuals and increase te number of te removed individuals greatly. Te results obtained sows also tat te effectiveness of vaccination campaign can reac 87%. References [1] K. Hattaf,. Yousfi- Matematical Model of te Influenza AH11 Infection- Advanced Studies in Biology, vol. 1, 2009 [2] L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, E. F. Miscenko, Te Matematical Teory of Optimal Processes, Wiley, ew York, [3] H.Raj josi, S. Lenart, M Y Li, L. Wang, Optimal control metods applied to disease models-contemporary Matematics, Volume 410, 2006 [4] A. B. Gumel, K. C. Patidar, and R. J. Spiteri, editors, Asymptotically Consistent on-standard Finite-Difference Metods for Solving Matematical Models Arising in Population Biology, R. E. Mickens and Worl Scientific, Singapore, [5] A. B. Gumel, P.. Sivakumar, and B. M. Saai, A matematical model for te dynamics of HIV-1 during te typical course of infection, Tird world congress of nonlinear analysts, 2001, 47: [6] K. Hattaf, M. Racik, S. Saadi, Y. Tabit et. Yousfi- Optimal Control of Tuberculosis wit Exogenous Reinfection- Applied Matematical Sciences, Vol. 3, 2009, no. 5, [7] J. Karrakcou, M. Racik, and S. Gourari, Optimal control and Infectiology: Application to an HIV/AIDS Model, Applied Matematics and Computation, 2006, 177: Received: Marc, 2012
Research Article Optimal Control of an SIR Model with Delay in State and Control Variables
ISRN Biomathematics Volume, Article ID 4549, 7 pages http://dx.doi.org/.55//4549 Research Article Optimal Control of an SIR Model with Delay in State and Control Variables Mohamed Elhia, Mostafa Rachik,
More informationBifurcation Analysis of a Vaccination Model of Tuberculosis Infection
merican Journal of pplied and Industrial Cemistry 2017; 1(1): 5-9 ttp://www.sciencepublisinggroup.com/j/ajaic doi: 10.11648/j.ajaic.20170101.12 Bifurcation nalysis of a Vaccination Model of Tuberculosis
More informationON THE GLOBAL STABILITY OF AN SIRS EPIDEMIC MODEL WITH DISTRIBUTED DELAYS. Yukihiko Nakata. Yoichi Enatsu. Yoshiaki Muroya
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: ttp://aimsciences.org pp. X XX ON THE GLOBAL STABILITY OF AN SIRS EPIDEMIC MODEL WITH DISTRIBUTED DELAYS Yukiiko Nakata Basque Center
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationOptimal Control of an SIR Epidemic Model with a Saturated Treatment
Appl. Math. Inf. Sci., No., 85-9 (26) 85 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/.8576/amis/7 Optimal Control of an SIR Epidemic Model with a Saturated Treatment
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationOptimal control of an epidemic model with a saturated incidence rate
448 Nonlinear Analysis: Modelling and Control, 2012, Vol. 17, No. 4, 448 459 Optimal control of an epidemic model with a saturated incidence rate Hassan Laarabi a, El Houssine Labriji a, Mostafa Rachik
More informationCervical Cancer and HIV Diseases Co-dynamics with Optimal Control and Cost Effectiveness
Pure and Applied Matematics Journal 017; 6(4: 14-136 ttp://www.scienceblisinggroup.com/j/pamj doi: 10.11648/j.pamj.0170604.14 ISS: 36-9790 (Print; ISS: 36-981 (Online Cervical Cancer and HIV Diseases Co-dynamics
More informationBioSystems 111 (2013) Contents lists available at SciVerse ScienceDirect. BioSystems
BioSystems 111 (213) 83 11 Contents lists available at SciVerse ScienceDirect BioSystems journa l o me pa g e: www.elsevier.com/locate/biosystems Optimal control strategies and cost-effectiveness analysis
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationA Mathematical Model of Malaria. and the Effectiveness of Drugs
Applied Matematical Sciences, Vol. 7,, no. 6, 79-95 HIKARI Ltd, www.m-ikari.com A Matematical Model of Malaria and te Effectiveness of Drugs Moammed Baba Abdullai Scool of Matematical Sciences Universiti
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationLocal and Global Stability of Host-Vector Disease Models
Local and Global Stability of Host-Vector Disease Models Marc 4, 2008 Cruz Vargas-De-León 1, Jorge Armando Castro Hernández Unidad Académica de Matemáticas, Universidad Autónoma de Guerrero, México and
More informationMathematical Modeling of Malaria
Matematical Modeling of Malaria - Metods for Simulation of Epidemics Patrik Joansson patrijo@student.calmers.se Jacob Leander jaclea@student.calmers.se Course: Matematical Modeling MVE160 Examiner: Alexei
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More information. If lim. x 2 x 1. f(x+h) f(x)
Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value
More informationAnalysis of an SEIR Epidemic Model with a General Feedback Vaccination Law
Proceedings of te World Congress on Engineering 5 Vol WCE 5 July - 3 5 London U.K. Analysis of an ER Epidemic Model wit a General Feedback Vaccination Law M. De la en. Alonso-Quesada A.beas and R. Nistal
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationPrinted Name: Section #: Instructor:
Printed Name: Section #: Instructor: Please do not ask questions during tis eam. If you consider a question to be ambiguous, state your assumptions in te margin and do te best you can to provide te correct
More informationThe Derivative as a Function
Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )
More informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More informationExam 1 Solutions. x(x 2) (x + 1)(x 2) = x
Eam Solutions Question (0%) Consider f() = 2 2 2 2. (a) By calculating relevant its, determine te equations of all vertical asymptotes of te grap of f(). If tere are none, say so. f() = ( 2) ( + )( 2)
More informationGlobal Analysis of a Mathematical Model of HCV Transmission among Injecting Drug Users and the Impact of Vaccination
Applied Mathematical Sciences, Vol. 8, 2014, no. 128, 6379-6388 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.48625 Global Analysis of a Mathematical Model of HCV Transmission among
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationMath 1210 Midterm 1 January 31st, 2014
Mat 110 Midterm 1 January 1st, 01 Tis exam consists of sections, A and B. Section A is conceptual, wereas section B is more computational. Te value of every question is indicated at te beginning of it.
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationOutline. MS121: IT Mathematics. Limits & Continuity Rates of Change & Tangents. Is there a limit to how fast a man can run?
Outline MS11: IT Matematics Limits & Continuity & 1 Limits: Atletics Perspective Jon Carroll Scool of Matematical Sciences Dublin City University 3 Atletics Atletics Outline Is tere a limit to ow fast
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationMATH CALCULUS I 2.1: Derivatives and Rates of Change
MATH 12002 - CALCULUS I 2.1: Derivatives and Rates of Cange Professor Donald L. Wite Department of Matematical Sciences Kent State University D.L. Wite (Kent State University) 1 / 1 Introduction Our main
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationKey words: HIV, prison system, epidemic model, equilibrium point, Newton s method, screening policy, quarantine policy.
COMMUICATIOS I IFORMATIO AD SYSTEMS c 27 International Press Vol. 7, o. 4, pp. 313-336, 27 1 O OPTIMAL SCREEIG AD QUARATIIG POLICY I A ETWORK OF PRISOS WAI-KI CHIG, YAG COG, ZHEG-JIA BAI, AD TUE-WAI G
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationarxiv: v1 [math.ap] 4 Aug 2017
New Two Step Laplace Adam-Basfort Metod for Integer an Non integer Order Partial Differential Equations arxiv:178.1417v1 [mat.ap] 4 Aug 17 Abstract Rodrigue Gnitcogna*, Abdon Atangana** *Department of
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More informationAPPLICATION OF OPTIMAL CONTROL TO THE EPIDEMIOLOGY OF MALARIA
Electronic Journal of Differential Equations, Vol. 212 (212, No. 81, pp. 1 22. ISSN: 172-6691. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu APPLICATION OF OPTIMAL
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More information1 Lecture 13: The derivative as a function.
1 Lecture 13: Te erivative as a function. 1.1 Outline Definition of te erivative as a function. efinitions of ifferentiability. Power rule, erivative te exponential function Derivative of a sum an a multiple
More information2.1 THE DEFINITION OF DERIVATIVE
2.1 Te Derivative Contemporary Calculus 2.1 THE DEFINITION OF DERIVATIVE 1 Te grapical idea of a slope of a tangent line is very useful, but for some uses we need a more algebraic definition of te derivative
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationGradient Descent etc.
1 Gradient Descent etc EE 13: Networked estimation and control Prof Kan) I DERIVATIVE Consider f : R R x fx) Te derivative is defined as d fx) = lim dx fx + ) fx) Te cain rule states tat if d d f gx) )
More informationNew families of estimators and test statistics in log-linear models
Journal of Multivariate Analysis 99 008 1590 1609 www.elsevier.com/locate/jmva ew families of estimators and test statistics in log-linear models irian Martín a,, Leandro Pardo b a Department of Statistics
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More information1 Limits and Continuity
1 Limits and Continuity 1.0 Tangent Lines, Velocities, Growt In tion 0.2, we estimated te slope of a line tangent to te grap of a function at a point. At te end of tion 0.3, we constructed a new function
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More information1 (10) 2 (10) 3 (10) 4 (10) 5 (10) 6 (10) Total (60)
First Name: OSU Number: Last Name: Signature: OKLAHOMA STATE UNIVERSITY Department of Matematics MATH 2144 (Calculus I) Instructor: Dr. Matias Sculze MIDTERM 1 September 17, 2008 Duration: 50 minutes No
More informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationestimate results from a recursive sceme tat generalizes te algoritms of Efron (967), Turnbull (976) and Li et al (997) by kernel smooting te data at e
A kernel density estimate for interval censored data Tierry Ducesne and James E Staord y Abstract In tis paper we propose a kernel density estimate for interval-censored data It retains te simplicity andintuitive
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationTHE IMPLICIT FUNCTION THEOREM
THE IMPLICIT FUNCTION THEOREM ALEXANDRU ALEMAN 1. Motivation and statement We want to understand a general situation wic occurs in almost any area wic uses matematics. Suppose we are given number of equations
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationEnergy-preserving variant of collocation methods 1
European Society of Computational Metods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Matematics (JNAIAM) vol.?, no.?, 21, pp. 1-12 ISSN 179 814 Energy-preserving
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationThe Complexity of Computing the MCD-Estimator
Te Complexity of Computing te MCD-Estimator Torsten Bernolt Lerstul Informatik 2 Universität Dortmund, Germany torstenbernolt@uni-dortmundde Paul Fiscer IMM, Danisc Tecnical University Kongens Lyngby,
More informationTopics in Generalized Differentiation
Topics in Generalized Differentiation J. Marsall As Abstract Te course will be built around tree topics: ) Prove te almost everywere equivalence of te L p n-t symmetric quantum derivative and te L p Peano
More informationMA119-A Applied Calculus for Business Fall Homework 4 Solutions Due 9/29/ :30AM
MA9-A Applied Calculus for Business 006 Fall Homework Solutions Due 9/9/006 0:0AM. #0 Find te it 5 0 + +.. #8 Find te it. #6 Find te it 5 0 + + = (0) 5 0 (0) + (0) + =.!! r + +. r s r + + = () + 0 () +
More informationMathematics 105 Calculus I. Exam 1. February 13, Solution Guide
Matematics 05 Calculus I Exam February, 009 Your Name: Solution Guide Tere are 6 total problems in tis exam. On eac problem, you must sow all your work, or oterwise torougly explain your conclusions. Tere
More informationChapter 4 Derivatives [ ] = ( ) ( )= + ( ) + + = ()= + ()+ Exercise 4.1. Review of Prerequisite Skills. 1. f. 6. d. 4. b. lim. x x. = lim = c.
Capter Derivatives Review of Prerequisite Skills. f. p p p 7 9 p p p Eercise.. i. ( a ) a ( b) a [ ] b a b ab b a. d. f. 9. c. + + ( ) ( + ) + ( + ) ( + ) ( + ) + + + + ( ) ( + ) + + ( ) ( ) ( + ) + 7
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationA domain decomposition convergence for elasticity equations
A domain decomposition convergence for elasticity equations A. Cakib a, A. Ellabib b, A. Nacaoui c, a Département de Matématiques Appliquées et Informatique FST de Beni-Mellal, Université Sultan My Slimane
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationBob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk
Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More information