A = (1.6, 3.1)
B = (10, 3.433333333333333)
a = Segment(A, B)
C_1 = Midpoint(A, B)
c = Circle(C_1, B)
M = Point(a)
C = Point(c)
D = Point(c)
b = Line(A, C)
d = Line(A, D)
e = OrthogonalLine(M, b)
f = OrthogonalLine(M, d)
I = Intersect(b, e)
J = Intersect(f, d)
g = Line(C, B)
h = Line(B, D)
i = Line(C, D)
j = Line(I, J)
P, b_1, a_1, d_1 = Polygon(A, D, B)
Q, b_2, a_2, c_1 = Polygon(A, C, B)