0;
            ones(m+1,1)*u' &lt;= u*ones(1,m+1)+(g_x*ones(1,m+1)).*<span class="keyword">...</span>
              (ones(m+1,1)*data(:,1)'-data(:,1)*ones(1,m+1))+<span class="keyword">...</span>
              (g_y*ones(1,m+1)).*(ones(m+1,1)*data(:,2)'-data(:,2)*ones(1,m+1));
            (u*ones(1,m+1)).*Pweak &gt;= (ones(m+1,1)*u').*Pweak;
    cvx_end
    bounds(k,1) = cvx_optval;
    fprintf( 1,<span class="string">'%g'</span>, round(cvx_optval) );

    <span class="comment">% Check for u(0.5,0.5) &gt;= u