Helper function for symdynN, used to simulate dynamics between disturbance events. Note that K is set to 1 for all species.
df_col(time, state, pars)rate of change for each species
time step - ignored, but required for consistency with ode function
vector of current states
parameter list, including matrix A with full interaction matrix (and growth rates along the diagonal), and a value Ifrac, which is the dispersal rate (D in the main text), and Ksim, which is the carrying capacity