svyCprod function

Computations for survey variances