On the gravitational field of static and stationary axial symmetric bodies with multi-polar structure

We give a physical interpretation to the multi-polar Erez-Rozen-Quevedo solution of the Einstein Equations in terms of bars. We find that each multi-pole correspond to the Newtonian potential of a bar with linear density proportional to a Legendre Polynomial. We use this fact to find an integral representation of the $\gamma$ function. These integral representations are used in the context of the inverse scattering method to find solutions associated to one or more rotating bodies each one with their own multi-polar structure.