function dy = MM_3CM_allp(p, t, y)

% exponential of parameter constrains optimization to
% nonnegative values

%UNITS IN ml!

%lamb = 1e1;
%dT = 0.01;
%k = 4.57e-5;
%dI = 0.4;
%NdI = 38;
%dV = 2.4;


dy = [p(1) - p(2)*y(1)*y(3) - p(3)*y(1) + (p(4)/(p(5) + y(3)))*y(1)*y(3); 
      p(2)*y(3)*y(1) - p(6)*y(2); 
      p(7)*y(2) - p(8)*y(3)];

% 3CM_MM - FULL parameter space
% dy = [p(1) - p(2)*y(1)*y(3) - p(3)*y(1) + (p(4)/(p(5) + y(3)))*y(1)*y(3); 
%     p(2)*y(3)*y(1) - p(6)*y(2); 
%     p(7)*p(6)*y(2) - p(8)*y(3)];


end
            
            
            