function dydt = selkov(t,y,p)

su=y(1);
sd=y(2);

global g v n a
   
f=su*sd^g/(su*sd^g+sd^g+1);

dydt = [v-f; a*f-n*sd];



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%