x(4);<span class="keyword">...</span>
      x(1)*x(4);<span class="keyword">...</span>
      x(1)*x(4)+1;<span class="keyword">...</span>
      x(1)^2+x(1)*x(2)+x(1)*x(3);<span class="keyword">...</span>
      0; 0], userdata);
  penm.objhess = @(x,Y,userdata) deal([2*x(4) x(4) x(4) 2*x(1)+x(2)+x(3) 0 0;<span class="keyword">...</span>
      x(4) 0 0 x(1) 0 0;<span class="keyword">...</span>
      x(4) 0 0 x(1) 0 0;<span class="keyword">...</span>
      2*x(1)+x(2)+x(3) x(1) x(1) 0 0 0;<span class="keyword"