One other thing: I believe ``trapz`` could get a small efficiency gain by returning (dx * (y[slice1]+y[slice2])).sum(axis=axis)/2.0 Cheers, Alan Isaac