Compute the matrix of joint inclusion probabilities from the quadratic form of a Horvitz-Thompson variance estimator.
Useful links