function z = vline(x,y,s)

% vline(x,y,s) returns a vertical line of type s at x=x from 
% (the largest value of the vector y+0.5) to (the smallest value of y-0.5)

z = plot([x,x],[min(y)-1,max(y)+1],s);
