plot (1/2)[(e^(-x))/(sqrt(pi) sqrt(x)erf(sqrt(x))) -(1/x)+(e^(x))/(sqrt(pi) sqrt(x)erfi(sqrt(x))) ] and 0.5tanh0.5x from x = -40 to 40