large, and it is necessary to introduce a new renormalization method, in which the counter terms remove the unnaturally large terms and make the terms in a natural size. In other words, we assume that fitting polynomial functions are represented as a natural power series at the low energy region, and to maintain such polynomial functions, large cancellations for the first and second terms with the $r_l$ and $P_l$ effective terms, respectively, are expected. So we include the three effective range parameters, $r_l$, $P_l$, and $Q_l$, for the $l=0$ channel, as the counter terms. The same tendency can be seen in the $l=1,2$ channels whereas one needs four effective range parameters for the $l=3$ channel. Thus, we employ the three effective range parameters, $r_l$, $P_l$, $Q_l$ for the $l=0,1,2$ channels and the four effective range parameters, $r_l$, $P_l$, $Q_l$, $R_l$ for the $l=3$ channel when fitting the parameters to the phase shift data below.[^4]

[**4. Numerical results**]{}

As discussed above, we fit the three effective range parameters, $r_l$, $P_l$ and $Q_l$ to the phase shift data for the $l=0,1,2$ channels and the four effective range parameters, $r_l$, $P_l$, $Q_l$, and $R_l$ to those for the $l=3$ channel, while $a_l$ are 