Use pytest from Astropy
Fixed special case for integration where xmin is equal to one of the values in the array (but not the first)
Use latest stable release of Astropy.