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