The integral equation for the inhomogeneous electron correlation function is established on the basis of the many-body wave function of metal surface, and the divergent terms in this equation caused by the long-range correlation are eliminated. Then we get the numerical results of the correlation function for Al surface in the first iteration.