See if this does what you want:
x = linspace(0, 10, 10);
y = randi(9, 1, 10);
ybars = [2 6];
figure(1)
hp = plot(x, y, 'bp');
hold on
patch([min(xlim) max(xlim) max(xlim) min(xlim)], [ybars(1) ybars(1), ybars(2) ybars(2)], [0.8 0.8 0.8])
plot(x, y, 'bp')
hold off
axis([0 10 0 10])
You will have to modifiy it to do what you want, but that should be straightforward. Note the repeated plot call to prevent the patch from hiding it.
Best Answer