diagnostics {deSolve}R Documentation

Print Diagnostic Characteristics of ODE Solvers

Description

Prints several diagnostics of the simulation to the screen, e.g. number of steps taken, the last step size,...

Usage

diagnostics(obj)

Arguments

obj is the output matrix as produced by one of the integration routines.

Details

When the integration output is saved as a data.frame, then the required attributes are lost and method diagnostics will not work anymore.

Value

none

Examples

## 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)

[Package deSolve version 1.3 Index]