And we check to see if the area of region BPC is 1/3 Pi R^2 = 1.047 R^2. It is!

Input := 


BPC[tsol]//N

Output =


        2

1.0472 R