Displaying 1 result from an estimated 1 matches for "fluxao".
Did you mean:
flux
2011 Aug 16
1
deSolve, extracting variable values from inside ode function
...2.389)
state <- c(QA=9.437)
MM <- function(t, state, parameters) {
with(as.list(c(state, parameters)),{
## concentration of substrate in pool A, quantity
## divided by the volume of the solvent pool
CA <- QA/vola
## output of substrate from pool A
fluxao <- vmax / (1 + (mmk/CA))
## rate of change of quantity of substrate
## in pool A
dA <- fluxoa - fluxao
## return the rate of change
list(dA)
})
}
times <- seq(0, 100, by = 5)
library(deSolve)
out1 <- ode(y = state, times = times, func = MM, parms = parameter...