Skip to content
Prev 282278 / 398498 Next

Brillouin index

It's a pretty simple formula, according to the sources I found.

Here's a function that looks right to me, but I have no independent
calculation with which to check it.
(no guarantees!)

Hb <- function(ns) {
  N <- sum(ns)
  (lfactorial(N) - sum(lfactorial(ns)))/N
}
[1] 1.028802