The correlated basis function theory is applied to the study of medium-heavy doubly closed shell nuclei with different wave functions for protons and neutrons and in the jj coupling scheme. State dependent correlations including tensor correlations are used. Realistic two-body interactions of Argonne and Urbana type, together with three-body interactions have been used to calculate ground state energies and density distributions of the 12 C , 16 O , 40 Ca , 48 Ca and 208 Pb nuclei.