  Recently, it has been recently shown that the linear response theory in
symmetric nuclear matter can be used as a tool to detect finite size
instabilities for different Skyrme functionals. In particular it has been shown
that there is a correlation between the density at which instabilities occur in
infinite matter and the instabilities in finite nuclei. In this article we
present a new fitting protocol that uses this correlation to add new additional
constraint in Symmetric Infinite Nuclear Matter in order to ensure the
stability of finite nuclei against matter fluctuation in all spin and isospin
channels. As an application, we give the parameters set for a new Skyrme
functional which includes central and spin-orbit parts and which is free from
instabilities by construction.
