function optimalMatchingDemoPlot(sources, detections, a)hold onp1 = size(sources,2);p2 = size(detections,2);for i=1:p1  text(sources(1,i), sources(2,i), sprintf('s%d', i));endfor i=1:p2  text(detections(1,i), detections(2,i), sprintf('d%d', i));endif nargin < 3, return; endfor i=1:p1  j = a(i);  if j==0 % i not matched to anything    continue  end  line([sources(1,i) detections(1,j)], [sources(2,i) detections(2,j)])endaxis_pct;