b A2*xrob]' lambda*eye(2) zeros(2,1); <span class="keyword">...</span>
         [A0*xrob-b]' zeros(1,2) t] &gt;= 0;
cvx_end

<span class="comment">% Generate Random Trials</span>
notrials=100000;
r = sqrt(rand(notrials,1));     <span class="comment">% random on [0,1] with pdf g(r) = 2r;</span>
theta = 2*pi*rand(notrials,1);  <span class="comment">% uniform on [0,2pi]</span>
v = [r.*cos(theta)  r.*sin(theta)];
ls_res = zeros(1,notrials);
rob2_res = zeros(1,notrials);
rob_res = zeros(1,notrials);
tych_res = z