function plerr(aa,vf)

hold on; for i=1:length(aa),
 plot([aa(i) aa(i)],[vf(i)*(1-.1) vf(i)*(1+.1)]); end; hold off

