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

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

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

The Dynamics of Two Harvesting Species with variable Effort Rate with the Optimum Harvest Policy

ON THE STABILITY OF DELAY POPULATION DYNAMICS RELATED WITH ALLEE EFFECTS. O. A. Gumus and H. Kose

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method

3. Renewal Limit Theorems

Contraction Mapping Principle Approach to Differential Equations

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration.

Think of the Relationship Between Time and Space Again

5.1-The Initial-Value Problems For Ordinary Differential Equations

Necessary and Sufficient Conditions for Asynchronous Exponential Growth in Age Structured Cell Populations with Quiescence

f t f a f x dx By Lin McMullin f x dx= f b f a. 2

Solutions to Problems from Chapter 2

A RECURSIVE PROCEDURE FOR COMMENSAL- HOST ECOLOGICAL MODEL WITH REPLENISHMENT RATE FOR BOTH THE SPECIES - A NUMERICAL APPROACH

1. Find a basis for the row space of each of the following matrices. Your basis should consist of rows of the original matrix.

Neural assembly binding in linguistic representation

GENERALIZATION OF SOME INEQUALITIES VIA RIEMANN-LIOUVILLE FRACTIONAL CALCULUS

An Integral Two Space-Variables Condition for Parabolic Equations

A Kalman filtering simulation

e t dt e t dt = lim e t dt T (1 e T ) = 1

A Simple Method to Solve Quartic Equations. Key words: Polynomials, Quartics, Equations of the Fourth Degree INTRODUCTION

RESPONSE UNDER A GENERAL PERIODIC FORCE. When the external force F(t) is periodic with periodτ = 2π

Probability, Estimators, and Stationarity

How to prove the Riemann Hypothesis

I = I = I for this case of symmetry about the x axis, we find from

1.0 Electrical Systems

The order of reaction is defined as the number of atoms or molecules whose concentration change during the chemical reaction.

On the Pseudo-Spectral Method of Solving Linear Ordinary Differential Equations

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points)

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

A Time Truncated Improved Group Sampling Plans for Rayleigh and Log - Logistic Distributions

Asymptotic relationship between trajectories of nominal and uncertain nonlinear systems on time scales

A LIMIT-POINT CRITERION FOR A SECOND-ORDER LINEAR DIFFERENTIAL OPERATOR IAN KNOWLES

15/03/1439. Lecture 4: Linear Time Invariant (LTI) systems

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

4.8 Improper Integrals

How to Prove the Riemann Hypothesis Author: Fayez Fok Al Adeh.

Hermite-Hadamard-Fejér type inequalities for convex functions via fractional integrals

New Inequalities in Fractional Integrals

Soliton Scattering on the External Potential in Weakly Nonlocal Nonlinear Media

Procedia Computer Science

A PREY-PREDATOR MODEL WITH AN ALTERNATIVE FOOD FOR THE PREDATOR AND OPTIMAL HARVESTING OF THE PREDATOR

Physics 2A HW #3 Solutions

LAPLACE TRANSFORM OVERCOMING PRINCIPLE DRAWBACKS IN APPLICATION OF THE VARIATIONAL ITERATION METHOD TO FRACTIONAL HEAT EQUATIONS

September 20 Homework Solutions

Temperature Rise of the Earth

Positive and negative solutions of a boundary value problem for a

Average & instantaneous velocity and acceleration Motion with constant acceleration

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment

Journal of Mathematical Analysis and Applications. Two normality criteria and the converse of the Bloch principle

Analytic solution of linear fractional differential equation with Jumarie derivative in term of Mittag-Leffler function

Minimum Squared Error

Minimum Squared Error

A new model for limit order book dynamics

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM

Bioeconomic modelling of a prey predator system using differential algebraic equations

0 for t < 0 1 for t > 0

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

Thermal neutron self-shielding factor in foils: a universal curve

..,..,.,

1. Introduction. 1 b b

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

UNIQUENESS OF AMMENSAL & COMMENSAL MATHEMATICAL MODELS WITH LIMITED RESOURCES IN FULLY WASHED-OUT STATE

Some Inequalities variations on a common theme Lecture I, UL 2007

Convergence of Singular Integral Operators in Weighted Lebesgue Spaces

CALDERON S REPRODUCING FORMULA FOR DUNKL CONVOLUTION

LAGRANGIAN AND HAMILTONIAN MECHANICS WITH FRACTIONAL DERIVATIVES

Scalar Conservation Laws

USING ITERATIVE LINEAR REGRESSION MODEL TO TIME SERIES MODELS

Tax Audit and Vertical Externalities

Physic 231 Lecture 4. Mi it ftd l t. Main points of today s lecture: Example: addition of velocities Trajectories of objects in 2 = =

Solitary Wave Solutions for the Boussinesq and Fisher Equations by the Modified Simple Equation Method

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008)

PARABOLA. moves such that PM. = e (constant > 0) (eccentricity) then locus of P is called a conic. or conic section.

ON NEW INEQUALITIES OF SIMPSON S TYPE FOR FUNCTIONS WHOSE SECOND DERIVATIVES ABSOLUTE VALUES ARE CONVEX

PHYSICS 1210 Exam 1 University of Wyoming 14 February points

Predict Global Earth Temperature using Linier Regression

1. Consider a PSA initially at rest in the beginning of the left-hand end of a long ISS corridor. Assume xo = 0 on the left end of the ISS corridor.

S Radio transmission and network access Exercise 1-2

A SPECIAL CASE OF ECOLOGICAL COMMENSALISM- PHASE PLANE ANALYSIS

ON THE OSTROWSKI-GRÜSS TYPE INEQUALITY FOR TWICE DIFFERENTIABLE FUNCTIONS

A Prey-Predator Model with an Alternative Food for the Predator and Optimal Harvesting of the Prey

Non-oscillation of perturbed half-linear differential equations with sums of periodic coefficients

MASS, STIFFNESS, AND DAMPING MATRICES FROM MEASURED MODAL PARAMETERS

IX.2 THE FOURIER TRANSFORM

LAPLACE TRANSFORMS. 1. Basic transforms

ON A NEW SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATION USING COMPLEX TRANSFORM IN THE UNIT DISK

A new model for solving fuzzy linear fractional programming problem with ranking function

Application on Inner Product Space with. Fixed Point Theorem in Probabilistic

2k 1. . And when n is odd number, ) The conclusion is when n is even number, an. ( 1) ( 2 1) ( k 0,1,2 L )

Chapter 2. Motion along a straight line. 9/9/2015 Physics 218

Chapter Direct Method of Interpolation

Fault-Tolerant Guaranteed Cost Control of Uncertain Networked Control Systems with Time-varying Delay

FURTHER GENERALIZATIONS. QI Feng. The value of the integral of f(x) over [a; b] can be estimated in a variety ofways. b a. 2(M m)

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

Hermite-Hadamard and Simpson Type Inequalities for Differentiable Quasi-Geometrically Convex Functions

graph of unit step function t

Green s Functions and Comparison Theorems for Differential Equations on Measure Chains

Transcription:

VOL. 6, NO. 4, APRIL 0 ISSN 89-6608 ARPN Jornl of Engineering nd Applied Sciences 006-0 Asin Reserch Pblishing Nework (ARPN). All righs reserved. www.rpnjornls.com A MATHEMATICAL MODEL OF FOUR SPECIES SYN-ECOSYMBIOSIS COMPRISING OF PREY-PREDATION, MUTUALISM AND COMMENSALISMS-I (FULLY WASHED OUT STATE) R. Srilh nd N. Ch. Pbhirmchryl Reserch Scholr, JNTUH, Kkplly, Hyderbd, Indi Professor (Red.) of Mhemics, NIT, Wrngl, Indi E-Mil: bsrilh8@gmil.com ABSTRACT This invesigion dels wih mhemicl model of for species (S, S, S nd S 4 ) Syn-Ecologicl sysem (Flly Wshed o Se). S is predor srviving on he prey S : he prey is commensl o he hos S which iself is in mlism wih he forh species S 4. S nd S 4 re nerl. The mhemicl model eqions chrcerizing he synecosysem consie se of for firs order non-liner copled differenil eqions. There re in ll sixeen eqilibrim poins. Crieri for he sympoic sbiliy of one of he sixeen eqilibrim poins: he flly wshed o se is esblished. The linerised eqions for he perrbions over he eqilibrim poin re nlyzed o esblish he crieri for sbiliy. The sysem is noiced o be loclly sble. Trjecories of he perrbions hve been illsred. Keywords: mhemicl model, species, syn-ecologicl sysem, mlism, commenslism, differenil eqions. INTRODUCTION Mhemicl modeling is n imporn inerdisciplinry civiy which involves he sdy of some specs of diverse disciplines. Biology, Epidemiodology, Physiology, Ecology, Immnology, Bio-economics, Geneics, Phrmocokineics re some of hose disciplines. This mhemicl modeling hs rised o he zenih in recen yers nd spred o ll brnches of life nd drew he enion of every one. Mhemicl modeling of ecosysems ws iniied by Lok [8] nd by Volerr [4]. The generl concep of modeling hs been presened in he reises of Meyer [9], Cshing [], Pl Colinvx [0], Freedmn [], Kpr [5, 6]. The ecologicl inercions cn be brodly clssified s prey-predion, compeiion, mlism nd so on. N.C. Srinivs [] sdied he compeiive eco-sysems of wo species nd hree species wih regrd o limied nd nlimied resorces. Ler, Lkshmi Nryn [7] hs invesiged he wo species prey-predor models. Recenly sbiliy nlysis of compeiive species ws invesiged by Archn Reddy []. Locl sbiliy nlysis for wospecies ecologicl mlism model hs been invesiged by B. Rvindr Reddy e l., [, ]. The presen invesigion is devoed o n nlyicl sdy of for species Syn-Ecologicl sysem. S is predor srviving on he prey S : he prey is commensl o he hos S which iself is in mlism wih he forh species S 4 ; S nd S 4 re nerl. Figre- shows he Schemic Skech of he sysem nder invesigion. The model eqions of he sysem consie se of for firs order non-liner ordinry differenil copled eqions. In ll he sixeen eqilibrim poins of he sysem re idenified nd he sbiliy nlysis is crried o only for he flly wshed o se. The linerized perrbed eqions over he eqilibrim ses re solved nd he rjecories illsred. Figre-. Schemic skech of he Syn Eco- sysem. BASIC EQUATIONS Noion doped N (): The Poplion of he Prey (S ) N (): The Poplion of he Predor (S ) N (): The Poplion of he Hos (S ) of he Prey (S ) nd ml o S 4 N 4 (): The Poplion of S 4 ml o S : Time insn,,, 4 : Nrl growh res of S, S, S, S 4,,, 44 : Self inhibiion coefficiens of S, S, S, S 4, : Inercion (Prey-Predor) coefficiens of S de o S nd S de o S : Coefficien for commensl for S de o he Hos S 4, 4 : Mlly inercion beween S nd S 4,,, 44 : Crrying cpciies of S, S, S, S 4 Frher he vribles N, N, N, N 4 re nonnegive nd he model prmeers,,, 4 ;,,, 44 ;,,, 4 re ssmed o be non-negive consns. 4

VOL. 6, NO. 4, APRIL 0 ISSN 89-6608 ARPN Jornl of Engineering nd Applied Sciences 006-0 Asin Reserch Pblishing Nework (ARPN). All righs reserved. www.rpnjornls.com The model eqions for he growh res of S, S, S, S 4 re d d d N N NN + NN. () N N + NN. () N N + NN. () 4 4 d 4 N N + NN. (4) 4 4 44 4 4 4 EQUILIBRIUM STATES The sysem nder invesigion hs sixeen eqilibrim ses defined by i 0, i,,, 4 d re given in he following ble. (5) S. No. Eqilibrim Ses Eqilibrim poin Flly Wshed o se N 0, N 0, N 0, N4 0 Only S 4 srvives 4 N 0, N 0, N 0, N4 44 Only he hos (S )of S srvives N 0, N 0, N, N 0 4 Only he predor S srvives 5 Only he prey S srvives 6 7 Prey (S ) nd predor (S ) wshed o Prey (S ) nd hos (S ) of S wshed o 8 Prey (S ) nd S 4 wshed o 9 Predor (S ) nd Hos (S ) of S wshed o 0 Predor (S ) nd S 4 wshed o Prey (S ) nd predor (S )srvives Only he prey (S ) wshed o Only he predor (S ) wshed o N 0, N, N 0, N 0 N, N 0, N 0, N 0 + + N 0, N 0, N, N 4 4 44 4 4 44 44 44 44 N 0, N, N 0, N 4 44 N 0, N, N, N 0 N, N 0, N 0, N 4 44 + N, N 0, N, N 0 + N, N, N 0, N 0 + + + + N 0, N, N, N 4 4 44 4 4 44 44 44 44 α + + N, N 0, N, N 4 4 44 4 4 α 44 44 44 44 where α ( 44 + 44 ) + ( 44 44) α ( ) 44 4 4 4

VOL. 6, NO. 4, APRIL 0 ISSN 89-6608 ARPN Jornl of Engineering nd Applied Sciences 006-0 Asin Reserch Pblishing Nework (ARPN). All righs reserved. www.rpnjornls.com 4 Only he Hos (S ) of S wshed o + N, N, N 0, N 4 + + 44 5 Only S 4 wshed o 6 The co-exisen se (or) Norml sedy se β β N, N, N, N 0 β β where β ( + ) β ( + ) β ( + ) + N N γ + γ γ + γ, N, 4 γ γ + +, N 4 4 44 4 4 4 44 44 44 44 where γ ( + )( ) 44 4 4 γ + 44 4 4 γ ( + )( ) 44 4 4 γ ( )( ) 4 4 44 4 The presen pper dels wih he flly wshed o se only. The sbiliy of he oher eqilibrim ses will be presened in he forh coming commnicions. STABILITY OF THE FULLY WASHED OUT EQUILIBRIUM STATE (Sl. No. in he bove Tble) To discss he sbiliy of eqilibrim poin N 0, N 0, N 0, N4 0 Le s consider smll deviions (), (), (), 4 () from he sedy se i.e., N ( ) N + ( ), i,,,4.. (6) i i i Where i () is smll perrbions in he species S i. Sbsiing (6) in (), (), (), (4) nd neglecing prodcs nd higher powers of,,, 4, we ge d d d d d d (7).. (8) (9) d d 4.. (0) 4 4 The chrcerisic eqion of which is ( λ )( λ )( λ )( λ ) 0.. () he roos,,, 4 of which re ll posiive. Hence he Flly Wshed o Se is nsble. The solions of he eqions (7), (8), (9), (0) re 0e () 0e () 0e. (4) 4 4 e. (5) Where 0, 0, 0, re he iniil vles of,,, 4 respecively. There wold rise in ll 576 cses depending pon he ordering of he mgnides of he growh res,,, 4 nd he iniil vles of he perrbions 0 (), 0 (), 0 (), () of he species S, S, S, S 4. Of hese 576 siions some ypicl vriions re illsred hrogh respecive solion crves h wold fcilie o mke some resonble observions. 44

VOL. 6, NO. 4, APRIL 0 ISSN 89-6608 ARPN Jornl of Engineering nd Applied Sciences 006-0 Asin Reserch Pblishing Nework (ARPN). All righs reserved. www.rpnjornls.com Cse (i): If 0 < 0 < 0 <, < < < 4 In his cse prey (S ) hs he les nrl birh re nd S 4 domines he prey (S ), predor (S ) nd he hos (S ) of S in nrl growh re s well s in is poplion srengh. Cse (iv): If 0<<0<0, <<<4 In his cse he hos (S ) of S hs he les nrl birh re. Iniilly he hos (S ) of S domines over S 4 nd he predor (S ) ill he imes insn * 4, * respecively. Therefer he dominnce is reversed. ; * * 0 4 4 0 0 Also he prey (S ) domines over S 4, Predor (S ) ill he ime insn * 4, * respecively nd herefer he dominnce is reversed. ; * * 0 4 4 0 0 Figre- Cse (ii): If 0 < 0 < 0 <, < < < 4 In his cse predor (S ) hs he les nrl birh re. Iniilly he predor (S ) domines over he prey (S ) ill he ime insn * 0 0 nd here-fer he prey (S ) domined he predor (S ). The * ime my be clled he dominnce ime of he predor (S ) over he prey (S ). Figre-4 Cse (v): If 0 < 0 < 0 <, < < 4 < In his cse he Hos (S ) of S hs he les nrl birh re. Iniilly i is domines over he predor (S ) ill he ime insn * 0 0 nd herefer he dominnce is reversed. Also S 4 domines over he prey (S ) ill he ime insn * 4 nd herefer he dominnce is 4 0 Figre- Cse (iii): If 0 < 0 < < 0, < < < 4 In his cse predor (S ) hs he les nrl birh re. Iniilly he predor (S ) domines over S 4, hos (S ) of S nd prey (S ) ill he ime insn,, respecively nd here fer he dominnce is * * * 4 reversed * * 0 * 0 4,, 4 0 0 0 Figre-5 45

VOL. 6, NO. 4, APRIL 0 ISSN 89-6608 ARPN Jornl of Engineering nd Applied Sciences 006-0 Asin Reserch Pblishing Nework (ARPN). All righs reserved. www.rpnjornls.com Cse (vi): If 0 < 0 < < 0, < 4 < < In his cse he prey (S ) hs he les nrl birh re. Iniilly he prey (S ) domines over is hos, S 4 nd Predor (S ) ill he ime insn *, * 4, * respecively nd herefer he dominnce is reversed. Also S 4 domines over he hos (S ) of S, nd he predor (S ) ill he ime insn * 4, * 4 nd herefer he dominnce is reversed. Similrly, he hos (S ) of S domines over he predor (S ) ill he ime insn * nd he dominnce ges reversed herefer. ; ; * 0 * * 0 4 0 4 0 0 Cse (viii): If 0 < 0 < 0 <, < < 4 < In his cse he predor (S ) hs he les nrl birh re. Iniilly he predor (S ) domines over he prey (S ), hos (S ) of S ill he ime insn *, * respecively nd herefer he dominnce is reversed. Also he prey (S ) domines over is hos ill he ime insn * nd herefer he dominnce is reversed. Similrly S 4 domines over he hos (S ) of S ill he ime insn * 4 he dominnce ges reversed fer. ; ; * * * 0 4 4 4 0 4 0 0 Figre-8 Figre-6 Cse (vii): If 0 < 0 < 0 <, < 4 < < In his cse he hos (S ) of S hs he les nrl birh re. Iniilly S 4 domines over boh he prey (S ) nd predor (S ) ill he ime insn * 4, * 4 respecively nd herefer he dominnce is reversed. Also he Prey (S ) domines over he Predor (S ) po he ime insn * nd he dominnce ges reversed fer. nd ; * * 4 4 4 0 4 0 * 0 0 ; * 0 * 0 0 0 ; * 0 * 4 0 4 0 Cse (ix): If 0 < 0 < < 0, < 4 < < In his cse he prey (S ) hs he les nrl birh re. Iniilly he prey (S ) domines over is hos, predor (S ) nd S 4 ill he ime insn *, *, * 4 respecively nd here fer he dominnce is reversed. Also S 4 domines over he predor (S ) nd he hos (S ) of S ill he imes insn * 4 nd * 4 nd herefer he dominnce is reversed. ; ; * 0 * 0 * 4 0 0 4 0 nd ; * * 4 4 4 0 4 0 Figre-7 46

VOL. 6, NO. 4, APRIL 0 ISSN 89-6608 ARPN Jornl of Engineering nd Applied Sciences 006-0 Asin Reserch Pblishing Nework (ARPN). All righs reserved. www.rpnjornls.com ; * * 0 4 4 0 0 ; * 0 * 4 0 4 0 nd * 4 4 0 Figre-9 Cse (x): If < 0 < 0 < 0, < < 4 < In his cse he prey (S ) hs he les nrl birh re. Iniilly he prey (S ) domines over S 4 ill he ime insn * 4 nd herefer he dominnce is reversed. Also he predor (S ) domines over S nd S 4 ill he ime insns *, * 4 respecively nd herefer he dominnce is reversed. * * 4 ; 4 nd 0 4 0 * 0 0 Figre- Cse (xii): If < 0 < 0 < 0, 4 < < < In his cse S 4 hs he les nrl birh re. Iniilly he prey (S ) domines over is hos nd predor (S ) ill he ime insn *, * respecively nd herefer he dominnce is reversed. * 0 * 0 ; 0 0 Figre-0 Cse (xi): If < 0 < 0 < 0, < < < 4 In his cse he predor (S ) hs he les nrl birh re. Iniilly he predor (S ) domines over S 4 nd prey (S ) ill he ime insn * 4, * respecively nd herefer he dominnce is reversed. Also he hos (S ) of S domines over he prey (S ) nd S 4 ill he ime insn *, * 4 respecively nd herefer he dominnce is reversed. Similrly he prey (S ) domines over S 4 ill he ime insn * 4 nd he dominnce ges reversed fer. Figre- TRAJECTORIES OF PERTURBATIONS The rjecories in he - plne given by ( ) ( ) nd re shown in Figre-. 0 0 47

VOL. 6, NO. 4, APRIL 0 ISSN 89-6608 ARPN Jornl of Engineering nd Applied Sciences 006-0 Asin Reserch Pblishing Nework (ARPN). All righs reserved. www.rpnjornls.com [4] George F. 974. Simmons: Differenil Eqions wih pplicions nd hisoricl noes. T McGrw- Hill, New Delhi. [5] Kpr J. N. 988. Mhemicl Modeling. Wiley - Esern. [6] Kpr J.N. 985. Mhemicl Models in Biology nd Medicine Affilied Es - Wes. Figre- Also he rjecories in he plne given by ( ) ( ) nd re shown in Figre-4. 0 0 Similrly he rjecories in he - 4, -, - 4, - 4 4 4 plnes re ( ) ( ),( ) ( ), 0 0 0 4 4 4 4 ( ) ( ),( ) ( ) respecively. 0 0 [7] Lkshmi Nryn K. 004. A Mhemicl sdy of Prey-Predor Ecologicl Models wih pril covers for he prey nd lernive food for he predor. PhD Thesis. J. N. T. Universiy. [8] Lok A. j. 95. Elemens of Physicl biology. Willims nd Wilkins, Blimore. [9] Meyer W.J. 985. Conceps of Mhemicl Modeling. McGrw - Hill. [0] Pl Colinvx. 986. Ecology. John Wiley nd Sons Inc., New York. [] Rvindr Reddy B., Lkshmi Nryn, K. nd Pbhirmchryl N.Ch. 009. A model of wo mlly inercing species wih limied resorces for boh he species. Inernionl J. of Engg. Reserch nd Ind. Appls. (II): 8-9. REFERENCES Figre-4 [] Rvindr Reddy B., Lkshmi Nryn K. nd Pbhirmchryl N.Ch. 00. A model of wo mlly inercing species wih limied resorces nd hrvesing of boh he species consn re. Inernionl J. of Mh. Sci nd Engg. Appls. (IJMSEA). 4(III): 97-06. [] Srinivs N.C. 99. Some Mhemicl specs of modeling in Bio Medicl Sciences. PhD Thesis. Kkiy Universiy. [4] Volerr V. 9. Leconssen l heorie mhemiqe de l leie po lvie. Ghier - Villrs, Pris. [] Archn Reddy R. 009. On he sbiliy of some mhemicl models in biosciences- inercing species. PhD Thesis. JNTU. [] Cshing J. M. 977. Inegro - differenil eqions nd Dely Models in Poplion Dynmics. Lecre Noes in Biomhemics. Springer- Verlg, Heidelberg. Vol. 0. [] Freedmn H. I. 980. Deerminisic Mhemicl Models in Poplion Ecology. Mrcel - Decker, New York. 48