function y = hsinc(t) % Hilbert transform of sinc y = (cos(pi*(t + (t == 0)*eps)) - 1) ./ pi ./ (t + (t == 0)*eps);