function v = z_hit(a,b,step)
%

v = a:step:b;
num_pts = length(v);
for p = 1:num_pts
    v(p) = quad(@(x)f_norm(x,v(p)),-100,0);
end
