In this paper we are concerned with energy-conserving methods for Poisson problems, which effectively solved by defining a suitable generalization of HBVMs, class Hamiltonian problems. The actual implementation the is fully discussed, particular emphasis on conservation Casimirs. Some numerical tests reported, in order to assess theoretical findings.