Critical points of natural gas mixtures containing linear alkane molecules (C-1-n-C-8), branch alkanes (i-C(4)i-C-5), nitrogen (N-2), and carbon dioxide (CO2) are calculated using a criterion based upon the Helmholtz free energy. We use the GERG equation of state with its mixing rules together with needed derivatives evaluated numerically. A multivariable iterative method is used to solve the system of non-linear equations that appear in this procedure. Results are checked for stability using mechanical, diffusional and global stability tests. Our results are compared with the literature values. Maximum error for the temperature and pressure are 1.301 and 11.087% respectively. The predictive capability of the GERG-2008 equation of state is compared with PR, SRK, PSRK, SPHCT, and PC-SAFT equations of state. For binary mixtures, all equations of state are capable of representing global phase behavior types Ito IV. Results show that the GERG-2008 equation of state with binary interaction parameters is adequate for the prediction of critical points. (C) 2017 Elsevier B.V. All rights reserved.