calculated using the PHH basis with the aforementioned values of K , M , and β but for the J = 3 / 2 + state. The J = 3 / 2 + state does not have any bound state, however six negative eigenvalues appear, all of them satisfying | E$_{n}$ | < | E$_{d}$ | . The positive energies E 0 n are also given in Table III. As in the previous case, the corresponding eigenvectors approximately describe the N − d scattering states, though asymptotically they go to zero. It should be observed that changing the values of K , M and β , the number of these states and the corresponding energies at which the eigenvalues appear change. They do not present the stability that a true bound state shows. If we call E D 0 the bound state energy calculated using a basis of dimension D , the variational principle establishes that E D 0 ≥ E$_{3}$$_{B}$ , with E$_{3}$$_{B}$ the energy corresponding to D → ∞ . When the value of D is sufficiently high, a further increase of the dimension of the basis will not give an appreciable improvement in E D 0 , showing a pattern of convergence of the type given in Table II. On the other hand, the eigenvalues E$_{n}$ are embedded in the continuum spectrum of H which starts at E$_{d}$ . Accordingly, increasing D these eigenvalues tend to E$_{d}$ and the number of them also increases. Similarly to what it has been done in the two-body case, we can consider these states to be approximate solutions of ( H − E$_{n}$ )Ψ$_{n}$ = 0 in the interaction region and use them as inputs in the integral relation to compute second order estimates of the phase-shifts. Defining k 2 n = (4 / 3) E 0 $_{n}$/ ( ℏ $^{2}$/m ), with Ψ$_{n}$ the corresponding eigenvector and | F$_{n}$ >, | ˜ G$_{n}$ > the asymptotic functions of Eq. (41) calculated at k$_{n}$ , the second order estimate for the phase shift at each energy is obtained as The second order estimates of the phase-shifts for the six cases given in Table III are shown in Figs. 4 and 5 as solid points in the effective range functions. We can observe an extremely good agreement with the scattering calculations. This method allows for the calculation of phase-shifts using bound state type functions, even in the case of charged particles. These results can be compared to the analysis of Ref. [27], in which N − d phase- shifts were obtained solving the Faddeev equations in configuration space. For the n − d case, the results presented here and those from Ref. [27] are in complete agreement. In the p − d case the results of Ref. [27] were obtained considering the Coulomb potential in s -wave, without including the correction obtained using the KVP as was done in Ref. [16]. 16 − 2 m ℏ 2 < Ψ$_{n}$ | H − E$_{n}$ | F$_{n}$ > = B$_{n}$ 2 m ℏ 2 < Ψ$_{n}$ | H − E$_{n}$ | ˜ G$_{n}$ > = A$_{n}$ [tan δ$_{n}$ ] 2 nd = B$_{n}$ A$_{n}$ . (46)