function y = p(t,tau)
    
    y = u(t,-tau/2) - u(t, tau/2);

end
