Internally redefine Heaviside function

Internally define the function h to be a sigmoid function similar in shape to the Heaviside function. Then, test the correct answer over an appropriately large domain using a lot of pre-specified sample points.