Finding normal forms of nonintegrable Hamiltonians is a problem encountered in many branches of physics. In this article, a symbolic program is presented which uses a succession of canonical transformations to construct the Birkhoff-Gustavson normal form of a nonseparable Hamiltonian.