[R] deSolve output

Hi, Try with the following ODE function. This should give you an extra column with the derivative of G in your THAAC matrix. degradation = function (t, state, parameters) {  with(as.list(c(state, parameters)),    {dG = (-a*(t+i)^b)*(G)    list(c(dG),dG=dG)    }) } Any additional variables that y

[R] deSolve output

Hi, I've solved a simple differential equation describing the degradation of amino acid carbon (THAA-C) using deSolve. Code is a follows: # Input of model parameters, a and b describes form of curve, i is apparent initial age of Org. C. parameters <- c(a = a, b = b, i=i) # Initial val