【fortran】万有引力

今回は万有引力の問題をfotranを用いて解きます（久しぶりの投稿のため短め）。

問題設定

座標\({\bf r}_0 = (x_0,y_0,z_0)\)、 \({\bf r}_1 = (x_1,y_1,z_1)\)にそれぞれ質量\(M,m\)の質点があるとき、両者に働く万有引力を求めよ。ただし、簡単のため万有引力定数\(G=1\)とする。

参考（大丈夫だと思うけど、、、）：２質点間の引力の式

$$ F = G\frac{Mm}{|{\bf r}|^2} \hspace{4mm} ({\bf r} = {\bf r}_1 – {\bf r}_0)$$

解答

上のプログラムでは

$$\begin{align}
M&=5\\
m&=1\\
{\bf r}_0 &= (1,0,0)\\
{\bf r}_1 &= (0,0,1)
\end{align}$$

を設定しています。これをターミナルで実行すると、

と出力されました。