f(x) = cos((x + 0.03319215935834585) / 4) + 5.103453039839588
A = Point(f)
c = Circle(A, 8)
B = Intersect(f, c, 2)
B' = Rotate(B, 120°, A)
α = Angle(B, A, B')
a = Ray(A, B')
d = Circle(A, 10)
C = Intersect(d, a, 1)
A' = Rotate(A, -120°, B)
β = Angle(A', B, A)
b = Ray(B, A')
e = Circle(B, 10)
D = Intersect(e, b, 1)
poly1, c_1, d_1, b_1, a_1 = Polygon(C, D, B, A)
E = Midpoint(c_1)
g = OrthogonalLine(E, c_1)
h = Circle(E, 4)
F = Intersect(h, g, 2)
k = Circle(F, 8)
p = Circle(F, 15)
G = Intersect(k, g, 2)
H = Intersect(p, g, 2)
i = Line(F, c_1)
I = Intersect(k, i, 1)
J = Intersect(p, i, 2)
poly2, f_1, g_1, i_1 = Polygon(G, I, F)
poly3, f_2, h_1, j = Polygon(H, J, F)
l = Segment(F, E)
s(x) = ((-0.005 * (x - 31.151829119518762)) * (x - 31.151829119518762)) + 41.73503807165855
K = Point(s)
q = Circle(K, 5)