A Multivariate Fay-Herriot model is used to aid the prediction of small area parameters of dependent variables with sample data aggregated to area level. The empirical best linear unbiased predictor of the parameter vector is used, and an approximation of the elements of the mean cross product error matrix is obtained by an extension of the results of Prasad and Rao (1990) to the multiparameter case. Three different bootstrap approximations of those elements are introduced, and a simulation study is developed in order to compare the efficiency of all presented approximations, including a comparison under lack of normality. Further, the number of replications needed for the bootstrap procedures to get stabilized are studied
