A Model of two mutually interacting Species with Mortality Rate for the Second Species

Similar documents
A PREY-PREDATOR MODEL WITH COVER FOR THE PREY AND AN ALTERNATIVE FOOD FOR THE PREDATOR AND CONSTANT HARVESTING OF BOTH THE SPECIES *

ON PREY-PREDATOR MODEL WITH HOLLING-TYPE II AND LESLIE-GOWER SCHEMES AHMED BUSERI ASHINE

Available online at Pelagia Research Library. Advances in Applied Science Research, 2011, 2 (3):

A MATHEMATICAL MODEL OF FOUR SPECIES SYN-ECOSYMBIOSIS COMPRISING OF PREY-PREDATION, MUTUALISM AND COMMENSALISMS-I (FULLY WASHED OUT STATE)

Solutions of Klein - Gordan equations, using Finite Fourier Sine Transform

Linear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Math 31S. Rumbos Fall Solutions to Assignment #16

Ordinary differential equations

MATH 3795 Lecture 18. Numerical Solution of Ordinary Differential Equations.

Heteroclinic cycles in coupled cell systems

Module 6: LINEAR TRANSFORMATIONS

INVESTIGATION OF MATHEMATICAL MODEL OF COMMUNICATION NETWORK WITH UNSTEADY FLOW OF REQUESTS

than 1. It means in particular that the function is decreasing and approaching the x-

A General Dynamic Inequality of Opial Type

MAC-solutions of the nonexistent solutions of mathematical physics

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

Lyapunov function for cosmological dynamical system

Conservation Law. Chapter Goal. 5.2 Theory

Modeling and Controller Design for the Air-to-Air Missile Uncertain System

TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

BIFURCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS

Review of Calculus, cont d

ADAPTIVE CONTROL AND SYNCHRONIZATION OF A HIGHLY CHAOTIC ATTRACTOR

Math 42 Chapter 7 Practice Problems Set B

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Kinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2)

A New Generalization of Lemma Gronwall-Bellman

AMATH 731: Applied Functional Analysis Fall Some basics of integral equations

MATHS NOTES. SUBJECT: Maths LEVEL: Higher TEACHER: Aidan Roantree. The Institute of Education Topics Covered: Powers and Logs

New Expansion and Infinite Series

Factors affecting the phonation threshold pressure and frequency

On a predator prey system with cross-diffusion representing the tendency of prey to keep away from its predators

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer

Improper Integrals, and Differential Equations

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

Construction and Selection of Single Sampling Quick Switching Variables System for given Control Limits Involving Minimum Sum of Risks

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

ANALYSIS OF A PREDATOR-PREY POPULATION MODEL

Consequently, the temperature must be the same at each point in the cross section at x. Let:

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

Generation of Lyapunov Functions by Neural Networks

Chaos in drive systems

Here are the graphs of some power functions with negative index y (x) =ax n = a n is a positive integer, and a 6= 0acoe±cient.

An iterative method for solving nonlinear functional equations

On the Linear Stability of Compound Capillary Jets

Deteriorating Inventory Model for Waiting. Time Partial Backlogging

A Modified ADM for Solving Systems of Linear Fredholm Integral Equations of the Second Kind

Introduction to the Calculus of Variations

1 2-D Second Order Equations: Separation of Variables

A Fuzzy Differential Approach to strong Allee Effect based on the Fuzzy Extension Principle

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Research Article On Existence and Uniqueness of Solutions of a Nonlinear Integral Equation

Mathematics Extension 1

Math 1B, lecture 4: Error bounds for numerical methods

A Bernstein polynomial approach for solution of nonlinear integral equations

Introduction to Electronic Circuits. DC Circuit Analysis: Transient Response of RC Circuits

International ejournals

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

Solution to Fredholm Fuzzy Integral Equations with Degenerate Kernel

Kinematics equations, some numbers

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Unit 1 Exponentials and Logarithms

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

Time delay in a basic model of the immune response

Introduction to Group Theory

Section 4: Integration ECO4112F 2011

The Wave Equation I. MA 436 Kurt Bryan

Arithmetic Mean Derivative Based Midpoint Rule

1.9 C 2 inner variations

du = C dy = 1 dy = dy W is invertible with inverse U, so that y = W(t) is exactly the same thing as t = U(y),

Songklanakarin Journal of Science and Technology SJST R1 Thongchan. A Modified Hyperbolic Secant Distribution

ME 501A Seminar in Engineering Analysis Page 1

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

Math 360: A primitive integral and elementary functions

CHM Physical Chemistry I Chapter 1 - Supplementary Material

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

Lesson 1: Quadratic Equations

3 Conservation Laws, Constitutive Relations, and Some Classical PDEs

Modelling Intermittant Androgen Deprivation Therapy. Alex Browning Dr Matthew Simpson Queensland University of Technology

DIRECT CURRENT CIRCUITS

Operations with Polynomials

Continuous Random Variables

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY

Math 6395: Hyperbolic Conservation Laws and Numerical Methods. Hyperbolic Conservation Laws are a type of PDEs arise from physics: fluid/gas dynamics.

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

THIELE CENTRE. Linear stochastic differential equations with anticipating initial conditions

Conservation Law. Chapter Goal. 6.2 Theory

Part I: Basic Concepts of Thermodynamics

5.3 Nonlinear stability of Rayleigh-Bénard convection

The Solution of Volterra Integral Equation of the Second Kind by Using the Elzaki Transform

Remark on boundary value problems arising in Ginzburg-Landau theory

Transcription:

Avilble online t www.pelgireserchlibrry.com Advnces in Applied Science Reserch, 0, 3 ():757-764 ISS: 0976-860 CODE (USA): AASRFC A Model of two mutully intercting Species with Mortlity Rte for the Second Species B. Rvindr Reddy JTUH College of Engineering, chuplly, Krimngr, Indi _ ABSTRACT The present pper concerns with model of two mutully intercting species with limited resources for first species nd unlimited resources for second species with mortlity rte. The model is chrcterized by coupled system of first order non-liner ordinry differentil equtions. In ll two equilibrium points re identified. If the deth rte of the second species is greter thn its birth rte, it is found tht there re three equilibrium points. The criteri for symptotic stbility hve been estblished for ll the equilibrium points. The co-existent stte is lwys stble. The solutions of the linerised bsic equtions re obtined nd their trends re illustrted. Key words: Equilibrium points, Mutulism, Coexistence stte, Stbility. _ ITRODUCTIO Mthemticl modeling is n importnt interdisciplinry ctivity which involves the study of some spects of diverse disciplines. Biology, Epidemiodology, Physiology, Ecology, Immunology, Bio-economics, Genetics, Phrmocokinetics re some of those disciplines. This mthemticl modeling hs rised to the zenith in recent yers nd spred to ll brnches of life nd drew the ttention of every one. Mthemticl modeling of ecosystems ws initited by Lotk [8] nd by Volterr [4]. The generl concept of modeling hs been presented in the tretises of Meyer [9], Cushing [], Pul Colinvux [0], Freedmn [3], Kpur [5, 6]. The ecologicl interctions cn be brodly clssified s Prey-Predtion, Competition, Mutulism nd so on..c. Srinivs [3] studied the competitive eco-systems of two species nd three species with regrd to limited nd unlimited resources. Lkshmi ryn [7] investigted the two species prey-predtor models nd stbility nlysis of competitive species ws investigted by Archn Reddy []. Locl stbility nlysis for two-species ecologicl mutulism model hs been presented by the present uthor et l. [, ]. Recently, stbility nlysis of three species ws crried out by Siv Reddy [8]. Further Srilth et l. [5, 6] nd Hri Prsd et.l [7] studied stbility nlysis of four species. The present investigtion is devoted to the nlyticl study of model of two mutully intercting species with mortlity rte for the second species. Before describing model, first we mke the following ssumptions: is the popultion of the first species,, the popultion of the second species, is the rte of nturl growth of the first species, is the rte of nturl growth of the second species, is the rte of decrese of the first species e to insufficient food, is the rte of increse of the first species e to interction with the second species, is the rte of increse of the second species e to interction with the first species. Further note tht the vribles, nd the model prmeters,,,, re non-negtive nd tht the rte of difference between the deth nd birth rtes is identified s the nturl growth rte with pproprite sign. The model equtions for two species mutulising re governed by system of non-liner ordinry differentil equtions. 757

B. Rvindr Reddy Adv. Appl. Sci. Res., 0, 3():757-764. Bsic Equtions: The bsic equtions re given by d d + = (.) = + (.) Here we come cross three equilibrium sttes: = 0; = 0, I. the stte in which both the species re wshed out. II. = ; 0 =. Here the first species ( ) survives while the second species ( ) is wshed out. III. = ; =. In this stte both the species co-exist nd this cn exist only when > 0. Equilibrium stte I (fully wshed out stte): To discuss the stbility of equilibrium stte u ( ) t from the stedy stte, i.e. we write = 0 ; = 0, we consider smll perturbtions u ( t) nd = + u ( t), (.3) = + u ( t) (.4) Substituting (.3) nd (.4) in (.) nd (.), we get = u u + u u After lineriztion, we get nd = u + u u u = (.5) = u (.6) The chrcteristic eqution is (λ - )(λ + ) =0 One root of this eqution is λ = which is positive nd the other root is λ = which is negtive. Hence the equilibrium stte is unstble. 758

B. Rvindr Reddy Adv. Appl. Sci. Res., 0, 3():757-764 The solutions of equtions (.5) nd (.6) re t u = u0 e (.7) t u = u e 0 where u 0, u 0 re the initil vlues of u nd u. The solution curves re illustrted in figures nd. (.8) Cse : u 0 > u 0 i.e. initilly the first species domintes the second species. We notice tht the first species is going wy from the equilibrium point while the second species pproches symptoticlly to the equilibrium point. Hence the stte is unstble. Cse : u 0 < u 0 i.e. initilly the second species domintes the first species. In this cse the second species out numbers the first species till the time, ln { u0 / u0} t = t* = ( + ) fter tht the first species out numbers the second species nd grows indefinitely while the second species pproches symptoticlly to the equilibrium point. Hence the stte is unstble. Further the trjectories in the ( u, u ) plne re given by u u 0 = u u 0 Equilibrium stte II ( exists while is wshed out): We hve = ; = 0 Substituting (.3) nd (.4) in (.) nd (.), we get u = u u + u u + u = u + uu + After lineriztion, we get u = u + (.9) = The chrcteristic eqution is ( λ + ) λ u { -[ ]} (.0) = 0 (.) One root of this eqution (.) is λ = which is negtive. 759

B. Rvindr Reddy Adv. Appl. Sci. Res., 0, 3():757-764 Cse A: When >, the other root of eqution (.) is λ = The trjectories re given by where u = γ - u0 e λ t + { u0γ - u 0 } e t λ t which is positive. Hence the equilibrium stte is unstble. (.) u = u0e (.3) γ = [ + ] The solution curves re illustrted in figures 3&4. u < u, we hve CASE : For 0 0 In this cse the second species is noted to be going wy from the equilibrium point while the first species would become extinct t the instnt + t * = 0 0 (λ ) ln u u γ u As such the stte is unstble. u > u, we hve 0 CASE : If 0 0 Here the first species out numbers the second species till the time, t = t * = u0 (λ + ) u0 ln λ + u0[ (λ + ) ] there fter the second species out numbers the first species. And lso the second species is noted to be going wy from the equilibrium point while the first species would become extinct t the instnt u t * = 0 u0 γ ln u 0 (λ + ) As such the stte is unstble. Cse B: When < One root of the eqution (.) is λ = which is negtive. λ = which is negtive nd the other root is As the roots of the eqution (.) re both negtive, the equilibrium stte is stble. The trjectories in this cse re the sme s in (.) nd (.3). 760

B. Rvindr Reddy Adv. Appl. Sci. Res., 0, 3():757-764 Tht is where u = γ λ u = u e t 0 - u0 e λ t + { u0γ - u 0 } e t γ = [ + ] The solution curves re illustrted in figures 5&6. u < u i.e. initilly the second species domintes the first species. CASE : 0 0 In this cse the second species lwys out numbers the first species. It is evident tht both the species converging symptoticlly to the equilibrium point. Hence the stte is stble. u > u i.e. initilly the first species domintes the second species. CASE : 0 0 Here the first species out numbers the second species till the time, t = t * = u0 (λ + ) u0 ln λ + u0[ (λ + ) ] there fter the second species out numbers the first species. As t both u & u pproch to the equilibrium point. Hence the stte is stble. Further the trjectories in the ( u, u ) plne re given by where ( q ) u q = = q cu + pu p = nd c is n rbitrry constnt. The solution curves re illustrted in figure 7. Equilibrium stte III (coexistence stte): We hve = ; ; > 0 = wherein Substituting (.3) nd (.4) in (.) nd (.), we get = u + u u u + u After lineriztion, we get nd = u + u u = u + u (.4) 76

B. Rvindr Reddy Adv. Appl. Sci. Res., 0, 3():757-764 u = (.5) The chrcteristic eqution is λ λ + = 0 One root of this eqution is positive nd the other root is negtive. Hence the equilibrium stte is unstble. The trjectories re given by u = u0λ + u0 λ λ u ( λ + ) + u λ λ u = 0 0 e λ t u0λ + u0 + e λ t λ λ t e λ + u ( λ + ) + u λ λ 0 0 e λ t The curves re illustrted in figures 8&9. Cse : u 0 > u 0 i.e. initilly the first species domintes the second species. In this cse, the first species is noted to be going wy from the equilibrium point while the second species pproches symptoticlly to the equilibrium point. Hence the stte is unstble. Cse : u 0 < u 0 i.e. initilly the second species domintes the first species. In this cse the second species out numbers the first species till the time, t = t * = λ λ ( b - λ ) u + ( b ) u ln 0 3 0 ( b λ ) u0 + ( 4 b ) u0 where b = ; b 3 = ; = λ + ; 4 = λ + fter tht the first species out numbers the second species nd grows indefinitely while the second species pproches symptoticlly to the equilibrium point. Hence the stte is unstble. Further the trjectories in the ( u, u ) plne re given by v ( )( v v ) ( u u v ) [ ] u d = v ( u vu ) where v nd v re roots of the qudrtic eqution c = nd d is n rbitrry constnt. v +bv+ c =0 with = ; ; b= 76

B. Rvindr Reddy Adv. Appl. Sci. Res., 0, 3():757-764 Trjectories: REFERECES [] Archn Reddy R, On the stbility of some mthemticl models in biosciences- intercting species, Ph.D thesis, 009, JTU. [] Cushing JM, Integro differentil equtions nd Dely Models in Popultion Dynmics, Lecture otes in Biomthemtics, Springer Verlg, Heidelberg, 977, Vol. 0. [3] Freedmn HI, Deterministic Mthemticl Models in Popultion Ecology, Mrcel Decker, ew York, 980. [4] George F Simmons, Differentil Equtions with pplictions nd historicl notes, Tt McGrw-Hill, ew Delhi, 974. [5] Kpur J, Mthemticl Modelling, Wiley Estern, 988. [6] Kpur J, Mthemticl Models in Biology nd Medicine Affilited Est West, 985. [7] Lkshmi ryn K, A Mthemticl study of Prey-Predtor Ecologicl Models with prtil covers for the prey nd lterntive food for the predtor, Ph.D thesis, J..T.University, 004. [8] Lotk AJ, Elements of Physicl biology, Willims nd Wilkins, Bltimore, 95. [9] Meyer WJ, Concepts of Mthemticl Modelling, Mc Grw Hill, 985. [0] Pul Colinvux, Ecology, John Wiley nd Sons Inc., ew York, 986. 763

B. Rvindr Reddy Adv. Appl. Sci. Res., 0, 3():757-764 [] Rvindr Reddy B, Lkshmi ryn K nd Pttbhirmchryulu Ch, Interntionl J.of Engg. Reserch & In. Appls, Vol., o.ii (009), 8-9. [] Rvindr Reddy B, Lkshmi ryn K nd Pttbhirmchryulu Ch, Interntionl J. of Mth. Sci & Engg. Appls. (IJMSEA), Vol. 4, o. III (August, 00), pp.97-06. [3] Srinivs C, Some Mthemticl spects of modeling in Bio Medicl Sciences, Ph.D thesis, Kktiy University, 99. [4] Volterr V, Leconssen l theorie mthemtique de l leitte pou lvie, Guthier Villrs, Pris, 93. [5] Srilth R, Rvindr Reddy B nd Pttbhirmchryulu Ch, Advnces in Applied Science Reserch, 0, (3): 5-65. [6] Srilth R nd Pttbhirmchryulu Ch, Advnces in Applied Science Reserch, 0, (3): 66-78. [7] Hri Prsd B nd Pttbhirmchryulu Ch, Advnces in Applied Science Reserch, 0, (5): 97-06. [8] Shiv Reddy K nd Pttbhi Rmchryulu Ch, Advnces in Applied Science Reserch, 0, (3): 08-8. 764