given the distance matrix the unbiased estimate for mu3 is computed
mu3.unbiased(B, b2ob = sum(tcrossprod(B) * B))