Numerical renormalization group study of the Kondo effect in superconductors