mputing the solution of the QP for the case where a^T*x0 &lt;=b...Done! 
Computing the solution of the QP for the case where a^T*x0 &gt; b...Done! 
-----------------------------------------------------------------
Verifying that p_C(x0) and x0_star are equal in the case where a^T*x0 &lt;=b
||p_C(x0) - x0_star|| = 6.6309e-10
Hence they are equal to working precision
Verifying that p_C(x1) and x1_star are equal in the case where a^T*x1 &gt; b
||p_C(x1) - x1_star|| = 2.0401e-10
Hence they are equal to working 