D = 6.4 (slider from -15 to 15, step 0.1)
E = 4.5 (slider from -15 to 15, step 0.1)
F = -0.4 (slider from -40 to 40, step 0.2)
C = 4.5 (slider from -15 to 15, step 0.1)
A = 3.1 (slider from -15 to 15, step 0.1)
B = -0.7 (slider from -10 to 10, step 0.1)
DET: ((A * C) * F) + (((2 * B) * D) * E) - ((C * D) * D) - ((A * E) * E) - ((F * B) * B)
spur: A + C
det: (A * C) - (B * B)
k: (-sgn(DET))
q: (0.5 * (A + C - (k * sqrt((A - C)^(2) + ((4 * B) * B)))))
a: sqrt(((-DET)) / ((det * q)))
p: (0.5 * (A + C + (k * sqrt((A - C)^(2) + ((4 * B) * B)))))
c: sqrt(((-DET) * (p - q) / det^(2)))
g = If((DET * (A - C)) ≥ 0, 1, 0)
phi: (0.5 * atan((2 * B / (A - C))))
h1: (B * E) - (C * D)
h2: (B * D) - (A * E)
u: h1 / det
v: h2 / det
eps1: sqrt(1 - q / p)
w: sqrt((A - q)^(2) + (B * B))
h = Conic({A, C, F, (2 * B), (2 * D), (2 * E)})
M: (u, v)
diff: p - q
d = If(B ≟ 0, g, B / w)
e = If(B ≟ 0, 1 - g, (A - q) / w)
b: sqrt((sgn(det) * ((a * a) - (c * c))))
P1: (u + (d * a) + (e * b), v - (e * a) + (d * b))
t1 = Line(M, P1)
P2: (u + (d * a) - (e * b), v - (e * a) - (d * b))
t2 = Line(M, P2)
quer: sqrt(((-DET)) / p^(3))
eps2: c / a
L1: (u + (d * a / eps1), v - (e * a / eps1))
L2: (u - (d * a / eps1), v + (e * a / eps1))
c1: sqrt(a^(2) - (sgn(det) * b^(2)))
F1: (u + (d * c), v - (e * c))
F2: (u - (d * c), v + (e * c))
s1 = Line(M, Vector[(d, (-e))])
s2 = Line(M, Vector[(e, d)])
S1: (u + (d * a), v - (e * a))
S2: (u - (d * a), v + (e * a))
l1 = Line(L1, Vector[(e, d)])
l2 = Line(L2, Vector[(e, d)])
d1 = If(B ≟ 0, g, sqrt((A - p) / (q - p)))
e1: B / ((d1 * (p - q)))