function barxd(x, y, varargin)

w = (x(1)-x(2))/3;

plot([x-w; x-w], [0*y; y], varargin{:});
hold on;
plot([x-w; x+w], [y; y], varargin{:});
plot([x+w; x+w], [y; y*0], varargin{:});

end