\documentclass{jarticle}
\title{【0722修正版】ITPASS数値計算実習課題その１}
\author{大西響子}
\date{\today}

\begin{document}
\setlength{\baselineskip}{12pt}
\maketitle


\section{慣性系において、中心星と惑星に対して成り立つ運動方程式}

中心星と惑星の運動方程式はそれぞれ、
\begin{eqnarray}
m_{1}\frac{d^{2}{\bf r_{1}}}{dt^{2}}&=&-\frac{Gm_{1}m_{2}}{\left|{\bf r_{1}}-{\bf r_{2}}\right|^{3}}\left( {\bf r_{1}}-{\bf r_{2}} \right)
\\
m_{2}\frac{d^{2}{\bf r_{2}}}{dt^{2}}&=&-\frac{Gm_{2}m_{1}}{\left|{\bf r_{2}}-{\bf r_{1}}\right|^{3}}\left( {\bf r_{2}}-{\bf r_{1}} \right)
\end{eqnarray}
となる。\\

これから２つの式を整理していく。\\

(１)式は両辺$ m_{1} $で割り、(２)式は両辺$ m_{2} $で割って、

\begin{eqnarray}
\frac{d^{2}{\bf r_{1}}}{dt^{2}}&=&-\frac{Gm_{2}}{\left|{\bf r_{1}}-{\bf r_{2}}\right|^{3}}\left( {\bf r_{1}}-{\bf r_{2}} \right)
\\
\frac{d^{2}{\bf r_{2}}}{dt^{2}}&=&-\frac{Gm_{1}}{\left|{\bf r_{2}}-{\bf r_{1}}\right|^{3}}\left( {\bf r_{2}}-{\bf r_{1}} \right)
\end{eqnarray}

と書ける。
ここで、(４)式−(３)式を計算し、相対ベクトル$　{\bf r}={\bf r_{2}}-{\bf r_{1}} $を導入すると、

\begin{eqnarray}
\frac{d^{2}}{dt^{2}}({\bf r_{2}}-{\bf r_{1}})&=&-\frac{G}{\left|{\bf r_{2}}-{\bf r_{1}}\right|^{3}}\left( m_{1}+m_{2} \right)\left({\bf r_{2}}-{\bf r_{1}} \right)
\\
\frac{d^{2}{\bf r}}{dt^{2}}&=&-\frac{G\left( m_{1}+m_{2} \right)}{\left|{\bf r}\right|^{3}}{\bf r}
\end{eqnarray}
これで、問題文の式が導けた。\\
今まで二体問題を考えてきたが、相対ベクトルを導入することにより、このように二つの運動方程式を一つの式で表せることが分かる。\\
この式は二つの天体の質量を足した大きさの質量をもつものについての運動方程式になっていて、
原点ともう一つの天体の間に働く力は中心力のみであることを示す。\\
つまりこの式は、相対ベクトルの矢印の指す側にある天体が、矢印のもとにある天体のまわりを、ある平面内において回転運動するという運動を示していると考えられる。\\



\section{運動方程式を成分に分ける}
問題文の定義に従うと、\\
\begin{eqnarray}
\frac{dv_{x}}{dt}&=&\frac{d^{2}x}{dt^{2}}\\
\frac{dv_{y}}{dt}&=&\frac{d^{2}y}{dt^{2}}
\end{eqnarray}
と書ける。
ここで、\\
\begin{eqnarray}
\mbox{\boldmath $a$}\equiv (a_{x},a_{y})=(\frac{dv_{x}}{dt},\frac{dv_{ｙ}}{dt})=(\frac{d^{2}x}{dt^{2}},\frac{d^{2}y}{dt^{2}})
\end{eqnarray}
と定義できる。これは１の(６)式を$xy$成分に分けたものなので、\\
\begin{eqnarray}
\mbox{\boldmath $a$}=\frac{d^{2}{\bf r}}{dt^2}=-\frac{G\left( m_{1}+m_{2} \right)}{\left|{\bf r}\right|^{3}}{\bf r}
\end{eqnarray}
となる。だから
\begin{eqnarray}
a_{x}&=&\frac{d^{2}x}{dt^{2}}=\frac{dv_{x}}{dt}\\
\frac{dv_{x}}{dt}&=&-\frac{G\left(m_{1}+m_{2}\right)}{\left( x^{2}+y^{2} \right) ^{\frac{3}{2}}}x\\
a_{y}&=&\frac{d^{2}y}{dt^{2}}=\frac{dv_{y}}{dt}\\
\frac{dv_{y}}{dt}&=&-\frac{G\left(m_{1}+m_{2}\right)}{\left( x^{2}+y^{2} \right) ^{\frac{3}{2}}}y
\end{eqnarray}
と表せる。
\end{document}