function t = ust_prog(a)
a = abs(a);
a = sort(a,'descend');
if size(a,1) < numa
    t = a(end);
else
    t = a(numa);
end

% % % t = ustal_prog(a);
% % % 
% % % %y = zeros(size(x));
% % % %[m,n]=size(x);
% % % y = zeros(n);
% % % for j = 1:n
% % %     vect = A(:,j);
% % %     absvect = abs(vect);
% % %     I = find(absvect >= t);
% % %     y(I,j) = vect(I);
% % % end

absc = abs(c);
a = sort(absc,'descend');
if size(a,2) < numa
    t = a(end);
else
    t = a(numa);
end 

y = zeros(1,size(c,2));

I = find(absc >= t);
y(I) = c(I);