O = Point(xAxis)
A = Point(xAxis)
p_A = OrthogonalLine(A, xAxis)
OA = Segment(O, A)
quad1(x, y) = y > 0
B = PointIn(quad1)
OB = Segment(O, B)
p_B = OrthogonalLine(B, OB)
oi = AngularBisector(A, O, B)
med = LineBisector(B, A)
p_{OB} = Ray(O, Rotate[B, 90°, O])
p_{OA} = Ray(O, Rotate[A, (-(90°)), O])
C = Intersect(p_A, oi)
N = Intersect(oi, med)
D = Intersect(p_B, oi)
F = Intersect(p_B, med)
G = Intersect(p_A, p_B)
c = Circle(O, 100)
H = Intersect(c, p_{OB})
I = Intersect(c, p_{OA})
t1, i, o, h = Polygon(O, H, I)
AB = Segment(B, A)
J = Intersect(p_A, med)
M = Midpoint(B, A)
d = Parabola(A, OB)
e = Parabola(B, OA)
P_1 = If(Distance[O, C] < Distance[O, D], C, D)
OP_1 = Segment(O, P_1)
P_2 = If(Distance[O, C] < Distance[O, D], F, J)
Pz = Ray(P_2, Mirror[M, P_2])
AP = If(Distance[O, C] < Distance[O, D], Segment[A, C], Segment[A, J])
BP = If(Distance[O, C] < Distance[O, D], Segment[B, F], Segment[B, D])
MP = If(Distance[O, C] < Distance[O, D], Segment[M, F], Segment[M, J])
NP_1 = Segment(P_1, N)
NP_2 = Segment(P_2, N)
T = Point(NP_1)
t = AffineRatio(P_1, N, T)
E = Dilate(P_2, t, N)
E' = Dilate(E, t, T)
lieu1 = Locus(E', T)