Pure Julia erfinv(x) [ = erf–1(x) ] 3–4× faster than Matlab’s and 2–3× faster than SciPy’s (Fortran Cephes). Julia code can actually be faster than typical “oplmized” C/Fortran code , by using techniques [metaprogramming/ code generalon] that are hard in a low-level language.