The Hopf algebra structure of Feynman graphs has been explored extensively in the last years. It proved to be valuable for the analytic computation of Feynman amplitudes by means of systematic parametric integration techniques and could lead to new non-perturbative results in the scope of Dyson-Schwinger equations. feyngen and feyncop were developed to provide input for the powerful new techniques. feyngen relies on the established nauty package by McKay for the fast generation of higher loop diagrams. feyncop can be used to calculate the coproduct of graphs on the Hopf algebra of Feynman graphs. This coproduct encodes the BPHZ algorithm necessary to evaluate the finite amplitude in a mathematically sound framework and fits well into the world of Dyson-Schwinger equations. In the framework of these equations certain identities can be obtained which were used to validate the two programs.