In this paper a method for non-linear robust stabilization based on solving a bilinear matrix inequality (BMI) feasibility problem is developed. Robustness against model uncertainty is handled. In different non-overlapping regions of the slate-space known as clusters the plant is assumed to be an element in a polytope whose vertices (local models) are affine systems. In the clusters containing the origin in their closure, the local models are restricted to being linear systems. The clusters cover the region of interest in the state-space. An affine state-feedback is associated with each cluster. By utilizing the affinity of the local models and the state-feedback, a set of linear matrix inequalities (LMIs) combined with a single non-convex BMI are obtained which, if feasible, guarantee quadratic stability of the origin of the closed loop. The feasibility problem is attacked by a branch-and-bound-based global approach. If the feasibility check is successful, the Lyapunov matrix and the piecewise affine state-feedback are given directly by the feasible solution. Control constraints are shown to be representable by LMIs or BMIs, and an application of the control design method to robustify constrained non-linear model predictive control is presented. In addition, the control design method is applied to a simple example.