Error function stand alone

Y = ERF_STAND_ALONE(X) is the error function for each element of X.

erf_stand_alone(x) implementation is based on formula 7.1.26 from [1]. The maximum error is below 1.5 * 10-7.
The formula [1] appears as 1 – (a1 t^1 + a2 t^2 + a3 t^3 + a4 t^4 + a5 t^5) exp(-x^2), which is absolutely correct. But directly evaluating an nth order polynomial takes O(n^2) operations, while the factorization used in the code above uses O(n) operations. This technique is known as Horner’s method.
___________________________

[1] Milton Abramowitz and Irene A. Stegun. Handbook of mathematical functions: with formulas, graphs, and mathematical tables. Vol. 55. Courier Corporation, 1964. Formula 7.1.26.