%constants and equations

%these are mostly from L&L chapter 3
L = 1e-9; %M , constant
rho = 10^10; %cells/L, constant
Nav = 6.022e23; %#/mol, constant
%all of these rates are also per cell
kf = 7.2e7; %M^-1 min^-1
kr = 3.4e-1; %min^-1
VR = 1.3e2; %#/min
keR = 3.0e-2; %min^-1
keC = 1e-1; %min^-1
kxR = 5.8e-2; %min^-1
kxC = 5.3e-2; %min^-1
kh = 1e-1; %min^-1
kfprime = 2e6; %M^-1 min^-1
krprime = 1.2e-2; %min^-1
kxL = 1e-2; %min^-1
kfp = 4e-14; %min^-1

dRsdt = VR - (kf*L*Rs - kr*Cs) - keR*Rs + kxR*Ri;
dCsdt = kf*L*Rs - kr*Cs - keC*Cs + keC*Ci;
dRidt = keR*Rs - kxR*Ri - kh*Ri - (kfprime*Li*Ri - krprime*Ci);
dCidt = keC*Cs - kxC*Ci - kh*Ci + kfprime*Li*Ri - krprime*Ci;
dLidt = krprime*Ci - kfprime*Li*Ri - kh*Li - kxL*Li + kfp*L;