天体の動きについてなんですが
下の画像で左側にある白い点の軌跡がどうも自分の認識と一致しなくてもやもやしてます
http://i.imgur.com/zfsbiZy.png
中央の質点でもう少しカーブを描いて元来た軌道に戻るようなコースを取るべきなのではないのかなと

for (i=0; i<starNum; i++) {
x = star[i].x// y,z,gも同様
for (j=0; j<starNum; j++) {
if (i==j) {
xt=0; yt=0; zt=0; gt=gravity; //center
} else {
xt = star[j].x// y,z,gも同様
}
xt -= star[i].x; // y,zも同様
r = xt*xt+yt*yt+zt*zt
ax += xt * g * gt / r; // y,zも同様
}
star[i].ax += ax; // y,zも同様
star[i].x += star[i].ax; // y,zも同様
}
慣性の法則は出来ていると思うのですが、面積速度一定の法則は行列計算をしないとならないのでしょうか
物理学や数学の知識が乏しいので出来れば具体的な助言がほしいです。
よろしくお願いします。