cosmo tail hdr

function [tail_values, idxs] = cosmo_tail(values, to_select)