int (r^3 sqrt(1-r/b))/(sqrt(1+k^2 r^2)) dr