m* ( (v+u)^2 ) + M*((u-s)^2), where v = u * ((M-m)/(M+m)), and s = (2*(m*u/(M+m)))