% LEP 9-3 : Batch Enzymatic Reactors
%

Vmax1=1.2;
KM=0.0233;
Co=0.1;
Et2=0.001;
Et1=5;
Vmax2=Vmax1*Et2/Et1;
X=linspace(0,1,200);
for i=1:200
t(i)=(KM/Vmax2)*log(1/(1-X(i)))+Co*X(i)/Vmax2;
end
plot(t(1:200),X(1:200),'LineWidth',2)
xlabel('time (s)')
ylabel('Conversion')
axis([0 1000 0 1])
title('Batch Enzymatic Reactors')
grid on

%