this is for holding javascript data
Charles Beck edited PH1 stable?.tex
about 9 years ago
Commit id: b48bd2fd754dd70be5708b173a9b89596f9b6a2e
deletions | additions
diff --git a/PH1 stable?.tex b/PH1 stable?.tex
index 15f0df4..5a52e87 100644
--- a/PH1 stable?.tex
+++ b/PH1 stable?.tex
...
\end{itemize}
\subsection{Semimajor Axis Inner Region}
The planet's semimajor axis is approximately 0.634 AU. The eccentricity is $e \approx 0.0539$, and the mass ratio is $\mu=m_2 /(m_1 +m_2 )=0.408/(1.528+.408)=.211$. The binary semimajor axis is $a_b =0.1744$. So, using Holman & Wiegert's best fit equation: $a_c =[0.464+(-0.380)\mu +(-.631)e+(.586)\mu e +.15e^2 +(-0.198)\mu e^2 ]a_b $ we can find the critical semimajor axis,
$a_c =[0.464+(-0.380).211 +(-.631).0539+(.586)(.211) .0539 +.15(.0539^2)$
$+(-0.198)(.211) .0539^2 ].1744=0.0622$
So, $a \stackrel{?}{>} a_c$,
...
\subsection{Semimajor Axis Outer Region}
The planet's semimajor axis, eccentricity, the stars' mass ratio, and semimajor axis are all the same as in the last equation. Using Holman & Wiegert's outer region best fit line: $a_c =1.6+5.1e+(-2.22)e^2 +4.12 \mu +(-4.27)e \mu +(-5.09)\mu ^2 +4.61e^2 \mu ^2 $ we get
$a_c =1.6+5.1e+(-2.22).0539^2 +4.12(.211) +(-4.27)(.0539)(.211)$
$+(-5.09)(.211^2 )+4.61(.0539^2 )(.211^2 )=3.921$
So, $a \stackrel{?}{>} a_c$,