for x1,x2 in zip(domain[:-1], domain[1:])

OT, to be sure, but: Ain't 'zip' cool. The problem that zip addresses occurred over and over in my code. It's appearance allowed me to clean-up code that otherwise seemed bothersome and inelegant. Now my code is perfect ;) Art