function s = singleton(x,y,a) n = find(x>=a,1,'first'); if(n<=1) s = y(n); else s = (y(n)+y(n-1))*0.5; end end