library(pacman)
p_load(deSolve)


derivs = function(t, y, parms)
    list(c(y*0.03,y*0.04))

times = (seq(from=0, to=10, by=1))

out = ode(y=c(y=1000,y=10), times = times, func = derivs, parms = NULL)

head(out, n=10)