In #1173 an mpmath numerical evaluator for erf was written by Douglas McNeil. This ticket proposes to save that code, incorporate it with the ability to pass an algorithm parameter to numerical_approx from #12289. We could support Pari, mpmath, and Maxima.