The purpose of the present work is, starting from a given assumed nucleon-nucleon interaction $V({\mathrm{r}}_{i}\ensuremath{-}{\mathrm{r}}_{j})$ and the independent-particle, central-field approximation, to deduce all the consequences of this model by means of the Hartree-Fock-Slater theory. The $V({\mathrm{r}}_{i}\ensuremath{-}{\mathrm{r}}_{j})$ is chosen to satisfy the saturation requirement and to be consistent with some of the properties of the two-nucleon system, but contains the strength factor ${V}_{0}$ as an undetermined parameter. This ${V}_{0}$ and the single-particle wave functions ${\ensuremath{\psi}}_{n}({\mathrm{r}}_{i})$ are determined by the variational principle together with the requirement that the total binding energy of the nucleus be equal to its empirical value. The binding energies of the individual nucleons in the various shells are themselves approximately given by the Fock equations which also lead to a central field which is different for the different shells. The average central field implied in the usual treatment of the shell model is, however, the same for a nucleon in any shell, and this must be identified with some approximate average field obtained by a procedure such as that suggested by Slater for an electron in an atom or in a metal. On the other hand, the central field in the sense of Hartree (i.e., obtained from the Fock theory by neglecting the exchange terms) would be a very poor approximation as the exchange terms are not negligible compared with the direct terms. It is emphasized in the present work that the application of the variational principle to the problem rids the shell model of the inconsistent procedure in the usual treatments in which two independent assumptions concerning $V({\mathrm{r}}_{i}\ensuremath{-}{\mathrm{r}}_{j})$ and the average central field $V(\mathrm{r})$ are made. A comparison of the result of the present program with the empirical facts will form a correct basis on which to judge the fairness or failure of the central-field approximation in the shell model.Numerical calculations have been carried out for the ${\mathrm{O}}^{16}$, ${\mathrm{O}}^{15}$, and ${\mathrm{O}}^{17}$ nuclei to illustrate the ideas of the self-consistent treatment of the independent-particle, central-field nuclear model. In a first approximation, the ${V}_{0}$ and the wave functions ${\ensuremath{\psi}}_{n}({\mathrm{r}}_{i})$ are determined by the modified Ritz method. The binding energies and the effective "central fields" for the individual nucleons in the various shells are then calculated from the Fock equations by means of these ${\ensuremath{\psi}}_{n}({\mathrm{r}}_{i})$. It is found that the parameter ${V}_{0}$ in ${\mathrm{O}}^{16}$, ${\mathrm{O}}^{17}$, and ${\mathrm{O}}^{15}$ comes out to be very nearly the same, that the Fock central potentials are considerably lower than the Hartree potentials; that the former exhibit a general shape of a wine-bottle with a diffuse boundary; and that the "binding energy of the last nucleon" does not come out very well from the Fock equations.