We prove that the homogeneously polyanalytic functions of total order m, defined by system equations $$\overline{D}^{(k_1,\ldots ,k_n)} f=0$$ with $$k_1+\cdots +k_n=m$$ , can be written as polynomials degree $$<m$$ in variables $$\overline{z_1},\ldots ,\overline{z_n}$$ some analytic coefficients. establish a weighted mean value property for such functions, using reproducing Jacobi polynomials. ...