Multivariate Conformal Prediction via Conformalized Gaussian Scoring