_pos(norm(xs1 - x1)) )
    a'*xs1 &lt;= b;
cvx_end
fprintf(1,<span class="string">'Done! \n'</span>);

<span class="comment">% Verification</span>
disp(<span class="string">'-----------------------------------------------------------------'</span>);
disp(<span class="string">'Verifying that p_C(x0) and x0_star are equal in the case where a^T*x0 &lt;=b'</span>);
disp([<span class="string">'||p_C(x0) - x0_star|| = '</span> num2str(norm(xs0 - pc_x0))]);
disp(<span class="string">'Hence they are equal to workin