function run_z_hit
%

v = z_hit(-10,10,0.1);
v = v/sum(v);
plot([-10:0.1:10],v);
