function UP_celldivision2
% Plott cell division rate depending on the ampicillin concentration
global KPLUSEINS;
global KMINUSEINS;
global KEXPRPROT;
global KPLUSZWEI;
global KMINUSZWEI;
global KATALYSE;
global KATALYSEZWEI;
global KEXPRTOR;
global KTRANSTOR;
global KTRANSAMP;
global KDEGREINS;
global KDEGRZWEI;
global AMPNULL;
KPLUSEINS =  8.e-6;
Kplus1 = KPLUSEINS*(6.022*1.e+23)*(9.e-16);
%KMINUSEINS = 4.e-4;
KD = KMINUSEINS/Kplus1;
KEXPRPROT =  0.6;
KDEGREINS =  17.e-4;
%
KPLUSZWEI =  4.e-5;
KMINUSZWEI = 2.e-4;
KATALYSE =   8;
KEXPRTOR =   1;
KDEGRZWEI =  2.e-3;
KTRANSTOR =  1.e-3;
%
KTRANSAMP = 5.e-3;
KATALYSEZWEI = 0.001;
AMPNULL = 8.e+2;
%
m = 6;
for n=0:0.2:m-1                  %parameter range
    %using parameter: increasing and solving ODEs in each loop cycle
    N = 1.e-8*((10^(n)));
    KMINUSEINS = N;
    KD=KMINUSEINS/Kplus1;
    q = fix((5*n)+1);
    %--- solving odes
    conzini1 = [3000 , 0 , 0]';
    [t1, csol1] = ode23s('conzt1', 4000, conzini1);
    l1=length(t1);
    %--- solving odes
    conzini2 = [csol1(l1,1) , csol1(l1,2) , csol1(l1,3) , 0 , 0 , 0]';
    [t2, csol2] = ode23s('conzt2', 4000, conzini2);
    l2=length(t2);
    %--- solving odes
    conzini3 = [csol2(l2,1) , csol2(l2,2) , csol2(l2,3) , csol2(l2,4) , csol2(l2,5) , csol2(l2,6), 0]';
    [tNR, csolNR] = ode23s('conzt3', (0:10:3000), conzini3);
    %---
    %calculate Ampicillin concentration in µg/mL out of N Molecules
    %using specific weight, perplam volume and Avogrado const.
    vek11 = ((csolNR(:,7).*((3.e-17)^(-1))).*((6.022*1.e23)^(-1))).*(349.41).*(1.e+3);
    %calculate cell doubling time rate out of Ampicillin concentration in
    %µg/mL by our self determined function
    vek12 = 20+(10.*2.^((2.*vek11)-1));
    %
    M(:,q)=vek12;
    T(:,q)=tNR;
    y(q)=log10(N);
    KDvek(q)=log10(KD);
end
% make a nice surface plot
figure
l=length(y);
t=T(:,l);
surf(KDvek,((t+8000)/60)',M,'EdgeColor','none')
xlabel('log10(k_D) in mol*L^-^1');
ylabel('Time in min');
zlabel('Cell doubling time in min');
title(['Cell doubling time']);
colormap jet;
colorbar;
camlight left;
lighting phong;
