Thomas--Fermi calculation of the interlayer force in graphite

A model of a graphite crystal is proposed in which planar layers of positive charge are considered instead of the point charges of nuclei. The interlayer electronic density is calculated integrating both the Thomas-Fermi and the Thomas-Fermi-Dirac equations. From these densities, the total energy of the electrons is calculated including corrections for inhomogeneity in the form of Weizs\"acker and Kirzhnits. The influence of the different corrections is studied with the result that the best method is to calculate the density from the Thomas-Fermi-Dirac equation and to take into account the inhomogeneity corrections in the form of Kirzhnits.