integrate ((sqrt((1-r^2)/(1+r^2)))*r*dr)