In this paper, a quantum dot mathematical model based on a twodimensional Schrödinger equation assuming the 1/r inter-electronic potential is revisited. Generally, it is argued that the solutions of this model obtained by solving a biconfluent Heun equation have some limitations. The known polynomial solutions are confronted with new numerical calculations based on the Numerov method. A good qualitative agreement between them emerges. The numerical method being more general gives rise to new solutions. In particular, we are now able to calculate the quantum dot eigenfunctions for a much larger spectrum of external harmonic frequencies as compared to previous results. Also the existence of bound state for such planar system, in the case = 0, is predicted and its respective eigenvalue is determined. Keywords Quantum dot model · Numerov numerical method · two-electron system · Schrödinger equation PACS PACS 81.07.Ta · 78.67.Hc · 36.10.-k 1 IntroductionModern technics in nanometer-scale semiconductor manufacturing enable the creation of quantum confinement of only a few electrons. These few-body systems are often called quantum dots [1]. They can be described by a model F. Caruso