p-homomorphisms between the first jet prolongation of a harmonic module and an arbitrary harmonic module are described using Casimir type computations