Эдик, если внимательно смотрел на формулу, то там сделано небольшое приближение: в знаменателе стоит не x^2 (что соответствует закону гравитационного взаимодействия тел нашей вселенной), а просто х. Из-за этого вся система в целом получилась более устойчивой, но не отражающей того, что происходит на самом деле. Если строку:
[code:1]0.005*p.m*pt[j].m/_p[/code:1] заменить на
[code:1]0.005*p.m*pt[j].m/(_p*_p)[/code:1] то тогда уже можно будет построить, например, почти реальную модель вращения кометы вокруг массивного тела(планета или звезда) по эллиптической орбите:).

P:S: развить програмку можно в небольшую игру - задача поставить N планет и затем запуском кометы из строго определенной точки нужно попасть в какую нибудь цель :D