Milling stability prediction using an improved complete discretization method