1 1.1]) ;
axis <span class="string">equal</span>;
title(<span class="string">'Quadratic placement problem'</span>);
<span class="comment">% print -deps placement-quadr.eps</span>

figure(2)
all = [free_sum; fixed];
bins = 0.05:0.1:1.95;
lengths = sqrt(sum((A*all).^2')');
[N2,hist2] = hist(lengths,bins);
bar(hist2,N2);
hold <span class="string">on</span>;
xx = linspace(0,2,1000); yy = (4/1.5^2)*xx.^2;
plot(xx,yy,<span class="string">'--'</span>);
axis([0 1.5 0 4.5]);
hold <span class="string">on</span>
plot(