# GRIND for MATLAB

Version 2.00

January 2017

MATLAB versions: R2008a to R2017a

GRIND for MATLAB is a modelling platform that supports ordinary differential equations,

difference equations, delay differential equations (DDE), vector models and

matrices (cellular automata). It has many tools for the analysis of models, including bifurcation analysis.

**New feature! ** GRIND can now do 1D and 2D bifurcation analysis, using specialized software packages. See the command "conteq". Two continuation tools are supported:

- MatCont (Dhooge, et al. 2003 ACM Transactions on Mathematical Software 29:141-164)
- COCO (Dankowicz, H. and Schilder, F. 2013 Recipes for Continuation. SIAM, Philadelphia. 584 pp.) .

- Download grind.zip (2.0 Mb) tested on newer MATLAB versions (>=R2008a)
- Contact
- Download grindr11.zip (1.02 Mb) GRIND 1.0, for MATLAB versions 5.3 - R2007b

**Description**

GRIND is a convenient DOS program for analyzing sets of differential equations made by Rob de Boer.

GRIND for MATLAB is a MATLAB version of GRIND, which also supports difference equations, cellular automata and time lags.

Use the command MODEL to create a model (which is saved to an ini file). Use the upper panel in the window to enter differential equations or difference equations. New: you can use the command VISMOD (visual model) to enter a simple model as a Forester diagram.

None of the MATLAB toolboxes are required, but for a few commands the Statistics and (especially) Symbolic Math toolboxes are recommended.

**Example:**

Logistic differential equation:

N'=N*r*(1-N/K)

Logistic difference equation:

N(t+1)=N(t)*r*(1-N(t)/K)

Note that parameters and state variables are case-sensitive. Parameter names can be all alphanumeric names except the following reserved words:

t time

pi pi=3.1416

inf Infinity

Inf Infinity

nan NaN (not-a-number)

NaN NaN (not-a-number)

eps Floating point relative accuracy eps=2.2204e-016

In the lower panel default values of the parameters and initial values of the state variables are entered.

Furthermore default commands can be entered here.**Example:**

N=0.01;

r=0.5;

ax x N [0 10];

(the semicolon is not required, but suppresses unnecessary output)**Author:**

Egbert van Nes

Wageningen University

Aquatic Ecology and Water Quality Management group

P.O. Box 47

6700 AA Wageningen

The Netherlands