# Michaelis–Mentenin model

K: 8$
R: 5$
M: 0.5$
F: 1$
v(e)=K*e0*e/(M+e)$
rj: rk([-F*e*s+(R+K)*c,-F*e*s+R*c,F*e*s-(R+K)*c,K*c ],[e,s,c,p ],[40,100,0.0,0.0 ],[ t,0,1,0.001 ])$
rje: makelist([ f[1], f[2 ]],f,rj)$
rjs: makelist([ f[1], f[3 ]],f,rj)$
rjc: makelist([ f[1], f[4 ]],f,rj)$
rjp: makelist([ f[1], f[5 ]],f,rj)$
wxplot2d([[discrete,rje ],[discrete,rjs ],[discrete,rjc ],[discrete,rjp ]],[legend, "e","s","c","p"]  );

# Kimotripsin

kill(K,R,M,F);

K: 0.14;
M: 1.5/100;

solve([M=(R+K)/F],[F,R]);

K: 0.14;
R: 0.001;
F: (200*R+28)/3;
M: (R+K)/F;

load(rkf45);

rj: rkf45([-F*e*s+(R+K)*c,-F*e*s+R*c,F*e*s-(R+K)*c,K*c ],[e,s,c,p ],[40,100,0.0,0.0 ],[t,0,0.01])$
rje: makelist([ f[1], f[2 ]],f,rj)$
rjs: makelist([ f[1], f[3 ]],f,rj)$
rjc: makelist([ f[1], f[4 ]],f,rj)$
rjp: makelist([ f[1], f[5 ]],f,rj)$
wxplot2d([[discrete,rje ],[discrete,rjs ],[discrete,rjc ],[discrete,rjp ]],[legend, "e","s","c","p"]);

rj: rkf45([-F*e*s+(R+K)*c,-F*e*s+R*c,F*e*s-(R+K)*c,K*c ],[e,s,c,p ],[40,100,0.0,0.0 ],[t,0,25])$
rje: makelist([ f[1], f[2 ]],f,rj)$
rjs: makelist([ f[1], f[3 ]],f,rj)$
rjc: makelist([ f[1], f[4 ]],f,rj)$
rjp: makelist([ f[1], f[5 ]],f,rj)$
wxplot2d([[discrete,rje ],[discrete,rjs ],[discrete,rjc ],[discrete,rjp ]],[legend, "e","s","c","p"]);