diagnostics {deSolve} | R Documentation |
Prints several diagnostics of the simulation to the screen, e.g. number of steps taken, the last step size,...
diagnostics(obj)
obj |
is the output matrix as produced by one of the integration routines. |
When the integration output is saved as a data.frame, then the required
attributes are lost and method diagnostics
will not work anymore.
none
## The famous Lorenz equations: chaos in the earth's atmosphere ## Lorenz 1963. J. Atmos. Sci. 20, 130-141. chaos <- function(t, state, parameters) { with(as.list(c(state)),{ dx <- -8/3*x+y*z dy <- -10*(y-z) dz <- -x*y+28*y-z list(c(dx, dy, dz)) }) } state <- c(x = 1, y = 1, z = 1) times <- seq(0, 50, 0.01) out <- vode(state, times, chaos, 0) pairs(out,pch=".") diagnostics(out)