integrate [1/r * sqrt((R^2 - r^2)/(r^2 - k^2))] dr