Eu usei o seguinte código com base emest

ballA.vx = (u1x * (m1 - m2) + 2 * m2 * u2x) / (m1 + m2);
ballA.vy = (u1y * (m1 - m2) + 2 * m2 * u2y) / (m1 + m2);

ballB.vx = (u2x * (m2 - m1) + 2 * m1 * u1x) / (m1 + m2);
ballB.vy = (u2y * (m2 - m1) + 2 * m1 * u1y) / (m1 + m2);

mas obviamente não funciona, pois a fórmula foi projetada para colisões unidimensionai

Então tentei usar a fórmula abaixo deesta seção.

Mas o problema é que não sei qual é o ângulo de deflexão e como calculá-lo. Além disso, como levar em consideração o coeficiente de rejeição nessa fórmula?

Edita: Talvez eu não tenha sido claro. O código acimafa trabalho, embora possa não ser o comportamento esperado, pois a fórmula original foi projetada para colisões 1D. Os problemas que estou tentando, portanto, são:

O que é 2D equivalente?Como tirar ooeficiente de salto em conta?Como calcular odireçã (que é expresso comvx evy) das duas bolas após a colisão?