% LEP 9-2 : Evaluation of Michaelis-Menten Parameters Vmax and KM
%

Vmax=1.33;
KM=0.0266;
S=linspace(0,0.25,200);
Sinv=linspace(0,500,200);
p=linspace(0,50,200);
for i=1:200
rurea(i)=Vmax*S(i)/(KM+S(i));
LB(i)=(1/Vmax)+(KM/Vmax)*Sinv(i);
EF(i)=Vmax-KM*p(i);
HW(i)=(KM/Vmax)+S(i)/Vmax;
end
tiledlayout(2,2)
set (gcf,'Position',[0,0.1,800,600])
nexttile
plot(S(1:200),rurea(1:200),'LineWidth',2)
xlabel('C_{urea} (kmol/m^3)')
ylabel('-r_{urea} (kmol/m^3.s)')
axis([0 0.25 0 5])
title('Michaelis-Menten Plot')
grid on
nexttile
plot(Sinv(1:200),LB(1:200),'LineWidth',2)
xlabel('1/C_{urea} (m^3/kmol)')
ylabel('1/-r_{urea} (m^3.s/kmol)')
axis([0 500 0 100])
title('Lineweaver-Burk Plot')
grid on
nexttile
plot(p(1:200),EF(1:200),'LineWidth',2)
xlabel('-r_{urea}/C_{urea} (s^{-1})')
ylabel('-r_{urea}(kmol/m^3.s)')
axis([0 50 0 5])
title('Eadie-Hofstee Plot')
grid on
nexttile
plot(S(1:200),HW(1:200),'LineWidth',2)
xlabel('C_{urea} (kmol/m^3)')
ylabel('C_{urea}/-r_{urea} (s)')
axis([0 0.25 0 2])
title('Hanes-Woolf Plot')
grid on

%