function [ dSdt ] = AlyssaGomesDE(t,S) % this function defines the enzyme differential equation global V K dSdt = zeros(size(S)); s=S(1); %substrate value dsdt = -V*s/(K+s); dpdt = V*s/(K+s); dSdt(1) = dsdt; dSdt(2) = dpdt; end