In this paper we analyze the effect of introducing a numerical integration in the piecewise linear finite element approximation of the Steklov eigenvalue problem. We obtain optimal order error estimates for the eigenfunctions when this numerical integration is used and we prove that, for singular eigenfunctions, the eigenvalues obtained using this reduced integration are better approximations than those obtained using exact integration when the mesh size is small enough.