0

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

0 Answers0