function [ dXdt ] = KaraDismukeDE(t,X) %this function defines the enzymatic DE global V K dXdt = zeros(size(X)); s = X(1); dsdt = -V*s/(K+s); dpdt = V*s/(K+s); dXdt(1) = dsdt; dXdt(2) = dpdt; end