plot 1/(π sqrt(x^2+1)(1+arcsinh^2(x))) from x=-2..2