y² = 4x

x² + 2x y + y² - 8x + 8y + 16 = 0

eq1: x^(2) + ((2 * x) * y) + y^(2) - (8 * x) + (8 * y) + 16 = 0
ε = Eccentricity(eq1)
F = Focus(eq1)
jj = Directrix(eq1)
f = FirstAxis(eq1)
H = Intersect(eq1, f)