f=fopen('netdiff.out');
A=fscanf(f,'%d');
fclose(f);
steiner=zeros(size(A,1)/2,1);
actual=zeros(size(A,1)/2,1);
for idx=1:size(A,1)/2,
    steiner(idx)=A(2*idx-1);
    actual(idx)=A(2*idx);
end
err=steiner-actual;
hist(err,100);