“…The numerical approximation of polyharmonic problems has been first addressed in the eighties by [18] and, more recently, in [13,40,45,42,36]. It is worth mentioning an increasing interest in the numerical approximation of models involving high-order differential operators, e.g., [46,47,41] in the context of sixth order Cahn-Hilliard equations. To the best of our knowledge, the conforming VEM proposed in this article is the first work addressing the approximation of arbitrary-order polyharmonic problems on polygonal meshes.…”