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 liqu...