Different thermodynamic properties (the vapour density, the liquid density and the saturation pressure) were calculated by the model based on the Nakamura-Breedveld-Prausnitz equation of state (NBP EOS). Since the original form of the NBP EOS often generates inaccurate results for liquids, it was modified to describe this phase better. The calculations were realized in the subcritical region. So far, the temperature-dependent NBP EOS parameters have been obtained by special correlations. Their constants were fitted to a lot of experimental data. In this paper the equation of state temperature-dependent parameters were obtained by a new method which was based at piecewise cubic Hermite interpolating polynomials (PCHIPs). In the proposed method some experimental data (called the key ones) were used, thus reducing the experimental effort. Seven substances were chosen for the test calculations. Each of them is common in industry. The calculation results were compared with the experimental data. The new method has made an accurate description of vapour-liquid equilibrium for the considered pure substances over a wide temperature range possible.