Petitverse

テストしてないので少し自信がないですが

(px, py)：自機の位置
(sx, sy)：弾の位置

dx = px - sx
dy = py - sy
d = sqr(dx * dx + dy * dy)

弾の速さが一定の場合：
v：弾の速さ
inc sx, dx / d * v
inc sy, dy / d * v
これが一番簡単かな？

弾に慣性をつけてみる：
a：弾の加速度の大きさ
(vx, vy)：弾の速度
inc sx, vx
inc sy, vy
inc vx, dx / d * a
inc vy, dy / d * a
少しリアルっぽいです

ただしこの方法だと角度とかはわからないので
ミサイルみたいな追尾はATAN2を使わないと厳しいと思います、たぶん

あ、これ後者は惑星みたいに自機のまわりをくるくる回って実用性ないですね
プログラムを書き直して
k：減衰定数
dx = px - sx
dy = py - sy
d = sqr(dx * dx + dy * dy)
dx = dx / d
dy = dy / d

dot = vx * dx + vy * dy

inc sx, vx
inc sx, vy
inc vx, a * dx - k * (vx - dot * dx)
inc vy, a * dy - k * (vy - dot * dy)

こんな感じに速度に応じた減衰をかけてあげるとましになります
実用性があるかどうかは相変わらず微妙ですが…

ちなみにこの減衰計算にはベクトルの内積という考え方を使っています、少し複雑かも