        strict feasibility of <tt>Y</tt> during the whole computation. The way to do it is similar to <tt>lbxbar</tt> and <tt>ubxbar</tt>. All such matrix variable bounds need to be listed in arrays penm.lbYbar and penm.ubYbar. Any named constraints which actually
                  don't exist (infinite bounds) are ignored.
               </li>
            </ul>
         </div><pre> penm.lbYbar = [1:penm.NY];
 penm.ubYbar = [ 3 ];</pre><div>
            <ul>
               <li>[optional] <b>box constraints 