function y = expfun(t) y = ustep(t) .* exp(-t);