% Created by Nick Fidalgo, Chris Gatesman, Phillip Less, Brett Peters, Samantha Sandwick, and Alex Shetler
function dydt = hill_M2(t,s,K_m,K_5,om,sig,gam,xi,phi,A_0,a)
%ODE45
B = s(1);
D = s(2);
E = s(3);
B_0 = 174;
E_0 = 45;
A = A_0*exp(-a*t);
dydt = zeros(3,1);
dydt(1) = -K_m*(A*A_0).^(om)/(sig*(K_5^(om) + (A*A_0).^om))*B; 
dydt(2) = K_m*(A*A_0).^(om)/(sig*(K_5^(om) + (A*A_0).^om))*B - D;
dydt(3) = (gam*D*B_0)/E_0 + (xi*gam*B*B_0)/(sig*E_0) - (phi*E)/(sig);
end