%
% unit step function
%
function u = stp_fn(t)
  eps = 1e-10;
  u = 0.5*(sign(t+eps)+1);
return;

