The geodesic flow of the BGPP metric is Liouville integrable

We prove that the geodesics equations corresponding to the BGPP metric are integrable in the Liouville sense. The $\mathrm{SO}(3,\mathbb{R})$ symmetry of the model allows to reduce the system from four to two degrees of freedom. Moreover, solutions of the reduced system and its degenerations can be given explicitly or reduced to a certain quadrature. In degenerated cases BGPP metric coincides with the Eguchi-Hanson metric and for this case the mentioned quadrature can be calculated explicitly in terms of elliptic integrals.