Phase equilibrium calculations at high pressures have been a continuous challenge for scientists and engineers. Traditionally, this task has been performed by solving a system of nonlinear algebraic equations originating from isofugacity equations. The reliability and accuracy of the solutions are strongly dependent on the initial guess, especially due to the fact that the phase equilibrium problems frequently have multiple roots. This work is focused on the application of a subdivision algorithm for thermodynamic calculations at high pressures. The subdivision algorithm consists in the application of successive subdivisions at a given initial interval (rectangle) of variables and a systematic test to verify the existence of roots in each subinterval. If the interval checked passes in the test, then it is retained; otherwise it is discharged. The algorithm was applied for vapor-liquid, solid-fluid and solid-vapor-liquid equilibrium as well as for phase stability calculations for binary and multicomponent systems. The results show that the proposed algorithm was capable of finding all roots of all high-pressure thermodynamic problems investigated, independent of the initial guess used.