search for: tderiv

...ctor of initial conditions (see example below). I know that this is my fault, but it might be possible/useful to "user-proof" lsoda() some more (try func() once and see if it returns a list of the appropriate form, check the derivative vector for length matching the IC vector, etc ... tderiv <- function(t,x,parms) { # return(list(c(x[1]*parms[1],x[2]*parms[2],parms[3]*(x[1]+x[2])),NULL)) return(c(x[1]*parms[1],x[2]*parms[2],parms[3]*(x[1]+x[2]))) } lsoda(rep(1,2),c(0:20),tderiv,c(0.01,-0.01,1.2)) --please do not edit the information below-- Version: platform = i686-pc-linux...