The equilibrium shape of the liquid bridge interface is analyzed theoretically and experimentally. Both axisymmetric and nonaxisymmetric perturbations are considered. The axisymmetric deviations are those related to volume effects, the difference between the radii of the disks, and the axial forces acting on the liquid bridge. The nonaxisymmetric deviations are those due to the eccentricity of the disk and the action of lateral forces. The theoretical study is performed using three different techniques: (i) an analytical expansion around the cylindrical solution, (ii) a finite difference scheme, and (iii) an approximate numerical approach valid only for slight nonaxisymmetric deviations. The results of the three methods are compared systematically. There is a very good agreement between the analytical and the numerical approaches for contours which are close to cylindrical, and the agreement extends to configurations with only moderate deviations from cylindrical. Experiments are performed using the so-called neutral buoyancy or plateau technique. Theoretical and experimental contours are compared considering a wide range of values for the parameters characterizing the perturbations. In general, the finite difference method provides reasonably accurate predictions even for large deviations of the liquid bridge contour from cylindrical.