%
%   Discrete Unit Pulse Function
%

function Dp = drect(n,L)
    Dp = n*0;
    Dp(find(n==0)-(L-1)/2:find(n==0)+(L-1)/2)=1;
end