A method of determination of thermodynamic properties of metastable and unstable fluid states by molecular simulation is presented. As illustrative examples spinodal and binodal states are calculated for the Lennard-Jones fluid, ethane, and n-butane using the optimized potentials for liquid simulations for the representation of the real molecules.