function y = r(t,to)

    y=(t-to).*(t>=to);
end
