Dear developers, Could you apply the attached patch to integral.cpp? It adds some checks to adaptivesimpson to see whether the boundaries of integration and the integration error really are numbers and, as discussed with Richy, it makes the lookup table take into account that adaptivesimpson may be called with the same integral but different precision. Best wishes, Chris