%
%  implement the unit step function 
%     u(t+t0) = 1 for t+t0>=0
%
function x = unit_step(t,t0);

x = 1*((t+t0)>=0);

return;
