The function calculates all possible submodels and find model with the minimal value of Cp according Gilmour
Useful links