We present an algorithm for automatic computation of generalized Zhang-Zhang (GZZ) polynomials of benzenoids. The Fortran 95 implementation of the presented algorithm has been used for the determination of GZZ polynomials for several basic classes of benzenoids. Correctness of the derived formulas has been demonstrated using standard decomposition recurrence relations.