Helper function for computing the posterior mean, posterior variance
Basic.integrand(u, y, k, tau, data.var)