function ne=find_sat(nb,I)
if nb+1 >length(I)
    ne=nb;
    return
else 
  for i=nb:length(I)
    if i==length(I)
    ne=i;
    return
    else
      if I(i+1)~=I(i)+1
      break
      end
    end
  end
  ne=i; 
end