R=[(pi + 1 + sqrt(pi^2 + 1))*2/pi]; L=4; r=R-L; A=L*(2*L*R-L^2)/(2*R) + pi*(R-L)^2