The R-script grind.R is a wrapper around the deSolve, rootSolve and FME R-libraries making it very easy to perform phase plane analysis and fit ODE models to data. Scripts and tutorials are provided on https://tbb.bio.uu.nl/rdb/grind.html Enjoy Rob de Boer ----------- History: 15-03-2024 - added the shading of nullclines to cube() - added tutorials in the form of webpages using R-markdown 21-01-2024 - added option silent to newton (and removed value) - simplified definition of option eps in plane 10-09-2023 - improved automatic scaling of axes in timeplot - when using plane(...,add=TRUE) plane remembers the previous axes - added option logpar to fit the log of the parameter value - added option addone to redefine 1 as 0 on logarithmic axes (like log1p). - no longer hide the control option when calling Pseudo in fit - updated the tutorial 17-03-2023 - added option value to newton to return the largest eigenvalue 09-09-2021 - arrest now uses the event handling in deSolve and stops better 28-04-2021 - corrected a minor bug in continue 07-10-2020 - allowed more arguments to be passed on to the plot functions 31-03-2020 - removed a bug appearing when bootstrapping with a single parameter 13-08-2019 - removed bug introduced by optimizing "after" - missing data is automatically removed before fitting and plotting 17-01-2019 - optimization of option "after" in run - optimization of nullcline algorithm (also in cube) 01-12-2017 - minor bug fixes 06-11-2017 - minor bug fixes 23-10-2017 - bug in the call of trajectories corrected - added options col and pch to plane() and run() for the color and symbol used for trajectories 10-10-2017 - added automatic scaling of ymin in fit() - add-on cube.R for 3D projection was added to directory 31-7-2017 - added the option delay to run() for solving delay differential equations 21-7-2017 - renamed who into free - continue() now turns at saddle-node bifurcations - improved the passing of ellipsis (...) to other functions - grind.R now sets: options(stringsAsFactors=FALSE) - changed all T and F in TRUE and FALSE 31-5-2017 - Added numeric events to arrest - Added elipses (...) to modCost - Added timePlot=, cost=, add= to options of fit - Corrected bug in call of timePlot in fit