Linear and quadratic prediction problems in finite populations have become of great interest to many authors recently. In the present paper, we mainly aim to extend the problem of quadratic prediction from a general linear model, of form y D X C e; e N. 0; 2V/, to a multivariate linear model, denoted by Y D XB C E; Vec.E/N.0;6 V/with Y D.yij/n q D. y1, ……yq/. Firstly, the optimal invariant quadratic unbiased (OIQU) predictor and the optimal invariant quadratic (potentially) biased (OIQB) predictor of Y0HY for any particular symmetric nonnegative definite matrix H satisfying HX D 0 are derived. Secondly, we consider predicting a0Y0HYb and tr.Y0HY/. The corresponding restricted OIQU predictor and restricted OIQB predictor for them are given. In addition, we also offer four concluding remarks. One concerns the generalization of predicting a0Y0HYb and tr.Y0HY/, and the others are concerned with three possible extensions from multivariate linear models to growth curve models, to restricted multivariate linear models, and to matrix elliptical linear models.