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