For the second harmonic generation at metal suriaces, the existing theoretial results obtained from hydrodynamics and the density functional theory are one order of magnitude deviated from the experimental data. We calculate the second harmonic generation at Al, Mg and Ag surfaces by using the correlated basis function theory. Since the electron correlation can be considered more property in this approach, our theoretical results are closer to the experimental data.