The Schrödinger equation with pseudo-Gaussian potential is investigated. The pseudo-Gaussian potential can be written as an infinite power series. Technically, by an ansatz to the wave-functions, exact solutions can be found by analytic approach . However, to calculate the solutions for each state, a condition that will stop the series has to be introduced. In this way the calculated energy values may suffer modifications by imposing the convergence of series. Our presentation, based on numerical methods, is to compare the results with those obtained in the analytic case and to determine if the results are stable under different stopping conditions.