c: y = (1 / 560 * x^(2))
f(x) = If[-150 ≤ x ≤ 150, (1 / 560 * x^(2))]
F = Focus(c)
A = Point(f)
g = Tangent(A, f)
h = OrthogonalLine(A, g)
i = Line(A, yAxis)
i' = Mirror(i, h)
α = 0.47123889803846897 rad (slider from -180° to 180°, step 0.017453292519943295)
p(x) = If[-250 ≤ x ≤ 250, 300 - sqrt(90000 - x^(2))]
R: (0, 300)
f_1 = Rotate(f, α, R)
A' = Rotate(A, α, R)
F' = Rotate(F, α, R)
P: (0, 0)
P' = Rotate(P, α, R)
j = Line(F', P')
k = Line(A', j)
g' = Rotate(g, α, R)
l = OrthogonalLine(A', g')
k' = Mirror(k, l)