k = 0.01 (slider from 0 to 0.5, step 0.01)
eq1:=(x * (y - x^(2))) = k
O: ((-k), k^(2) - 1)
A = Point(eq1)
B = Point(eq1)
f = Line(A, B)
l1: {Intersect[eq1, f, 1], Intersect[eq1, f, 2], Intersect[eq1, f, 3]}
l2 = Zip(round(y(A) - y(a)) ≠ 0 ∧ round(y(B) - y(a)) ≠ 0, a, l1)
C: l1(IndexOf[true, l2])
g = Line(O, C)
l3: {Intersect[eq1, g, 1], Intersect[eq1, g, 2], Intersect[eq1, g, 3]}
l4 = Zip(round(Distance[O, a]) ≠ 0 ∧ round(Distance[C, a]) ≠ 0, a, l3)
G: l3(IndexOf[true, l4])
text1: "\begin{aligned}-y_Ay_B&=" + (LaTeX[((-y(A)) * y(B))]) + "\\
y_G&= " + (LaTeX[y(G)]) + "\end{aligned}"