function y=heaviside(t)
% Esim:
% t=[-2,-.0001 .0001 2]; plot(t,heaviside(t));grid;ylim([-.2,1.2]);shg
y=t>0;

