It is probably amplified background numerical noise due to machine-number precision. The z
values are or the order of 10^-16
. Compare:
Plot3D[(t/(2 e)) (HeavisideTheta[t + e] -
HeavisideTheta[t - e]), {t, -3, 3}, {e, -2, 2},
AxesLabel -> {t, e}, PlotRange -> All]
The automatic value for PlotRange
comes out abnormally small in this case.