We investigate the quasinormal modes (QNMs) of the scalar field coupled to the Einstein's tensor in the non-commutative geometry inspired black hole spacetime. It is found that the lapse function of the non-commutative black hole metric can be represented by a Kummer's confluent hypergeometric function, which can effectively solve the problem that the numerical results of the QNMs are sensitive to the model parameters and make the QNMs values more reliable. We make a careful analysis of the scalar QNM frequencies by using several numerical methods, and find that the numerical results obtained by the new WKB method (the Padé approximation) and the Mashhoon method (Pöschl-Teller potential method) are quite different from those obtained by the AIM method (the asymptotic iterative method) and time-domain integration method when the parameters θ and η are large. The most obvious difference is that the numerical results obtained by the AIM method and the time-domain integration method appear a critical η c with an increase of η, which leads to the dynamical unstability. After carefully analyzing the numeral results, we conclude that the numerical results obtained by the AIM method and the time-domain integration method are closer to the theoretical values than those obtained by the WKB method and the Mashhoon method, when the θ and η are large. Moreover, through a numerical fitting, we obtain that the functional relationship between the threshold η c and the non-commutative parameter θ satisfies η c = aθ b + c for a fixed l approximately.