N(2,kept,const): 18 states, reaction groups {'c1': 6, 'c2': 6, 'a': 3, 'delta': 9, 'doff': 4, 'don': 2, 'c3': 6}
  random search (3000 vectors, rates in [1e-4,1e4], L in [1e-8,1e8]): min slope 5.9595e-13  (12s)
  40 Nelder-Mead runs: min slope 1.9287e-15 at log10 params [  1.96 -27.88  -1.78 -10.4   -6.61   5.66  -6.16]
N(3,kept,const): 26 states, reaction groups {'c1': 9, 'c2': 9, 'a': 4, 'delta': 13, 'doff': 6, 'don': 3, 'c3': 9}
  random search (3000 vectors, rates in [1e-4,1e4], L in [1e-8,1e8]): min slope 3.9779e-13  (16s)
  40 Nelder-Mead runs: min slope 9.6433e-16 at log10 params [ 1.1  -8.18  3.42  0.65 -7.68 -0.02 -0.26]
