In the present work, the neutron skin thickness of neutron-rich closed-shell nuclei 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb, has been calculated using self-consistent Hartree-Fock (HF) with 20 types of Skyrme interactions: SKM*, SQMC650, SV-K218, SQMC700, SkO, NRAPR, KDE0, T31, Skχm∗, v075, SkSC10, RATP, SAMi, SK255, SkD, GS2, SIV, QMC2 and SIII of deferent nuclear matter properties. Neutron skin thickness (Rskin) can be obtained through the proton and neutron density distribution which is the difference between the neutron and proton root-mean-square radii, and depends on the properties of the nuclear matter. We examined the sensitivity of neutron skin thickness (Rskin) to nuclear matter properties, such as symmetry energy density J (MeV) and the symmetry energy coefficient L (MeV), as well as we deduced the value of neutron skin thickness (Rskin) as function of isospin asymmetric ratio I=((N-Z))⁄A. The present work was supported by statistical calculations and determining and discussion the Pearson linear correlation coefficient.