g precision'</span>);
disp(<span class="string">'Verifying that p_C(x1) and x1_star are equal in the case where a^T*x1 &gt; b'</span>);
disp([<span class="string">'||p_C(x1) - x1_star|| = '</span> num2str(norm(xs1 - pc_x1))]);
disp(<span class="string">'Hence they are equal to working precision'</span>);
</pre>
<a id="output"></a>
<pre class="codeoutput">
Computing the analytical solution for the case where a^T*x0 &lt;=b...Done! 
Computing the analytical solution for the case where a^T*x0 &gt; b...Done! 
Co