function enzy_subst(k1,k2,k3,s0,e0,tf)
%
% Utilizamos la funcion ODE45 para resolver las ecuaciones
% del modelo Michaelis-Menten de enzima-substrato.
%
[t,y]=ode45('en_sust',[0,tf],[s0,e0,0,0]',[],k1,k2,k3);
%
% s representa la concentracion del substrato.
%
s=y(:,1);
%
% e representa la concetracion de enzima.
%
e=y(:,2);
%
% c representa la concetracion del compuesto substrato-enzima.
%
c=y(:,3);
%
% p representa la concetracion del producto.
%
p=y(:,4);
figure(1);clf
plot(t,s,t,e,t,c,t,p)
legend('Substrato','Enzima','SE','Producto')
