Feedback linearization suffers from a number of restrictions which have limited its use in Model-based Predictive Control. Some of these restrictions do not apply to the case of bilinear systems, but problems with input constraints and unstable zero dynamics persist. The present paper overcomes these difficulties by means of an interpolation strategy. Involved in this interpolation is a stabilizing trajectory which is computed through the use of invariant feasible sets (defined for the biline...