disp(sprintf('thin_cap: %d', thin_cap));
disp(sprintf('wild errors (> 500mm): %f %%', ...
	sum(error_ary>=500)/length(error_ary)*100 ...
	));
tame = error_ary(find(error_ary<500));
disp(sprintf('tame errors (< 500mm): mean %f std %f',...
	mean(tame), std(tame) ...
	));
