function doublehist(xx,y1,y2)
width=(xx(end)-xx(1))/(length(xx)-1);
space=.45;
shift=space*width/2;
bar(xx-shift,y1,space,'r')
hold on
bar(xx+shift,y2,space,'b')
hold off 

