impls_fn(t, delta)

%	A function generating a rectangular approximation for
%	the unit impulse function; width is delta
%
function imp = impls_fn(t, delta)
imp = zeros(size(t));
imp = (stp_fn(t+delta/2)-stp_fn(t-delta/2))/delta;