clear all

% LB plot of basic Michaelis-Menten speed vs substrate behavior.  
%
% Executes properly on MATLAB version 9.4.0.813654 (R2018a).
% Executes mostly properly on OCTAVE version 4.2.1.
% Joshua C. Vaughan 2/27/2023

Km = 5; % mM
Vmax = 4; % mM/sec
S = linspace(0,20,2^6); % mM

% Michaelis-Menten basic kinetics
v = Vmax*S./(Km+S);

% create LB 1/v vs 1/[S] plot

ax_lw = 2;
pl_lw = 2;
ax_fs = 24;
la_fs = 24;
ti_fs = 24;
fontname = 'Arial';

figure(1); clf; set(gcf,'color',[1 1 1])
plot(1./S,1./v,'k-','linewidth',2);
%hold on; plot(1./[S(2) S(end)],1./(Vmax*[1 1]),'--','linewidth',2,'color',0.5*[1 1 1]); hold off

set(gca,'tickdir','out')
set(gca,'linewidth',2)
xlabel('1/[S], (1/mM)','fontsize',la_fs)
ylabel('1/{\itv}, (sec/mM)','fontsize',la_fs)
set(gca,'fontsize',ax_fs)
ylim([0 3]); xlim([0 2])
title('Lineweaver-Burk Plot, {\itK}_M=5 mM','fontsize',ti_fs)