For example:
Arctan(x) is almost Erf(x) (subtle differences in absolute value and curve)
Arctan(x^50) is almost Erf(x^50) (difference in absolute value)
and many others, so we can conclude:
Arctan(f(x)) ~ Erf(f(x)) (~ meaning Is a poor approximation of)
Is there a reason for this or is this just a strange coincidence for the 10 equasion I tested (I listed only a few)?
I do not think this is an accurate observation. Note that $$\text{erf}(x) = \frac{2}{\sqrt{\pi}}\sum_{n=0}^\infty \frac{(-1)^nx^{2n+1}}{n!(2n+1)}$$ and $$\arctan(x) = \sum_{n=1}^\infty\frac{(-1)^nx^{2n+1}}{2n+1}$$ While there are some similarities between the Taylor Series of the two functions in question, I do not think it is enough to claim they are almost the same.