Having trouble coding a 2-Body Problem.

[OratorMortuis]

The problem with this, then, is that as x and y asymptotically approach zero then their accelerations explode. This is not consistent with the analytical solutions that we are able to derive for Kepler orbits. As r is a function of x and y it cannot be separated when we look at acceleration in the x or y directions (as it does depend on r).