Displaying 1 result from an estimated 1 matches for "dbadt".
2008 Nov 21
1
lsoda warning "too much accuracy requested"
...e.
Many thanks
Colleen
___________
My code is:
library(odesolve)
SI <- 80
model <- function(t, x, parms) {
H <- x[1]
BA <- x[2]
N <- x[3]
with(as.list(parms), {
dHdt <- (b/c)*(((a**c)*((H)**(1-c))-H))
dBAdt <- -(BA*b)*(c0+(c1*SI)-log(BA))/(log(1-((H/a)**c)))
dNdt <- N*alpha*(((log(1-((H/a)**c)))/b)**beta) - (gamma*BA)
list(c(dHdt, dBAdt, dNdt))
})
}
times <- seq(0, 40, 1)
parms <- c(a=(SI*1.258621)-1.32759, b=0.1, c=0.4, c0=4.6012,...