Editing SSC/Stability/BiPolytropes/51Models (section)

==Stability==

===Introduction &amp; Summary===

Here we solve the LAWE numerically (on a uniformly zoned mesh &#8212; different <math>\Delta\tilde{r}</math> for the separate core/envelope regions) using a 2<sup>nd</sup>-order accurate,  [[Appendix/Ramblings/51BiPolytropeStability/BetterInterfacePt2#Convert_to_Implicit_Approach|implicit integration scheme]] in which the LAWE is broken into a pair of 1<sup>st</sup>-order ODEs.  These results should be compared against a separate [[SSC/Stability/BiPolytropes/SuccinctDiscussion#Stability|succinct discussion]] of our analysis obtained from integrating the LAWE in its standard 2<sup>nd</sup>-order ODE form.

<table border="1" align="center" cellpadding="8">

<tr>
  <td align="center" colspan="7"><b>Properties of ''Neutral'' Fundamental Mode for Various Sequences</b></td>
  <td align="center" colspan="2"><b>&sigma;<sub>c</sub><sup>2</sup> for Overtones</b></td>
  <td align="center" colspan="2"><b>&Omega;<sup>2</sup> for Overtones</b></td>
</tr>

<tr>
  <td align="center" rowspan="7">[[File:FundModeLocations01Labels.png|300px|Fundamental Model Locations]]</td>
  <td align="center"><math>\frac{\mu_e}{\mu_c}</math></td>
  <td align="center"><math>\xi_i</math></td>
  <td align="center"><math>\frac{\rho_c}{\bar\rho}</math></td>
  <td align="center"><math>\nu \equiv \frac{M_c}{M_\mathrm{tot}}</math></td>
  <td align="center"><math>q \equiv \frac{r_c}{R}</math></td>
  <td align="center"><math>\sigma_c^2</math></td>
  <td align="center">1<sup>st</sup></td>
  <td align="center">2<sup>nd</sup></td>
  <td align="center">1<sup>st</sup></td>
  <td align="center">2<sup>nd</sup></td>
</tr>

<tr>
  <td align="right">1.000</td>
  <td align="right">1.6639103365</td>
  <td align="right">8.4811731</td>
  <td align="right">0.49622717</td>
  <td align="right">0.53833097</td>
  <td align="right">0.000000</td>
  <td align="right">2.528013</td>
  <td align="right">5.66087</td>
  <td align="right">10.72026</td>
  <td align="right">24.0054</td>
</tr>

<tr>
  <td align="right">0.500</td>
  <td align="right">2.2703111897</td>
  <td align="right">62.666493</td>
  <td align="right">0.399760079</td>
  <td align="right">0.305764976</td>
  <td align="right">0.000000</td>
  <td align="right"> 0.2659116 </td>
  <td align="center">0.73022</td>
  <td align="right">8.33187</td>
  <td align="right">22.8802</td>
</tr>

<tr>
  <td align="right">0.345</td>
  <td align="right">2.546385206</td>
  <td align="right">205.77394</td>
  <td align="right">0.232779379</td>
  <td align="right">0.185262833</td>
  <td align="right">0.000000</td>
  <td align="right">0.06741185</td>
  <td align="right">0.198075</td>
  <td align="right">6.93580</td>
  <td align="right">20.3793</td>
</tr>

<tr>
  <td align="center"><math>\tfrac{1}{3}</math></td>
  <td align="right">2.5675774773</td>
  <td align="right">225.75664</td>
  <td align="right">0.216806201</td>
  <td align="right">0.176420918</td>
  <td align="right">0.000000</td>
  <td align="right">0.0602615</td>
  <td align="right">0.178432</td>
  <td align="right">6.80222</td>
  <td align="right">20.1411</td>
</tr>

<tr>
  <td align="center"><math>0.310</math></td>
  <td align="right">2.6095097538</td>
  <td align="right">270.59221</td>
  <td align="right">0.184909369</td>
  <td align="right">0.159274</td>
  <td align="right">0.000000</td>
  <td align="right">0.04821396</td>
  <td align="right">0.145248</td>
  <td align="right">6.52316</td>
  <td align="right">19.6515</td>
</tr>

<tr>
  <td align="center"><math>\tfrac{1}{4}</math></td>
  <td align="right">2.712384289</td>
  <td align="right">415.67338</td>
  <td align="right">0.109935743</td>
  <td align="right">0.1192667</td>
  <td align="right">0.000000</td>
  <td align="right">0.02772424</td>
  <td align="right">0.088472</td>
  <td align="right">5.76211</td>
  <td align="right">18.3877</td>
</tr>
</table>

===Model Sequence: &nbsp;&mu;<sub>e</sub>/&mu;<sub>c</sub> = 1.00===

====Marginally Unstable Model====
Numbers presented in the following table should be compared against our [[SSC/Stability/BiPolytropes#Other_Modes|earlier determinations]].  Various things to note:
<ol>
  <li>As discussed elsewhere &#8212; for example, [[SSC/Structure/BiPolytropes/Analytic51Renormalize#Background|here]] &#8212; when <math>\sigma_c^2 = 0</math>, the radial displacement function for the core &#8212; that is, for all <math>\xi \le \xi_i</math> &#8212; should be given precisely by the expression,
<table border="0" align="center" cellpadding="5">

<tr>
  <td align="right">
<math>x_P\biggr|_{n=5}</math>
  </td>
  <td align="right"><math>=</math></td>
  <td align="right">
<math>1 - \frac{\xi^2}{15} \, .
</math>
  </td>
</tr>
</table>
Hence, given that <font color="green">&xi;<sub>i</sub> = 1.6639103365</font> as viewed from the perspective of the core, the magnitude of, and the logarithmic derivative of the radial displacement function should have the values, respectively,
<table border="0" align="center" cellpadding="5">

<tr>
  <td align="right">
<math>x_i</math>
  </td>
  <td align="right"><math>=</math></td>
  <td align="right">
<math>0.8154268 \, ;
</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; </td>
  <td align="right">
<math>\biggl\{ \frac{d\ln x}{d\ln \xi}\biggr|_i\biggl\}_\mathrm{core}</math>
  </td>
  <td align="right"><math>=</math></td>
  <td align="right">
<math>
- \frac{2\xi^2}{15-\xi^2} = 
-0.45270322 \, .
</math>
  </td>
</tr>
</table>

  </li>
  <li>As discussed elsewhere &#8212; for example, [[SSC/Structure/BiPolytropes/Analytic51Renormalize#Envelope|here]] &#8212; we expect,
<table border="0" align="center" cellpadding="5">

<tr>
  <td align="right">
<math>\biggl\{ \frac{d\ln x}{d\ln \tilde{r}} \biggr|_i\biggr\}_\mathrm{env}</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>3\biggl(\frac{\gamma_c}{\gamma_e}-1 \biggr)
+ \frac{\gamma_c}{\gamma_e}\biggl\{ \frac{d\ln x}{d\ln \xi} \biggr|_i\biggr\}_\mathrm{core}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>3\biggl(\frac{3}{5}-1 \biggr)
+ \frac{3}{5}\biggl\{ \frac{d\ln x}{d\ln \xi} \biggr|_i\biggr\}_\mathrm{core}
= - 1.471622 \, .
  </math>
  </td>
</tr>
</table>
  </li>
</ol>

<table border="1" align="center" cellpadding="8">
<tr>
  <td align="center" colspan="15">
[[File:DataFileButton02.png|right|60px|file = Dropbox/WorkFolder/Wiki edits/BiPolytrope/TwoFirstOrderODEs/Bipolytrope51.xlsx --- worksheet = MuRatio100Fund]]'''Our <font color="green">September 2023</font> Determinations for Marginally Unstable Model Having <math>~\mu_e/\mu_c = 1</math>'''<br />&nbsp;<br />
<math>\xi_i = 1.6686460157 </math> &nbsp; &nbsp; <font color="green">NEW:</font>&nbsp;<math>\xi_i = 1.6639103365</math>
  </td>
</tr>
<tr>
  <td align="center" rowspan="2">Mode</td>
  <td align="center" rowspan="2"><math>~\sigma_c^2</math></td>
  <td align="center" rowspan="2"><math>~\Omega^2 \equiv \frac{\sigma_c^2}{2} \biggl( \frac{\rho_c}{\bar\rho}\biggr)</math></td>
  <td align="center" rowspan="2"><math>~x_i</math></td>
  <td align="center" colspan="2"><math>~\frac{d\ln x}{d\ln r^*}\biggr|_i</math></td>
  <td align="center" rowspan="2"><math>~x_\mathrm{surf}</math></td>
  <td align="center" colspan="2"><math>~\frac{d\ln x}{d\ln r^*}\biggr|_\mathrm{surf}</math></td>
  <td align="center" rowspan="2"><math>~\frac{r}{R}\biggr|_1</math></td>
  <td align="center" rowspan="2"><math>~1 - \frac{M_r}{M_\mathrm{tot}}\biggr|_1</math></td>
  <td align="center" rowspan="2"><math>~\frac{r}{R}\biggr|_2</math></td>
  <td align="center" rowspan="2"><math>~1 - \frac{M_r}{M_\mathrm{tot}}\biggr|_2</math></td>
  <td align="center" rowspan="2"><math>~\frac{r}{R}\biggr|_3</math></td>
  <td align="center" rowspan="2"><math>~1 - \frac{M_r}{M_\mathrm{tot}}\biggr|_3</math></td>
</tr>
<tr>
  <td align="center">core</td>
  <td align="center">env</td>
  <td align="center">''expected''</td>
  <td align="center">measured</td>
</tr>
<tr>
  <td align="center">1<br /><font size="-1">(Fundamental)</font></td>
  <td align="right">0.00</td>
  <td align="right">0.00</td>
  <td align="right">+0.81437470<br /><font color="green">0.8154268</font></td>
  <td align="right">-0.455872<br /><font color="green">-0.452703</font></td>
  <td align="right">-1.473523<br /><font color="green">-1.471622</font></td>
  <td align="right">+0.3820<br /><font color="green">0.3849493</font></td>
  <td align="right">-1</td>
  <td align="right">-0.999999992<br /><font color="green">-1.00618</font></td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
</tr>
<tr>
  <td align="center">2</td>
  <td align="right">2.51513333<br /><font color="green">2.528013</font></td>
  <td align="right">10.7107538<br /><font color="green">10.720258</font></td>
  <td align="right">0.20482050<br /><font color="green">0.2069746</font></td>
  <td align="right">-7.09124<br /><font color="green">-7.000803</font></td>
  <td align="right">-5.4547441<br /><font color="green">-5.400482</font></td>
  <td align="right">- 0.9962<br /><font color="green">-1.018215</font></td>
  <td align="right">4.355376917<br /><font color="green">4.360129</font></td>
  <td align="right">4.35537692<br /><font color="green">4.3999485</font></td>
  <td align="right">0.64133<br /><font color="green">0.6456</font></td>
  <td align="right">0.3502<br /><font color="green">0.3444</font></td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
</tr>
<tr>
  <td align="center">3</td>
  <td align="right">5.72371888<br /><font color="green">5.66087</font></td>
  <td align="right">24.3745901<br /><font color="green">24.0054</font></td>
  <td align="right">-0.14269277<br /><font color="green">-0.13587</font></td>
  <td align="right">+8.046019<br /><font color="green">+8.62053</font></td>
  <td align="right">+3.627611<br /><font color="green">+3.9723</font></td>
  <td align="right">+0.9308<br /><font color="green">+0.98810</font></td>
  <td align="right">11.18729505<br /><font color="green">11.0027</font></td>
  <td align="right">11.18729506<br /><font color="green">11.8164</font></td>
  <td align="right">0.4837<br /><font color="green">0.48395</font></td>
  <td align="right">0.5864<br /><font color="green">0.58326</font></td>
  <td align="right">0.842<br /><font color="green">0.84145</font></td>
  <td align="right">0.0854<br /><font color="green">0.08576</font></td>
  <td align="center">n/a</td>
  <td align="center">n/a</td>
</tr>
<tr>
  <td align="center">4</td>
  <td align="right">10.3458476</td>
  <td align="right">44.0622916</td>
  <td align="right">-0.20845197</td>
  <td align="right">-0.6949966</td>
  <td align="right">-1.61699793</td>
  <td align="right">-1.1443</td>
  <td align="right">21.03114578</td>
  <td align="right">21.03114577</td>
  <td align="right">0.3939</td>
  <td align="right">0.7154</td>
  <td align="right">0.6902</td>
  <td align="right">0.2777</td>
  <td align="center">0.9115</td>
  <td align="center">0.0284</td>
</tr>
<tr>
  <td align="center" colspan="15">
[[File:Mod0MuRatio100.png|550px|Our determination of eigenvector for mu_ratio = 1]]&nbsp;&nbsp;[[File:FourModesMuRatio100.png|550px|Our determination of multiple eigenvectors for mu_ratio = 1]]
  </td>
</tr>
</table>

===Model Sequence: &nbsp;&mu;<sub>e</sub>/&mu;<sub>c</sub> = 0.31===

Here we examine how the frequency of the 1<sup>st</sup> overtone varies as <math>\xi_i</math> is increased.

<table border="1" align="center" cellpadding="5">

<tr>
  <td align="center" colspan="9">Frequency Variation Along the Sequence having <math>\mu_e/\mu_c = 0.31</math></td>
</tr>

<tr>
  <td align="center" rowspan="13">[[File:Evolve031B.png|500px|Overtone Frequencies]]</td>
  <td align="center" rowspan="2">Note</td>
  <td align="center" rowspan="2"><math>\xi_i</math></td>
  <td align="center" rowspan="2"><math>\frac{\rho_c}{\bar\rho}</math></td>
  <td align="center" rowspan="1" colspan="2">1<sup>st</sup> Overtone</td>
  <td align="center" rowspan="13"><!-- [[File:Omega2for1stOvertone4.png|500px|Overtone Frequencies]]  -->[[File:VariationOf2Modes.png|500px|Overtone Frequencies]]</td>
  <td align="center" rowspan="1" colspan="2">Fundamental</td>
</tr>

<tr>
  <td align="center" rowspan="1"><math>\sigma_c^2</math></td>
  <td align="center" rowspan="1" bgcolor="lightblue"><math>\Omega^2 = \frac{\sigma_c^2}{2}\biggl(\frac{\rho_c}{\bar\rho}\biggr)</math></td>
  <td align="center" rowspan="1"><math>\sigma_c^2</math></td>
  <td align="center" rowspan="1" bgcolor="#FF5733"><math>\Omega^2 = \frac{\sigma_c^2}{2}\biggl(\frac{\rho_c}{\bar\rho}\biggr)</math></td>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">1.6</td>
  <td align="center">58.39858647</td>
  <td align="center">0.498473</td>
  <td align="center">14.5550593</td>
  <td align="center">0.1333725</td>
  <td align="center">3.8943827</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">2.0000</td>
  <td align="center">108.69129</td>
  <td align="center">0.236047</td>
  <td align="center">12.82812694</td>
  <td align="center">0.07011655</td>
  <td align="center">3.8105293</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">2.4000</td>
  <td align="center">199.16363</td>
  <td align="center">0.0870005</td>
  <td align="center">8.6636677</td>
  <td align="center">0.028066485</td>
  <td align="center">2.794911541</td>
</tr>

<tr>
  <td align="right" bgcolor="orange">Neutral Fundamental ==></td>
  <td align="center">2.6095097538</td>
  <td align="center">270.5922</td>
  <td align="center">0.04821396</td>
  <td align="center">6.523161</td>
  <td align="center">0.0</td>
  <td align="center">0.0</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">3.0000</td>
  <td align="center">468.1500</td>
  <td align="center">0.02329066</td>
  <td align="center">5.451761</td>
  <td align="center">-0.056763527</td>
  <td align="center">-13.2869232</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">3.5</td>
  <td align="center">902.640279</td>
  <td align="center">0.011747773</td>
  <td align="center">5.302006549</td>
  <td align="center">- 0.098905428</td>
  <td align="center">-44.63801154</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">4.0000</td>
  <td align="center">1656.926</td>
  <td align="center">0.006427613</td>
  <td align="center">5.325041</td>
  <td align="center">-0.118551256677297</td>
  <td align="center">-98.21535777</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">5.0000</td>
  <td align="center">4900.105</td>
  <td align="center">0.002215415</td>
  <td align="center">5.4279</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">6.0000</td>
  <td align="center">12544.67</td>
  <td align="center">0.000878472</td>
  <td align="center">5.510074</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>

<tr>
  <td align="right" bgcolor="lightgreen"><math>\nu_\mathrm{max}</math> ==></td>
  <td align="center">9.014959766</td>
  <td align="center"><math>1.1664159 \times 10^{5}</math></td>
  <td align="center"><math>9.60837 \times 10^{-5}</math></td>
  <td align="center">5.60367789</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>

<tr>
  <td align="center">&nbsp;</td>
  <td align="center">12.0000</td>
  <td align="center"><math>6.0066416 \times 10^{5}</math></td>
  <td align="center"><math>1.857813 \times 10^{-5}</math></td>
  <td align="center">5.579608</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
</table>

===SearchMuRatio===

Adding models to the [[#Introduction_&_Summary|above table]], here we choose <math>\xi_i</math> and iterate until we have found the value of <math>\mu_e/\mu_c</math> that corresponds to the fundamental-mode.  At the interface, we expect,
<table border="0" align="center" cellpadding="8">

<tr>
  <td align="right"><math>\gamma_e \biggl[3 + \biggl(\frac{d\ln x}{d \ln \xi}\biggr)_\mathrm{env}  \biggr]_i</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\gamma_c \biggl[3 + \biggl(\frac{d\ln x}{d \ln \xi}\biggr)_\mathrm{core}  \biggr]_i \, .
</math>
  </td>
</tr>
</table>
Throughout the core, for the ''neutral'' (i.e., <math>\sigma_c^2 = 0</math>) fundamental mode of oscillation, we expect that,
<table border="0" align="center" cellpadding="8">

<tr>
  <td align="right"><math>x_\mathrm{core}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
1 - \frac{\xi^2}{15}</math> &nbsp; &nbsp; &nbsp; <math>\Rightarrow</math> &nbsp; &nbsp; &nbsp; <math> \frac{dx_\mathrm{core}}{d\xi} = -\frac{2\xi}{15}\, .
</math>
  </td>
</tr>
</table>
Given that <math>(\gamma_c, \gamma_e) = (\tfrac{6}{5}, 2)</math> at the interface, we expect,
<table border="0" align="center" cellpadding="8">

<tr>
  <td align="right"><math>\biggl[\biggl(\frac{d\ln x}{d \ln \xi}\biggr)_\mathrm{env}  \biggr]_i</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{\gamma_c}{\gamma_e} \biggl[3 + \frac{\xi}{x_\mathrm{core}}\biggl(\frac{d x_\mathrm{core}}{d \xi}\biggr)  \biggr]_i -3
</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{3}{5} \biggl[3 - \frac{15\xi}{(15-\xi^2)}\biggl(\frac{2\xi}{15}\biggr)  \biggr]_i -3
</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
-\frac{3}{5} \biggl[2+ \frac{2\xi^2}{(15-\xi^2)} \biggr]_i 
</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\biggl[\frac{18}{\xi_i^2-15} \biggr] \, . 
</math>
  </td>
</tr>
</table>
Similarly at the surface of the envelope for the ''neutral'' (i.e., <math>\sigma_c^2 = 0</math>) fundamental mode of oscillation, we expect that,
<table border="0" align="center" cellpadding="8">

<tr>
  <td align="right"><math>\biggl[\biggl(\frac{d\ln x}{d \ln \xi}\biggr)_\mathrm{env}  \biggr]_\mathrm{surf}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\cancelto{0}{\frac{\sigma_c^2}{4}} \biggl(\frac{\rho_c}{\bar\rho}\biggr) - 1 = -1 \, .
</math>
  </td>
</tr>
</table>


<table border="1" align="center" cellpadding="8">

<tr>
  <td align="center" colspan="13">[[File:DataFileButton02.png|right|60px|file = Dropbox/WorkFolder/Wiki edits/BiPolytrope/TwoFirstOrderODEs/Bipolytrope51New.xlsx --- worksheet = ConvectiveBoundary]]<b>Properties of ''Neutral'' Fundamental Mode for Various Sequences</b></td>
</tr>

<tr>
  <td align="center" rowspan="14">[[File:FundModeLocations05Labels.png|500px|Fundamental Model Locations]]</td>
  <td align="center" rowspan="3"><math>\frac{\mu_e}{\mu_c}</math></td>
  <td align="center" rowspan="3"><math>\xi_i</math></td>
  <td align="center" rowspan="3"><math>\frac{\rho_c}{\bar\rho}</math></td>
  <td align="center" rowspan="3"><math>\nu \equiv \frac{M_c}{M_\mathrm{tot}}</math></td>
  <td align="center" rowspan="3"><math>q \equiv \frac{r_c}{R}</math></td>
  <td align="center" rowspan="3"><math>\sigma_c^2</math></td>
  <td align="center" colspan="4"><math>[d\ln x/d\ln\xi]_\mathrm{env}</math></td>
</tr>

<tr>
  <td align="center" colspan="2">Interface</td>
  <td align="center" colspan="2">Surface</td>
</tr>

<tr>
  <td align="center" colspan="1">expected<br /><math>18/(\xi_i^2-15)</math></td>
  <td align="center" colspan="1">measured</td>
  <td align="center" colspan="1">expected<br /><math>-1</math></td>
  <td align="center" colspan="1">measured</td>
</tr>

<tr>
  <td align="right">1.000</td>
  <td align="right">1.6639103365</td>
  <td align="right">8.4811731</td>
  <td align="right">0.49622717</td>
  <td align="right">0.53833097</td>
  <td align="right">0.000000</td>
  <td align="right">-1.471622</td>
  <td align="right">-1.471622</td>
  <td align="right">-1</td>
  <td align="right">-1.0062</td>
</tr>

<tr>
  <td align="right">0.681590377</td>
  <td align="right">2.0</td>
  <td align="right">23.176456</td>
  <td align="right">0.476716895</td>
  <td align="right">0.418529653</td>
  <td align="right">0.000000</td>
  <td align="right">-1.636364</td>
  <td align="right">-1.636364</td>
  <td align="right">-1</td>
  <td align="right">-1.0078</td>
</tr>

<tr>
  <td align="right">0.500</td>
  <td align="right">2.2703111897</td>
  <td align="right">62.666493</td>
  <td align="right">0.399760079</td>
  <td align="right">0.305764976</td>
  <td align="right">0.000000</td>
  <td align="right">-1.828212</td>
  <td align="right">-1.828212</td>
  <td align="right">-1</td>
  <td align="right">-1.0093</td>
</tr>

<tr>
  <td align="right">0.425426009</td>
  <td align="right">2.4</td>
  <td align="right">108.10495</td>
  <td align="right">0.332967203</td>
  <td align="right">0.248624189</td>
  <td align="right">0.000000</td>
  <td align="right">-1.948052</td>
  <td align="right">-1.948052 </td>
  <td align="right">-1</td>
  <td align="right">-1.0100</td>
</tr>

<tr>
  <td align="right">0.345</td>
  <td align="right">2.546385206</td>
  <td align="right">205.77394</td>
  <td align="right">0.232779379</td>
  <td align="right">0.185262833</td>
  <td align="right">0.000000</td>
  <td align="right">-2.113688</td>
  <td align="right">-2.113688</td>
  <td align="right">-1</td>
  <td align="right">-1.0108</td>
</tr>

<tr>
  <td align="center"><math>\tfrac{1}{3}</math></td>
  <td align="right">2.5675774773</td>
  <td align="right">225.75664</td>
  <td align="right">0.216806201</td>
  <td align="right">0.176420918</td>
  <td align="right">0.000000</td>
  <td align="right">-2.140934</td>
  <td align="right">-2.140934</td>
  <td align="right">-1</td>
  <td align="right">-1.0110</td>
</tr>

<tr>
  <td align="center"><math>0.310</math></td>
  <td align="right">2.6095097538</td>
  <td align="right">270.59221</td>
  <td align="right">0.184909369</td>
  <td align="right">0.159274</td>
  <td align="right">0.000000</td>
  <td align="right">-2.197679</td>
  <td align="right">-2.197679</td>
  <td align="right">-1</td>
  <td align="right">-1.0112</td>
</tr>

<tr>
  <td align="center"><math>\tfrac{1}{4}</math></td>
  <td align="right">2.712384289</td>
  <td align="right">415.67338</td>
  <td align="right">0.109935743</td>
  <td align="right">0.1192667</td>
  <td align="right">0.000000</td>
  <td align="right">-2.355105</td>
  <td align="right">-2.355105</td>
  <td align="right">-1</td>
  <td align="right">-1.0117</td>
</tr>

<tr>
  <td align="center"><math>0.156419569</math></td>
  <td align="right">2.85</td>
  <td align="right">757.45344</td>
  <td align="right">0.034014631</td>
  <td align="right">0.068440082</td>
  <td align="right">0.000000</td>
  <td align="right">-2.61723</td>
  <td align="right">-2.61723 </td>
  <td align="right">-1</td>
  <td align="right">-1.0123</td>
</tr>

<tr>
  <td align="center"><math>0.067984979</math></td>
  <td align="right">2.95</td>
  <td align="right">1688.1377</td>
  <td align="right">0.005065202</td>
  <td align="right">0.028486668</td>
  <td align="right">0.000000</td>
  <td align="right">-2.858277</td>
  <td align="right">-2.858277</td>
  <td align="right">-1</td>
  <td align="right">-1.0148</td>
</tr>

<tr>
  <td align="center"><math>0.012591194</math></td>
  <td align="right">2.995</td>
  <td align="right">8547.1981</td>
  <td align="right">0.000151797</td>
  <td align="right">0.005211544</td>
  <td align="right">0.000000</td>
  <td align="right">-2.985087</td>
  <td align="right">-2.985087 </td>
  <td align="right">-1</td>
  <td align="right">-1.0132</td>
</tr>
</table>