A PREY-PREDATOR MODEL WITH AN ALTERNATIVE FOOD FOR THE PREDATOR AND OPTIMAL HARVESTING OF THE PREDATOR
|
|
- Marlene Blake
- 5 years ago
- Views:
Transcription
1 VOL. 7, O. 8, AUGUST ISS ARP Journal of Engineering and Applied Sciences 6- Asian Research Publishing ework (ARP). All righs reserved. A PREY-PREDATOR MODEL WITH A ALTERATIVE FOOD FOR THE PREDATOR AD OPTIMAL HARVESTIG OF THE PREDATOR K. Madhusudhan Reddy and K. Lakshmi arayan Vardhaman College of Engineering, Hyderabad, India SLC S Insiue of Engineering and Technology, Hyderabad, Piglipur, Hyderabad, India kmsreddy.vce@gmail.com ASTRACT The presen paper deals wih a prey - predaor model comprising an alernaive food for he predaor in addiion o he prey, and he predaor is harvesed under opimal condiions. The model is characerized by a pair of firs order nonlinear ordinary differenial equaions. All he four equilibrium poins of he model are idenified and he crieria for he sabiliy are discussed. The possibiliy of exisence of bioeconomic equilibrium is discussed. The opimal harves policy is sudied wih he help of Ponryagin s maximum principle []. Finally, some numerical examples are discussed. Keywords: prey; predaor, equilibrium poins, sabiliy, bionomic equilibrium, opimal harvesing, normal seady sae.. ITRODUCTIO Ecology relaes o he sudy of living beings in relaion o heir living syles. Research in he area of heoreical ecology was iniiaed by Loka [] and by Volerra [3]. Since hen many mahemaicians and ecologiss conribued o he growh of his area of knowledge as repored in he reaises of Paul Colinvaux [4], Freedman [5], Kapur [6, 7] ec. Harvesing of mulispecies fisheries is an imporan area of sudy in fishery modeling. The issues and echniques relaed o his field of sudy and he problem of combined harvesing of wo ecologically independen populaions obeying he logisic law of growh are discussed in deail by Clark [8, 9]. Chaudhuri [, ] formulaed an opimal conrol problem for he combined harvesing of wo compeing species. Models on he combined harvesing of a wospecies prey-predaor fishery have been discussed by Chaudhuri and Saha Ray []. iological and bionomic equilibria of a mulispecies fishery model wih opimal harvesing policy is discussed in deail by Kar and Chaudhari [3]. Recenly Archana Reddy [4] discussed he sabiliy analysis of wo ineracing species wih harvesing of boh species. Lakshmi arayan, Paabhiramacharyulu [5, 6] and Shiva Reddy [7] e al., have discussed differen prey-predaor models in deail. Mos of he mahemaical models on he harvesing of a mulispecies fishery have so far assumed ha he species are affeced by harvesing only. A populaion model proposed by Kar and Chaudhuri, (c.f. Harvesing in a wo-prey one-predaor fishery: ioeconomic model, AZIAM J.45 (4), ) moivaed he presen invesigaion.. MATHEMATICAL MODEL The model equaions for a wo species preypredaor sysem are given by he following sysem of nonlinear ordinary differenial equaions employing he following noaion: and are he populaions of he prey and predaor wih naural growh raes a and a respecively, is rae of decrease of he prey due o insufficien food, is rae of decrease of he prey due o inhibiion by he predaor, is rae of increase of he predaor due o successful aacks on he prey, is rae of decrease of he predaor due o insufficien food oher han he prey; q is he cach abiliy co-efficien of he predaor, E is he harvesing effor and qe is he cach-rae funcion based on he CPUE (cach-per-uni-effor) hypohesis. Furher boh he variables and are non-negaive and he model parameers a, a,,,,, q, E a q E are assumed o be non-negaive and ( ) consans. d d a () ( ) a q E + () 3. EQUILIRIUM STATES The sysem under invesigaion has four equilibrium saes defined by: a) The fully washed ou sae wih he equilibrium poin ; (3) b) The sae in which, only he predaor survives given ( a qe) by ; (4) c) The sae in which, only he prey survives given a by ; (5) d) The co-exisen sae (normal seady sae) given by: 3
2 VOL. 7, O. 8, AUGUST ISS ARP Journal of Engineering and Applied Sciences 6- Asian Research Publishing ework (ARP). All righs reserved. a ( a q E) ( a q E) + a ; + + This sae would exi only when a ( a q E) (6) > (7) 4. STAILITY OF THE EQUILIRIUM STATES To invesigae he sabiliy of he equilibrium saes we consider small perurbaions u, u in and over and respecively, so ha + u ; + u (8) y subsiuing (8) in () and () and neglecing second and higher powers of he perurbaions u, u, we ge he equaions of he perurbed sae: du AU (9) can be noed o be negaive. Hence he co-exisen equilibrium sae is sable. The soluions curves are: u u( λ + ) u λ λ u( λ + ) u λ λ e λ + u ( λ + ) + u u λ λ u( λ + ) + u λ λ e λ (3) e λ + e λ (4) where a () A ( a qe) + The characerisic equaion for he sysem is: [ λi ] de A () The equilibrium sae is sable only when he roos of he equaion () are negaive in case hey are real or have negaive real pars in case hey are complex. The equilibrium saes I, II, and III are found o unsable, so we resriced our sudy o he normal seady sae only. 4. Sabiliy of he normal seady sae In his case he characerisic equaion is λ + ( + ) λ+ [ ] () + Since he sum of he roos of () is negaive and he produc of he roos is posiive, he roos of which Where λ, λ, are he roos of he equaion (). 5. GLOAL STAILITY Theorem: The Equilibrium sae (, ) is globally asympoically sable. Le us consider he following Liapunov s funcion V(, ) ln + l ln (5) where l is posiive consan,o be chosen laer Differeniaing V w.r., we ge dv d d + l Subsiuing () and () in (6), we ge { ( ) ( )( ) } l ( )( ) ( ) { } dv + (6) Choosing l and wih some algebraic manipulaion yields dv ( ) ( ) < (7) Therefore, he equilibrium sae (, ) is globally asympoically sable. 6. IOOMIC EQUILIRIUM The erm bionomic equilibrium is an amalgamaion of he conceps of biological equilibrium as well as economic equilibrium. The economic equilibrium is said o be achieved when he oal revenue obained by selling he harvesed biomass equals he oal cos for he effor devoed o harvesing. Le c fishing cos per uni effor of he predaor, p price per uni biomass of he predaor. The ne economic revenue for he predaor a any ime is given by: 4
3 VOL. 7, O. 8, AUGUST ISS ARP Journal of Engineering and Applied Sciences 6- Asian Research Publishing ework (ARP). All righs reserved. ( ) R p q c E (8) ( ) The biological equilibrium is ( ),( ),( E) ( ),( ),( E), where are he posiive soluions of a (9) ( ) a q E + () and( ) pq c E () From (), we have c { pq ( ) c}( E) ( ) () p q c a From (9) and (), we ge ( ) pq (3) From (), () and (3), we ge ( E) ( ) a (4) q pq I is clear ha ( E) > if ( ) > a (5) pq Thus he bionomic equilibrium (( ),( ),( E) ) exiss, if inequaliy (5) holds. 7. OPTIMAL HARVESTIG POLICY The presen value J of a coninuous ime-sream of revenues is given by: ( ) (6) J e p q c E Where denoes he insananeous annual rae of discoun. Our problem is o maximize J subjec o he sae equaions () and () and conrol consrains E ( E) max by invoking Ponryagin s maximum principle. The Hamilonian for he problem is given by: ( ) λ( ) ( a qe) H e p q c E+ a + λ + Where λ, λ, are he adjoin variables. (7) Le us assume ha he conrol consrains are no binding i.e., he opimal soluion does no occur a( E ) max. A ( E) max we have a singular conrol. y Ponryagin s maximum principle, dλ dλ ; ; E e ( pq c) λq λ e p E q (8) c dλ + dλ { λ( a ) λ( ) } ( λ λ E) (9) dλ dλ { e pqe λ( ) λ ( a qe ) } ( λ λ e pqe) + (3) From (8) and (9), we ge c Where p q Whose soluion is given by dλ λ e λ ( + ) e (3) dλ From (3) and (3), we ge λ e Where pqe ( + ) Whose soluion is given by λ e (3) ( + ) From (8) and (3), we ge a singular pah p q ( + ) Thus (33) can be wrien as: F ( ) p q ( + ) (33) There exis a unique posiive roo ( ) of F ( ) in he inerval < < k, if he 5
4 VOL. 7, O. 8, AUGUST ISS ARP Journal of Engineering and Applied Sciences 6- Asian Research Publishing ework (ARP). All righs reserved. following hold F () <, F (k ) >, F ( ) ) > for >. ( ), we ge ( ) For and ( ) ( ) a pq c E a q pq (34) (35) Hence once he opimal equilibrium, is deermined, he opimal harvesing (( ) ( ) ) effor ( E) can be deermined. From (8), (3) and (3), we found ha λ, λ do no vary wih ime in opimal equilibrium. Hence hey remain bounded as. From (33), we also noe ha Figure-. Shows he rajecory of he prey and predaor populaions beginning wih 45 and 35. Le a 8;.;.5; a ;.;.4; q.& E 5. p q ( + ) as Thus, he ne economic revenue of he predaor R. This implies ha if he discoun rae increases, hen he ne economic revenue decreases and even may end o zero if he discoun rae end o infiniy. Thus i has been concluded ha high ineres rae will cause high inflaion rae. This conclusion was also drawn by Clark [9] in he combined harvesing of wo ecologically independen populaions and by Chaudhuri [] in he combined harvesing of wo compeing species. 8. UMERICAL EXAMPLES Figure-3. Shows he variaion of he populaions agains he ime. Le a 3;.;.; a.5;.3;.4; q.& E. Figure-. Shows he variaion of he populaions agains he ime. Figure-4. Shows he rajecory of he prey and predaor populaions beginning wih 5 and 3. 6
5 VOL. 7, O. 8, AUGUST ISS ARP Journal of Engineering and Applied Sciences 6- Asian Research Publishing ework (ARP). All righs reserved. COCLUSIOS A prey - predaor model wih an alernaive food for he predaor, and he predaor is harvesed under opimal condiions. The model is characerized by a pair of firs order non-linear ordinary differenial equaions. All he four equilibrium poins of he model are idenified and he crieria for he sabiliy were discussed by using Lyponuv s mehod. The bio-economic equilibrium poin was idenified and he opimal harves policy was sudied, and some numerical examples were discussed. Finally i was observed ha he model was sabilized by harvesing. REFERECES [] Ponryagin L.S., olyanskii V.S., Gamkerlidge R.. and Mischenko E.F. 96. The mahemaical heory of opimal process, Wiley, ew York, USA. [] Loka A.J. 95. Elemens of Physical iology, Williams and Wilkins, alimore, ew York, USA. [3] Volerra V. 93. Leçons sur la héorie mahémaique de la lue pour la vie, Gauhier-Villars, Paris, France. [4] Paul Colinvaux A Ecology. John Wiley and Sons Inc., ew York, USA. [3] Kar T.K. and ChaudhuriK. S. 4. Harvesing in a wo-prey one-predaor fishery: A bioeconomic model. AZIAM J. 45: [4] ArchanaReddy R., Paabhiramacharyulu.Ch and KrishnaGandhi. 7. A sabiliy analysis of wo compeiive ineracing species wih harvesing of boh he species a a consan rae. In. J. of scienific compuing. () Jan- June. pp [5] Lakshmiarayan K. 5. A mahemaical sudy of a prey-predaor ecological model wih a parial cover for he prey and alernaive food for he predaor, Ph.D hesis, JTU. [6] Lakshmi arayan K. and Paabhiramacharyulu. Ch. 8. A Prey-Predaor Model wih an Alernaive Food for he Predaor, Harvesing of oh he Species and wih a Gesaion Period for Ineracion. In. J. Open Problems Comp. Mah. : [7] Shiva Reddy K., Lakshmi arayan K. and Paabhiramacharyulu. Ch.. A Three Species Ecosysem Consising of a Prey and Two Predaors. Inernaional J. of Mah. Sci. and Engg. Appls. 4(IV): [5] Freedman H. I. 98. Deerminisic Mahemaical Models in Populaion Ecology, Decker, ew York, USA. [6] Kapur J Mahemaical modeling in biology and Medicine. Affiliaed Eas Wes Press. [7] Kapur J Mahemaical modeling, wiley, easer. [8] Clark C.W ioeconomic modelling and fisheries managemen. Wiley, ew York, USA. [9] Clark C.W Mahemaical bioeconomics: he opimal managemen of renewable resources. Wiley, ew York, USA. [] Chaudhuri K.S A bioeconomic model of harvesing a mulispecies fishery. Ecol. Model. 3: [] Chaudhuri K.S Dynamic opimizaion of combined harvesing of a wo species fishery. Ecol. Model. 4: 7-5. [] Chaudhuri K.S. and Saha Ray S On he combined harvesing of a prey predaor sysem. J. iol. Sis. pp
A Prey-Predator Model with an Alternative Food for the Predator and Optimal Harvesting of the Prey
Available online a www.pelagiaresearchlibrary.com Advances in Applied Science Research, 0, (4):45-459 A Prey-Predaor Model wih an Alernaive Food for he Predaor and Opimal Harvesing of he Prey K. Madhusudhan
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationA Prey-Predator Model with a Reserved Area
Nonlinear Analysis: Modelling and Conrol, 2007, Vol. 12, No. 4, 479 494 A Prey-Predaor Model wih a Reserved Area B. Dubey Mahemaics Group, Birla Insiue of Technology and Science Pilani - 333 031, Rajashan,
More informationMath 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.
Mah 36. Rumbos Spring 1 1 Soluions o Assignmen #6 1. Suppose he growh of a populaion is governed by he differenial equaion where k is a posiive consan. d d = k (a Explain why his model predics ha he populaion
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationThis document was generated at 7:34 PM, 07/27/09 Copyright 2009 Richard T. Woodward
his documen was generaed a 7:34 PM, 07/27/09 Copyrigh 2009 Richard. Woodward 15. Bang-bang and mos rapid approach problems AGEC 637 - Summer 2009 here are some problems for which he opimal pah does no
More informationDIFFUSION ANALYSIS OF A PREY PREDATOR FISHERY MODEL WITH HARVESTING OF PREY SPECIES
VOL., NO., DECEMBE 7 ISSN 89-668 APN Journal of Engineering and Applied Sciences 6-7 Asian esearch Publishing Nework (APN). All righs reserved. DIFFUSION ANALYSIS OF A PEY PEDATO FISHEY MODEL WITH HAVESTING
More informationAnn. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:
Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS
More informationADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations
ADVANCED MATHEMATICS FOR ECONOMICS - /3 Shee 3: Di erenial equaions Check ha x() =± p ln(c( + )), where C is a posiive consan, is soluion of he ODE x () = Solve he following di erenial equaions: (a) x
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationBifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays
Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationLecture Notes 3: Quantitative Analysis in DSGE Models: New Keynesian Model
Lecure Noes 3: Quaniaive Analysis in DSGE Models: New Keynesian Model Zhiwei Xu, Email: xuzhiwei@sju.edu.cn The moneary policy plays lile role in he basic moneary model wihou price sickiness. We now urn
More informationMath 2214 Solution Test 1A Spring 2016
Mah 14 Soluion Tes 1A Spring 016 sec Problem 1: Wha is he larges -inerval for which ( 4) = has a guaraneed + unique soluion for iniial value (-1) = 3 according o he Exisence Uniqueness Theorem? Soluion
More informationOn a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration
Journal of Agriculure and Life Sciences Vol., No. ; June 4 On a Discree-In-Time Order Level Invenory Model for Iems wih Random Deerioraion Dr Biswaranjan Mandal Associae Professor of Mahemaics Acharya
More informationIn Game Theory and Applications, Volume X WITH MIGRATION: RESERVED TERRITORY APPROACH. Vladimir V. Mazalov. Anna N. Rettieva
In Game Theory and Applicaions, Volume X Edied by L.A. Perosjan and V.V. Mazalov, pp.97{10 ISBN 1-59033-34-X c 004 Nova Science Publishers, Inc. A FISHERY GAME MODEL WITH MIGRATION: RESERVED TERRITORY
More informationProblem Set on Differential Equations
Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationThe motions of the celt on a horizontal plane with viscous friction
The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and
More informationThis document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC
his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,
More informationProperties Of Solutions To A Generalized Liénard Equation With Forcing Term
Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick
More informationEssential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems
Essenial Microeconomics -- 6.5: OPIMAL CONROL Consider he following class of opimizaion problems Max{ U( k, x) + U+ ( k+ ) k+ k F( k, x)}. { x, k+ } = In he language of conrol heory, he vecor k is he vecor
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationA Dynamic Model of Economic Fluctuations
CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model
More informationMacroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3
Macroeconomic Theory Ph.D. Qualifying Examinaion Fall 2005 Comprehensive Examinaion UCLA Dep. of Economics You have 4 hours o complee he exam. There are hree pars o he exam. Answer all pars. Each par has
More informationAn Inventory Model for Time Dependent Weibull Deterioration with Partial Backlogging
American Journal of Operaional Research 0, (): -5 OI: 0.593/j.ajor.000.0 An Invenory Model for Time ependen Weibull eerioraion wih Parial Backlogging Umakana Mishra,, Chaianya Kumar Tripahy eparmen of
More informationProblem Set #3: AK models
Universiy of Warwick EC9A2 Advanced Macroeconomic Analysis Problem Se #3: AK models Jorge F. Chavez December 3, 2012 Problem 1 Consider a compeiive economy, in which he level of echnology, which is exernal
More informationOnline Appendix to Solution Methods for Models with Rare Disasters
Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,
More informationANSWERS TO EVEN NUMBERED EXERCISES IN CHAPTER 6 SECTION 6.1: LIFE CYCLE CONSUMPTION AND WEALTH T 1. . Let ct. ) is a strictly concave function of c
John Riley December 00 S O EVEN NUMBERED EXERCISES IN CHAPER 6 SECION 6: LIFE CYCLE CONSUMPION AND WEALH Eercise 6-: Opimal saving wih more han one commodiy A consumer has a period uiliy funcion δ u (
More informationOscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,
More informationUSP. Surplus-Production Models
USP Surplus-Producion Models 2 Overview Purpose of slides: Inroducion o he producion model Overview of differen mehods of fiing Go over some criique of he mehod Source: Haddon 2001, Chaper 10 Hilborn and
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More information4. The multiple use forestry maximum principle. This principle will be derived as in Heaps (1984) by considering perturbations
4. The muliple use foresry maximum principle This principle will be derived as in Heaps (1984) by considering perurbaions H(; ) of a logging plan H() in A where H(; 0) = H() and H(; ) A is di ereniable
More informationt + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that
ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so
More informationSliding Mode Controller for Unstable Systems
S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.
More information1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population
Scienific Comuing I Module : Poulaion Modelling Coninuous Models Michael Bader Par I ODE Models Lehrsuhl Informaik V Winer 7/ Discree vs. Coniuous Models d d = F,,...) ) =? discree model: coninuous model:
More informationarxiv: v1 [math.gm] 4 Nov 2018
Unpredicable Soluions of Linear Differenial Equaions Mara Akhme 1,, Mehme Onur Fen 2, Madina Tleubergenova 3,4, Akylbek Zhamanshin 3,4 1 Deparmen of Mahemaics, Middle Eas Technical Universiy, 06800, Ankara,
More informationMA 366 Review - Test # 1
MA 366 Review - Tes # 1 Fall 5 () Resuls from Calculus: differeniaion formulas, implici differeniaion, Chain Rule; inegraion formulas, inegraion b pars, parial fracions, oher inegraion echniques. (1) Order
More informationKeywords: competition models; density-dependence; ecology; population dynamics; predation models; stochastic models UNESCO EOLSS
POPULATIO MODELS Michael B. Bonsall Deparmen of Zoology, Universiy of Oxford, Oxford, UK Keywords: compeiion models; densiy-dependence; ecology; populaion dynamics; predaion models; sochasic models Conens.
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More information= ( ) ) or a system of differential equations with continuous parametrization (T = R
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationResearch Article Hopf Bifurcation of a Differential-Algebraic Bioeconomic Model with Time Delay
Journal of Applied Mahemaics Volume 01, Aricle ID 768364, 15 pages doi:10.1155/01/768364 Research Aricle Hopf Bifurcaion of a Differenial-Algebraic Bioeconomic Model wih Time Delay Xiaojian Zhou, 1, Xin
More informationTheory of! Partial Differential Equations-I!
hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationSolutions Problem Set 3 Macro II (14.452)
Soluions Problem Se 3 Macro II (14.452) Francisco A. Gallego 04/27/2005 1 Q heory of invesmen in coninuous ime and no uncerainy Consider he in nie horizon model of a rm facing adjusmen coss o invesmen.
More informationAppendix 14.1 The optimal control problem and its solution using
1 Appendix 14.1 he opimal conrol problem and is soluion using he maximum principle NOE: Many occurrences of f, x, u, and in his file (in equaions or as whole words in ex) are purposefully in bold in order
More informationBoundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms
Advances in Dynamical Sysems and Applicaions. ISSN 0973-531 Volume Number 1 007, pp. 107 11 Research India Publicaions hp://www.ripublicaion.com/adsa.hm Boundedness and Exponenial Asympoic Sabiliy in Dynamical
More information1 Differential Equation Investigations using Customizable
Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationUNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS
UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS Exam: ECON4325 Moneary Policy Dae of exam: Tuesday, May 24, 206 Grades are given: June 4, 206 Time for exam: 2.30 p.m. 5.30 p.m. The problem se covers 5 pages
More informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More information8. Basic RL and RC Circuits
8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics
More informationAn Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.
1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More informationdi Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.
di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationpe pt dt = e pt Probabilty of death given survival till t : pe pt = p Expected life at t : pe(s t)p ds = e (s t)p t =
BLANCHARD Probabiliy of Deah: π () = pe p ; Probabily of living ill : Ω () = pe p d = e p Probabily of deah given survival ill : pe p = p e p Expeced life a : (s ) pe (s )p ds = p 1 Populaion normalized
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationdy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page
Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,
More informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationASTR415: Problem Set #5
ASTR45: Problem Se #5 Curran D. Muhlberger Universi of Marland (Daed: April 25, 27) Three ssems of coupled differenial equaions were sudied using inegraors based on Euler s mehod, a fourh-order Runge-Kua
More informationReview - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y
Review - Quiz # 1 (1) Solving Special Tpes of Firs Order Equaions I. Separable Equaions (SE). d = f() g() Mehod of Soluion : 1 g() d = f() (The soluions ma be given implicil b he above formula. Remember,
More informationMath Final Exam Solutions
Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,
More informationMA Study Guide #1
MA 66 Su Guide #1 (1) Special Tpes of Firs Order Equaions I. Firs Order Linear Equaion (FOL): + p() = g() Soluion : = 1 µ() [ ] µ()g() + C, where µ() = e p() II. Separable Equaion (SEP): dx = h(x) g()
More informationTheory of! Partial Differential Equations!
hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationElementary Differential Equations and Boundary Value Problems
Elemenar Differenial Equaions and Boundar Value Problems Boce. & DiPrima 9 h Ediion Chaper 1: Inroducion 1006003 คณ ตศาสตร ว ศวกรรม 3 สาขาว ชาว ศวกรรมคอมพ วเตอร ป การศ กษา 1/2555 ผศ.ดร.อร ญญา ผศ.ดร.สมศ
More information( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationPOSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION
Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial
More informationIntermediate Differential Equations Review and Basic Ideas
Inermediae Differenial Equaions Review and Basic Ideas John A. Burns Cener for Opimal Design And Conrol Inerdisciplinary Cener forappliedmahemaics Virginia Polyechnic Insiue and Sae Universiy Blacksburg,
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationCash Flow Valuation Mode Lin Discrete Time
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationChapter 7 Response of First-order RL and RC Circuits
Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial
More informationSolutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore
Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s
More informationLet ( α, β be the eigenvector associated with the eigenvalue λ i
ENGI 940 4.05 - Sabiliy Analysis (Linear) Page 4.5 Le ( α, be he eigenvecor associaed wih he eigenvalue λ i of he coefficien i i) marix A Le c, c be arbirary consans. a b c d Case of real, disinc, negaive
More informationImproved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method
Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics
More informationOptimal Fishery. Modelling Population Dynamics with Optimal Harvest Strategies. Degree Project in Engineering Physics, First Level, SA114X
Opimal Fishery Modelling Populaion Dynamics wih Opimal Harves Sraegies Degree Projec in Engineering Physics, Firs Level, SA4X Deparmen of Mahemaics, division of Numerical Analysis Mas Barkman Supervisor:
More informationSZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1
SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision
More informationComparing Theoretical and Practical Solution of the First Order First Degree Ordinary Differential Equation of Population Model
Open Access Journal of Mahemaical and Theoreical Physics Comparing Theoreical and Pracical Soluion of he Firs Order Firs Degree Ordinary Differenial Equaion of Populaion Model Absrac Populaion dynamics
More informationA Study of Inventory System with Ramp Type Demand Rate and Shortage in The Light Of Inflation I
Inernaional Journal of Mahemaics rends and echnology Volume 7 Number Jan 5 A Sudy of Invenory Sysem wih Ramp ype emand Rae and Shorage in he Ligh Of Inflaion I Sangeea Gupa, R.K. Srivasava, A.K. Singh
More informationENV 6015 Solution to Mixing Problem Set
EN 65 Soluion o ixing Problem Se. A slug of dye ( ) is injeced ino a single ank wih coninuous mixing. The flow in and ou of he ank is.5 gpm. The ank volume is 5 gallons. When will he dye concenraion equal
More information1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.
. Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN
Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.
More informationSTABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS
Elecronic Journal of Differenial Equaions, Vol. 217 217, No. 118, pp. 1 14. ISSN: 172-6691. URL: hp://ejde.mah.xsae.edu or hp://ejde.mah.un.edu STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS
More informationProduction Inventory Model with Different Deterioration Rates Under Shortages and Linear Demand
Inernaional Refereed Journal of Engineering and Science (IRJES) ISSN (Online) 39-83X, (Prin) 39-8 Volume 5, Issue 3 (March 6), PP.-7 Producion Invenory Model wih Differen Deerioraion Raes Under Shorages
More informationEconomics 8105 Macroeconomic Theory Recitation 6
Economics 8105 Macroeconomic Theory Reciaion 6 Conor Ryan Ocober 11h, 2016 Ouline: Opimal Taxaion wih Governmen Invesmen 1 Governmen Expendiure in Producion In hese noes we will examine a model in which
More informationSignal and System (Chapter 3. Continuous-Time Systems)
Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b
More informationMath 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm
Mah Spring 6 Differenial Equaions Final Exam Due Dae: Tuesday, June 6, 5pm Your name (please prin): Insrucions: This is an open book, open noes exam. You are free o use a calculaor or compuer o check your
More information