  While discrete harmonic functions have been objects of interest for quite
some time, this is not the case for discrete polyharmonic functions, as appear
for instance in the asymptotics of path counting problems. In this article, a
novel method to compute all discrete polyharmonic functions in the quarter
plane for non-singular models with small steps and zero drift is proposed. In
case of a finite group, an alternative method using decoupling functions is
given, which often leads to a basis consisting of rational functions. In a
similar manner one can obtain polyharmonic functions in the continuous setting,
and convergence between the discrete and continuous cases is proven. Lastly,
using a concrete example it is shown why the decoupling approach seems not to
work in the infinite group case.
