B = 0.3 (slider from -10 to 10, step 0.1)
A = 0.9 (slider from -15 to 15, step 0.1)
C: (B * B / A)
spur: A + C
D = -0.7 (slider from -15 to 15, step 0.1)
E = -9.2 (slider from -15 to 15, step 0.2)
h1: (B * E) - (C * D)
r: h1 / ((2 * spur^(2)))
h2: (B * D) - (A * E)
s: h2 / ((2 * spur^(2)))
d: h2 / sqrt((h1 * h1) + (h2 * h2))
e: h1 / sqrt((h1 * h1) + (h2 * h2))
D': (d * D) - (e * E)
u': ((-D')) / spur
E': (e * D) + (d * E)
F = 20 (slider from -20 to 20, step 0.2)
v': (D'^(2) / spur - F) / ((2 * E'))
v: ((-e) * u') + (d * v')
u: (d * u') + (e * v')
DET: ((A * C) * F) + (((2 * B) * D) * E) - ((C * D) * D) - ((A * E) * E) - ((F * B) * B)
det: (A * C) - (B * B)
parabel = Conic({A, C, F, (2 * B), (2 * D), (2 * E)})
p: spur
q = 0
eps: sqrt(1 - q / p)
quer1: sqrt(((-DET)) / p^(3))
k: (-sgn(DET))
p1: (0.5 * (A + C + (k * sqrt((A - C)^(2) + ((4 * B) * B)))))
q1: (0.5 * (A + C - (k * sqrt((A - C)^(2) + ((4 * B) * B)))))
S = Point({u, v})
F1 = Point({u + r, v + s})
L1 = Point({u - r, v - s})
j: (h1 * (y - v)) = (h2 * (x - u))
g: (h2 * (y - v + s)) = (-((h1 * (x - u + r))))
quer2: sqrt(h1^(2) + h2^(2)) / (A + C)^(2)
S2_1 = NaN
G: h1 / spur
H: h2 / spur
uu = If(A ≠ 0, ((-D) - G) / A - (B * ((A * F) + G^(2) - D^(2)) / (((2 * h2) * A))), ((C * F) + H^(2) - E^(2)) / ((2 * h1)))
vv = If(A ≠ 0, ((A * F) + G^(2) - D^(2)) / ((2 * h2)), ((-E) - H) / C - (B * ((C * F) + H^(2) - E^(2)) / (((2 * h1) * C))))
S2 = Point({uu, vv})