diff --git a/main_shared_calibration_2022_08_18.m b/main_shared_calibration_2022_08_18.m
index 1641e3d..ca3409e 100644
--- a/main_shared_calibration_2022_08_18.m
+++ b/main_shared_calibration_2022_08_18.m
@@ -37,8 +37,8 @@
 % GuessBOUNDS = array for boundaries of initial parameter guesses -----
 GuessBOUNDS(:,01) = [1 1]*mean(cat(1,IVs{1}{:})) ;% V iv [log(#/mL)] (co 1)
 GuessBOUNDS(:,02) = [1 1]*mean(cat(1,IVs{3}{:})) ;% V iv [log(#/mL)] (co 2)
-GuessBOUNDS(:,03) = [1 1]*mean(cat(1,IVs{2}{:})) ;% S+N iv [#/µL] (co 1)
-GuessBOUNDS(:,04) = [1 1]*mean(cat(1,IVs{4}{:})) ;% S+N iv [#/µL] (co 2)
+GuessBOUNDS(:,03) = [1 1]*mean(cat(1,IVs{2}{:})) ;% S+N iv [#/µL] (co 1)
+GuessBOUNDS(:,04) = [1 1]*mean(cat(1,IVs{4}{:})) ;% S+N iv [#/µL] (co 2)
 GuessBOUNDS(:,05) = [.03 .3]   ;% initial frequency: S/(S+N) (co 1)
 GuessBOUNDS(:,06) = [.01 1]    ;% normalized S/(S+N) (co 2)
 
@@ -49,7 +49,7 @@
 GuessBOUNDS(:,11) = [1 10] ;% normalized Sa reversion rate constant []
 GuessBOUNDS(:,12) = [1 10] ;% normalized Na reversion rate constant []
 
-GuessBOUNDS(:,13) = [300 3e3]  ;% 50% S+N proliferation saturation [#/µL]
+GuessBOUNDS(:,13) = [300 3e3]  ;% 50% S+N proliferation saturation [#/µL]
 GuessBOUNDS(:,14) = [0.1 1]    ;% S0/N0 normalized prolif rate constant [/d]
 GuessBOUNDS(:,15) = [1 4]      ;% Sa proliferation rate constant [/d]
 GuessBOUNDS(:,16) = [0.5 2]    ;% Na proliferation rate constant [/d]
@@ -87,7 +87,9 @@
 % costFunction = handle for calling cost function -------------------------
 costFunction = @(P) costFun_NLL(P,modelFunction,X,Y,'Xevaluate',Xeval) ;
 
+% tic
 % costFunction(GuessBOUNDS(1,:)) ;% uncomment to debug your costFunction
+% disp(toc)
 
 % filePrefix = string to prefix to file names -----------------------------
 filePrefix = 'N803_shared' ;
diff --git a/main_shared_noN_calibration_2022_11_02.m b/main_shared_noN_calibration_2022_11_02.m
index 9453dcc..293cee6 100644
--- a/main_shared_noN_calibration_2022_11_02.m
+++ b/main_shared_noN_calibration_2022_11_02.m
@@ -84,7 +84,9 @@
 % costFunction = handle for calling cost function -------------------------
 costFunction = @(P) costFun_NLL(P,modelFunction,X,Y,'Xevaluate',Xeval) ;
 
+% tic
 % costFunction(GuessBOUNDS(1,:)) ;% uncomment to debug your costFunction
+% disp(toc)
 
 % filePrefix = string to prefix to file names -----------------------------
 filePrefix = 'N803_shared_noN' ;
diff --git a/main_singles_calibration_2022_09_25.m b/main_singles_calibration_2022_09_25.m
index c4eab89..2fdc22a 100644
--- a/main_singles_calibration_2022_09_25.m
+++ b/main_singles_calibration_2022_09_25.m
@@ -95,7 +95,9 @@
 % costFunction = handle for calling cost function -------------------------
 costFunction = @(P) costFun_NLL(P,modelFunction,X,Y,'Xevaluate',Xeval) ;
 
+% tic
 % costFunction(GuessBOUNDS(1,:)) ;% uncomment to debug your costFunction
+% disp(toc)
 
 % filePrefix = string to prefix to file names -----------------------------
 filePrefix = ['N803_cohort_' num2str(c)];