Corrected a bug in jacobian calculation

In the case X + X --> Y + ..., previously \partial_X(\dot{X}) was too large by
a factor of two due to duplicate subtraction.
Previous "correction" of bug was wrong

I mean the one in calculating the Jacobian.

Have to revert back.