  A novel application of the Pade approximation is proposed in which the Pade
approximant is used as an interpolation for the small and large coupling
behaviors of a physical system, resulting in a prediction of the behavior of
the system at intermediate couplings. This method is applied to quarkonium
systems and reasonable values for the c and b quark masses are obtained.
