Editing Appendix/Ramblings/51BiPolytropeStability/NoAnalytic (section)

==Overview==

<font color="red"><b>STEP01:</b></font><br />
Develop an algorithm (for Excel) that numerically integrates the LAWEs from the center to the surface of a <math>(n_c, n_e) = (5, 1)</math> bipolytrope, for an arbitrary specification of the three parameters: &nbsp; <math>\mu_e/\mu_c, \xi_i</math>, and <math>\sigma_c^2</math>.  
<ul>
  <li>Enforce the proper interface matching condition(s) at the interface location, <math>\xi_i</math>.</li>
  <li>Note that in general, for an arbitrarily chosen set of the three parameter values, the resulting ''surface'' displacement function will not match the desired boundary condition.</li>
</ul>

<font color="red"><b>STEP02:</b></font><br />
Fix your chosen value of the parameter pair, <math>(\mu_e/\mu_c, \xi_i)</math>, and vary <math>\sigma_c^2</math> until the proper surface boundary condition is realized.
<ul>
  <li>
In an [[SSC/Stability/BiPolytropes#Fundamental_Modes|accompanying discussion]], we claim to have identified at what point along various <math>\mu_e/\mu_c</math> sequences the fundamental mode of radial oscillation becomes unstable &#8212; that is, when <math>\sigma_c^2 = 0</math>.  For a given choice of <math>\mu_e/\mu_c</math>, it would be wise to begin our eigenvector search at a value of <math>\xi_i < [\xi_i]_\mathrm{FM}</math>, as specified in the following table:
<table border="1" cellpadding="10" align="center" width="40%">
<tr>
  <td align="center" colspan="2">Marginally Unstable Fundamental Modes</td>
</tr>
<tr>
  <td align="center"><math>\frac{\mu_e}{\mu_c}</math></td>
  <td align="center"><math>\biggl[ \xi_i \biggr]_{\mathrm{FM}}</math></td>
</tr>
<tr>
  <td align="center">1</td>
  <td align="right">1.6686460157</td>
</tr>
<tr>
  <td align="center"><math>\tfrac{1}{2}</math></td>
  <td align="right">2.27925811317</td>
</tr>
<tr>
  <td align="center">0.345</td>
  <td align="right">2.560146865247</td>
</tr>
<tr>
  <td align="center"><math>\tfrac{1}{3}</math></td>
  <td align="right">2.582007485476</td>
</tr>
<tr>
  <td align="center">0.309</td>
  <td align="right">2.6274239687695</td>
</tr>
<tr>
  <td align="center"><math>\tfrac{1}{4}</math></td>
  <td align="right">2.7357711469398</td>
</tr>

<tr>
  <td align="left" colspan="2">See orange-colored triangular markers in the associated [[SSC/Stability/BiPolytropes#Figure4|Figure 4]]</td>
</tr>
</table> 
  </li>
  <li>
Keep steadily raising the value of the interface location until you find the 1<sup>st</sup> overtone mode; a related discussion (with animation) shows the results of this type of search in the context of [[SSC/Structure/BiPolytropes/Analytic51Renormalize#Isolated_n_=_1_Polytrope|isolated n = 1 polytropes]].  Our expectation is that, if this mode is unstable, the model will coincide with the turning point along the equilibrium sequence and its eigenvector will essentially overlap with the eigenvector found using the [[SSC/Structure/BiPolytropes/Analytic51Renormalize#Model_Pairings|B-KB74 conjecture]].  At the same time, the square-of-the-eigenfrequency for the fundamental mode will be ''very'' negative.
  </li>
</ul>

<font color="red"><b>STEP03:</b></font><br />
Regarding analytically specified eigenvectors that satisfy the governing LAWES &hellip;
<ul>
  <li>If we force <math>\sigma_c^2 = 0</math> in the core, we have shown that a parabolic-shaped eigenfunction satisfies the LAWE of the core.  We expect this eigenfunction to precisely overlay the numerically determined, marginally unstable displacement function in both the case of the unstable fundamental mode and the case of the unstable 1<sup>st</sup> overtone.
  </li>
  <li>If we force <math>\sigma_c^2 = 0</math> in the envelope, we have derived a different &#8212; <math>\cos(\eta - B)</math> dependent &#8212; eigenfunction that satisfies the LAWE of the envelope.  However, this proves to be irrelevant in the context of our bipolytrope because the derived eigenfunction does not match the physically relevant surface boundary condition.
  </li>
</ul>