Fast and stable computation of the barycentric representation of rational interpolants