% grossbergReqns.m, m-function by ralph, 14 oct 2003
function dydt = grossberg(t,y,p,q);
% exemplary dyn scheme (2D flow with one param fixed)
b = 1.0;
dydt = [-y*(1+q) + b*p*q];
end;