A program to compute Birkhoff normal forms of symplectic maps in R4