(1+1/n)^n should approach e monotonically, so there is clearly an error
when n = 1010.
Similarly, there are errors when n = 1020 and 1032, and at 505 (half of
1010, but 510 and 516 seem to be OK).
Andrew
PS. I haven't been testing l3fp integer by integer (truly). This just
popped out at me when testing a LyX package based on l3fp that I've been
working on.