11]*xhub, <span class="string">'-'</span>, <span class="keyword">...</span>
     [-11; 11], [1 -11; 1 11]*xls, <span class="string">'--'</span>);
axis([-11 11 -20 25])
title(<span class="string">'Least-square fit vs robust least-squares fit (Huber-penalized)'</span>);
xlabel(<span class="string">'x'</span>);
ylabel(<span class="string">'y'</span>);
legend(<span class="string">'Data points'</span>,<span class="string">'Huber penalty'</span>,<span class="string">'Regular LS'</span>,<span class="string">'Locat