 
% *Define box constraints on variables* 
%
% *  *box constraints* (lower/upper bounds) on the variables, |lbx <= x <=
% ubx|. If |lbx = -Inf|, it is not considered; similarly for |ubx = Inf|. So,
% e.g., by setting |lbx=-Inf|, |ubx=Inf| the variable is unbounded. Length of
% each array should be |Nx|. If one or both are missing, 'x' is considered as
% unconstrained from that side. Any |lbx>ubx| is considered as error. If
% lbx==ubx the variable will be fixed and won't have further effect. 
%
%  