% grossberg2Reqns.m, m-function by ralph, 15 oct 2003
function dydt = grossberg(t,y,W,P);
% grossberg shunt model, 2nd layer 
% args = (time, state scalar, row vector weight, col vector pattern)
dydt = -y + (1-y)*(10*y^2/(1+y^2) + W*P) - y*10*y^2/(1+y^2);
% end;