xx = arange(start, end, step) returns numbers in the closed interval [start, end] with a step size of step.
This is similar to numpy's arange, but care is taken that the end point is included if step divides the interval evenly. Also, the number zero is protected from numerical imprecision.