A generalized predictive controller for discrete-time MIMO bilinear processes is derived based on a new bilinear CARIMA model with integral action. Using this model a multi-step ahead optimal predictor is derived. Two alternative solution methods-rigorous and short-cut-of the minimization problem of a long-range cost-function are introduced and used to calculate control input. Several simulation results show that the proposed control methods have robustness to the limited prior knowledge of the process.