function [ reproduction_rate ] = doubling_ecoli_time( molL )
%calculates the growth rate depending on the amp concentration in the
%periplasm
%   Input: amp concentration in periplasm in mol/l
%   Output: growth rate

%computes the concentration of amp in ug/ml 
%   molar mass of amp: 349.41 g/mol
ugMl = (molL/1000)*349410000;

%calculates how long it takes until one E. coli cell divides in dependence of the amp
%concentration
doubling = 20+10*2^(1*ugMl+1);
%calulates the reproduction rate per hour
reproduction_rate = 1/(doubling/60);
end

