function y=shrink(th,x)
% shrink - Shrink values towards zero
%
% shrink :: nonneg~'threshold', [Size] -> [Size].

y=max(0,x-th);
