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.