Integrate[(1+(1-(b/x)^1/2)^2)^1/2,x]0..b