global k1 km1 k2 e0
format long e
x0 = [1;0.2;0;0];
% x0 is representative of the variables s0, e0, es0/c0, and p0

t = 0:0.1:100;

k1 = 2;
km1 = 0;
k2 = 10;

K = (km1+k2)/k1;

e0 = x0(2);

[tt,xx] = ode45('catalysis',t,x0);


y0 = x0(1);

[ts,ys] = ode45('michaelis_menten',t,y0);
cs      = e0*ys./(ys+K);
figure
plot(tt,xx),legend('susbtrate','enzyme','complex','product'),title('Concentration per time');
saveas (gcf, 'Concentration', 'jpeg')

figure
plot(tt,[xx(:,1),xx(:,3),ys,cs]),legend('susbtrate','complex','MM substrate','MM complex'),title('MM vs Regular');
saveas (gcf, 'MM comparison', 'jpeg')

close all