multiple harmonic function