function plot2dbar(xmean, ymean, xstd, ystd)

plot(xmean, ymean, 'r.');
hold on;
plot([xmean, xmean], [ymean-ystd, ymean+ystd], 'b-');
plot([xmean-xstd, xmean+xstd], [ymean, ymean], 'b-');

end