An Extensive Study of Approximating the Periodic. Solutions of the Prey Predator System

Similar documents
Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

A Prey-Predator Model with an Alternative Food for the Predator, Harvesting of Both the Species and with A Gestation Period for Interaction

10. The Discrete-Time Fourier Transform (DTFT)

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

LINEAR DELAY DIFFERENTIAL EQUATION WITH A POSITIVE AND A NEGATIVE TERM

Homotopy perturbation technique

A Propagating Wave Packet Group Velocity Dispersion

PROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

Rational Approximation for the one-dimensional Bratu Equation

2008 AP Calculus BC Multiple Choice Exam

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

Hydrogen Atom and One Electron Ions

Exercise 1. Sketch the graph of the following function. (x 2

DIFFERENTIAL EQUATION

5.80 Small-Molecule Spectroscopy and Dynamics

Finite element discretization of Laplace and Poisson equations

The Matrix Exponential

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation

ON A SECOND ORDER RATIONAL DIFFERENCE EQUATION

The Matrix Exponential

Inference Methods for Stochastic Volatility Models

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

Introduction to Condensed Matter Physics

Legendre Wavelets for Systems of Fredholm Integral Equations of the Second Kind

Limiting value of higher Mahler measure

On the irreducibility of some polynomials in two variables

Content Skills Assessments Lessons. Identify, classify, and apply properties of negative and positive angles.

[1] (20 points) Find the general solutions of y y 2y = sin(t) + e t. Solution: y(t) = y c (t) + y p (t). Complementary Solutions: y

4. (5a + b) 7 & x 1 = (3x 1)log 10 4 = log (M1) [4] d = 3 [4] T 2 = 5 + = 16 or or 16.

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

Einstein Equations for Tetrad Fields

Addition of angular momentum

cycle that does not cross any edges (including its own), then it has at least

Derangements and Applications

Recall that by Theorems 10.3 and 10.4 together provide us the estimate o(n2 ), S(q) q 9, q=1

1 General boundary conditions in diffusion

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Differential Equations

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

Chapter 1. Chapter 10. Chapter 2. Chapter 11. Chapter 3. Chapter 12. Chapter 4. Chapter 13. Chapter 5. Chapter 14. Chapter 6. Chapter 7.

Where k is either given or determined from the data and c is an arbitrary constant.

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

Lecture 2: Discrete-Time Signals & Systems. Reza Mohammadkhani, Digital Signal Processing, 2015 University of Kurdistan eng.uok.ac.

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

4037 ADDITIONAL MATHEMATICS

Problem Set 6 Solutions

(1) Then we could wave our hands over this and it would become:

Evaluating Reliability Systems by Using Weibull & New Weibull Extension Distributions Mushtak A.K. Shiker

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

MA 262, Spring 2018, Final exam Version 01 (Green)

Setting up the problem. Constructing a Model. Stability Analysis of ODE. Cheyne-Stokes. A Specific Case of the ODE 7/20/2011

(Upside-Down o Direct Rotation) β - Numbers

Title: Vibrational structure of electronic transition

Einstein Rosen inflationary Universe in general relativity

Southern Taiwan University

EXST Regression Techniques Page 1

2.3 Matrix Formulation

Section 11.6: Directional Derivatives and the Gradient Vector

Engineering Differential Equations Practice Final Exam Solutions Fall 2011

INTEGRATION BY PARTS

Supplementary Materials

BSc Engineering Sciences A. Y. 2017/18 Written exam of the course Mathematical Analysis 2 August 30, x n, ) n 2

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Complex Powers and Logs (5A) Young Won Lim 10/17/13

Construction of asymmetric orthogonal arrays of strength three via a replacement method

Note If the candidate believes that e x = 0 solves to x = 0 or gives an extra solution of x = 0, then withhold the final accuracy mark.

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Higher order derivatives

2F1120 Spektrala transformer för Media Solutions to Steiglitz, Chapter 1

Calculus II (MAC )

Optics and Non-Linear Optics I Non-linear Optics Tutorial Sheet November 2007

Basic Polyhedral theory

Hardy-Littlewood Conjecture and Exceptional real Zero. JinHua Fei. ChangLing Company of Electronic Technology Baoji Shannxi P.R.

nd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1).

Random Process Part 1

DISTRIBUTION OF DIFFERENCE BETWEEN INVERSES OF CONSECUTIVE INTEGERS MODULO P

Deift/Zhou Steepest descent, Part I

ABEL TYPE THEOREMS FOR THE WAVELET TRANSFORM THROUGH THE QUASIASYMPTOTIC BOUNDEDNESS

Sec 2.3 Modeling with First Order Equations

Properties of Phase Space Wavefunctions and Eigenvalue Equation of Momentum Dispersion Operator

The Equitable Dominating Graph

1973 AP Calculus AB: Section I

ON THE DISTRIBUTION OF THE ELLIPTIC SUBSET SUM GENERATOR OF PSEUDORANDOM NUMBERS

That is, we start with a general matrix: And end with a simpler matrix:

Mathematics. Complex Number rectangular form. Quadratic equation. Quadratic equation. Complex number Functions: sinusoids. Differentiation Integration

SAFE HANDS & IIT-ian's PACE EDT-15 (JEE) SOLUTIONS

EEO 401 Digital Signal Processing Prof. Mark Fowler

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

Addition of angular momentum

Chapter 10. The singular integral Introducing S(n) and J(n)

( ) Differential Equations. Unit-7. Exact Differential Equations: M d x + N d y = 0. Verify the condition

Calculus II Solutions review final problems

First order differential equation Linear equation; Method of integrating factors

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

On spanning trees and cycles of multicolored point sets with few intersections

surface of a dielectric-metal interface. It is commonly used today for discovering the ways in

On the optimality of a general production lot size inventory model with variable parameters

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

VTU NOTES QUESTION PAPERS NEWS RESULTS FORUMS

Transcription:

pplid athmatical Scincs Vol. 00 no. 5 5 - n xtnsiv Study of pproximating th Priodic Solutions of th Pry Prdator Systm D. Vnu Gopala Rao * ailing addrss: Plot No.59 Sctor-.V.P.Colony Visahapatnam 50 07 Dsignation: Dpt. of pplid athmatics GI Collg of Scinc GI Univrsity Visahapatnam.P. India * corrsponding author. -mail: dvgrgitam@gmail.com P. Sri Hari Krishna Dpt of pplid athmatics GI Collg of Scinc GI Univrsity Visahapatnam.P. India bstract In this papr an approximation to th priodic solutions of th gnral Lota Voltrra pry prdator systm is obtaind using th Poincar Lindst mthod. In [] th proprtis of gnral solutions wr obtaind without considration of ral world complxitis. Hnc w analyz th systm aftr modifying it to bttr approximat th ral lif scnario and prov that th wor don by [] is a particular cas of our modl. athmatics subjct classifications: B5 D0 799 Kywords: Poincar - Lindst prturbation Lota Voltrra Pry Prdator. Introduction h mathmatical study of pry prdator systm in population dynamics has bn th subjct of svral rcnt paprs starting with th wor of Lota Voltrra. In th study of non linar systm of diffrntial quations such as th Lota Voltrra

5 D. Vnu Gopala Rao and P. Sri Hari Krishna quations analytical quations ar usually unnown. In this cas in ordr to analyz th bhaviour of th systm w usually rsort to intgral tchniqus such as prturbation tchniqu. Prturbation tchniqu dpnds on th xistnc of small or larg paramtrs in th non linar problms. Vrma Wilson and Burnsid [0] [] [] obtaind xact solutions of dx X a by d dy Y C X - d. d with various assumptions. K.N.urty and D.V.G.Rao [7] obtaind approximat analytical solutions of th quations dx a - b X - cy d X dy - a c X Y. d whr a b c a b ar positiv constants with out maing any of th abov assumptions on th constants.. Grozdanovsi and J.J. Shphrd [] considrd a simpl modl of pry prdator systm. and obtaind approximat priodic solutions. In th ral world to analyz th bhaviour of th solutions such simpl modl is of not much us. In som situations ach spcis may hav a ngativ fd bac rsulting in stabl coxistnc of both spcis. For our analysis w introduc ngativ fd bac for th pry population and dvlop approximat priodic solutions of th gnral Lota Voltrra pry prdator systm. using th Poincar Lindst mthod and analyz th bhaviour of th solutions nar th non trivial critical point. Hr it is shown that th rsults of [] ar particular cas of th prsnt wor. h phas plan rprsntation of th rsults is also givn. h gnral form of Lota Voltrra systm is considrd introducing ngativ fd bac for th pry population givn in. Introducing a chang in th systm by dfining th variabls c cc x t X yt Y a a c a b a ba t a. a c a For a stabl and positiv quilibrium point of th systm w rquir th inquality. a c < which implis 0 < <. a b hn th systm bcoms d xt xt x t x t yt

pproximating th priodic solutions 5 d y t yt xt yt. h non zro critical point of th systm is obtaind at xt yt. Nar th nontrivial quilibrium w introduc nw variabls that masur our distanc from. x t ut yt νt whr is small.5 s a rsult du t u t u t ν t u t ν t dv t u t u t ν t.. Priodic natur of th scond ordr quations. h systm of quations by taing 0 or ignoring th non linar trms in. du t u t ν t du t u t. combining th quations and writing thm as a scond ordr quation d ut dut ut 0. Lt ut b any solution of. considr a transformation ut mt yt whr m t yt ar twic continuously diffrntiabl functions. h computations ofu u and thir substitutions in. lads to m y m m y m m mt yt 0 hn quating th cofficint of y to zro it is sn that mt xp t hrfor yt satisfis a diffrntial quation y yt 0. whr is a constant and 0 < <. So it is concludd that if ut is a solution of. thn yt ut xp t is a solution of.. Similarly if yt is a solution of. thn ut yt xp t is a solution of. i... is oscillatory if and only if. is oscillatory. Hnc th priodic natur of. is provd.

5 D. Vnu Gopala Rao and P. Sri Hari Krishna. h Poincar Lindst mthod Introducing a tim scal wt whr w dpnds on w0. h non linar wdu systm. bcoms u u ν - u ν. d wd ν u u ν.. d and solutions can b xpctd as u u 0 u u.. ν ν 0 ν ν... w K w w.. Whr u i u i and ν i ν c ar priodic in. Now substituting ths xprssions into. and collcting li powrs of w obtain a squnc of diffrntial quations. Collcting th trms non containing w gt u o u o ν o ν o u o. quating th cofficint of on bothsids u wu o u uo ν u0ν 0 ν wν o u u. 0ν 0 Collcting th cofficint of on both sids w u o w u u u u0u ν u0ν uν 0 ν wν wν 0 u u0ν uν 0.5 h systm. has solutions u cos. 0 ν 0 sin cos Whr ar arbitary constants. By giving valus to initial conditions ar obtaind. On substituting th solutions. into. th following diffrntial quations ar obtaind. u u ν cos sin

pproximating th priodic solutions 55 w sin ν cos u [ sin cos ] w cos.7 Hr th following simpl rsult is usful for th systm of quations to hav a priodic solutions. dx For th systm y sin t B cos t highr harmonics dy x C sin t D cos t highr harmonics. to hav priodic solutions it is ncssary and sufficint that D 0 and B C 0. pplying th abov rsult. th systm.7 will hav priodic solutions only if w choos w 0 giving whr u [ Pcos Q sin R] [ cos U sin S] v.9 [ ] P 5 Q [ ] R S [ QK P ] U [ Q P ] Hr.9 is a particular solution of.7 and no complmntary solution is includd. Sinc th gnral solution of. will involv two arbitrary constants. Substituting.9 in to.5 and noting that w 0 w obtain

5 D. Vnu Gopala Rao and P. Sri Hari Krishna K u u θ ν [ cos θ sin cos θ sin θ ] W sin θ cos ν u [ cos θ sin θ w cos θ sin θ ] cos θ.0 whr [ S P R KQ] [ U Q P R ] [ P KQ] [ U Q P] [ QK { R P }] [ R U P U Q ] [ QK P ] [ P U U Q ] pplying th rsult. to.0 it obsrvd that to gt π priodic solutions w must choos w G. K whr G [ P K QK R KU ] and so w obtain th following particular solutions u sin KG K

pproximating th priodic solutions 57 cos } K KG K K { } [ sin N N K { } N N cos ] cos 9 ν sin. Whr ] 9 [ KQ U P P U KQ L } 9 { } { U P Q Q P U L 0 } { L L 0 } { L L [ ] N [ N Substituting u i v i & w i i in x u y ν w obtain cos { cos P t x [ G R Q { } sin G { }sin }] } Cos

5 D. Vnu Gopala Rao and P. Sri Hari Krishna [ } sin { N N } { Cos N N. a } sin { Cos t y } sin { U Cos S }cos [{ 9 }sin {.b G W. Whr is vry small amplitud of th dviations of x y from h rsults of. giv an approximation to th priod of ths orbits ] [ ω G Π Π.5 givs a highr ordr approximations to th priod for small amplitud. Put 0 U O S O R Q p thn L L G N N

pproximating th priodic solutions 59 W in.a.b and. thn w gt approximat priodic solutions of simpl modl q.. ar of th form. a x a cos sin Cos a a a Sm cos sin a Cos o a a y a sin sin cos a sin cos O a a W whr a 0 a Hr it is shown that th solutions of th modl. bcom a particular cas of th gnral Lota Voltrra quations with ngativ fd bac for th pry population... RSULS Using th approximations.a and.b w plot th approximat solutions of x and y for various initial conditions. Figur shows th rsults for Є 0.05 and 0. for. Figur shows th rsults for Є 0.5 and 0. for. Figur shows th rsults for Є 0.05 and 0.5 for. Figurs - show that th orbits associatd with our modl ar curvs that spiral in toward th critical point. If th valu of th amplitud Є is incrasd w obsrv distortion in th path of th spiral as shown in Fig whn compard to that in Fig. Similarly with incras in valu of at constant amplitud w obsrv a mard diffrnc in th path travrsd by th spiral.

0 D. Vnu Gopala Rao and P. Sri Hari Krishna.0.0.0 Y 0.9 0.9 0.9 0.9 0.9 0.9.0.0.0 X Fig : Plot of th approximat priodic solution for amplitud Є 0.05 0..... Y. 0. 0. 0. 0..... X Fig : Plot of th approximat priodic solution for amplitud Є 0. 5 0.

pproximating th priodic solutions..0.0.0 Y.0 0.9 0.9 0.97 0.9 0.99.0.0.0.0.05.0.07 X Fig : Plot of th approximat priodic solution for amplitud Є 0.05 0.5.0.0.0 X 0.9 0.9 0.9 0 0 0 0 0 50 0 t Fig : Plot of x and t for th approximat solution shown in Figur

D. Vnu Gopala Rao and P. Sri Hari Krishna.0.0.0 Y 0.9 0.9 0.9 0 0 0 0 0 50 0 t Fig 5: Plot of y and t for th approximat solution shown in Figur.0.0.0.0 XY 0.9 0.9 0.9 0 0 0 0 0 50 0 t Fig : Plot of th product xy and t for th approximat solution shown in Figur

pproximating th priodic solutions.0.0.0 XY 0.9 0.9 0.9 0 0 0 0 0 50 0 t Fig 7: Plot of xsolid lin and ydottd lin for th approximat priodic solution obtaind in Figur Figurs -7 compar th approximations for x y and xy obtaind from.a and.b. In Figur 7 th x approximation is qual in amplitud with th y approximation and hnc rsultd in th propr spiral shown in Figur. 5. Conclusions h approximat priodic solutions of th gnral Lota Voltrra pry prdator systm wr obtaind using th Poincar Lindst mthod aftr introducing ngativ fdbac in th pry population as th additional complxity thrby giving a bttr approximation to th ral lif scnario. W also provd that th wor don by [] is a particular cas of our modl. Rfrncs [] R.R. Burnsid not on xact solutions of th pry prdator quations. Bull. ath Biol. 9 9-97.

D. Vnu Gopala Rao and P. Sri Hari Krishna [] W.J. Cunningham Introduction to nonlinar nalysis. c Graw Hill Nw Yor 95. [] R. Grimshaw Nonlinar ordinary diffrntial quations. CRC prss US 99. []. Grozdanovsi and J.J. Shphrd pproximating th priodic solutions of Lota - Voltrra quations. ustralian mathmatical soc. 00. nziam J. 9 C- 007 PP C-C57 00. [5].J. Lota lmnts of physical biology. Williams & Wilins Batismor aryland US 95. [] J.D. urray athmatical Biology n introduction. Springr Vrlog Brlin Hidlbrg 00. [7] K.N. urty and D.V.G.Rao pproximat analytical solutions of th gnral Lota Voltrra quations. Journal of mathmatical analysis and applications 5-5 97. [].G. Pilov n introduction to mathmatical cology. Wily Intrscinc Nw Yor 99. [9] F. Roth h priods of th Voltrra Lota systms. J.Rsn. ngrw ath 55 9-95. [0] V.S. Varma xact solutions for a spcial pry-prdator compting spcis systm. Bull. ath. Biol 9 977 9-. [] V. Voltrra Variations and fluctuations of th numbr of individuals in animal spcis living togthr. ranslatd from 9 dition by R.N.Chapman. nimal cology rno Nw Yor US 9. [] J. Waldvogl h priod in th Voltrra - Lota prdator pry modl. SI J. Numrical nalysis 0-7 9. [] J. Waldvogl h priod in th Voltrra - Lota systm is monotonic. J. ath nalysis and pplications - 79-9. C C55. [].J. Wilson In Varma s pry-prdator problms. Bull. ath. Biol 599-00 90. Rcivd: pril 00