global k1, k2,k3;
k1=100;k2=0.1;k3=100;

%initial conditions
S0=1; E0=1; ES0=0; P0=0; tend=100;
y0=[S0 E0 ES0 P0];

[t,y]=ode45(@MMmodel,[0 tend],y0);
%result

S=y(:,1);E=y(:,2);ES=y(:,3);P=y(:,4);

subplot(211)
plot(t,S,'k',t,P,'y');

subplot(212)
plot(t,ES,'k',t,E,'y');
