samSige function

Internal function to sample variances for eta

samSige(v1, v2, eta, mu, e)