Based on the Gross-Pitaevskii(G-P) energy functional and variational method, we present a new method to find the ground state wave function for Bose-condensed gas in a harmonic trap at zero temperature. With this method we are able to find the analytic expression of ground state wave function and explore the relevant quantities, such as energy, chemical potential and so on. These results agree well with previous ground state numerical solutions of the G-P equation given by Dalfovo et al. Under the large number of particles N limit for repulsive interaction, our results and Thomas-Fermi approximation model approach to reaching unanimity. This new method is simple compared to other methods used to solve numerically the G-P equation, and one can obtained the analytic and reliable results.