integrate BesselH(0,x)^2*dx/x