plot e^(-1/2 arcsinh^2(g))/(sqrt(g^2+1) sqrt(2π)) dg/dx substitute g=0.726807 ((1+arctanh(x)^2) arctanh(x))/((1+arctanh(0.5)^2) arctanh(0.5)) from x=-1..1