Recursive computation of the Frobenius-Padé approximants: theory, stability and conditioning