This is not correct. The discrete form of the Heaviside step function has
the value 1 for x == 0.
Yeah, I was looking at it and wondering if I'd misremembered the definition. Assuming you're implementing the discrete Heaviside function, H[0] = 1 as Sturla notes.