Conformal prediction for multi-dimensional time series by ellipsoidal sets