Functions for drawing samples according to Sunter's sampling design, and for computing first and second order inclusion probabilities