For the first time, Fermi HyperNetted Chain (FHNC) techniques have been applied to describe the ground-state of medium and heavy doubly closed shell nuclei, with fully realistic nuclear interactions, including both two- and three-body forces, and operator dependent correlation functions. Calculations for the 12C, 16O,40Ca, 48Ca and 208Pb nuclei, have been done by using Argonne V8' two-body potential together with Urbana IX three-body force. These calculations reach an accuracy comparable to that of the best nuclear matter variational calculations. We have also investigated the effects produced by the short range correlations (SRC) on some ground state quantities related to observables.