Talk:Lotka–Volterra equations

message: Deleted a link to a computer game clan - Some guy.
Help! Can anyone give me the equations for Volterra-Lotka equations for three species.
What I'm talking about is
Fox eats Rabbit
Rabbit eats Carrot
Carrot eats sunlight (infinite supply of sunlight)
Sure, I can give lots of them. Let C, R, F denote mass of Carrots, number of Rabbits, and number of Foxes. For each of C, R, and F I will give several equations, increading in realisticness. All parameters are >= 0.
Carrot equation:
C1.    C' = a - bR.  This assumes the number of carrots is fixed ( they do not get killed, or propagate),  a is the total growth due to sunlight, and b is how much is eaten by a rabbit.
C2.    C' = aC - bR.  This assumes a = p - d is the difference of the propagation rate and the 'old age' death rate for carrots. This implies unlimited exponential growth in the abscence of rabbits.
C3.    C' = a ( 1 - C/maxC ) C - bR.  This is called a 'logistic' equation, it assumes a 'maximum mass of carrots' called MaxC. This is appropriate in the likely event that you are dealing with a finite area for your ecology experiment. This ODE gives an 'S shaped curve', rapidly growing from 0, and slowly approaching 1. This is the best equation to use, unless you are looking for an easy equation. This ignores the effect that if carrots get low, the rabbits will have a hard time finding them. if you want to deal with that effect, see the F' equations below for ideas. In many situations, the foxes keep the rabbits from eating too many carrots.
Rabbit equations. Most Lotka-Volterra models use an absurd simple Rabbit - fox interaction term of R*F, which would indicate each fox eats a fixed proportion of all the rabbits. The only good thing about that is it leads to the simplest nonlinear term.
R1.     R' = cR - d RF.  This is the standard lame equation ( ignoring for the moment rabbit food).
R2.     R' = R [ -e + f(C/R) ] / [ 1 + g(C/R) ] = R [ -eR + fC ] / [ R + gC ].  This is a pretty good equation for the rabbit food aspect, ignoring the foxes for the moment. Here -e indicates how fast rabbits starve if the ratio C/R goes to zero, and f/g indicates how fast rabbit population grows when there are plenty of carrots to go around. For a fixed f/g, you can fiddle with f and g to better fit whatever data you have.
R3.     R' = R [ -eR + fC ] / [ R + gC ] - hF.  Here we assume the foxes have no trouble finding rabbits, and each one eats h rabbits per unit time.
R4.     R' = R [ -eR + fC ] / [ R + gC ] - hRF / [ kF + R ].  Here each fox eats h rabbits if there are plenty to go around, and 0 as the R/F ratio goes to zero. k is a free parameter to help fit whatever curve you thing is best.
F1      F' = F [ - m + pR ]. Traditional lame equation.
F2.     F' = F [ -mF + pR ] / [ F + nR ].   Here -m is the death rate due to starvation when there are no rabbits around.  p/n is the birth - (old age) death rate for well fed rabbits. p and n can be adjusted to fit whatever. For example, if you have a guess for what the R/F ratio is that keeps the Fox population constant, choose p so that m/p = stable(R/F).
equation summary: A good balance between reality and simplicity is
C' = a ( 1 - C/maxC ) C - bR
R' = R [ -eR + fC ] / [ R + gC ] - hRF / [ kF + R ]
F' = F [ -mF + pR ] / [ F + nR ]
Before working with this set of equations, simplify by 'scaling away parameters'. You can scale four paramiters (C, R, F, and time) which will allow you to eliminate 4 of the 11 parameters. The simplist is to measure C in units of maxC, so that maxC becomes 1. There are many ways to simplify the remaining equations. For example, letting replacing R with gR and F with (g/k)F removes the parameters in the denominators of the R' equation. Scaling time by a eliminates the a. This leads to ( the parameters will now have different values ):
C' = C ( 1 - C ) - b R
R' = R [ -eR + fC ] / [ R + C ] - hRF / [ F + R ]
F' = F [ -mF + pR ] / [ F + nR ]
Thus you have a three dimensional ODE with seven parameters.
The first step in the analysis is to find the (nontrivial) critical points, which are values of (C,R,F) where all right hand sides are zero.
This will require
R = C(1-C)/b,
F = pR/m,
[-eR + fC ] / [ R+C ] = hF / [F+R ].
subbing the second eqn into the third:
-eR + fC = [ hp / (m+p) ] [ R + C ].
subbing in the first equation:
(m+p) [ (bf-e)C + eC^2 ] = hp [ ( 1+b)C - C^2 ].
Dividing thru by C:
(m+p) [ (bf-e) + eC ] = hp [ ( 1+b) - C ]
which leads to
C = [ ( e - bf ) ( m + p ) + ( 1 + b ) hp ] / [ ( m + p ) e + hp ].
From this you can find R, and then F.
Obviously the parameters ( b,e,f,h,m,p ) will have to be such that C, R, and F are > 0.
The parameter n does not affect the critical point. Interesting.
From this point on you will probably need computer assistance.
Next you classify this critical point. Find the three eigenvalues of the frechet derivative (this is a three by three matrix of partial derivatives) evaluated at the CP. If you are looking for a cyclical population situation, you probably need a complex conjugate pair with positive real part, or maybe a pure imaginary complex pair, and the third eigenvalue a positive real.
Finding the regions in (b,e,f,h,m,p) space where these two events happen would be nice, but likely involve a lot of UGLY work, maybe Monte Carlo simulation. (BTW be sure to use a GOOD random number generator. Contact me for advice.  It might be a good time to find some realistic value ranges for  b, e, f, h, m, n, and p.
Ralph Kelsey
--132.235.14.20 08:08, 5 Nov 2004 (UTC)
replies: text: You should be holding  this discussion at Competitive Lotka–Volterra equations. Cuzkatzimhut (talk) 15:13, 24 June 2015 (UTC)
thread_title: Old discussions

message: The
Lotka-Volterra equations are not necessarily the ones given in the
article now, which form the Lotka-Volterra predator-prey model.
There is also the well known Lotka-Volterra competition model. They
have in common that the per capita growth rates are linear functions of
population densities or sizes. Which of these two models deserves the
name Lotka-Volterra sec is unclear to me, and if nobody justifies
the choice for the predator-prey model, I think it should be changed.
On
the other hand, maybe mathematicians refer to the predator-prey model
with "the Lotka-Volterra equations", while ecologists call it the
Lotka-Volterra predator-prey model.
Also, the name of the article is a singular, while there are clearly two equations.
OpenScience (talk) 12:13, 6 February 2010 (UTC)
thread_title: Are these really "the" Lotka-Volterra equations?

message: Is there a discrete version? I would expect it to be able to produce  chaotic results.--Henrygb 16:35, 17 Mar 2005 (UTC)
thread_title: Discrete version

message: There is only approximate solution revealed in this article. But what about strong analytic solution? I know (thanks for Maple) that it is expressed in terms of LambertW function. But it may be interesting if this solution would be referred in article not in so vaguely manner. 216.66.24.58 07:17, 20 November 2005‎ 
:I disagree the numerical solutions are "only approximate". The "analytic" Lambert W function solution is a hardly superior implicit renaming, and it does not shed that much light, unless one understands Lambert   functions well, in which case it is self-evident. 
:Absorbing the superfluous parameters into the variables as described in the text, retaining only the essential ratio , the solved parametric equation for constant  is  already provided in the article, 
:-y e^{-y}= - K (e^{x}/x)^{\gamma/\alpha}, so that    y=-W( - K (e^{x}/x)^{\gamma/\alpha}). 
:Without mastery of the nasty properties of the Lambert  for negative real arguments, this is all but useless here, and not suitable for inclusion. Cuzkatzimhut (talk) 18:57, 23 June 2015 (UTC)
thread_title:  Analytic solution 

message: In the analysis of the second EP, the terminology and tools used are good for linear systems, not for nonlinear systems such as the LV system.
~~
That's correct and I noticed it too.  All the complex eigen values reveal is that the fixed point is periodic, but does not reveal if it is a center or a stable spiral.
You have to eliminate the time dependance by dividing the dy/dt equation by dx/dt equation, which by the chain rule gives a single seperable equation for dy/dx.  When you seperate this and integrate you are left with an equation in x and y which you can relate to a third variable (magnitude) giving level sets of the origional equations.  Since f(x) and f(y) approach 0 as x&y approach infinity, the level curve is closed.  This translates into the closed orbits.
Vortmax 22:46, 2 December 2006 (UTC)
thread_title:  Little problem with second EP 

message: The article has the text:
In this equation, δxy represents the growth of the predator population. 
This is misleading in (at least) two respects.
Most simply, it seems to be saying that the Delta-x-y term represents the totality of change in y. This is clearly in error, since there is a second term as well. However, this problem is certainly just a matter of carelessness. Presumably the writer meant that the Delta-x-y term represents the rate of births of predators.
More troubling is question, why the birth rate is proportional to meetings between predator and prey. Wouldn't it be more logical to say that births among predators are proportional to meetings between male and female predators, i.e. proportional to y2? However, this would seem to be only appropriate for modelling in cases where predator-predator encounters are generally infrequent (in comparison with the period of gestation). For instance, the predators may be an endangered species. If predator-predator encounters are very frequent, we could ignore this term. 
Or maybe the x-y term is meant to show the dependency on food supply. (This would make sense, else why introduce the relationships predator and prey). In this case, the xy proportionality of the term would seem to say that predators on average have some characteristic number of progeny and that this characteristic number is proportional to the amount of available food. (But since we are talking of an average birth rate for the entire predator population, shouldn't that be the amount of available food per predator, ie y/x?) However, this explanation only appears (approximately) valid for moderate ranges of populations: As I've already mentioned we ought to include a y2 factor if y is small. Similarly, we ought to model a cap on the birth rate in case of very abundant prey. 
In any case, I think I've succeeded in making my point that the applicability for these equations needs to be addressed carefully and thorough if intelligent newcomers to the subject (I include myself) are to be able to learn from the article.
--Philopedia 19:32, 4 January 2007 (UTC)
thread_title:  Confusing; needs explication 

message: The initial conditions are 80 baboons and 40 cheetah, but
what are the parameters α, β, γ and δ ?
--141.53.37.95 10:59, 13 August 2007 (UTC)
replies: text: What about the black one? 98.245.92.246 (talk) 08:31, 30 April 2010 (UTC)
thread_title: Example: baboon (prey) and cheetah (predator)

message: A. Lotka originally proposed this equation in 1910 in chemistry and extended it in 1920 to biological interactions, most definitely not as late as 1925(cited two different ways). This article could do with at lot more citations and not just some vague ref list at the end. Kae1is (talk) 00:42, 9 March 2010 (UTC)
thread_title:  Proposal of Equation and citations 

message: In case anyone wants to know this: The Lotka-Volterra equations are analytically solvable.
Here are some papers:
.
At some point someone should put these information into the article --Svebert (talk) 20:41, 24 October 2011 (UTC)
replies: text: Hi,you provide two references, for the exact solutions. The first reference has a broken link, and the second one links to an article by 'Martin J. Gander' that gives an analytical solution without derivation , and cites: A. Steiner and M. Arrigoni, \Die Losung gewisser Rauber-Beute-Systeme", Studia Biophysica vol. 123 (1988) No. 2. I cannot find this last reference by searching online. Can you please provide another reference? Ami-de-la-Terre (talk) 07:19, 12 October 2021 (UTC)
thread_title:  Exact solutions 

message: An opinion from Wikiproject Mathematics suggests that information from the above article be merged into this article.  I don't know how  to do that without losing attribution, so I am leaving a note here in the hopes that an interested editor will take on this task.  —Anne Delong (talk) 23:18, 11 November 2013 (UTC)
replies: text: The Atto-Fox material is now in the history of this redirect: Atto-fox problem.  —Anne Delong (talk) 02:53, 26 June 2014 (UTC)
thread_title:  [[Wikipedia talk:Articles for creation/Atto-fox problem]] 

message: I can't believe this title said "Lotka–Volterra equation", with no final "s".  I changed it.  This is not about a type of equation of which any particular instance is a "Lotka–Volterra equation", nor about just one equation.  In either of those situations, the singular title would be required under WP:MOS.  It's about a particular pair of equations that constitute a system of equations. Michael Hardy (talk) 20:54, 29 September 2015 (UTC)
thread_title:  plural!! 

message: Hello fellow Wikipedians,
I have just modified one external link on Lotka–Volterra equations. Please take a moment to review my edit. If you have any questions, or need the bot to ignore the links, or the page altogether, please visit this simple FaQ for additional information. I made the following changes:
*Added archive https://web.archive.org/web/20100531204042/http://entomology.wsu.edu/profiles/06BerrymanWeb/Berryman%2892%29Origins.pdf to http://entomology.wsu.edu/profiles/06BerrymanWeb/Berryman%2892%29Origins.pdf
When you have finished reviewing my changes, you may follow the instructions on the template below to fix any issues with the URLs.
Cheers.—InternetArchiveBot (Report bug) 10:55, 11 December 2017 (UTC)
thread_title:  External links modified 

message: I thought it would be a good idea to update the "baboon/cheetah" example to use svg figures. I can make them and propose the changes, but wanted to float the idea here first. --Eitanlees (talk) 13:10, 12 June 2018 (UTC)
replies: text: The existing figures are poor quality, better ones would be great.  I don't see any reason at all to seek preapproval.  --JBL (talk) 18:26, 12 June 2018 (UTC)
thread_title:  Updating Figures 

message: I reorganised the first section to bring various important parts mentioned downstream together, and made the text more relevant to a non-mathematical audience. I  have added a sub-section on the biological relevance of the model by providing some examples how the behaviour of the model relates to biology (in part by replacing previous text and figure). I also have replaced the words "baboon" and "cheetah" with "rabbit" and "fox". This is in the first place to keep to one example: the baboons didn't add anything (there is nothing baboon specific in the example), and now the predator ties in with the atto-fox. Also because the cheetah is not a main predator of the baboon. I hope this improves the quality of the article. --Pips15 (talk) 14:46, 20 April 2023 (UTC)
thread_title: Reorganised top section

message: I find it confusing that this page is "Lotka-Volterra equations" but only considers predator-prey, while the  competitive and  generalized equations have separate pages. I'd say a single page for all, or otherwise simply a retitling of this page to specific "Predator-prey" (and perhaps a clearer cross-reference in the intro) would resolve. Kyle MoJo (talk) 05:12, 17 February 2024 (UTC)
thread_title:  Proposal: Retitle as "Predator-Prey" or merge with other Lotka-Volterra equations