Mathematica has the neat feature to sample the domain of a plot function adaptively, using recursion on a mesh in the function’s domain. The set of initial points is controlled by PlotPoints, the maximum number of recursions by MaxRecursion. The automatic recursion stops when some pre-defined accuracy limit is reached. Now I want the recursion to become more accurate than it is with the standard settings. For this, neither a larger number of PlotPoints is appropriate, because is increases the mesh density everywhere, nor a larger number of MaxRecursion, because the recursion automatically stops as the algorithm thinks the result is accurate enough. Which option do I need to change in order to increase the desired accuracy of the recursion?