PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 40, Number 7, July 0, Pages 453 463 S 000-9939(0)8-X Article electronically published on November, 0 MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS TIZIANA GIORGI AND ROERT SMITS (Communicated by Matthew J. Gursky) Abstract. We derive a mean value property for p-harmonic functions in two dimensions, <p<, which holds asymptotically in the viscosity sense. The formula coincides with the classical mean value property for harmonic functions, when p =, and is a consequence of a representation for the Game p-laplacian obtained via p-averaging.. Introduction A recent article by Manfredi et al. [6] (see also [9]) characterizes p-harmonic functions via a weak asymptotic formula which holds in a suitably defined viscosity sense. Inspired by their results and by our recent work [4], where we present a numerical algorithm for the Game p-laplace operator based on the idea of p-average, we derive a generalization in a viscosity sense to two-dimensional p-harmonic functions, <p<, of the classical mean value property for harmonic functions. The variational p-laplace operator is defined, for <p<, as () Δ p u div ( u p u ), while the Game p-laplacian, recently introduced in [7] to model a stochastic game called Tug of war with noise,readsas () Δ G p u p u p div ( u p u ). A function u C 0 (Ω), with Ω R a smooth domain, is called p-harmonic in ΩifitisaviscositysolutionofΔ p u = 0 (see Definition.). The focus of this paper is in providing a representation of p-harmonic functions that for the case p = reproduces the mean value property. Nevertheless, it will be clear that our main interest is the Game p-laplacian and that our approach sheds light on the local properties of the solution of the Game p-laplace operator. The representation formula derived was suggested to us by the numerical approximation we propose in an upcoming paper [4]. An insight on the local properties of the Game p-laplacian suggests that the value of a solution at a given point is related to the p-average on small balls centered at that point. The numerical solution that we construct, in the case of dimension n =, using this idea satisfies a discrete analogue of our proposed generalized mean value formula. We derive the following main results. Our first theorem finds an expansion for C functions in terms of the Game p-laplacian: Received by the editors November, 00 and, in revised form, February 6, 0. 00 Mathematics Subject Classification. Primary 35J9, 35D40, 35J60, 35J70. Funding for the first author was provided by National Science Foundation Grant #DMS- 0604843. c 0 American Mathematical Society Reverts to public domain 8 years from publication 453
454 TIZIANA GIORGI AND ROERT SMITS Theorem.. Let φ C (Ω), whereω R is a smooth domain, and let x 0 Ω. If φ(x 0 ) 0, then for any ɛ>0 such that ɛ (x 0 ) Ω we have φ(x 0 )= ɛ (x 0 ) φ(x p ɛ 0) (x x 0 ) p dx p + ΔG p φ(x 0 )+o(ɛ ). Here ɛ (x 0 ) denotes the ball of radius ɛ and center x 0. We then use this representation to derive our weak mean value formula for p-harmonic functions. Theorem.. Let u be a continuous function in Ω R,andletx 0 Ω. For any ɛ>0 such that ɛ (x 0 ) Ω we have that ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + o(ɛ ) 0) (x x 0 ) p dx holds in the viscosity sense if and only if u is p-harmonic; that is, u is a viscosity solution of Δ p 0. We present detailed proofs for the case of smooth domains Ω R, but from our treatment it will be clear how to obtain generalizations to dimensions n>. The paper is organized as follows. In Section we recall some definitions and background results. In Section 3 we derive the representation formula for C functions and discuss why this is the correct local way of describing p-harmonic functions. In Section 4 we prove Theorem.. To conclude, in Section 5 we derive a similar result for the non-homogeneous Game p-laplacian.. p-laplacian and Game p-laplacian Our representation formula for smooth functions is based on the so-called Game p-laplacian introduced by Peres and Sheffield [7], and its proof is based on the characterization of Δ G p as a convex combination of two limiting operators. When p =, traditionally the -Laplacian is given by Δ u u u u, x i,j i x j x i x j while the Game -Laplacian is its -homogeneous renormalized version: Δ G u u u u u. x i,j i x j x i x j For p =,wecansetp = in () and obtain Δ u div ( u u ), while for the Game -Laplacian we follow [7] and define it in terms of the Laplace operator and the Game -Laplacian: (3) Δ G u Δ u Δ G u. If u is a smooth function, by expanding the derivatives, one obtains Δ G p u = p Δ u + p p u i,j u x i u x j u x i x j,
MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS 455 which allows us to think of Δ G p that is, as the convex combination of the two limiting cases, (4) Δ G p = p ΔG + q ΔG, with q the conjugate exponent of p. Furthermore, the Game -Laplacian and the Game p-laplacian for u 0 can then be rewritten as the second derivative in the orthogonal direction of u andinthedirectionof u, respectively. That is, (5) Δ G u = u D u u, u, and (6) Δ G u = u D u u, u, where D u denotes the Hessian matrix. In the homogeneous case, solutions to the Game p-laplacian agree with the ones of the p-laplacian. Also note that the Game -Laplacian is the limit as p of the Game p-laplacian, a fact which is not true for the p-laplacian. The fundamental difference between the classical p-laplacian and the Game p-laplacian is that the former can be obtained as the Euler-Lagrange equation of an energy functional. Additionally, while for <p< both operators are degenerate, for u = 0theGamep-Laplacian and the -Laplacian are singular, so suitable definitions of viscosity solutions need to be given. We recall that Juutinen and the coauthors in [5] show that for the p-laplacian when <p< the notions of viscosity solution and weak solution are equivalent. Therefore, we will work with viscosity solutions for both operators. We consider the definition of viscosity solution for the p-laplacian provided in [6]: Definition.. Let Ω R be a smooth domain, and <p<. Then: (i) We say that an upper semi-continuous function u is a viscosity subsolution of Δ p u = 0 in Ω if for any φ C such that u φ has a strict local maximum at x Ω, we have (7) (p )Δ G φ(x) Δ φ(x) 0 whenever φ(x) 0. (ii) We say that a lower semi-continuous function u is a viscosity supersolution of Δ p u = 0 in Ω if for any φ C such that u φ has a strict local minimum at x Ω, we have (8) (p )Δ G φ(x) Δ φ(x) 0 whenever φ(x) 0. (iii) If u is both a subsolution and a supersolution everywhere in Ω, we say that u is a viscosity solution of Δ p u =0inΩ. Various equivalent definitions of viscosity solutions for the Game p-laplacian operator can be given and are found in the literature. The most suitable for our treatment is the one obtained by following the definition in the classical paper of arles and Souganidis []. Definition.. Consider a smooth domain Ω R,andlet<p<. Iff is a continuous function, we say that an upper semi-continuous function u (respectively, lower semi-continuous) is a viscosity subsolution (respectively, supersolution) of Δ G p u = f(x) inω
456 TIZIANA GIORGI AND ROERT SMITS if for any φ C such that u φ has a local maximum (respectively, local minimum) at x Ω, we have: (i) Δ G p φ(x) f(x) if φ(x) 0 (respectively, Δ G p φ(x) f(x)); (ii) λ p λ f(x) if φ(x) =0andp q (respectively, λ q λ p f(x)); λ q λ p f(x) (respectively, λ if φ(x) =0and<p< p λ q f(x)); here p + q =,andλ λ are the eigenvalues of D φ(x). Remark.. Part (ii) of the definition of viscosity subsolution (supersolution) is implied by the condition (ii) Δ G φ(x) f(x) whenever φ(x) = 0 (respectively, Δ G φ(x) f(x)). This is a consequence of the fact that p λ q λ Δ G φ(x) q λ p λ, if p, and q λ p λ Δ G φ(x) p λ q λ, if <p<. Recall that Δ G φ = Δ φ. Uniqueness for viscosity solutions of non-linear operators that are singular at isolated points typically does not depend on the particular value one assigns to these points as long as this is chosen in a consistent manner (see for example Section 9 in [3]). Additionally, our numerical results in [4] show that the numerical approximation converges to solutions that verify (ii). Therefore, we will use the following definition for a viscosity solution of the Game p-laplacian: Definition.3. We say that a function u is a viscosity solution of Δ G p u = f for <p< if u is a subsolution and supersolution according to (ii) of Definition. and (ii) in Remark.. 3. Representation formula Theorem 3.. Let Ω R be a smooth domain. Given φ C (Ω) and x 0 Ω for which φ(x 0 ) 0, we have that for any ɛ>0 such that ɛ (x 0 ) Ω it holds that φ(x 0 )= ɛ (x 0 ) φ(x p ɛ 0) (x x 0 ) p dx p + ΔG p φ(x 0 )+o(ɛ ). Here ɛ (x 0 ) denotes the ball of radius ɛ and center x 0. Proof. Take x =(x,x ) R and denote by e =(, 0) the unit director of the x -axis. Assume φ C (Ω), x 0 Ωand φ(x 0 ) 0. Without loss of generality, we can assume x 0 =0and φ(x 0 )= φ(x 0 ) e. Equation (6) then gives (9) Δ G φ(0) = φ(0),
MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS 457 while (5) yields (0) Δ G φ(0) = φ(0). For any ɛ>0 such that ɛ ɛ (0) Ω, if <p<, we integrate over ɛ to obtain ɛ φ(0) x p D φ(0)x, x dx which by (9) and (0) implies = φ(0) p x p D φ(0)x, x dx ɛ = φ(0) p x p D φ(0)x, x dx, ɛ {x >0} φ(0) x p D φ(0)x, x dx ɛ = φ(0) [Δ p G φ(0) x p dx ɛ {x >0} + φ(0) x p x x dx ɛ {x >0} ] () +Δ G φ(0) x p x dx. ɛ {x >0} Also, note that () φ(0) x p dx = φ(0) p x p dx. ɛ ɛ {x >0} To compute the integrals in () and (), we use polar coordinates and derive ɛ π π x p dx = (r cos θ) p dθ r dr = ɛp+ (cos θ) p dθ, ɛ {x >0} 0 π p + π as well as ɛ π x p x dx = (r cos θ) p (r sin θ) dθ r dr ɛ {x >0} = ɛp+ p + = ɛp+ p + π π [ π π 0 π (cos θ) p ( cos θ) dθ while by symmetry we see that x p x x dx = ɛ {x >0} = ɛp+ p + π π π ] (cos θ) p dθ (cos θ) p dθ π, ɛ 0 π π (cos θ) p sin θdθ=0. (r cos θ) p (r sin θ) dθ r dr
458 TIZIANA GIORGI AND ROERT SMITS For the last integral that we need, we find that ɛ π π (3) x p dx = (r cos θ) p dθ r dr = ɛp ɛ {x >0} 0 π p π We next substitute the above integrals in () and () to gather (4) ɛ φ(0) x p D φ(0)x, x dx ɛ φ(0) x p dx where we used the elementary equality π π π π (cos θ) p dθ (cos θ) p dθ = q = ɛ p p + ΔG p φ(0), for any <p<, and p + q =. (cos θ) p dθ. However, now the theorem is a consequence of the Taylor expansion. In fact, given x ɛ,sinceφ C, we know that φ(x) =φ(0) + φ(0) x + D φ(0)x, x + o( x ), as x 0. Therefore, φ(0) x p φ(x) dx = φ(0) φ(0) x p dx ɛ ɛ + (5) φ(0) x p D φ(0)x, x dx + o(ɛ +p ), ɛ since by symmetry one has φ(0) x p φ(0) xdx=0, ɛ and, using () and (3), we find 0 (6) ɛ +p φ(0) x p o( x ) dx ɛ o(ɛ ) ɛ p π π (cos θ) p dθ 0 as ɛ 0. Dividing by the coefficient of φ(0) in (5) and by using (4) and (6), we obtain φ(0) = ɛ φ(0) x p φ(x) dx ɛ p ɛ φ(0) x p dx p + ΔG p φ(0) + o(ɛ ). This proves the statement of the theorem. An analogous expansion for C functions in terms of surface integrals can also be derived. Its precise expression is given in the statement of Proposition 3. below. Proposition 3.. Let Ω R be a smooth domain. Given φ C (Ω) and x 0 Ω for which φ(x 0 ) 0, we have that for any ɛ>0 such that ɛ (x 0 ) Ω it holds that (7) φ(x 0 )= ɛ (x 0 ) φ(x ɛ 0) (x x 0 ) p dx ΔG p φ(x 0 )+o(ɛ ).
MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS 459 Proof. Equation (7) is obtained by following the proof of Theorem 3. step by step, with straightforward modifications. At a point x 0 where the gradient of the function is zero, when considering the Game p-laplacian in Section 5 we will use the following classical formulas (in the spirit of the works [, 8]), which are obtained by integrating over ɛ (x 0 )and ɛ (x 0 ) the Taylor expansion. Lemma 3.3. Let φ C (Ω) R. At each point x 0 where φ(x 0 )=0, given any ɛ>0 such that ɛ (x 0 ) Ω, it holds that φ(x 0 )= φ(x) dx ɛ ɛ (x 0 ) ɛ (x 0 ) 4 ΔG φ(x 0 )+o(ɛ ) and φ(x 0 )= φ(x) dx ɛ ɛ (x 0 ) ɛ (x 0 ) ΔG φ(x 0 )+o(ɛ ). 4. Weak mean value property Manfredi et al. in [6] introduce a definition of asymptotic equality in the viscosity sense, which we use to specify in which sense we claim that a p-harmonic function verifies a mean value property. Definition 4.. Let Ω R be a smooth domain, and let x 0 Ω. We say that u C 0 (Ω) verifies the equality ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + o(ɛ ), 0) (x x 0 ) p dx in the viscosity sense, if the following conditions hold: (i) For any φ C for which u φ has a strict local maximum at x 0 Ω, there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, (8) φ(x 0 ) ɛ (x 0 ) φ(x + o(ɛ ), 0) (x x 0 ) p dx whenever φ(x) 0. (ii) For any φ C for which u φ has a strict local minimum at x 0 Ω, there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, (9) φ(x 0 ) ɛ (x 0 ) φ(x + o(ɛ ), 0) (x x 0 ) p dx whenever φ(x) 0. We are now ready to prove our representation formula. Theorem 4.. Let u be a continuous function in Ω R,andletx 0 Ω. For any ɛ>0 such that ɛ (x 0 ) Ω we have that (0) ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + o(ɛ ) 0) (x x 0 ) p dx
460 TIZIANA GIORGI AND ROERT SMITS holds in the viscosity sense if and only if u is p-harmonic, that is, if u is a viscosity solution of () Δ p 0. Proof. Assume (0) holds; we need to show that u is p-harmonic in Ω. We will show that u is a subsolution of () according to Definition.; the proof that u is a supersolution is similar. Let φ be such that u φ has a strict local maximum at x 0 Ω, and assume that φ(x 0 ) 0. Then by equation (8) in Definition 4. (8) there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, φ(x 0 ) ɛ (x 0 ) φ(x + o(ɛ ). 0) (x x 0 ) p dx On the other hand, since φ C by Theorem 3., as long as ɛ (x 0 ) Ωwehave φ(x 0 )= ɛ (x 0 ) φ(x p ɛ 0) (x x 0 ) p dx p + ΔG p φ(x 0 )+o(ɛ ). Thus for ɛ<ɛ 0 it holds that p ɛ p + ΔG p φ(x 0 )+o(ɛ ) o(ɛ ), an inequality which we divide by ɛ and let ɛ gotozerotoobtain Δ G p φ(x 0 ) 0. Finally, since φ is a C function, we can use (4) and (3) to conclude that () 0 Δ G p φ = ( Δ φ Δ G p φ ) q ΔG φ = p Δ φ p Δ G p φ, which gives that u is a subsolution of Δ p 0. Assume next that u is a viscosity supersolution. We are going to show that (9) holds. Let φ C be such that u φ has a strict local minimum at x 0 Ωand φ(x 0 ) 0. y (8) it follows that (p )Δ G φ(x 0 ) Δ φ(x 0 ) 0, which as in () gives Δ G p φ(x 0 ) 0. Hence, by Theorem 3. for any ɛ<ɛ 0,where ɛ 0 is picked so that ɛ0 (x 0 ) Ω, we obtain φ(x 0 ) ɛ (x 0 ) φ(x + o(ɛ ), 0) (x x 0 ) p dx which is exactly (9). In a similar way, if u is a subsolution, one can show that (8) holds, which proves that if u is p-harmonic, then it verifies (0) in the viscosity sense. This concludes the proof of our theorem. With the due modifications in Definition 4., one can also derive a representation for p-harmonic functions which uses only the values on the surface of the ball.
MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS 46 Theorem 4.. Let u be a continuous function in Ω R,andletx 0 Ω. For any ɛ>0 such that ɛ (x 0 ) Ω we have that ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + o(ɛ ) 0) (x x 0 ) p dx holds in the viscosity sense if and only if u is p-harmonic, that is, if u is a viscosity solution of Δ p 0. Proof. The result follows as in the proof of Theorem 4., thanks to Proposition 3.. 5. Representation for the Game p-laplacian The definition of asymptotic equality in the viscosity sense for the case of the Game p-laplacian needs to be adapted to a different definition of viscosity solution since one needs to account for the fact the the operator is singular. Definition 5.. Let Ω R be a smooth domain, and let x 0 Ω. If f is a continuous function, we say that u C 0 (Ω) verifies the equality ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + ɛ p 0) (x x 0 ) p dx p + f(x 0)+o(ɛ ), in the viscosity sense, if the following conditions hold: (i) For any φ C (Ω) for which u φ has a local maximum at x 0 Ω, there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, (3) φ(x 0 ) ɛ (x 0 ) φ(x + ɛ p 0) (x x 0 ) p dx p + f(x 0)+o(ɛ ), whenever φ(x) 0;and (4) φ(x 0 ) φ(x) dx + ɛ ɛ (x 0 ) ɛ (x 0 ) 4 f(x 0)+o(ɛ ), whenever φ(x) =0. (ii) For any φ C (Ω) for which u φ has a local minimum at x 0 Ω, there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, (5) φ(x 0 ) ɛ (x 0 ) φ(x + ɛ p 0) (x x 0 ) p dx p + f(x 0)+o(ɛ ), whenever φ(x) 0;and (6) φ(x 0 ) φ(x) dx + ɛ ɛ (x 0 ) ɛ (x 0 ) 4 f(x 0)+o(ɛ ), whenever φ(x) =0. We can then prove the following result. Theorem 5.. Let f and u be continuous functions in Ω R,andletx 0 Ω. We have that (7) ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + ɛ p 0) (x x 0 ) p dx p + f(x 0)+o(ɛ )
46 TIZIANA GIORGI AND ROERT SMITS holds in the viscosity sense in Ω if and only if u is a viscosity solution of Δ G p f(x 0 ) in Ω. Proof. Assume that (7) holds; we need to show that u is a viscosity solution. Recalling Definition.3 of viscosity solution it should be clear that part (i) of the definition of a subsolution can be obtained with almost the same proof as in Theorem 4.. If instead φ C (Ω) is such that u φ has a local maximum at x 0 Ωand φ(x 0 ) = 0, by (4) we know there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, φ(x 0 ) φ(x) dx + ɛ ɛ (x 0 ) ɛ (x 0 ) 4 f(x 0)+o(ɛ ), and by Lemma 3.3, we conclude as in Theorem 4. that Δ G φ(x 0 ) f(x 0 ). The case u supersolution is analogous. Assume next that u is a viscosity supersolution of Δ G p φ = f in Ω. We need to show that (5) and (6) are satisfied. Again (5) follows as in Theorem 4., assume then that φ C (Ω)issuchthat u φ has a local minimum at x 0 Ωand φ(x 0 ) = 0. y part (ii) of the definition of a supersolution this implies Δ G φ(x 0 ) f(x 0 ), but by Lemma 3.3 for any ɛ<ɛ 0,whereɛ 0 is picked so that ɛ0 (x 0 ) Ω, we derive φ(x 0 ) ɛ (x 0 ) φ(x 0) (x x 0 ) p dx + ɛ 4 f(x 0)+o(ɛ ), which is exactly (5). In a similar fashion, if u is a subsolution one can prove that (3) and (4) hold. We conclude this section with the corresponding theorem with surface integrals, similar to what was done in Theorem 4. for p-harmonic functions. We start by giving the appropriate definition of asymptotic equality in the viscosity sense. Definition 5.. Let Ω R be a smooth domain, and let x 0 Ω. If f is a continuous function, we say that u C 0 (Ω) verifies the equality ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + ɛ 0) (x x 0 ) p dx f(x 0)+o(ɛ ), in the viscosity sense, if the following conditions hold: (i) For any φ C (Ω) for which u φ has a local maximum at x 0 Ω, there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, φ(x 0 ) ɛ (x 0 ) φ(x + ɛ 0) (x x 0 ) p dx f(x 0)+o(ɛ ), whenever φ(x) 0;and φ(x 0 ) φ(x) dx + ɛ ɛ (x 0 ) f(x 0)+o(ɛ ), whenever φ(x) =0. ɛ (x 0 )
MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS 463 (ii) For any φ C (Ω) for which u φ has a local minimum at x 0 Ω, there exists an ɛ 0 > 0 such that for every ɛ<ɛ 0, φ(x 0 ) ɛ (x 0 ) φ(x + ɛ 0) (x x 0 ) p dx f(x 0)+o(ɛ ), whenever φ(x) 0;and φ(x 0 ) φ(x) dx + ɛ ɛ (x 0 ) ɛ (x 0 ) f(x 0)+o(ɛ ), whenever φ(x) =0. We leave to the reader the proof of the last theorem, since it is a simple modification of the proof of Theorem 5. above. Theorem 5.. Let f and u be continuous functions in Ω R,andletx 0 Ω. We have that ɛ (x 0 ) u(x 0) (x x 0 ) p u(x) dx ɛ (x 0 ) u(x + ɛ 0) (x x 0 ) p dx f(x 0)+o(ɛ ) holds in the viscosity sense in Ω if and only if u is a viscosity solution of Δ G p f(x 0 ) in Ω. References [] G. arles and P.E. Souganidis. Convergence of approximation schemes for fully nonlinear second order equations. Asymptotic Anal. 4 (99), 7-83. MR5933 (9d:3537) [] W laschke. Ein Mittelwertsatz und eine kennzeichnende Eigenschaft des logarithmischen Potentials. Leipz. er. 68 (96), 37. [3] M. G. Crandall, H. Ishii and P.-L. Lions. User s guide to viscosity solutions of second order partial differential equations. ull. Amer. Math. Soc. (N.S.) 7 (99), no., 67. MR8699 (9j:35050) [4] M. Falcone, S. Finzi Vita, T. Giorgi and R. Smits. A semi-lagrangian scheme for the Game p-laplacian via p-averaging. Submitted. [5] P. Juutinen, P. Lindqvist and J. Manfredi. On the equivalence of viscosity solutions and weak solutions for a quasi-linear equation. SIAM J. Math. Anal. 33 (00), no. 3, 699 77. MR8747 (00m:3505) [6] J. Manfredi, M. Parviainen and J. D. Rossi. An asymptotic mean value characterization for p-harmonic functions. Proc. Amer. Math. Soc. 38 (00), no. 3, 88-889. MR566554 (00k:3500) [7] Y. Peres and S. Sheffield. Tug-of-war with noise: a game-theoretic view of the p-laplacian. Duke Math. J. 45 (008), 9-0. MR459 (00i:3500) [8] I. Privaloff. Sur les fonctions harmoniques. Rec. Math. Moscou (Mat. Sbornik) 3 (95), 464-47. [9] P. Wang. A formula for smooth -harmonic functions. PanAmerican Mathematical Journal 6 (006), no., 57-65. MR86538 Department of Mathematical Sciences, New Mexico State University, Las Cruces, New Mexico 88003-800 E-mail address: tgiorgi@nmsu.edu Department of Mathematical Sciences, New Mexico State University, Las Cruces, New Mexico 88003-800 E-mail address: rsmits@nmsu.edu