Harmonicity of a function via harmonicity of its spherical means

It is proved that harmonic functions are characterized by harmonicity of their spherical means, for which purpose the iterated spherical means are used. The similar characterization of solutions to the modified Helmholtz equation (panharmonic functions) is given. Another description of harmonic functions is the pointwise equality of a function and its iterated mean over an admissible pair of spheres.