We propose a numerical algorithm based on a discrete variable representation and shifted inverse iterations and apply it to for the analysis of the bound states of edge dislocation modelled by a quantum dipole in a plane. The good agreement with results of recent papers of Amore et al [J. Phys. B 45, 235004 (2012)] was obtained. The error estimates of the previous results of low-lying states energies of other authors were not known due to limitations of the variational approaches and this paper fills this gap presenting calculated low-lying bound states energies by non-variational technique. The probability densities of low-lying states were calculated.