“…Clearly, there is a price for this because total derivatives of the f function are involved in the difference equation defining the method, and thus, a suitable smoothness requirement for f is necessary. Multiderivative methods have been considered often in the past for the numerical treatment of ODEs, for example also in the context of boundary value methods [2], and in the last years, there has been a renewed interest in this topic, also considering its application to the numerical solution of differential algebraic equations; see, e.g., [3][4][5][6][7][8]. Here, we consider the numerical solution of Hamiltonian problems which in canonical form can be written as follows:…”