gradient function for ODE runs

gradfun(t, y, parms, M)

Arguments

t: time vector
y: state vector
parms: parameter vector
M: rate matrix