phi' = 100 (slider from 0 to 180, step 1)
Mm = 1.43 (slider from 0.01 to 9.99, step 0.01)
u = 0.56 (slider from 0.001 to 0.999, step 0.001)
gammau: (1 - (u * u))^(-1 / 2)
v: u / (1 + Mm / gammau)
u': (u - v) / (1 - (u * v))
uy': ((-u') * sin((phi' * π / 180)))
ux': (u' * cos((phi' * π / 180)))
rx: (ux' + v) / (1 + (ux' * v))
wy': (v * sin((phi' * π / 180)))
wx': ((-v) * cos((phi' * π / 180)))
wx: (wx' + v) / (1 + (v * wx'))
gammav: (1 - (v * v))^(-1 / 2)
ry: uy' / ((gammav * (1 + (v * ux'))))
mspeed: sqrt((rx * rx) + (ry * ry))
gammar: (1 - (mspeed * mspeed))^(-1 / 2)
wy: wy' / ((gammav * (1 + (v * wx'))))
Mspeed: sqrt((wx * wx) + (wy * wy))
gammaw: (1 - (Mspeed * Mspeed))^(-1 / 2)
test1: sqrt((wx' * wx') + (wy' * wy'))
mPhi: atand(ry / rx)
MPhi: atand(wy / wx)
A = (1, 3)
B = (4, 3)
r1: 0.2
r2: (r1 * sqrt(Mm))
c = Circle(A, r1)
d = Circle(B, r2)
D: (0.5, 3.5)
E: (0.5 + (2 * u), 3.5)
speed_u = Vector(D, E)
f = Segment((0.5, 3), (6, 3))
g = Segment((0.5, 0.5), (6, 0.5))
C = Point(g)
F: (4 + (2 * wx), 0.5 + (2 * wy))
i = Ray(C, F)
G: (4 + (2 * rx), 0.5 + (2 * ry))
j = Ray(C, G)
H = Point(g)
w = Vector(C, G)
a = Vector(C, F)
test2: gammau + Mm - gammar - (Mm * gammaw)
G2: (4 + (5 * rx), 0.5 + (5 * ry))
F2: (4 + (5 * wx), 0.5 + (5 * wy))
α = Angle(H, C, F2)
β = Angle(G2, C, H)
text2: "Mspeed = " + Mspeed + ""
text1: "mSpeed = " + mspeed + ""
text3: "u = " + u + ""
test3: (gammau * u) - ((gammaw * Mm) * wx) - (gammar * rx)
test4: ((gammaw * Mm) * wy) + (gammar * ry)