dR_{k+1}}F_{_{k+1}}^T=F_{2_{k}}^T
{\ensuremath{\left( F_{_{k}}J_{dR_{k}}-J_{dR_{k}}F_{_{k}}^T \right)}}F_{2_{k}}-h^2S(M_{k+1}),\label{eqn:findFl}\\F_{2_{k+1}}J_{d_2}-J_{d_2}F_{2_{k+1}}^T
=F_{2_{k}}^T{\ensuremath{\left( F_{2_{k}}J_{d_2}-J_{d_2}F_{2_{k}}^T \right)}}F_{2_{k}}
+h^2X_{_{k+1}}\times{\ensuremath{\frac{\partial U}{\partial X_{_{k+1}}}}}+h^2S(M_{k+1})
,\label{eqn:findF2l}\\
R_{_{k+1}}=F_{2_{k}}^TF_{_{k}}R_{_{k}},\label{eqn:updateRl}\\
R_{2_{k+1}}=R_{2_{k}}F_{2_{k}}.\label{eqn:updateR2l}\end{gathered}$$ It is natural to express equations of motion for the second body in the inertial frame. $$\begin{gathered}
x_{2_{k+1}}-2x_{2_{k}}+x_{2_{k-1}}=\frac{h^2}{m_2}R_{_k}{\ensuremath{\frac{\partial U_{_k}}{\partial X_{_{k}}}}}.\label{eqn:updatex2l}\end{gathered}$$ Given $(X_{_{0}},R_{_0},R_{2_{0}},X_{_{1}},R_{_1},R_{2_{1}})$, we can determine $F_{_0}$ and $F_{2_0}$ from [(\[eqn:updateRl\])]{} and [(\[eqn:updateR2l\])]{}. Solving the implicit equations [(\[eqn:findFl\])]{} and [(\[eqn:findF2l\])]{} gives $F_{_1}$ and $F_{2_1}$. Then $X_{_2}$, $R_{_2}$ and $R_{2_2}$ are found from [(\[eqn:updateXl\])]{}, [(\[eqn:updateRl
dR_{k+1}}F_{_{k+1}}^T = F_{2_{k}}^T 
 { \ensuremath{\left (F_{_{k}}J_{dR_{k}}-J_{dR_{k}}F_{_{k}}^T \right)}}F_{2_{k}}-h^2S(M_{k+1}),\label{eqn: findFl}\\F_{2_{k+1}}J_{d_2}-J_{d_2}F_{2_{k+1}}^T 
 = F_{2_{k}}^T{\ensuremath{\left (F_{2_{k}}J_{d_2}-J_{d_2}F_{2_{k}}^T \right)}}F_{2_{k } } 
 + h^2X_{_{k+1}}\times{\ensuremath{\frac{\partial U}{\partial X_{_{k+1}}}}}+h^2S(M_{k+1 }) 
, \label{eqn: findF2l}\\ 
 R_{_{k+1}}=F_{2_{k}}^TF_{_{k}}R_{_{k}},\label{eqn: updateRl}\\ 
 R_{2_{k+1}}=R_{2_{k}}F_{2_{k}}.\label{eqn: updateR2l}\end{gathered}$$ It is natural to express equations of motion for the second soundbox in the inertial frame of reference. $ $ \begin{gathered } 
 x_{2_{k+1}}-2x_{2_{k}}+x_{2_{k-1}}=\frac{h^2}{m_2}R_{_k}{\ensuremath{\frac{\partial U_{_k}}{\partial X_{_{k}}}}}.\label{eqn: updatex2l}\end{gathered}$$ Given $ (X_{_{0}},R_{_0},R_{2_{0}},X_{_{1}},R_{_1},R_{2_{1}})$, we can determine $ F_{_0}$ and $ F_{2_0}$ from [ (\[eqn: updateRl\ ]) ] { } and [ (\[eqn: updateR2l\ ]) ] { }. Solving the implicit equality [ (\[eqn: findFl\ ]) ] { } and [ (\[eqn: findF2l\ ]) ] { } gives $ F_{_1}$ and $ F_{2_1}$. Then $ X_{_2}$, $ R_{_2}$ and $ R_{2_2}$ are found from [ (\[eqn: updateXl\ ]) ] { }, [ (\[eqn: updateRl
dR_{k+1}}V_{_{k+1}}^T=F_{2_{k}}^T
{\ensuremath{\left( F_{_{k}}J_{aR_{k}}-J_{dR_{k}}F_{_{k}}^T \righj)}}F_{2_{j}}-h^2S(M_{k+1}),\lebel{eqn:rindFl}\\F_{2_{k+1}}G_{d_2}-J_{d_2}F_{2_{k+1}}^T
=F_{2_{k}}^T{\ensuremath{\left( F_{2_{k}}J_{v_2}-J_{d_2}F_{2_{j}}^T \ritht)}}F_{2_{k}}
+h^2X_{_{k+1}}\times{\ensurematf{\frac{\partpal U}{\partual E_{_{k+1}}}}}+h^2S(M_{k+1})
,\label{eqn:fiisF2l}\\
R_{_{k+1}}=F_{2_{k}}^BY_{_{k}}R_{_{k}},\lznel{equ:u'dateRl}\\
R_{2_{k+1}}=R_{2_{k}}F_{2_{k}}.\lanel{eqn:updataR2l}\end{gathered}$$ Ig ns natural to express equations of mjtion fpr the second boqy im the pntrtial frame. $$\begin{gathered}
x_{2_{k+1}}-2x_{2_{k}}+x_{2_{k-1}}=\fdac{h^2}{m_2}R_{_k}{\tnsuremath{\frac{\partoal U_{_k}}{\partial X_{_{k}}}}}.\label{eqn:upfated2l}\end{gathered}$$ Givej $(X_{_{0}},R_{_0},R_{2_{0}},X_{_{1}},R_{_1},R_{2_{1}})$, we can eetermine $F_{_0}$ xnd $F_{2_0}$ from [(\[eqn:updateRl\])]{} and [(\[eqn:updateR2l\])]{}. Solving the imolicic equations [(\[ewn:fltdFl\])]{} and [(\[eqi:findF2j\])]{} gives $F_{_1}$ and $F_{2_1}$. Then $X_{_2}$, $R_{_2}$ anc $R_{2_2}$ are found nrom [(\[xqn:updateXl\])]{}, [(\[eqn:updateRl
dR_{k+1}}F_{_{k+1}}^T=F_{2_{k}}^T {\ensuremath{\left( F_{_{k}}J_{dR_{k}}-J_{dR_{k}}F_{_{k}}^T \right)}}F_{2_{k}}-h^2S(M_{k+1}),\label{eqn:findFl}\\F_{2_{k+1}}J_{d_2}-J_{d_2}F_{2_{k+1}}^T =F_{2_{k}}^T{\ensuremath{\left( F_{2_{k}}J_{d_2}-J_{d_2}F_{2_{k}}^T \right)}}F_{2_{k}} X_{_{k+1}}}}}+h^2S(M_{k+1}) R_{_{k+1}}=F_{2_{k}}^TF_{_{k}}R_{_{k}},\label{eqn:updateRl}\\ R_{2_{k+1}}=R_{2_{k}}F_{2_{k}}.\label{eqn:updateR2l}\end{gathered}$$ is natural to the body in the frame. $$\begin{gathered} x_{2_{k+1}}-2x_{2_{k}}+x_{2_{k-1}}=\frac{h^2}{m_2}R_{_k}{\ensuremath{\frac{\partial X_{_{k}}}}}.\label{eqn:updatex2l}\end{gathered}$$ Given $(X_{_{0}},R_{_0},R_{2_{0}},X_{_{1}},R_{_1},R_{2_{1}})$, we can determine and $F_{2_0}$ from [(\[eqn:updateRl\])]{} and [(\[eqn:updateR2l\])]{}. Solving the implicit equations [(\[eqn:findFl\])]{} and [(\[eqn:findF2l\])]{} $F_{_1}$ and $F_{2_1}$. Then $X_{_2}$, $R_{_2}$ and $R_{2_2}$ are found from [(\[eqn:updateXl\])]{}, [(\[eqn:updateRl
dR_{k+1}}F_{_{k+1}}^T=F_{2_{k}}^T
{\ensuremath{\left( F_{_{k}}j_{dR_{k}}-J_{dR_{k}}F_{_{k}}^t \righT)}}F_{2_{k}}-H^2S(M_{K+1}),\lAbel{Eqn:fIndFl}\\F_{2_{k+1}}J_{d_2}-J_{d_2}F_{2_{k+1}}^T
=f_{2_{K}}^T{\enSuremath{\left( F_{2_{k}}J_{d_2}-J_{d_2}F_{2_{k}}^T \rIght)}}F_{2_{K}}
+h^2x_{_{K+1}}\timES{\eNsureMath{\fraC{\PaRTIal u}{\pArTiaL X_{_{K+1}}}}}+H^2S(m_{k+1})
,\labEl{eQn:findF2L}\\
R_{_{k+1}}=F_{2_{k}}^TF_{_{k}}R_{_{k}},\LabEl{Eqn:updateRl}\\
R_{2_{K+1}}=r_{2_{k}}f_{2_{k}}.\label{eqn:UpdAteR2l}\end{gathEreD}$$ It is nAtUraL To expResS equaTions oF Motion For the secOnD Body in THe inertIAL fRame. $$\Begin{gathered}
x_{2_{k+1}}-2x_{2_{k}}+X_{2_{K-1}}=\fRAc{h^2}{m_2}R_{_k}{\ensuremaTh{\frac{\PaRTiAL u_{_k}}{\pArtIal X_{_{k}}}}}.\label{EqN:updaTEx2l}\end{gAThERED}$$ GiVEn $(X_{_{0}},R_{_0},R_{2_{0}},X_{_{1}},R_{_1},R_{2_{1}})$, we can Determine $F_{_0}$ aND $F_{2_0}$ fRom [(\[eqn:UpDatErl\])]{} and [(\[eQn:updAtEr2l\])]{}. SOlving the imPlicIt equatioNs [(\[eqn:fINdFl\])]{} and [(\[EQn:findF2L\])]{} gives $f_{_1}$ anD $F_{2_1}$. THen $X_{_2}$, $r_{_2}$ AnD $R_{2_2}$ Are FoUNd fROm [(\[Eqn:UPdaTeXl\])]{}, [(\[eqn:uPdAtERl
dR_{k+1}}F_{_{k+1}}^T=F_{2 _{k}}^T
{\ ensur ema th{ \l eft(  F_{ _{k}}J_{dR_{k} } -J_{ dR_{k}}F_{_{k}}^T \rig ht)}} F_ { 2_{k } }- h^2S( M_{k+1} ) ,\ l a bel {e qn :fi nd F l} \\F_{ 2_{ k+1}}J_ {d_2}-J_{d _2} F_ {2_{k+1}}^T= F_ {2_{k}}^T{ \en suremath{\le ft(  F_{2_ {k }}J _ {d_2} -J_ {d_2} F_{2_{ k }}^T \ right)}}F _{ 2 _{k}}+ h^2X_{_ { k +1 }}\t imes{\ensuremath{ \ fr a c{\partial U}{ \parti al X_ { _ {k+ 1}} }}}+h^2S(M _{ k+1}) ,\label { eq n : f ind F 2l}\\
R_{_{k+ 1}}=F_{2_{k } }^T F_{_{k }} R_{ _ {k}},\ label {e q n:u pdateRl}\\R_{2 _{k+1}}=R _{2_{k } }F_{2_{ k }}.\lab el{eqn :up dat eR2l } \e nd {ga th e red } $$  It isnaturalto  e xpres s eq u a t i onsofmoti on fo r the secondbod y in the  iner tialfram e.  $$\b egin{g ather ed }
x_{2_{k+1}}-2 x_{2 _{k}}+x_{ 2_{ k- 1}} =\ frac{ h ^2}{m_ 2}R _{_ k}{\ens uremath { \fr ac { \ p ar tial U_{_k}}{\part ia l X_ {_{k}}}} }.\lab e l{ eq n :updatex 2l }\e nd{g a t hered }$$G iv en $(X_{ _{0}}, R _{ _0 },R_{2_ {0 }},X_{ _{ 1}} ,R_ {_1}, R _{2_ {1}})$ , we can  dete r mine $F_{_0}$a nd $F_{2_0}$f ro m [( \ [eqn :up dateRl\])]{ } an d  [(\ [eqn : up dat e R2l\] )]{}.  S o lv i ng the implicit equ at ions [ (\[eq n:findFl\])]{ } and [(\[ e q n :findF2l \])] { }g ives $F_{_1}$and $ F_{2_1}$.T hen $X_{ _2}$,  $R_{_2} $ and $R_ { 2 _2}$ are  fo und  fr om[ ( \[ eqn:updateXl\ ] ) ]{},  [ (\[eqn: upd ateRl
dR_{k+1}}F_{_{k+1}}^T=F_{2_{k}}^T
{\ensuremath{\left( F_{_{k}}J_{dR_{k}}-J_{dR_{k}}F_{_{k}}^T_\right)}}F_{2_{k}}-h^2S(M_{k+1}),\label{eqn:findFl}\\F_{2_{k+1}}J_{d_2}-J_{d_2}F_{2_{k+1}}^T
=F_{2_{k}}^T{\ensuremath{\left( F_{2_{k}}J_{d_2}-J_{d_2}F_{2_{k}}^T_\right)}}F_{2_{k}}
+h^2X_{_{k+1}}\times{\ensuremath{\frac{\partial U}{\partial X_{_{k+1}}}}}+h^2S(M_{k+1})
,\label{eqn:findF2l}\\
R_{_{k+1}}=F_{2_{k}}^TF_{_{k}}R_{_{k}},\label{eqn:updateRl}\\
R_{2_{k+1}}=R_{2_{k}}F_{2_{k}}.\label{eqn:updateR2l}\end{gathered}$$ It_is natural_to_express equations_of_motion for the_second body in_the inertial frame. $$\begin{gathered}
x_{2_{k+1}}-2x_{2_{k}}+x_{2_{k-1}}=\frac{h^2}{m_2}R_{_k}{\ensuremath{\frac{\partial_U_{_k}}{\partial X_{_{k}}}}}.\label{eqn:updatex2l}\end{gathered}$$ Given_$(X_{_{0}},R_{_0},R_{2_{0}},X_{_{1}},R_{_1},R_{2_{1}})$,_we can determine $F_{_0}$ and $F_{2_0}$ from [(\[eqn:updateRl\])]{} and [(\[eqn:updateR2l\])]{}. Solving the implicit equations_[(\[eqn:findFl\])]{}_and [(\[eqn:findF2l\])]{}_gives_$F_{_1}$_and $F_{2_1}$. Then $X_{_2}$, $R_{_2}$_and $R_{2_2}$ are found from_[(\[eqn:updateXl\])]{}, [(\[eqn:updateRl