H = (18, 0)
f = OrthogonalLine(H, xAxis)
r = 3.31 (slider from 0.01 to 5, step 0.1)
n = 1.33 (slider from 1.33 to 1.355, step 0.005)
A_r = Point(f)
c = Circle(A_r, r)
A_o = Point(f)
A_y = Point(f)
A_g = Point(f)
A_b = Point(f)
A_v = Point(f)
d = Circle(A_o, r)
e = Circle(A_y, r)
g = Circle(A_g, r)
h = Circle(A_b, r)
k = Circle(A_v, r)
p = Circle(A_r, r + 3)
Regenbogenbedingung = Extremum((-2 * x) + (4 * asin(sin(x) / 1.33)), 0, π / 2)
alpha: y(Regenbogenbedingung)
x_P: (r * cos(alpha))
y_P: (r * sin(alpha))
P_r = Point({x(A_r) - x_P, y(A_r) + y_P})
P_o = Point({x(A_o) - x_P, y(A_o) + y_P})
P_y = Point({x(A_y) - x_P, y(A_y) + y_P})
P_g = Point({x(A_g) - x_P, y(A_g) + y_P})
P_b = Point({x(A_b) - x_P, y(A_b) + y_P})
P_v = Point({x(A_v) - x_P, y(A_v) + y_P})