search for: thaac

...eters)), #differential equation to be solved {dG = (-a*(t+i)^b)*(G) list(c(dG)) }) } # Makes sequence of times to be solved for. Start at t0 and proceed to iteno in steps of 1 day times = seq(t0, iteno, by = 1) #calls solver, to solve diffrential equation and saves results as variable THAAC THAAC = ode (y = state, times = times, func = degradation, parms = parameters) This all works fine, and output THAAC contains carbon concentrations of G at all times t. However, i would like the model to give me the rates, i.e. dG for all times. And possibly also the value of -a*(t+i)^b for all...