Fitting Multidimensional Data with Multivariate Gaussian Uncertainties