Changes version 1.10.6 (Thomas) 
================================
 o change declaration of variable dimensions from (1) to (*) in legacy
   Fortran code to pass automatic bounds check
 o remove the Jacobian examples from ?ode because banddown=0 can
   lead to problems on some systems; examples will come back in the
   next release
 o fixed bug in the "iteration" solver
 o small documentation updates

Changes version 1.10.5 (Karline, Thomas)
================================
 o extended subset.deSolve with argument arr, when TRUE returns an
   array for >2-D output
 o fixed the R compiler notes
 o plot.ode.2D now has an mtext argument, via the ..., to label multiple 
    figures in margin... CHECK - see ode.2D
 o subset can also be a vector with indices in addition to logical
 o image with legend = TRUE changed size of plot in different layouts - now solved
   (by adding par(mar = par("mar")) )
 o new method to output warnings and error messages 
 o add data type check for external outputs in rk_util.c
 o add interface for compiled dede models
 o emphasize consistent order of states in y and return value of func
 o changes of Fortran error messages (to be continued)

Changes version 1.10-4 (Thomas, Karline)
================================
 o allow reverted time vector for fixed step solvers
   - todo: find solution for dense output methods, and Livermore solvers
 o all solvers now have default atol = 1e-6; before this daspk and vode had 1e-8.
 o multiple warnings from daspk if num steps = 500 toggled off.
 o added input argument "nind" to daspk, to make it compatible with radau.
   this also changes the way the variables are weighed, 
   hence this differs from the original daspk 2.0 code.
 o improved warning printing in daspk and vode 
 o extended sparse Jacobian input in lsodes. (2-D and 3-D sparsity
   with mapping var and arbitrary sparsity in ian/jan format).

Changes version 1.10-3 (Karline)
================================
 o rwork and iwork in lsodes from Fortran -> C (to remove compiler warnings)
 o roots + events:  now certain roots can stop simulation + fixed bug in radau root
 o improved events\roots help file
 o diagnostics(out) gave error in case method=iteration (no rstate) now fixed
 o the package authors agreed to assign the maintainer role to T.P.,
   but the order of authorship and credits remain unchanged.

Changes version 1.10-2 (Karline)
================================
 o remove NAs from forcing functions - when used in DLL (file forcings.R)
 o new argument "restructure" in ode.1D, for use with implicit solvers not in deSolve
 o removed requirement to have eventfunc in compiled code when func is in compiled code
 o subsetting on summary.deSolve
 
Changes version 1.10-1 (Thomas)
===============================
 o remove several redundant variables from C code
 o add NEWS file

Changes version 1.10 (Karline, Thomas)
======================================
 o compiled code using mass in daspk
 o cleanEventTimes

Changes version 1.9+ (Karline)
==============================
 o roots, events, lags in radau
 o roots in lsodes
 o lags in daspk
 o ode (method = "iteration")

Changes version 1.9 (Karline, Thomas)
=====================================
 o summary.deSolve
 o subset.deSolve
 o plotting deSolve objects improved:
   - plot more than one output in same figures (scenarios), 
   - add observations
 o vignette improved
 o fixed bug in 'timesteps'

Changes version 1.8.1 (Thomas, Woody, Karline)
==============================================
 o fixed compiler warnings using valgrind
 o fixed compiler warning C-code

Changes version 1.8 (Thomas)
============================
 o Dormand-Prince 8(7) coefficients use now common 
   instead of decimal fractions


Changes version 1.8 (Karline)
=============================
 o Runge-Kuttas:
   - extra output: number of failed steps (see also 2)
   - number of function evaluations + 1 for initial condition
   - dense output for cash-karp
   - dopri8(7) added
   - radau added!! implicit runge kutta, solves also DAE up to index 3!
 o other:
   - image function for ode.2-D added.
   - changed warning printing in FORTRAN code
   - common interface for radau and daspk:
     both can solve systems written as M*dy = f(x,y).
     daspk can also solve systems written as 0 = g(x,y,dy) (=default for daspk)