A new variational method proposed by Neveu is applied to the Gross-Neveu model. Chiral condensate is computed perturbatively at zero and finite temperature by the use of effective potential. At both temperatures the expected results in the mean field analysis, nonzero condensate at T=0 and restoration of γ5 symmetry at T≠0, are recovered in an efficient manner by calculating only a few lowest order diagrams.