MaxValue[|erf(x)-tanh(2/sqrt(pi)*(x+0.0894343861*x^3))|,x]