function [FeOutB, FeOutV] = CellIron(EntTot)

format long

%finds concentrations of FeEnt in bacterial/viral infections,
%for a given EntTot concentration

%Prepare stock: find FeEnt from Ent, Fe and Che
FeTot = 180e-6; %intial Fe concentration units:M
CheTot = 180*3e-6; %inital Che concentration units:M
[~,~,~,initFeEnt,~] = FeEntChe(FeTot,EntTot,CheTot);

%FeEnt Binding: find [FeEnt] after Lip binding
BindTime = 0:1e0:1e3;
LipFreeB = 20e-9; % Lip concentration in blood (bacterial) units:M
LipFreeV = 4e-9; % Lip concentration in blood (viral) units:M
[~, FeEntB, ~] = LipoBinding(LipFreeB, initFeEnt, BindTime);
[~, FeEntV, ~] = LipoBinding(LipFreeV, initFeEnt, BindTime);

% equilibrium [FeEnt] after Lip binding
FeOutB = FeEntB(end);
FeOutV = FeEntV(end);

end