Computing the Rao distance for Gamma distributions

This paper deals with the calculation of the Rao distance between Gamma distributions, which is the Riemannian distance induced by the Fisher information matrix on the Gamma statistical model. In this case no closed form expression of the Rao distance is available and a numerical approach is thus necessary to compute the distance. A computer program based on a simple shooting algorithm is presented and discussed.