ラジアンにし忘れたのかなと思っています。
deg2radを使用しました
nt = 2^11; % FFT points
Tmax = 2^8; % window size
T0 = 25; % Pulse width, ps
dtau = (Tmax*T0)/nt; % step size in tau
tau = (-nt/2:nt/2-1)*dtau; % time array
tau_r = deg2rad(tau);
plot(tau_r, sech(tau_r))