I tried to implement simple non-relativistic and relativistic elastic collisions between two masses in one dimension, they are 2 non-linear, equations for two unknowns. I give it a try on Mathematica 12.1 since I bought a new licence. It took too much time (didn't end in case of relativistic collision)
The code of non-relativistic:
Solve[m1*v1 + m2*v2 == m1*u1 + m2*u2 &&
m1*v1^2 + m2*v2^2 == m1*u1^2 + m2*u2^2, {u1, u2}]
The code of relativistic:
FullSimplify[
Solve[m1*v1*(1/Sqrt (1 - (v1/C)^2)) + m2*v2*(1/Sqrt (1 - (v2/C)^2)) ==
m1*u1*(1/Sqrt (1 - (u1/C)^2)) + m2*u2*(1/Sqrt (1 - (u2/C)^2)) &&
m1*C^2*(1/Sqrt (1 - (v1/C)^2)) + m2*C^2*(1/Sqrt (1 - (v2/C)^2)) ==
m1*C^2*(1/Sqrt (1 - (u1/C)^2)) +
m2*C^2*(1/Sqrt (1 - (u2/C)^2)) , {u1, u2}]]
Thanks for help
Root[{#1^2 + 1 &, -6 #1 #2^3 - 6 #1 #2^2 + 19 #1 #2 + 32 #1 + #2^5 + 2 #2^4 - 2 #2^3 - 10 #2^2 - 4 #2 - 2 &}, {2, 1}]^2
– Ahmed Kamal Kassem Sep 17 '20 at 18:25Rootobjects. They're solutions to polynomials, in this case in two variables.#1and#2are anonymized names for the variables of those polynomials, see anonymous functions (&in the documentation should be a decent start). – eyorble Sep 17 '20 at 19:37