function w=rectangleValue(t, A,Tc, tau)
% function w=rectangleValue(t, A, Tc, tau)
% 
% input: A(=sqrt(Pulse Energy)), Tc(Time shift): [s]
%          tau(Impulse Width) [s]
% ouput: t(time) :[s], w(wave form)
%

if t>Tc-tau/2 & t<=Tc+tau/2
    w=A/sqrt(tau);
else
    w=0;
end

return;