A program to compute Birkhoff normal forms of symplectic maps in

A program to compute Birkhoff normal forms for symplectic maps in R4 is described. We consider the case of an elliptic fixed point, which is the most relevant for applications. We compute the normal forms both in the nonresonant and resonant case and we provide the interpolating Hamiltonian and the normalizing transformation.