function y = sech2pulse(t,p,tau)
%   y = sech2pulse(t,p,tau)
% Unit-height Sech2-shaped Pulse, where t is the independent variable, p is
% the position of the maximum on the t axis, and tau is the width of the
% pulse. See https://www.rp-photonics.com/sech2_shaped_pulses.html
   % Example:
   % x=-10:.1:10;
   % p=0;
   % tau=2;
   % y=sech2pulse(x,p, tau);
   % plot(x,y,'.')
   % grid
   %
   y = sech((t-p) / tau).^2;
end
