Displaying 1 result from an estimated 1 matches for "lotvol".
Did you mean:
logvol
2004 Sep 06
1
A naive lsoda question....
Hello,
I am an R newbie, trying to use lsoda to solve standard
Lotka-Volterra competition equations. My question is: how do I
pass a parameter that varies with time, like say, phix <- 0.7 +
runif(tmax) in the example below.
# defining function
lotvol <- function(t,n,p){
x <- n[1]; y <- n[2]
rx <- p["rx"]; ry <- p["ry"]
Kx <- p["Kx"]; Ky <- p["Ky"]
phix <- p["phix"]; phiy <- p["phiy"]
dx.dt <- rx*x*(1 - x/Kx) - phix*x*y
dy.dt <- r...