Matlab function of design equation for reactor 1:

function V = steadystate_cstr_volume(f,S,P) 
%calculates volume of CSTR in steady state -> design equation
% f = flow
% S = substrate concentration
% P = product concentration that will leave the reactor
% we assume that the inlet water does not contain any product

kcat = 18765; %[1/min]  AMO+HAO
kM = 1.34; % [mM] AMO+HAO
E = 8.457e-4; % estimation that cells take up 40% of the volume 1000 enz/cell

V = f*P/(((kcat*E*S)/(kM+S))*0.4);