  Correlated Basis Function theory and Fermi Hypernetted Chain technique are
extended to study medium-heavy, doubly closed shell nuclei in j-j coupling
scheme, with different single particle wave functions for protons and neutrons
and isospin dependent two-body correlations. Central semirealistic interactions
are used. Ground state energies, one-body densities, distribution functions and
momentum distributions are calculated for 12C, 16O, 40Ca, 48Ca and 208Pb
nuclei. The values of the ground state energies provided by isospin dependent
correlations are lower than those obtained with isospin independent
correlations. In finite nuclear systems, the two--body Euler equations provide
correlation functions variationally more effective than those obtained with the
same technique in infinite nuclear matter.
