          (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,2) = cvx_optval;
    fprintf( 1,<span class="string">' %g\n'</span>, round(cvx_optval) );

<span class="keyword">end</span>

figure(2);
hold <span class="string">off</span>

<span class="comment">% plot data pt and contour line through it</span>
val = 1.1*sqrt(0.5)+ 0.8*sqrt(.5);   <span class="comment">% value at center</span>
t = linspace((