With \( R \) you can then evaluate the phase shifts by noting that $$ R(k_{N+1},k_{N+1})=R(k_0,k_0)=-\frac{tan\delta}{mk_0}, $$ where \( \delta \) are the phase shifts.