Deutsch: 2*stetig diff'bare Hutfunktion b(t) = (t-t_0+\epsilon)^3*(t_0-t+\epsilon)^3 für t\in(t_0-\epsion,t_0+\epsilon)

Varianten dieses Bildes könnten als Approximationen der Dirac-Delta-Distribution dienen. Hier der Quelltext:

eps = 0.1;
x0 = 1;
f(x) = ( abs(x-x0) < eps ? (x-x0+eps)**3 * (x0-x+eps)**3 : 0 );
set samples 1000;
set xrange [x0-2*eps:x0+2*eps];
set yrange [-0.1*eps**6 : 1.1*eps**6];
set xlabel "t";
set ylabel "b(t)";
set notitle;
plot f(x) notitle;
set term png;
set output "Variationsrechnung_hut_polynom6Ordnung.png";
repl;
set term x11;