The Percus–Yevick (PY) and mean spherical model (MSM) approximation integral equations are used to calculate the radial distribution functions of liquid sodium, rubidium, and aluminium. The resulting radial distribution functions are compared with molecular dynamics (MD) simulation results. The agreement is satisfactory.