plot(expr, (x, 0, 100), adaptive=True, modules="mpmath")