%?Use?this?program?to?reproduce?Fig.?5.2?of?text
close?all
clear?all
eps?=?0.000001;
taup?=?3;
[x]?=?single_pulse_ambg?（taup）;
taux?=?linspace（-tauptaup?size（x1））;
fdy?=?linspace（-5/taup+eps5/taup-eps?size（