function output=unitstep(time)
output=0.5*(ones(size(time))+sign(time));
