We present here the momentum dependent local-ansatz wavefunction approach (MLA) with the best choice of the self-consistent variational parameters to describe the ground-state properties of the correlated electron systems in solids. With use of the self-consistent variational scheme, we performed the numerical calculations for the half-filled band as well as non-half-filled band Hubbard model on the hypercubic lattice in infinite dimensions. The ground-state energy in the MLA is lower than those of the local-ansatz approach (LA) and the Gutzwiller approach (GA) in weak and intermediate Coulomb interaction regimes. The double occupation number is suppressed as compared with the LA. We observe the distinct momentum dependence of the momentum distribution functions which is qualitatively different from those of the LA and the GA.