ct feasibility. To do so, it is  
% necessary to name all such box constraints in arrays penm.lbxbar, 
% penm.ubxbar. Any named box constraints which actually don't exist
% (infinite bounds) are ignored.
%
%   penm.lbxbar = [ 3, 4 ];
%   penm.ubxbar = [1:penm.Nx];
% 
% * *lower, upper bounds on elements of |Y| variable* If used, it is a cell
% array of |m| matrices of the same size as matrix variables
% |Y{1},...,Y{m}|. Use empty matrix when constraints are not defined.
% If not present, the elem