when we use the transformation u = e/h, then the density is not correctly normed
it should be ``1/h * norm.pdf(e/h)``

I don't see from the code whether normalizing the kernel to be a proper pdf is necessary or not, or if it is already done.