Among different ways to obtain the thermodynamic properties of a fluid in the liquid phase, there is the possibility to solve, numerically, a system of two differential equations where density and specific heat capacity are the unknown variables expressed as functions of temperature and pressure. By means of general methods, like Predictor-Corrector or Runge-Kutta, this system can be solved integrating along isothermal paths, once the initial conditions are known along an isobar and the speed of sound values are measured over the p-T region of interest. In this work a new perturbative method. called recursive equation method (REM), based on the analytic recursive determination of the coefficients describing density and specific heat capacity, is proposed. The main advantage of REM is the possibility to also get the value of the uncertainty associated with the solutions, obtained by means of standard methods of error analysis. Another improvement introduced by this algorithm is the possibility to solve the system of equations following arbitrary integration paths. In particular, a comparison between the results of density and heat capacity, obtained by integrating along isentropic and isothermal paths, is presented, both for water and for acetone. (C) 2009 Elsevier Ltd. All rights reserved.