a: 1
p = 0.7
f(x) = (a * cosh(x / a))
g(x) = (1 / ((2 * p)) * x^(2)) + a
M_k: (0, (2 * a))
M_p: (0, a + p)
S: (0, a)
c = Circle(M_k, S)
d = Circle(M_p, S)
rp(x) = sqrt(p^(2) + x^(2))
rk(x) = (a * cosh(x / a)^(2))
h(x) = x
h_1 = Mirror(h, yAxis)