maybe histogram without dice rolling for identical domains?

Further tension models:
  Exponential rising tension (ERT)
    f = E0 / alpha ( e^(alpha x) - 1 )
  Enhanced WLC (with stretching)
    f = kB T/p { 1/4 [1/(1-x/L-f/K)² - 1] + x/L - f/K }
  Enhanced FJC (with stretching)
    x/L = [coth(Lk f/ kB T)- kB T/Lk f](1+f/K)
from Leuba & Zlatanova.  Biology at the Single-Molecule Level.  p13