close to the boundary)
% * if option setting xinit_mod = 1, in which case the point is adjusted
%   always so that it lays within box constraints
%
%
% *Constraints*
%
% * Linear and nonlinear function constraints
%
%   penm.NgNLN = 2; 
%   penm.NgLIN = 1;
% 
% * Linear and nonlinear matrix constraints
%
%   penm.NANLN = 0; 
%   penm.NALIN = 1;
% 
% * Lower/upper bounds ... must be defined as finite at least on one
%   side; equality constraints defined by equal bounds 
%
%   penm.lbg = 