<tt>Y{i}</tt>, the starting point <tt>Yinit{i}</tt> must be a symmetric matrix of the same non-zero pattern as <tt>Y{i}</tt>.
         </p><pre> penm.xinit=rand(penm.Nx,1);
 penm.Yinit{1} = 3.*eye(n);</pre><p>The starting point might be automatically modified within solve() * if barrier (for example, penm.lbxbar) is used and the
            point is not   feasible w.r.t. the box constraints (or very close to the boundary) * if option setting xinit_mod = 1, in which
            case the point is adjusted   a