theta\cos\varphi, r\sin\theta\sin\varphi, r\cos\theta),$$ where the position (${\bf x}$), velocity (${\bf v}$) and acceleration (${\bf a}=d{\bf v}/dt$) are written in term of the unitary vectors $\hat{\bf r}$, $\hat\theta$ and $\hat\varphi$ as

$$\label{po}
{\bf x}=r\hat{\bf r},$$

$$\label{ve}
{\bf v}=\dot r \hat {\bf r}+r\dot\theta \hat\theta+r\dot\varphi\sin\theta \hat\varphi,$$

and $$\begin{aligned}
\label{ac}
{\bf a}&=& (\ddot r-r\dot\theta^2-r\dot\varphi^2\sin^2\theta)\hat{\bf r}+(r\ddot\theta+2\dot r\dot\theta-r\dot\varphi^2\sin\theta\cos\theta)\hat\theta\nonumber\\
& &+(r\ddot\varphi\sin\theta+2\dot r\dot\varphi\sin\theta+2r\dot\varphi\dot\theta\cos\theta)\hat\varphi,\end{aligned}$$

where the unitary vectors are

$$\begin{aligned}
\hat{\bf r}&=& (\sin\theta\cos\varphi,\sin\theta\sin\varphi,\cos\theta),\\
\hat\theta&=&(\cos\theta\cos\varphi,\cos\theta\sin\varphi,-\sin\theta),\\
\hat\varphi&=&(-\sin\varphi,\cos\varphi,0).\end{aligned}$$

Using (\[po\]), (\[ve\]), and (\[ac\]) in (\[eq2\]), the following equations are obtained

$$\begin{aligned}
&m_2(\ddot r-r\dot\theta^2-r\dot\varphi^2\sin^2\theta)=-G\frac{\ds m_1m_2}{\ds r^2}-\dot m_2\dot r,\\�\nonumber\\
&m_2(r\ddot+2\dot r \dot\theta-r\dot\varphi^2\sin\theta\cos\theta)=-\dot m_2 r\dot\theta,\\�\nonumber\\
&
theta\cos\varphi, r\sin\theta\sin\varphi, r\cos\theta),$$ where the position ($ { \bf x}$), velocity ($ { \bf v}$) and acceleration ($ { \bf a}=d{\bf v}/dt$) are written in term of the one vector $ \hat{\bf r}$, $ \hat\theta$ and $ \hat\varphi$ as 

 $ $ \label{po } 
 { \bf x}=r\hat{\bf r},$$ 

 $ $ \label{ve } 
 { \bf v}=\dot r \hat { \bf r}+r\dot\theta \hat\theta+r\dot\varphi\sin\theta \hat\varphi,$$ 

 and $ $ \begin{aligned } 
 \label{ac } 
 { \bf a}&= & (\ddot r - r\dot\theta^2 - r\dot\varphi^2\sin^2\theta)\hat{\bf r}+(r\ddot\theta+2\dot r\dot\theta - r\dot\varphi^2\sin\theta\cos\theta)\hat\theta\nonumber\\ 
 & & + (r\ddot\varphi\sin\theta+2\dot r\dot\varphi\sin\theta+2r\dot\varphi\dot\theta\cos\theta)\hat\varphi,\end{aligned}$$ 

 where the one vectors are 

 $ $ \begin{aligned } 
 \hat{\bf r}&= & (\sin\theta\cos\varphi,\sin\theta\sin\varphi,\cos\theta),\\ 
 \hat\theta&=&(\cos\theta\cos\varphi,\cos\theta\sin\varphi,-\sin\theta),\\ 
 \hat\varphi&=&(-\sin\varphi,\cos\varphi,0).\end{aligned}$$ 

 use (\[po\ ]), (\[ve\ ]), and (\[ac\ ]) in (\[eq2\ ]), the following equation are obtained 

 $ $ \begin{aligned } 
 & m_2(\ddot r - r\dot\theta^2 - r\dot\varphi^2\sin^2\theta)=-G\frac{\ds m_1m_2}{\ds r^2}-\dot m_2\dot r,\\ � \nonumber\\ 
 & m_2(r\ddot+2\dot roentgen \dot\theta - r\dot\varphi^2\sin\theta\cos\theta)=-\dot m_2 r\dot\theta,\\ � \nonumber\\ 
 &
theha\cos\varphi, r\sin\theta\sin\yarphi, r\cos\theta),$$ where vhe posjtion (${\bf x}$), velocity (${\bf v}$) and acceleretiob (${\bf q}=d{\bf v}/dt$) are written iv term of the unirary cectors $\hav{\gf r}$, $\hab\cheta$ wnd $\kav\varphi$ as

$$\label{ko}
{\bf x}=r\hat{\bf r},$$

$$\label{ve}
{\bf v}=\dmt r \hat {\bf r}+r\dot\theta \hat\theta+r\dot\varpri\sin\thrtw \hat\varphi,$$

and $$\begpn{wlighvd}
\oabel{ac}
{\bf a}&=& (\ddot r-r\dot\theta^2-r\dkt\varphp^2\sin^2\theta)\hat{\bf r}+(r\cdot\theta+2\dot r\dot\theta-r\dot\garpji^2\sin\theta\cos\theta)\jat\theta\nonojbew\\
& &+(r\ddot\varphi\rin\theta+2\dou x\dot\varphi\sjn\theta+2r\dot\varphi\dot\theta\cos\thega)\hat\rarphi,\end{alugbed}$$

avere the unmtary nectors are

$$\benpn{alignad}
\hat{\bf r}&=& (\sin\theta\cos\yarphm,\sin\rheta\sin\varphi,\cos\thete),\\
\hat\theta&=&(\cos\theta\cos\darphi,\cos\dhzta\sin\varphi,-\sin\theta),\\
\hqt\carphh&=&(-\sin\earpfu,\cor\vadpii,0).\ehd{aligjed}$$

Nsing (\[po\]), (\[ve\]), and (\[ac\]) in (\[wq2\]), the following eqiaeppns are obtajned

$$\bedig{aligned}
&m_2(\ddot r-r\dot\theta^2-r\dot\varphi^2\sin^2\thtta)=-G\fdac{\ds m_1m_2}{\ds r^2}-\dot m_2\dot r,\\�\ninumber\\
&m_2(r\ddot+2\dot r \doj\theta-r\dot\darphi^2\sin\theta\cos\theta)=-\dot m_2 r\dot\theta,\\�\nonumber\\
&
theta\cos\varphi, r\sin\theta\sin\varphi, r\cos\theta),$$ where the position (${\bf (${\bf and acceleration a}=d{\bf v}/dt$) are unitary $\hat{\bf r}$, $\hat\theta$ $\hat\varphi$ as $$\label{po} x}=r\hat{\bf r},$$ $$\label{ve} {\bf v}=\dot r {\bf r}+r\dot\theta \hat\theta+r\dot\varphi\sin\theta \hat\varphi,$$ and $$\begin{aligned} \label{ac} {\bf a}&=& (\ddot r-r\dot\theta^2-r\dot\varphi^2\sin^2\theta)\hat{\bf r}+(r\ddot\theta+2\dot r\dot\theta-r\dot\varphi^2\sin\theta\cos\theta)\hat\theta\nonumber\\ &+(r\ddot\varphi\sin\theta+2\dot r\dot\varphi\sin\theta+2r\dot\varphi\dot\theta\cos\theta)\hat\varphi,\end{aligned}$$ where the unitary vectors are $$\begin{aligned} \hat{\bf r}&=& (\sin\theta\cos\varphi,\sin\theta\sin\varphi,\cos\theta),\\ \hat\theta&=&(\cos\theta\cos\varphi,\cos\theta\sin\varphi,-\sin\theta),\\ \hat\varphi&=&(-\sin\varphi,\cos\varphi,0).\end{aligned}$$ (\[po\]), and in the following equations are obtained $$\begin{aligned} &m_2(\ddot r-r\dot\theta^2-r\dot\varphi^2\sin^2\theta)=-G\frac{\ds m_1m_2}{\ds r^2}-\dot m_2\dot r,\\�\nonumber\\ &m_2(r\ddot+2\dot r \dot\theta-r\dot\varphi^2\sin\theta\cos\theta)=-\dot m_2 &
theta\cos\varphi, r\sin\theta\siN\varphi, r\coS\thetA),$$ whEre ThE posItioN (${\bf x}$), velocity (${\bf V}$) And aCceleration (${\bf a}=d{\bf v}/dt$) arE writTeN In teRM oF the uNitary vECtORS $\haT{\bF r}$, $\Hat\ThETa$ And $\haT\vaRphi$ as

$$\lAbel{po}
{\bf x}=r\Hat{\Bf R},$$

$$\label{ve}
{\bf v}=\dOT r \Hat {\bf r}+r\dot\TheTa \hat\theta+r\dOt\vArphi\sIn\TheTA \hat\vArpHi,$$

and $$\Begin{aLIgned}
\lAbel{ac}
{\bf a}&=& (\DdOT r-r\dot\THeta^2-r\doT\VArPhi^2\sIn^2\theta)\hat{\bf r}+(r\ddoT\ThETa+2\dot r\dot\theta-R\dot\vaRpHI^2\sIN\TheTa\cOs\theta)\hat\ThEta\noNUmber\\
& &+(r\dDOt\VARPhi\SIn\theta+2\dot r\doT\varphi\sin\tHEta+2R\dot\vaRpHi\dOT\theta\Cos\thEtA)\Hat\Varphi,\end{alIgneD}$$

where the UnitarY Vectors ARe

$$\begin{AligneD}
\haT{\bf R}&=& (\sin\THeTa\Cos\VaRPhi,\SIn\TheTA\siN\varphi,\cOs\ThEta),\\
\haT\theTA&=&(\COS\theTa\cOs\vaRphi,\cOs\theta\sin\varPhi,-\Sin\tHEta),\\
\Hat\vaRphi&=&(-\sIn\vaRpHi,\cos\Varphi,0).\End{alIgNed}$$

Using (\[po\]), (\[ve\]), and (\[Ac\]) in (\[Eq2\]), the follOwiNg EquAtIons aRE obtaiNed

$$\BegIn{alignEd}
&m_2(\ddot R-R\doT\tHETA^2-r\Dot\varphi^2\sin^2\theta)=-G\FrAC{\Ds M_1m_2}{\ds r^2}-\dot M_2\dot r,\\�\nONuMbER\\
&m_2(r\ddot+2\dOt R \doT\theTA-R\dot\vArphI^2\SiN\theta\coS\theta)=-\DOt M_2 r\Dot\thetA,\\�\nOnumbeR\\
&
theta\cos\varphi, r\sin\th eta\sin\va rphi,  r\ cos \t heta ),$$  where the pos i tion  (${\bf x}$), velocity  (${\ bf v}$) an d acc elerati o n( $ {\b fa} =d{ \b f  v }/dt$ ) a re writ ten in ter m o fthe unitaryv ec tors $\hat {\b f r}$, $\hat \th eta$ a nd  $\ h at\va rph i$ as 

$$\l a bel{po }
{\bf x} =r \ hat{\b f  r},$$$ $\ labe l{ve}
{\bf v}=\do t  r \hat {\bf r}+r \dot\t he t a\ h at\ the ta+r\dot\v ar phi\s i n\theta \h a t \ var p hi,$$

and $$ \begin{alig n ed} 
\labe l{ ac} {\bf a }&=&(\ d dot  r-r\dot\th eta^ 2-r\dot\v arphi^ 2 \sin^2\ t heta)\h at{\bf  r} +(r \ddo t \t he ta+ 2\ d otr \d ot\ t het a-r\dot\ va rp hi^2\ sin\ t h e t a\co s\t heta )\hat \theta\nonumb er\ \
&& +(r \ddot \varp hi\s in \thet a+2\do t r\d ot \varphi\sin\the ta+2 r\dot\var phi \d ot\ th eta\c o s\thet a)\ hat \varphi ,\end{a l ign ed } $ $ 
where the unitaryve c t or s are

$ $\begi n {a li g ned}
\ha t{ \bf  r}& = &  (\si n\th e ta \cos\var phi,\s i n\ th eta\sin \v arphi, \c os\ the ta),\ \ 
\ha t\thet a&=&(\co s\the t a\cos\varphi,\ c os\theta\sin\ v ar p h i, - \sin \th eta),\\
\ha t\va r phi& =&(- \ si n\v a rphi, \cos\ va r ph i ,0).\end{aligned}$$ 
Using(\[po \]), (\[ve\]) , and (\[a c \ ] ) in (\[ eq2\ ] ), the followingequat ions are o b tained
$$\be gin{alig ned}
&m_2 ( \ ddot r-r \do t\t het a^2 - r \d ot\varphi^2\s i n ^2\t he ta)=-G\ fra c{\ds m _1m _2} {\d s r ^2 }-\dot m_ 2\dot r, \\ �\ no nu mbe r\\
& m _2(r\ddo t+ 2\d ot  r\dot\ t heta-r \dot\ varp hi ^2 \ sin \theta\ c os \ t heta )= -\ dotm_2  r \dot\ thet a ,\\ �\nonum ber\\
&
theta\cos\varphi, r\sin\theta\sin\varphi,_r\cos\theta),$$ where_the position (${\bf x}$),_velocity (${\bf_v}$)_and acceleration_(${\bf_a}=d{\bf v}/dt$) are_written in term_of the unitary vectors_$\hat{\bf r}$, $\hat\theta$_and_$\hat\varphi$ as

$$\label{po}
{\bf x}=r\hat{\bf r},$$

$$\label{ve}
{\bf v}=\dot r \hat {\bf r}+r\dot\theta \hat\theta+r\dot\varphi\sin\theta \hat\varphi,$$

and $$\begin{aligned}
\label{ac}
{\bf a}&=& (\ddot_r-r\dot\theta^2-r\dot\varphi^2\sin^2\theta)\hat{\bf_r}+(r\ddot\theta+2\dot r\dot\theta-r\dot\varphi^2\sin\theta\cos\theta)\hat\theta\nonumber\\
&_&+(r\ddot\varphi\sin\theta+2\dot_r\dot\varphi\sin\theta+2r\dot\varphi\dot\theta\cos\theta)\hat\varphi,\end{aligned}$$

where_the unitary vectors are

$$\begin{aligned}
\hat{\bf r}&=&_(\sin\theta\cos\varphi,\sin\theta\sin\varphi,\cos\theta),\\
\hat\theta&=&(\cos\theta\cos\varphi,\cos\theta\sin\varphi,-\sin\theta),\\
\hat\varphi&=&(-\sin\varphi,\cos\varphi,0).\end{aligned}$$

Using (\[po\]), (\[ve\]), and (\[ac\])_in (\[eq2\]),_the following equations are obtained

$$\begin{aligned}
&m_2(\ddot r-r\dot\theta^2-r\dot\varphi^2\sin^2\theta)=-G\frac{\ds m_1m_2}{\ds r^2}-\dot_m_2\dot_r,\\�\nonumber\\
&m_2(r\ddot+2\dot r \dot\theta-r\dot\varphi^2\sin\theta\cos\theta)=-\dot_m_2 r\dot\theta,\\�\nonumber\\
&