On the error in multivariate polynomial interpolation