The dependence of the interaction force between two undulating phospholipid bilayers and of the rootmean-square fluctuation of their separation distances on the average separation can be determined once the distribution of the intermembrane separation is known as a function of the applied pressure. However, most of the present theories for interacting membranes start by assuming that the distance distribution is symmetric, a hypothesis invalidated by Monte Carlo simulations. Here we present an approach to calculate the distribution of the intermembrane separation for any arbitrary interaction potential and applied pressure. The procedure is applied to a realistic interaction potential between neutral lipid bilayers in water, involving the hydration repulsion and van der Waals attraction. A comparison with existing, experiments is provided.