In this paper, at first we have formulated the lowest order constrained variational method for the relativistic case of an interacting fermion system at finite temperature. Then we have used this formalism to calculate some thermodynamic properties of liquid in the relativistic regime. The results show that the difference between total energies of relativistic and non-relativistic cases of liquid decreases by increasing the density. On the other hand, at densities smaller than , which is close to the theoretical saturation point of the system, this difference increases as the temperature increases. We have found that the relativistic calculations lead to the thermodynamic characteristics for liquid which are in a better agreement with the experimental data with respect to those of non-relativistic calculations.
