function [ dXdt ] = KarinaAlvarezDE( t,X )
global K V

dXdt=zeros(size(X));

s=X(1); %substrate

dsdt=-V*s/(K+s);
dpdt=V*s/(K+s);

dXdt(1)=dsdt;
dXdt(2)=dpdt;


end