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