We take a new approach to harmonic polynomials via differentiation.Surprisingly powerful results about harmonic functions can be obtained simply by differentiating the function |x|2~" and observing the patterns that emerge.This is one of our main themes and is the route we take to Theorem 1.7, which leads to a new proof of a harmonic decomposition theorem for homogeneous polynomials (Corollary 1.8) and a new proof of the identity in Corollary 1.10.We then discuss a fast algorithm for computing the Poisson integral of any polynomial.(Note: The algorithm involves differentiation, but no integration.)We show how this algorithm can be used for many other Dirichlettype problems with polynomial data.Finally, we show how Lemma 1.4 leads to the identity in (3.2), yielding a new simple proof that the Kelvin transform preserves harmonic functions. Derivatives of |x|2_"Unless otherwise stated, we work in R" , n > 2 ; the function |x|2~" is then harmonic and nonconstant on R"\{0}.(When n = 2 we need to replace |x|2~" with log|x| ; the minor modifications needed in this case are discussed in Section 4.)Letting Dj denote the partial derivative with respect to the 7th coordinate variable, we list here some standard differentiation formulas that will be useful later:The first two formulas are valid on R"\{0} for every real t, while the last formula holds on any open set where u and v are twice continuously differentiable (and real valued); as usual, A denotes the Laplacian and V denotes the gradient.