The incompressibility of both nuclear matter and finite nuclei is estimated by the monopole compression modes in nuclei in the framework of a nonrelativistic Hartree–Fock–Bogoliyubov method and the coherent density fluctuation model. The monopole states originate from vibrations of the nuclear density. The calculations in the model for the incompressibility in finite nuclei are based on the Brueckner energy–density functional for nuclear matter. Results for the energies of the breathing vibrational states and finite nuclei incompressibilities are obtained for various nuclei and their values are compared with recent experimental data. The evolution of the isoscalar giant monopole resonance (ISGMR) along Ni, Sn, and Pb isotopic chains is discussed. This approach can be applied to analyses of neutron stars properties, such as incompressibility, symmetry energy, slope parameter, and other astrophysical quantities, as well as for modelling dynamical behaviors within stellar environments.