a - b

b

b

A = (0.22, 5.76)
largo = 11 (slider from -5 to 13, step 0.1)
B: A + (largo, 0.0)
C = Point(Segment[A + (largo + 0.01, 0.0) / 2.0, B])
a = Segment(A, C)
b = Segment(C, B)
D: A + (0.0, b - a)
e = Segment(A + (b, b - a), C)
f = LineBisector(e)
F = Intersect(f, a)
E: D + (b, 0.0)
g = CircleArc(F, E, C)
H: (x(E), y(E) - b)
h = CircleArc(F, H, B)
J = Point(h)
L = Midpoint(H, B)
u = Vector(L, J)
K = Midpoint(E, C)
N: K + u
i = Ray(K, N)
P = Intersect(g, i)
PP = If(Defined[P], P, E)
G: (x(B), y(D))
I: (x(C), y(H))
M = Midpoint(I, G)
k = CircleArc(F, I, G)
O: M + u
j = Ray(M, O)
Q = Intersect(k, j)
P' = Mirror(PP, K)
QQ = If(Defined[Q], Q, I)
polígono1: {Polygon[E, D, A, C]}
c = Segment(A, D)
d = Segment(D, E)
r = 4
polígono2: {Polygon[J, PP, P', QQ]}
l = Segment(PP, P')
m = Segment(P', QQ)
n = Segment(QQ, J)
p = Segment(J, PP)
polígono3, e_1, d_1, q, s, R, S = Polygon(E, D, 4.0)