integral sqrt((1+h^2 / c^2)(1+a^2 h^2 / (c^4+c^2 h^2))) dh