The method of large-N expansion has been applied to an attractive radial Gaussian potential to obtain its bound state energy levels. It has been shown that the shifted 1/N expansion technique provides a better approximation than the method of unshifted 1/N expansion and yields energy values which are in good agreement with the accurate numerical and analytic results.