function dx = glureg(t,x,par)
% Glucolysis Regulation
% t = time
% x(1) = normalized ATP
% x(2) = normalized ADP
% par = Parameter set vector [4x1]
my = par(1);
v = par(2);
alfa = par(3);
eta = par(4);
f = x(1)*x(2)^my/(x(2)^my*x(1) + x(2)^my + 1);
dx(1,:) = v - f;
dx(2,:) = alfa*f - eta*x(2);