Here we shall introduce the concept of harmonic balls in subdomains of Rn, through a mean value property for a sub-class of harmonic functions on such domains. In the complex plane, and for analytic functions, a similar concept fails to exist due to the fact that analytic functions can not have prescribed data on the boundary. Nevertheless, a two-phase version of the problem does exists, and gi...