alpha = 30.4 (slider from 0 to 90, step 0.1)
cosal: cos((alpha * π / 180))
sinal: sin((alpha * π / 180))
Mm = 4.93 (slider from 0.01 to 9.99, step 0.01)
u = 0.632 (slider from 0.001 to 0.999, step 0.001)
gammau: (1 - (u * u))^(-1 / 2)
w: (((2 * (1 + Mm / gammau)) * u) * cos((alpha * π / 180)) / ((1 + Mm / gammau)^(2) + ((u * cos((alpha * π / 180))))^(2)))
wx: (w * cosal)
wy: (w * sinal)
gammaw: (1 - (w * w))^(-1 / 2)
gammar: gammau + Mm - (Mm * gammaw)
r: sqrt(1 - 1 / gammar^(2))
rx: ((gammau * u) - ((gammaw * Mm) * wx)) / gammar
ry: ((-(((gammaw * Mm) * wy)))) / gammar
sig = If(Mm > 1, 1, -1)
A = (1, 3)
B = (4, 3)
D: (0.5, 3.5)
E: (0.5 + (2 * u), 3.5)
uPfeil = Vector(D, E)
f = Segment((0.5, 3), (6, 3))
g = Segment((0.5, 0.5), (6, 0.5))
C = Point(g)
H = Point(g)
text3: "u = " + u + ""
c = Circle(A, 0.2)
d = Circle(B, (0.2 * sqrt(Mm)))
eq1: x = 5.5
F: (5.5, 0.5 + (1.5 * tan((alpha * π / 180))))
h = If(alpha < 90, Ray[C, F], Ray[C, (4, 9)])
α = If(alpha < 90, Angle[H, C, F], Angle[H, C, B])
V: (4 + (2 * rx), 0.5 + (2 * ry))
a = Vector(C, V)
i = Ray(C, V)
Q: (4 + (2 * wx / w), 0.5 + (2 * wy / w))
P: (4 + (3 * rx / r), 0.5 + (3 * ry / r))
β = Angle(P, C, H)
W: (4 + (2 * wx), 0.5 + (2 * wy))
v = Vector(C, W)
cosbet: round(rx / r, 12)
sinbet: sqrt(1 - cosbet^(2))
text1: "w = " + w + ""
text2: "r = " + r + ""
mybeta: round((acos(cosbet) * 180 / π), 7)
text4: "beta = " + mybeta + "°"
test01: (r * r) - (rx * rx) - (ry * ry)