TODO

   1 maybe histogram without dice rolling for identical domains?
   2
   3 Further tension models:
   4   Exponential rising tension (ERT)
   5     f = E0 / alpha ( e^(alpha x) - 1 )
   6   Enhanced WLC (with stretching)
   7     f = kB T/p { 1/4 [1/(1-x/L-f/K)² - 1] + x/L - f/K }
   8   Enhanced FJC (with stretching)
   9     x/L = [coth(Lk f/ kB T)- kB T/Lk f](1+f/K)
  10 from Leuba & Zlatanova.  Biology at the Single-Molecule Level.  p13
  11