Walsh functions, scrambled $(0, m, s)$-nets, and negative covariance: applying symbolic computation to quasi-Monte Carlo integration.