function translateModel(kc, Km, ics)
% function translateModel(kc, Km, ics)

% This work is licensed under the Creative Commons Attribution-Noncommercial-No Derivative 
% Works 3.0 United States License. To view a copy of this license, visit 
% http://creativecommons.org/licenses/by-nc-nd/3.0/us/ or send a letter to Creative Commons,
% 171 Second Street, Suite 300, San Francisco, California, 94105, USA.

nX = length(ics);

[e1 e2 s p] = jdeal(1:nX);

%% build model
[w{1}, dwdx{1}] = add_MM(s, e1, p, kc, Km, nX);
[w{2}, dwdx{2}] = add_MM(p, e2, s, kc, Km, nX);
[f dfdx] = flux2rhs(w, dwdx, nX);

%% simulate model
xSol = ode15s(f, [0 50], ics, odeset('Jacobian', dfdx));
plot(xSol.x, xSol.y)
end
