data(lv3)
plot(sim(lv3))
solver(lv3) <- "lsoda"
plot(sim(lv3)) # wrong! automatic time step overlooks internal inputs
plot(sim(lv3, hmax=1)) # integration with correct maximum time step
Run the code above in your browser using DataLab