Subtractive renormalization of the next-to-leading order NN interaction ⁎

The subtracted kernel method is implemented recursively to solve scattering equations for the 1S0 phase-shifts considering the leading and the next-to-leading order NN interaction.