A practical method based on Mayer's cluster expansion to calculate critical values for a quark-gluon plasma (QGP) phase transition to hadrons is represented. It can be applied to a high-temperature QGP for clustering of quarks to mesons and baryons. The potential used is the Cornell potential, i.e., a potential containing both confining and gluon exchange terms. Debye screening effects are included. An equation of state (EOS) for hadron production is found by analytical methods, which is valid near the critical point. The example of the formation of $J/\ensuremath{\psi}$ and \ensuremath{\Upsilon} is recalculated. It is shown that in the range of temperatures available by today's accelerators, the latter particles are suppressed. This is further confirmation for heavy quarkonia suppression and, hence, for a signature of a QGP. The EOS presented here also shows that in future colliders there will be no heavy quarkonia production by the mechanism of phase transition. Hence, if there will be heavy quarkonia production, it must be based on some other mechanisms, perhaps on the basis of some recently suggested possibilities.