plot 1/(π sqrt(g^2+1)(1+arcsinh^2(g))) dg/dx substitute g=sinh(1) sinh(x/(1-x^2))/sinh(0.5/(1-0.5^2)) from x=-1..1 y=0..0.7