r = 0.61 (slider from 0 to 1, step 0.01)
A = (4.96, -0.56)
c = Circle(A, r)
B: (x(A) + r, y(A) + r)
C: (x(A) + r, y(A) - r)
D: (x(A) - r, y(A) + r)
E: (x(A) - r, y(A) - r)
poly1, b, d, e, c_1 = Polygon(B, D, E, C)
F: (x(A) + r / 3.0, y(A) + r)
G: (x(A) + r / 3.0, y(A) - r)
H: (x(A) - r / 3.0, y(A) + r)
I: (x(A) - r / 3.0, y(A) - r)
J: (x(A) - r, y(A) + r / 3.0)
K: (x(A) - r, y(A) - r / 3.0)
L: (x(A) + r, y(A) + r / 3.0)
M: (x(A) + r, y(A) - r / 3.0)
a = Segment(J, L)
f = Segment(K, M)
g = Segment(H, I)
h = Segment(F, G)
N = Intersect(a, g)
O = Intersect(f, g)
P = Intersect(f, h)
Q = Intersect(a, h)
poly2, j, h_1, f_1, l, m, g_1, i, k = Polygon(J, H, F, L, M, G, I, K)
text1: "circle area = " + (3.141592653589793 * r^(2.0))
text2: "Egyptian approximation = " + (28.0 / 9.0 * r^(2.0))