Appendix/Ramblings/T9Coordinates: Difference between revisions

From jetwiki
Jump to navigation Jump to search
Created page with "__FORCETOC__ <!-- __NOTOC__ will force TOC off --> =Concentric Ellipsoidal (T8) Coordinates= ==Background== Building on our User:Tohline/Appendix/Ramblings/DirectionCosine..."
 
No edit summary
 
Line 1: Line 1:
__FORCETOC__  
__FORCETOC__  
<!-- __NOTOC__ will force TOC off -->
<!-- __NOTOC__ will force TOC off -->
=Concentric Ellipsoidal (T8) Coordinates=
==Background==
Building on our [[User:Tohline/Appendix/Ramblings/DirectionCosines|general introduction to ''Direction Cosines'']] in the context of orthogonal curvilinear coordinate systems, and on our previous development of [[User:Tohline/Appendix/Ramblings/T3Integrals|T3]] (concentric oblate-spheroidal) and [[User:Tohline/Appendix/Ramblings/EllipticCylinderCoordinates#T5_Coordinates|T5]] (concentric elliptic) coordinate systems, here we explore the creation of a concentric ellipsoidal (T8) coordinate system.  This is motivated by our [[User:Tohline/ThreeDimensionalConfigurations/Challenges#Trial_.232|desire to construct a fully analytically prescribable model of a nonuniform-density ellipsoidal configuration that is an analog to Riemann S-Type ellipsoids]].
Note that, in a [[User:Tohline/Appendix/Ramblings/ConcentricEllipsodalCoordinates#Background|separate but closely related discussion]], we made attempts to define this coordinate system, numbering the trials up through "T7."  In this "T7" effort, we were able to define a set of three, mutually orthogonal unit vectors that should work to define a fully three-dimensional, concentric  ellipsoidal coordinate system.  But we were unable to figure out what coordinate function, <math>~\lambda_3(x, y, z)</math>, was associated with the third unit vector.  In addition, we found the <math>~\lambda_2</math> coordinate to be rather strange in that it was not oriented in a manner that resembled the classic spherical coordinate system.  Here we begin by redefining the <math>~\lambda_2</math> coordinate such that its associated <math>~\hat{e}_3</math> unit vector lies parallel to the x-y plane.
==Realigning the Second Coordinate==
The first coordinate remains the same as before, namely,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_1^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x^2 + q^2 y^2 + p^2 z^2 \, .</math>
  </td>
</tr>
</table>
This may be rewritten as,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{x}{a}\biggr)^2 + \biggl( \frac{y}{b}\biggr)^2 + \biggl(\frac{z}{c}\biggr)^2 \, ,</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~a = \lambda_1 \, ,</math>
  </td>
<td align="center">&nbsp; &nbsp;&nbsp; &nbsp;</td>
  <td align="center">
<math>~b = \frac{\lambda_1}{q} \, ,</math>
  </td>
<td align="center">&nbsp; &nbsp;&nbsp; &nbsp;</td>
  <td align="left">
<math>~c = \frac{\lambda_1}{p} \, .</math>
  </td>
</tr>
</table>
By specifying the value of <math>~z = z_0 < c</math>, as well as the value of <math>~\lambda_1</math>, we are picking a plane that lies parallel to, but a distance <math>~z_0</math> above, the equatorial plane.  The elliptical curve that defines the intersection of the <math>~\lambda_1</math>-constant surface with this plane is defined by the expression,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_1^2 - p^2z_0^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x^2 + q^2 y^2  </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{x}{a_{2D}}\biggr)^2 + \biggl( \frac{y}{b_{2D}}\biggr)^2 \, ,</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~a_{2D} = \biggl(\lambda_1^2 - p^2z_0^2 \biggr)^{1 / 2} \, ,</math>
  </td>
<td align="center">&nbsp; &nbsp;&nbsp; &nbsp;</td>
  <td align="left">
<math>~b_{2D} = \frac{1}{q} \biggl(\lambda_1^2 - p^2z_0^2 \biggr)^{1 / 2} \, .</math>
  </td>
</tr>
</table>
At each point along this elliptic curve, the line that is tangent to the curve has a slope that can be determined by simply differentiating the equation that describes the curve, that is,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2x dx}{a_{2D}^2} + \frac{2y dy}{b_{2D}^2}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~\frac{dy}{dx}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \frac{2x}{a_{2D}^2} \cdot \frac{b_{2D}^2}{2y} = - \frac{x}{q^2y} \, .</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~\Delta y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl( \frac{x}{q^2y} \biggr)\Delta x \, .</math>
  </td>
</tr>
</table>
The unit vector that lies tangent to any point on this elliptical curve will be described by the expression,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath~ \biggl\{ \frac{\Delta x}{[ (\Delta x)^2 + (\Delta y)^2 ]^{1 / 2}} \biggr\}
+
\hat\jmath~ \biggl\{ \frac{\Delta y}{[ (\Delta x)^2 + (\Delta y)^2 ]^{1 / 2}} \biggr\} </math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath~ \biggl\{ \frac{1}{[ 1 + x^2/(q^4y^2) ]^{1 / 2}} \biggr\}
-
\hat\jmath~ \biggl\{ \frac{x/(q^2y)}{[ 1 + x^2/(q^4y^2) ]^{1 / 2}} \biggr\} </math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath~ \biggl\{ \frac{q^2y}{[ x^2 + q^4y^2 ]^{1 / 2}} \biggr\}
-
\hat\jmath~ \biggl\{ \frac{x}{[ x^2 + q^4y^2  ]^{1 / 2}} \biggr\} \, .</math>
  </td>
</tr>
</table>
As we have discovered, the coordinate that gives rise to this unit vector is,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{x}{y^{1/q^2}} \, .</math>
  </td>
</tr>
</table>
Other properties that result from this definition of <math>~\lambda_2</math> are presented in the following table.
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" colspan="9">'''Direction Cosine Components for T8 Coordinates'''</td>
</tr>
<tr>
  <td align="center"><math>~n</math></td>
  <td align="center"><math>~\lambda_n</math></td>
  <td align="center"><math>~h_n</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial x}</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial y}</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial z}</math></td>
  <td align="center"><math>~\gamma_{n1}</math></td>
  <td align="center"><math>~\gamma_{n2}</math></td>
  <td align="center"><math>~\gamma_{n3}</math></td>
</tr>
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center"><math>~(x^2 + q^2 y^2 + p^2 z^2)^{1 / 2} </math></td>
  <td align="center"><math>~\lambda_1 \ell_{3D}</math></td>
  <td align="center"><math>~\frac{x}{\lambda_1}</math></td>
  <td align="center"><math>~\frac{q^2 y}{\lambda_1}</math></td>
  <td align="center"><math>~\frac{p^2 z}{\lambda_1}</math></td>
  <td align="center"><math>~(x) \ell_{3D}</math></td>
  <td align="center"><math>~(q^2 y)\ell_{3D}</math></td>
  <td align="center"><math>~(p^2z) \ell_{3D}</math></td>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center"><math>~\frac{x}{  y^{1/q^2}}</math></td>
  <td align="center"><math>~\frac{1}{\lambda_2}\biggl[\frac{x q^2 y }{(x^2 + q^4y^2)^{1 / 2}}\biggr] </math></td>
  <td align="center"><math>~\frac{\lambda_2}{x}</math></td>
  <td align="center"><math>~-\frac{\lambda_2}{q^2 y}</math></td>
  <td align="center"><math>~0</math></td>
  <td align="center"><math>~\frac{q^2 y }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~- \frac{x }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~0</math></td>
</tr>
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
<tr>
  <td align="left" colspan="9">
<table border="0" cellpadding="8" align="center">
<tr>
  <td align="right">
<math>~\ell_{3D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2 + p^4 z^2 )^{- 1 / 2} </math>
  </td>
</tr>
</table>
  </td>
</tr>
</table>
The associated unit vector is, then,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath~\biggl[ \frac{q^2 y }{(x^2 + q^4y^2)^{1 / 2}}  \biggr] - \hat\jmath~\biggl[ \frac{x }{(x^2 + q^4y^2)^{1 / 2}}  \biggr] \, .
</math>
  </td>
</tr>
</table>
It is easy to see that <math>~\hat{e}_2 \cdot \hat{e}_2 = 1</math>.  We also see that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x \ell_{3D}\biggl[ \frac{q^2 y }{(x^2 + q^4y^2)^{1 / 2}}  \biggr]
-
q^2y \ell_{3D} \biggl[ \frac{x }{(x^2 + q^4y^2)^{1 / 2}}  \biggr] = 0 \, ,
</math>
  </td>
</tr>
</table>
so it is clear that these first two unit vectors are orthogonal to one another. 
==Search for the Third Coordinate==
===Cross Product of First Two Unit Vectors===
The cross-product of these two unit vectors should give the third, namely,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_3 = \hat{e}_1 \times \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath~\biggl[ {e}_{1y}{e}_{2z} - {e}_{1z}{e}_{2y} \biggr]
+
\hat\jmath~\biggl[ {e}_{1z}{e}_{2x} - {e}_{1x}{e}_{2z} \biggr]
+
\hat{k}~\biggl[ {e}_{1x}{e}_{2y} - {e}_{1y}{e}_{2x} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath~\biggl[  \frac{x p^2z\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} \biggr]
+
\hat\jmath~\biggl[ \frac{q^2 y p^2z \ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}}  \biggr]
-
\hat{k}~\biggl[ \frac{x^2\ell_{3D} }{(x^2 + q^4y^2)^{1 / 2}}  ~+~ \frac{q^4 y^2\ell_{3D} }{(x^2 + q^4y^2)^{1 / 2}}  \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} \biggl\{~
\hat\imath~( x p^2z )
+
\hat\jmath~(q^2 y p^2z  )
-
\hat{k}~(x^2  + q^4y^2)
~\biggr\} \, .
</math>
  </td>
</tr>
</table>
Inserting these component expressions into the last row of the T8 Direction Cosine table gives &hellip;
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" colspan="9">'''Direction Cosine Components for T8 Coordinates'''</td>
</tr>
<tr>
  <td align="center"><math>~n</math></td>
  <td align="center"><math>~\lambda_n</math></td>
  <td align="center"><math>~h_n</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial x}</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial y}</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial z}</math></td>
  <td align="center"><math>~\gamma_{n1}</math></td>
  <td align="center"><math>~\gamma_{n2}</math></td>
  <td align="center"><math>~\gamma_{n3}</math></td>
</tr>
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center"><math>~(x^2 + q^2 y^2 + p^2 z^2)^{1 / 2} </math></td>
  <td align="center"><math>~\lambda_1 \ell_{3D}</math></td>
  <td align="center"><math>~\frac{x}{\lambda_1}</math></td>
  <td align="center"><math>~\frac{q^2 y}{\lambda_1}</math></td>
  <td align="center"><math>~\frac{p^2 z}{\lambda_1}</math></td>
  <td align="center"><math>~(x) \ell_{3D}</math></td>
  <td align="center"><math>~(q^2 y)\ell_{3D}</math></td>
  <td align="center"><math>~(p^2z) \ell_{3D}</math></td>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center"><math>~\frac{x}{  y^{1/q^2}}</math></td>
  <td align="center"><math>~\frac{1}{\lambda_2}\biggl[\frac{x q^2 y }{(x^2 + q^4y^2)^{1 / 2}}\biggr] </math></td>
  <td align="center"><math>~\frac{\lambda_2}{x}</math></td>
  <td align="center"><math>~-\frac{\lambda_2}{q^2 y}</math></td>
  <td align="center"><math>~0</math></td>
  <td align="center"><math>~\frac{q^2 y }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~- \frac{x }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~0</math></td>
</tr>
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center"><math>~\frac{x p^2 z\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~\frac{q^2 y p^2 z\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~-\frac{(x^2 + q^4 y^2)\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
</tr>
<tr>
  <td align="left" colspan="9">
<table border="0" cellpadding="8" align="center">
<tr>
  <td align="right">
<math>~\ell_{3D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2 + p^4 z^2 )^{- 1 / 2} \, ,</math>
  </td>
</tr>
</table>
  </td>
</tr>
</table>
===Associated h<sub>3</sub> Scale Factor===
<table border="1" align="right" cellpadding="10"><tr><td align="center">
[[File:EUREKA 21Jan2021 sm.png|350px|Whiteboard EUREKA moment]]</td></tr></table>
After working through various scenarios on my whiteboard today (21 January 2021), I propose that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{xp^2z}{(x^2 + q^4y^2)} \, ;</math>
  </td>
<td align="center>&nbsp; &nbsp; &nbsp; &nbsp;</td>
  <td align="right">
<math>~\frac{\partial \lambda_3}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{q^2 y p^2z}{(x^2 + q^4y^2)} \, ;</math>
  </td>
<td align="center>&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp;</td>
  <td align="right">
<math>~\frac{\partial \lambda_3}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-1 \, .</math>
  </td>
</tr>
</table>
This means that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_3^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\sum_{i=1}^3 \biggl( \frac{\partial \lambda_3}{\partial x_i}\biggr)^2</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{xp^2z}{(x^2 + q^4y^2)} \biggr]^2
+
\biggl[ \frac{q^2 y p^2z}{(x^2 + q^4y^2)} \biggr]^2
+
\biggl[ -1 \biggr]^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{p^4z^2(x^2 + q^4y^2)}{(x^2 + q^4y^2)^2}
+
1
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{(x^2 + q^4y^2 +p^4z^2)}{(x^2 + q^4y^2)}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{\ell_{3D}^2 (x^2 + q^4y^2)}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~ h_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\ell_{3D} (x^2 + q^4y^2)^{1 / 2} \, .
</math>
  </td>
</tr>
</table>
This seems to work well because, when combined with the three separate expressions for <math>~\partial \lambda_3/\partial x_i</math>,  this single expression for <math>~h_3</math> generates all three components of the third unit vector, that is, all three direction cosines, <math>~\gamma_{3i}</math>.  All of the elements of this new "EUREKA moment" result have been entered into the following table.
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" colspan="9">'''Direction Cosine Components for T8 Coordinates'''</td>
</tr>
<tr>
  <td align="center"><math>~n</math></td>
  <td align="center"><math>~\lambda_n</math></td>
  <td align="center"><math>~h_n</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial x}</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial y}</math></td>
  <td align="center"><math>~\frac{\partial \lambda_n}{\partial z}</math></td>
  <td align="center"><math>~\gamma_{n1}</math></td>
  <td align="center"><math>~\gamma_{n2}</math></td>
  <td align="center"><math>~\gamma_{n3}</math></td>
</tr>
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center"><math>~(x^2 + q^2 y^2 + p^2 z^2)^{1 / 2} </math></td>
  <td align="center"><math>~\lambda_1 \ell_{3D}</math></td>
  <td align="center"><math>~\frac{x}{\lambda_1}</math></td>
  <td align="center"><math>~\frac{q^2 y}{\lambda_1}</math></td>
  <td align="center"><math>~\frac{p^2 z}{\lambda_1}</math></td>
  <td align="center"><math>~(x) \ell_{3D}</math></td>
  <td align="center"><math>~(q^2 y)\ell_{3D}</math></td>
  <td align="center"><math>~(p^2z) \ell_{3D}</math></td>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center"><math>~\frac{x}{  y^{1/q^2}}</math></td>
  <td align="center"><math>~\frac{1}{\lambda_2}\biggl[\frac{x q^2 y }{(x^2 + q^4y^2)^{1 / 2}}\biggr] </math></td>
  <td align="center"><math>~\frac{\lambda_2}{x}</math></td>
  <td align="center"><math>~-\frac{\lambda_2}{q^2 y}</math></td>
  <td align="center"><math>~0</math></td>
  <td align="center"><math>~\frac{q^2 y }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~- \frac{x }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~0</math></td>
</tr>
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center">---</td>
  <td align="center"><math>~\ell_{3D}(x^2 + q^4 y^2)^{1 / 2}</math></td>
  <td align="center"><math>~\frac{xp^2z}{(x^2 + q^4y^2)} </math></td>
  <td align="center"><math>~\frac{q^2 y p^2z}{(x^2 + q^4y^2)}</math></td>
  <td align="center"><math>~-1</math></td>
  <td align="center"><math>~\frac{x p^2 z\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~\frac{q^2 y p^2 z\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~-\frac{(x^2 + q^4 y^2)\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
</tr>
<tr>
  <td align="left" colspan="9">
<table border="0" cellpadding="8" align="center">
<tr>
  <td align="right">
<math>~\ell_{3D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2 + p^4 z^2 )^{- 1 / 2} </math>
  </td>
</tr>
</table>
  </td>
</tr>
</table>
===What is the Third Coordinate Function, &lambda;<sub>3</sub>===
The remaining [https://en.wikipedia.org/wiki/The_$64,000_Question $64,000 question] is, "What is the actual expression for <math>~\lambda_3(x, y, z)</math>&nbsp; ? &nbsp;" 
Notice that the (partial) derivatives of <math>~\lambda_3</math> with respect to <math>~x</math> and <math>~y</math> may be rewritten, respectively, in the form
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\biggl( \frac{q^2 y}{p^2z} \biggr) \frac{\partial \lambda_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{q^2 y}{x(1 + \eta^2)}
=
\frac{\eta}{(1 + \eta^2)} \, ,
</math>
&nbsp; &nbsp; &nbsp; and,
  </td>
</tr>
<tr>
  <td align="right">
<math>~\biggl( \frac{x}{p^2z} \biggr) \frac{\partial \lambda_3}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{q^2y}{x(1 + \eta^2)}
=
\frac{\eta}{(1 + \eta^2)} \, ,
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\eta</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{q^2y}{x} ~~~~\Rightarrow~~ \frac{\partial \ln\eta}{\partial \ln x} = -1 \, , ~~~~\frac{\partial \ln\eta}{\partial \ln y} = +1 \, .
</math>
  </td>
</tr>
</table>
Then, after searching through the [[User:Tohline/Appendix/References#CRC|CRC Mathematical Handbook's]] pages of familiar derivative expressions, we appreciate that
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{d}{dx_i} \biggl[\frac{1}{\cosh\gamma}\biggr]</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{\tanh\gamma}{\cosh\gamma} \biggr] \frac{d\gamma}{dx_i} \, .</math>
  </td>
</tr>
</table>
Hence, it will be useful to adopt the mapping,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\eta</math>
  </td>
  <td align="center">
<math>~~~\rightarrow~~~</math>
  </td>
  <td align="left">
<math>~\sinh \gamma \, ,</math>
  </td>
</tr>
</table>
because the right-hand side of both partial-derivative expressions becomes,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\eta}{(1+\eta^2)}</math>
  </td>
  <td align="center">
<math>~~~\rightarrow~~~</math>
  </td>
  <td align="left">
<math>~\frac{\sinh\gamma}{\cosh^2\gamma} = \frac{\tanh\gamma}{\cosh\gamma} \, .</math>
  </td>
</tr>
</table>
====Guess A====
In particular, this suggests that we set,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{A}{\cosh\gamma} \, ,</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\gamma</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\sinh^{-1}\eta = \pm \cosh^{-1}[\eta^2 + 1]^{1 / 2} \, .</math>
  </td>
</tr>
</table>
In other words,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~A[\eta^2 + 1]^{-1 / 2} \, .</math>
  </td>
</tr>
</table>
Let's check the first and second partial derivatives.
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \frac{A}{2} \biggl[ \frac{2\eta}{(\eta^2 + 1)^{3 / 2}} \biggr] \frac{\partial \eta}{\partial x}
</math>
  </td>
</tr>
</table>
====Guess B====
What if,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{2}\ln(1+\eta^2) \, .</math>
  </td>
</tr>
</table>
Then,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{d\lambda_3}{d\eta}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\eta}{(1+\eta^2)} \, ,
</math>
  </td>
</tr>
</table>
in which case we find,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial\lambda_3}{\partial x_i} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{d\lambda_3}{d\eta} \cdot \frac{\partial\eta}{\partial x_i}</math>
  </td>
</tr>
</table>
which means,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial\lambda_3}{\partial x} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\eta}{(1+\eta^2)}  \biggl[ - \frac{q^2 y}{x^2} \biggr]
=
- \frac{x}{(x^2 + q^4y^2)} \biggl[ \frac{q^4 y^2}{x^2} \biggr] \, .
</math>
  </td>
</tr>
</table>
====Guess C====
What if,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{2} \ln(1+\eta^{-2}) \, .</math>
  </td>
</tr>
</table>
Then,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{d\lambda_3}{d\eta}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-
\frac{\eta^{-3}}{(1+\eta^{-2})}
=
-\frac{ \eta^{-1}}{(1+\eta^{2})}
</math>
  </td>
</tr>
</table>
in which case we find,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial\lambda_3}{\partial x_i} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{d\lambda_3}{d\eta} \cdot \frac{\partial\eta}{\partial x_i}</math>
  </td>
</tr>
</table>
which means,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial\lambda_3}{\partial x} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{\eta^{-1}}{(1+\eta^2)}  \biggl[ - \frac{q^2 y}{x^2} \biggr]
=
\frac{x}{(x^2 + q^4y^2)}  \, ;
</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial\lambda_3}{\partial y} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{\eta^{-1}}{(1+\eta^2)}  \biggl[ \frac{q^2 }{x} \biggr]
=
-\frac{x^2}{y(x^2+ q^4y^2)} 
</math>
  </td>
</tr>
</table>
==Inverting Coordinate Relations==
===In a Plane Perpendicular to the Z-Axis===
====General Case====
At a fixed value of <math>~z = z_0</math>, let's invert the <math>~\lambda_1(x, y)</math> and <math>~\lambda_2(x, y)</math> relations to obtain expressions for <math>~x(\lambda_1, \lambda_2)</math> and <math>~y(\lambda_1, \lambda_2)</math>.  Perhaps this will help us determine what the third coordinate expression should be. 
We start with,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_1</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^2y^2 + p^2 z_0^2)^{1 / 2} \, ;</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\lambda_2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{x}{y^{1/q^2}} \, .</math>
  </td>
</tr>
</table>
This means that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\ln y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~q^2(\ln x - \ln\lambda_2 ) \, ,</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~q^2 y^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(\lambda_1^2 - p^2z_0^2) - x^2</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~2\ln y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\ln\biggl\{ \frac{1}{q^2} \biggl[ (\lambda_1^2 - p^2z_0^2) - x^2 \biggr] \biggr\} \, .</math>
  </td>
</tr>
</table>
Together, this gives,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\ln\biggl\{ \frac{1}{q^2} \biggl[ (\lambda_1^2 - p^2z_0^2) - x^2 \biggr] \biggr\} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2q^2(\ln x - \ln\lambda_2 ) </math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\ln\biggl[\frac{x}{\lambda_2} \biggr]^{2q^2}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~(\lambda_1^2 - p^2z_0^2) - x^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~q^2\biggl[\frac{x}{\lambda_2} \biggr]^{2q^2} </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~ x^2 +q^2\biggl[\frac{x}{\lambda_2} \biggr]^{2q^2} +  (p^2z_0^2 - \lambda_1^2) </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, .</math>
  </td>
</tr>
</table>
====What if Axisymmetric (q<sup>2</sup> = 1)====
In an axisymmetric configuration, <math>~q^2 = 1</math> and <math>~(\lambda_1^2 - p^2z_0^2) = \varpi^2</math>, so this general expression for <math>~x</math> becomes,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x^2 + \biggl[\frac{x}{\lambda_2} \biggr]^{2} -\varpi^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x^2\biggl[\frac{1 + \lambda_2^2}{\lambda_2^2}\biggr]  -\varpi^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~x</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\varpi \biggl[\frac{\lambda_2^2}{1 + \lambda_2^2}\biggr]^{1 / 2} \, .
</math>
  </td>
</tr>
</table>
Given that, for axisymmetric systems,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~x</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\varpi \cos\varphi \, ,</math>
  </td>
</tr>
</table>
we conclude that when <math>~q^2 = 1</math>,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\biggl[\frac{\lambda_2^2}{1 + \lambda_2^2}\biggr]^{1 / 2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\cos\varphi </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~\lambda_2^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\cot^2\varphi </math>
  </td>
</tr>
</table>
====What if q<sup>2</sup> = 2====
For example, if we choose <math>~q^2 = 2</math>, we have a quadratic expression for <math>~x^2</math>, namely,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~0 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x^2 + 2\biggl[\frac{x}{\lambda_2} \biggr]^{4} +  (p^2z_0^2 - \lambda_1^2) </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 0 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x^{4} + \frac{1}{2} \lambda_2^4~ x^2 +  \frac{1}{2}\lambda_2^4(p^2z_0^2 - \lambda_1^2) </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 2x^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \frac{\lambda_2^4}{2} \pm \biggl[ \frac{\lambda_2^8}{4} - 2\lambda_2^4(p^2z_0^2 - \lambda_1^2) \biggr]^{1 / 2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\lambda_2^4}{2} \biggl\{ \biggl[ 1 - \frac{8(p^2z_0^2 - \lambda_1^2)}{\lambda_2^4} \biggr]^{1 / 2} - 1 \biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ x^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\lambda_2^4}{4} \biggl\{ \biggl[ 1 - \frac{8(p^2z_0^2 - \lambda_1^2)}{\lambda_2^4} \biggr]^{1 / 2} - 1 \biggr\} \, .
</math>
  </td>
</tr>
</table>
Given that, for <math>~q^2 = 2</math>, one of the two defining expression means, <math>~\lambda_2 = x/\sqrt{y}</math>, we also have,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~y </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\lambda_2^2}{4} \biggl\{ \biggl[ 1 - \frac{8(p^2z_0^2 - \lambda_1^2)}{\lambda_2^4} \biggr]^{1 / 2} - 1 \biggr\} \, .
</math>
  </td>
</tr>
</table>
===New 2<sup>nd</sup> Coordinate===
Apparently it will be cleaner to define a new "2<sup>nd</sup> coordinate," <math>~\kappa_2</math>, such that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\kappa_2^{2q^2}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
q^2\biggl[\frac{1}{\lambda_2} \biggr]^{2q^2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~\kappa_2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
q^{1/q^2}\biggl[\frac{1}{\lambda_2} \biggr]
=
q^{1/q^2}\biggl[\frac{y^{1/q^2}}{x} \biggr]
=
\frac{(qy)^{1/q^2}}{x}  \, .
</math>
  </td>
</tr>
</table>
(With this new definition, <math>~\kappa_2 \sim \tan\varphi</math>; it is exactly this when <math>~q^2 = 1</math>.)  Then we can rewrite the last expression from the [[#General_Case|above general case]] as,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~x^2 +(\kappa_2 x)^{2q^2} -  (\lambda_1^2 - p^2z_0^2 ) </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, .</math>
  </td>
</tr>
</table>
When <math>~q=1</math> (the axisymmetric case), this gives,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~x^2(1+\kappa_2^2)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ (\lambda_1^2 - p^2z_0^2 ) </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~\frac{x^2}{(\lambda_1^2 - p^2z_0^2 ) }</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{1}{(1+\kappa_2^2)} \, ,</math>
  </td>
</tr>
</table>
which means that <math>~\kappa_2 = \tan\varphi</math>.  And, for the case of <math>~q^2 = 2</math>, after making the substitution,
<div align="center"><math>~\lambda_2 \rightarrow (q)^{1/q^2}\kappa_2^{-1} = \frac{2^{1 / 4}}{\kappa_2} \, ,</math></div>
we find,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~x^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\lambda_2^4}{4} \biggl\{ \biggl[ 1 - \frac{8(p^2z_0^2 - \lambda_1^2)}{\lambda_2^4} \biggr]^{1 / 2} - 1 \biggr\} \, .
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~ x^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{2\kappa_2^4} \biggl\{ \biggl[ 1 + 4\kappa_2^4(\lambda_1^2 - p^2z_0^2 ) \biggr]^{1 / 2} - 1 \biggr\} \, ,
</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~y = \biggl(\frac{\kappa_2^2}{2^{1 / 2}} \biggr) x^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{\kappa_2^2}{2^{1 / 2}} \biggr)
\frac{1}{2\kappa_2^4} \biggl\{ \biggl[ 1 + 4\kappa_2^4(\lambda_1^2 - p^2z_0^2 ) \biggr]^{1 / 2} - 1 \biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{2^{3 / 2}\kappa_2^2} \biggl\{ \biggl[ 1 + 4\kappa_2^4(\lambda_1^2 - p^2z_0^2 ) \biggr]^{1 / 2} - 1 \biggr\} \, .
</math>
  </td>
</tr>
</table>
=Angle Between Unit Vectors=
We begin by restating that the coordinate, scale factor, and unit vector associated with the normal to our ellipsoidal surface are,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_1</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^2y^2 + p^2z^2)^{1 / 2} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~h_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\lambda_1 \ell_{3D} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath (x\ell_{3D} )
+
\hat\jmath (q^2 y\ell_{3D} )
+
\hat{k} (p^2z\ell_{3D} ) \, .
</math>
  </td>
</tr>
</table>
In the [[#TableKappa8|Table below titled, "Direction Cosines Components for &kappa;8 Coordinates"]], there are two fully-formed unit vectors that are each orthogonal to the (first) unit vector that is normal to the ellipsoid's surface.  Here we will refer to the coordinates of these two fully-formed unit vectors as,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{xy^{1/q^2}}{z^{2/p^2}}</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp; </td>
  <td align="right">
<math>~\kappa_2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{(qy)^{1/q^2}}{x} \, .</math>
  </td>
</tr>
</table>
The associated scale factors and unit vectors are given by the following expressions:
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_{\lambda_2}^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{\partial\lambda_2}{\partial x} \biggr)^2
+
\biggl(\frac{\partial\lambda_2}{\partial y} \biggr)^2
+
\biggl(\frac{\partial\lambda_2}{\partial z} \biggr)^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{\lambda_2}{x} \biggr)^2
+
\biggl(\frac{\lambda_2}{q^2y} \biggr)^2
+
\biggl(- \frac{2\lambda_2}{p^2z} \biggr)^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \lambda_2^2 \biggl[
\frac{q^4y^2p^4z^2 + x^2p^4z^2 + 4x^2 q^4y^2}{x^2 q^4y^2 p^4z^2}
\biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ h_{\lambda_2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{x q^2y p^2z}{\lambda_2 \mathcal{D}}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \hat{e}_{\lambda_2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl[ h_{\lambda_2} \biggl(\frac{\partial \lambda_2}{\partial x}\biggr)\biggr]
+
\hat\jmath \biggl[ h_{\lambda_2} \biggl(\frac{\partial \lambda_2}{\partial y}\biggr)\biggr]
+
\hat{k} \biggl[ h_{\lambda_2} \biggl(\frac{\partial \lambda_2}{\partial z}\biggr)\biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl[ \frac{q^2y p^2z}{\mathcal{D}} \biggr]
+
\hat\jmath \biggl[ \frac{xp^2z}{\mathcal{D}} \biggr]
-
\hat{k} \biggl[ \frac{2x q^2y }{\mathcal{D}} \biggr] \, .
</math>
  </td>
</tr>
</table>
<table border="1" width="80%" align="center" cellpadding="8"><tr><td align="left">
With regard to orthogonality, note that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_{\lambda_2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[
\hat\imath (x\ell_{3D} )
+
\hat\jmath (q^2 y\ell_{3D} )
+
\hat{k} (p^2z\ell_{3D} )
\biggr] \cdot
\biggl\{
\hat\imath \biggl[ \frac{q^2y p^2z}{\mathcal{D}} \biggr]
+
\hat\jmath \biggl[ \frac{xp^2z}{\mathcal{D}} \biggr]
-
\hat{k} \biggl[ \frac{2x q^2y }{\mathcal{D}} \biggr]
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(x\ell_{3D} )\biggl[ \frac{q^2y p^2z}{\mathcal{D}} \biggr]
+
(q^2 y\ell_{3D} )\biggl[ \frac{xp^2z}{\mathcal{D}} \biggr]
-
(p^2z\ell_{3D} )\biggl[ \frac{2x q^2y }{\mathcal{D}} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
0 \, .</math>
  </td>
</tr>
</table>
</td></tr></table>
And,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_{\kappa_2}^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{\partial\kappa_2}{\partial x} \biggr)^2
+
\biggl(\frac{\partial\kappa_2}{\partial y} \biggr)^2
+
\biggl(\frac{\partial\kappa_2}{\partial z} \biggr)^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(-\frac{\kappa_2}{x} \biggr)^2
+
\biggl(\frac{\kappa_2}{q^2 y} \biggr)^2
=
\frac{\kappa_2^2}{x^2q^4y^2}\biggl[x^2 + q^4y^2\biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~h_{\kappa_2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{xq^2y}{\kappa_2 (x^2 + q^4y^2)^{1 / 2}}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~\hat{e}_{\kappa_2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-~\hat\imath \biggl[ \frac{q^2y}{(x^2 + q^4y^2)^{1 / 2}}\biggr]
+
\hat\jmath \biggl[ \frac{x}{(x^2 + q^4y^2)^{1 / 2}} \biggr]
+
\hat{k} \biggl[ 0 \biggr] \, .
</math>
  </td>
</tr>
</table>
<table border="1" width="80%" align="center" cellpadding="8"><tr><td align="left">
Again, note that with regard to orthogonality,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_{\kappa_2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[
\hat\imath (x\ell_{3D} )
+
\hat\jmath (q^2 y\ell_{3D} )
+
\hat{k} (p^2z\ell_{3D} )
\biggr] \cdot
\biggl\{
-\hat\imath \biggl[ \frac{q^2y}{(x^2 + q^4y^2)^{1 / 2}}\biggr]
+
\hat\jmath \biggl[ \frac{x}{(x^2 + q^4y^2)^{1 / 2}} \biggr]
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-~(x\ell_{3D} )\biggl[ \frac{q^2y}{(x^2 + q^4y^2)^{1 / 2}} \biggr]
+
(q^2 y\ell_{3D} )\biggl[ \frac{x}{(x^2 + q^4y^2)^{1 / 2}}  \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
0 \, .</math>
  </td>
</tr>
</table>
</td></tr></table>
From this pair of orthogonality checks, we appreciate that both unit vectors always lie in the plane that is tangent to the surface of our ellipsoid.  Next, let's determine the angle, <math>~\alpha</math>, between these two unit vectors as measured in the relevant tangent-plane.
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\cos\alpha \equiv \hat{e}_{\lambda_2} \cdot \hat{e}_{\kappa_2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{
\hat\imath \biggl[ \frac{q^2y p^2z}{\mathcal{D}} \biggr]
+
\hat\jmath \biggl[ \frac{xp^2z}{\mathcal{D}} \biggr]
-
\hat{k} \biggl[ \frac{2x q^2y }{\mathcal{D}} \biggr]
\biggr\}
\biggl\{
-\hat\imath \biggl[ \frac{q^2y}{(x^2 + q^4y^2)^{1 / 2}}\biggr]
+
\hat\jmath \biggl[ \frac{x}{(x^2 + q^4y^2)^{1 / 2}} \biggr]
+
\hat{k} \biggl[ 0 \biggr]
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\biggl[ \frac{q^2y}{(x^2 + q^4y^2)^{1 / 2}}\biggr]\biggl[ \frac{q^2y p^2z}{\mathcal{D}} \biggr]
+
\biggl[ \frac{x}{(x^2 + q^4y^2)^{1 / 2}} \biggr]\biggl[ \frac{xp^2z}{\mathcal{D}} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{p^2 z}{\mathcal{D} (x^2 + q^4y^2)^{1 / 2}} \biggr] \biggl[ x^2 - q^4y^2 \biggr] \, .
</math>
  </td>
</tr>
</table>
----
Let's again visit the unit vector that we know lies in the tangent-plane and is always orthogonal to <math>~\lambda_2</math>, namely,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_{\lambda_3}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-~\hat\imath \biggl[ x(2q^4 y^2 + p^4z^2 )\biggr]\frac{\ell_{3D}}{\mathcal{D}}
+
\hat\jmath \biggl[q^2y(p^4z^2 + 2x^2) \biggr]\frac{\ell_{3D}}{\mathcal{D}}
+
\hat{k} \biggl[p^2z (x^2 - q^4y^2) \biggr]\frac{\ell_{3D}}{\mathcal{D}}
\, .
</math>
  </td>
</tr>
</table>
We acknowledge that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_{\lambda_2} \cdot \hat{e}_{\lambda_3}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{
\hat\imath \biggl[ \frac{q^2y p^2z}{\mathcal{D}} \biggr]
+
\hat\jmath \biggl[ \frac{xp^2z}{\mathcal{D}} \biggr]
-
\hat{k} \biggl[ \frac{2x q^2y }{\mathcal{D}} \biggr]
\biggr\}
\biggl\{
-~\hat\imath \biggl[ x(2q^4 y^2 + p^4z^2 )\biggr]\frac{\ell_{3D}}{\mathcal{D}}
+
\hat\jmath \biggl[q^2y(p^4z^2 + 2x^2) \biggr]\frac{\ell_{3D}}{\mathcal{D}}
+
\hat{k} \biggl[p^2z (x^2 - q^4y^2) \biggr]\frac{\ell_{3D}}{\mathcal{D}}
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-~\biggl[ xq^2yp^2z(2q^4 y^2 + p^4z^2 )\biggr]\frac{\ell_{3D}}{\mathcal{D}^2}
+
\biggl[x q^2y p^2z(p^4z^2 + 2x^2) \biggr]\frac{\ell_{3D}}{\mathcal{D}^2}
-
\biggl[2xq^2y p^2z (x^2 - q^4y^2) \biggr]\frac{\ell_{3D}}{\mathcal{D}^2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[
-(2q^4 y^2 + p^4z^2 )
+
(p^4z^2 + 2x^2)
-
2 (x^2 - q^4y^2) 
\biggr] \frac{(x q^2y p^2z) \ell_{3D}}{\mathcal{D}^2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
0 \, .
</math>
  </td>
</tr>
</table>
=Kappa (&kappa;8) Coordinates=
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\kappa_1</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^2y^2 + p^2 z^2)^{1 / 2} \, ;</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\kappa_3</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\tan^{-1}\biggl[ \frac{(qy)^{1/q^2}}{x} \biggr] \, .</math>
  </td>
</tr>
</table>
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial\kappa_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[1 + \frac{(qy)^{2/q^2}}{x^2} \biggr]^{-1} \biggl[- \frac{(qy)^{1/q^2}}{x^2}  \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{\sin^2\kappa_3}{(qy)^{1/q^2}} = -\frac{\sin^2\kappa_3}{x\tan\kappa_3}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{\sin\kappa_3 \cos\kappa_3}{x}\, .
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial\kappa_3}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[1 + \frac{(qy)^{2/q^2}}{x^2} \biggr]^{-1} \biggl[\frac{q^{1/q^2} y^{1/q^2}}{q^2 x y}  \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[1 + \frac{(qy)^{2/q^2}}{x^2} \biggr]^{-1} \biggl[\frac{(qy)^{1/q^2}}{x }  \biggr] \frac{1}{q^2y}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{q^2y}\biggl[\frac{\tan\kappa_3}{1 + \tan^2\kappa_3 }\biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
+ \frac{\sin\kappa_3 \cos\kappa_3}{q^2y} \, .
</math>
  </td>
</tr>
</table>
Therefore,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_3^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[-\frac{\sin\kappa_3 \cos\kappa_3}{x} \biggr]^2
+
\biggl[\frac{\sin\kappa_3 \cos\kappa_3}{q^2y} \biggr]^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2)
\biggl[\frac{\sin\kappa_3 \cos\kappa_3}{xq^2y} \biggr]^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~ h_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2)^{-1 / 2}
\biggl[\frac{xq^2y}{\sin\kappa_3 \cos\kappa_3} \biggr] \, .
</math>
  </td>
</tr>
</table>
<span id="TableKappa8">&nbsp;</span>
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" colspan="9">'''Direction Cosine Components for &kappa;8 Coordinates'''</td>
</tr>
<tr>
  <td align="center"><math>~n</math></td>
  <td align="center"><math>~\kappa_n</math></td>
  <td align="center"><math>~h_n</math></td>
  <td align="center"><math>~\frac{\partial \kappa_n}{\partial x}</math></td>
  <td align="center"><math>~\frac{\partial \kappa_n}{\partial y}</math></td>
  <td align="center"><math>~\frac{\partial \kappa_n}{\partial z}</math></td>
  <td align="center"><math>~\gamma_{n1}</math></td>
  <td align="center"><math>~\gamma_{n2}</math></td>
  <td align="center"><math>~\gamma_{n3}</math></td>
</tr>
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center"><math>~(x^2 + q^2 y^2 + p^2 z^2)^{1 / 2} </math></td>
  <td align="center"><math>~\kappa_1 \ell_{3D}</math></td>
  <td align="center"><math>~\frac{x}{\kappa_1}</math></td>
  <td align="center"><math>~\frac{q^2 y}{\kappa_1}</math></td>
  <td align="center"><math>~\frac{p^2 z}{\kappa_1}</math></td>
  <td align="center"><math>~(x) \ell_{3D}</math></td>
  <td align="center"><math>~(q^2 y)\ell_{3D}</math></td>
  <td align="center"><math>~(p^2z) \ell_{3D}</math></td>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center">---</td>
  <td align="center"><math>~\ell_{3D}(x^2 + q^4 y^2)^{1 / 2}</math></td>
  <td align="center"><math>~\frac{xp^2z}{(x^2 + q^4y^2)} </math></td>
  <td align="center"><math>~\frac{q^2 y p^2z}{(x^2 + q^4y^2)}</math></td>
  <td align="center"><math>~-1</math></td>
  <td align="center"><math>~\frac{x p^2 z\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~\frac{q^2 y p^2 z\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~-\frac{(x^2 + q^4 y^2)\ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}} </math></td>
</tr>
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center"><math>~\tan^{-1}\biggl[ \frac{(qy)^{1/q^2}}{x} \biggr]</math></td>
  <td align="center"><math>~\frac{1}{\sin\kappa_3 \cos\kappa_3}\biggl[\frac{x q^2 y }{(x^2 + q^4y^2)^{1 / 2}}\biggr] </math></td>
  <td align="center"><math>~-\frac{\sin\kappa_3 \cos\kappa_3}{x}</math></td>
  <td align="center"><math>~\frac{\sin\kappa_3 \cos\kappa_3}{q^2 y}</math></td>
  <td align="center"><math>~0</math></td>
  <td align="center"><math>~- \frac{q^2 y }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~\frac{x }{(x^2 + q^4y^2)^{1 / 2}} </math></td>
  <td align="center"><math>~0</math></td>
</tr>
<tr>
  <td align="left" colspan="9">
<table border="0" cellpadding="8" align="center">
<tr>
  <td align="right">
<math>~\ell_{3D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2 + p^4 z^2 )^{- 1 / 2} </math>
  </td>
</tr>
</table>
  </td>
</tr>
<tr>
  <td align="center"><math>~4</math></td>
  <td align="center"><math>~\frac{x y^{1/q^2}}{ z^{2/p^2}}</math></td>
  <td align="center"><math>~\frac{1}{\kappa_4}\biggl[\frac{x q^2 y p^2 z}{ \mathcal{D}}\biggr] </math></td>
  <td align="center"><math>~\frac{\kappa_4}{x}</math></td>
  <td align="center"><math>~\frac{\kappa_4}{q^2 y}</math></td>
  <td align="center"><math>~-\frac{2\kappa_4}{p^2 z}</math></td>
  <td align="center"><math>~\frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(\frac{1}{x}\biggr)</math></td>
  <td align="center"><math>~ \frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(\frac{1}{q^2y}\biggr)</math></td>
  <td align="center"><math>~\frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(-\frac{2}{p^2z}\biggr)</math></td>
</tr>
<tr>
  <td align="center"><math>~5</math></td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center"><math>~-\frac{\ell_{3D}}{\mathcal{D}}\biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]</math></td>
  <td align="center"><math>~\frac{\ell_{3D}}{\mathcal{D}}\biggl[ q^2 y(p^4z^2 + 2x^2 )  \biggl]</math></td>
  <td align="center"><math>~\frac{\ell_{3D}}{\mathcal{D}}\biggl[ p^2z( x^2 - q^4y^2 ) \biggl]</math></td>
</tr>
<tr>
  <td align="left" colspan="9">
<table border="0" cellpadding="8" align="center">
<tr>
  <td align="right">
<math>~\mathcal{D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)^{1 / 2} \, .</math>
  </td>
</tr>
</table>
Also note &hellip;
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~AB \equiv \biggl( \frac{\mathcal{D}}{\ell_{3D}} \biggr)^{2} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]^2
+
\biggl[ q^2 y(p^4z^2 + 2x^2 )  \biggl]^2
+
\biggl[ p^2z( x^2 - q^4y^2 ) \biggl]^2 \, ,
</math>
  </td>
</tr>
</table>
and the partial derivatives of <math>~A</math> and <math>~B</math> are detailed in [[User:Tohline/Appendix/Ramblings/ConcentricEllipsodalCoordinates#ABderivatives|an accompanying discussion]].
  </td>
</tr>
</table>
The direction-cosines of the second unit vector &#8212; as has already been inserted into the "&kappa;8 coordinates" table &#8212; should be obtainable from the first and third unit vectors via the cross product,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_2 = \hat{e}_3 \times \hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl[ e_{3y}e_{1z} - e_{3z} e_{1y} \biggr]
+
\hat\jmath \biggl[ e_{3z}e_{1x} - e_{3x} e_{1z} \biggr]
+
\hat{k} \biggl[ e_{3x}e_{1y} - e_{3y} e_{1x} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl[ \frac{x (p^2z) \ell_{3D} }{(x^2 + q^4y^2)^{1 / 2}} \biggr]
+
\hat\jmath \biggl[  \frac{q^2 y  (p^2z) \ell_{3D}}{(x^2 + q^4y^2)^{1 / 2}}\biggr]
-
\hat{k} \biggl[ \frac{(x^2 + q^4y^2) \ell_{3D} }{(x^2 + q^4y^2)^{1 / 2}} \biggr] \, .
</math>
  </td>
</tr>
</table>
The other boxes in the n = 2 row have been drawn from [[User:Tohline/Appendix/Ramblings/ConcentricEllipsodalT8Coordinates#Associated_h3_Scale_Factor|our accompanying EUREKA! moment]] and the n = 3 row of the table that details "Direction Cosine Components for T8 Coordinates."
==Attempt 1==
Let's try &hellip;
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\kappa_2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\tan^{-1}\biggl[ \frac{x(qy)^{1/q^2}}{(pz)^{1/p^2}} \biggr] \, ,</math>
  </td>
</tr>
</table>
which leads to,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial\kappa_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[1 + \frac{x^2(qy)^{2/q^2}}{(pz)^{2/p^2}} \biggr]^{-1} \biggl[ \frac{(qy)^{1/q^2}}{(pz)^{1/p^2}} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{x}\biggl[1+\tan^2\kappa_2\biggr]^{-1} \tan\kappa_2 = \frac{\sin\kappa_2 \cos\kappa_2}{x} \, ;
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial\kappa_3}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[1 + \frac{x^2(qy)^{2/q^2}}{(pz)^{2/p^2}} \biggr]^{-1} \biggl[ \frac{xq^{1/q^2}(y)^{1/q^2}}{q^2 y(pz)^{1/p^2}} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{q^2y}\biggl[1+\tan^2\kappa_2\biggr]^{-1} \tan\kappa_2 = \frac{\sin\kappa_2 \cos\kappa_2}{q^2y} \, ;
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial\kappa_3}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \frac{\sin\kappa_2 \cos\kappa_2}{p^2z} \, .
</math>
  </td>
</tr>
</table>
Hence,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_2^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\sin^2\kappa_2 \cos^2\kappa_2 \biggl[ \frac{1}{x^2} + \frac{1}{q^4y^2} + \frac{1}{p^4z^2} \biggr]</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\sin^2\kappa_2 \cos^2\kappa_2 \biggl[ \frac{x^2 + q^4 y^2 + p^4 z^2}{x^2 q^4y^2 p^4 z^2}  \biggr]
=
\biggl[ \frac{\sin^2\kappa_2 \cos^2\kappa_2 }{x^2 q^4y^2 p^4 z^2 \ell_{3D}^2}  \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~ h_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x q^2y p^2 z \ell_{3D}}{\sin\kappa_2 \cos\kappa_2 }  \biggr] \, .
</math>
  </td>
</tr>
</table>
The three direction-cosines are, then,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\gamma_{21} = h_2 \biggl(\frac{\partial \kappa_2}{\partial x}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x q^2y p^2 z \ell_{3D}}{\sin\kappa_2 \cos\kappa_2 }  \biggr] \frac{\sin\kappa_2 \cos\kappa_2}{x}
=
q^2y p^2z \ell_{3D} \, .
</math>
  </td>
</tr>
</table>


=T9 Coordinates=
=T9 Coordinates=

Latest revision as of 18:41, 23 July 2021


T9 Coordinates

Establish 1st and 3rd Coordinates

λ1

(x2+q2y2+p2z2)1/2,

λ3

[q2yxq2]1/(q21)=q2/(q21)y1/(q21)xq2/(1q2).

λ3x

=

(q21q2)λ3x,

λ3y

=

(1q21)λ3y,

h32

=

[(q21q2)λ3x]2+[(1q21)λ3y]2=λ32(q21)2[x2+q4y2x2y2]

h3

=

(q21)λ3[xy(x2+q4y2)1/2].

Hence, the three direction-cosines are,

γ31=h3(λ3x)

=

(q21)λ3[xy(x2+q4y2)1/2](q21q2)λ3x=[q2y(x2+q4y2)1/2],

γ32=h3(λ3y)

=

(q21)λ3[xy(x2+q4y2)1/2](1q21)λ3y=[x(x2+q4y2)1/2].

And the position vector is given by the expression,

x

=

e^1(γ11x+γ12y+γ13z)+e^2(γ21x+γ22y+γ23z)+e^3(γ31x+γ32y+γ33z)

 

=

e^1[x23D+q2y23D+p2z23D]+e^2(γ21x+γ22y+γ23z)+e^3[q2xy(x2+q4y2)1/2+xy(x2+q4y2)1/2]

 

=

e^1[λ123D]+e^2(γ21x+γ22y+γ23z)e^3[(q21)xy(x2+q4y2)1/2].

Guess 2nd Coordinate

Unspecified Coefficients

λ2

(ax2+by2+ez2)1/2(ax2+by2)1/2

 

=

[(ax2+by2)(ax2+by2+ez2)]1/2

 

=

[1+ez2(ax2+by2)]1/2.

λ2x

=

12[1+ez2(ax2+by2)]3/2x[ez2(ax2+by2)1]

 

=

ez22[1+ez2(ax2+by2)]3/2(1)[(ax2+by2)2]2ax

 

=

axez2(ax2+by2)2[1+ez2(ax2+by2)]3/2

 

=

axez2λ23(ax2+by2)2,

λ2y

=

12[1+ez2(ax2+by2)]3/2y[ez2(ax2+by2)1]

 

=

ez22[1+ez2(ax2+by2)]3/2(1)[(ax2+by2)2]2by

 

=

byez2(ax2+by2)2[1+ez2(ax2+by2)]3/2

 

=

byez2λ23(ax2+by2)2,

λ2z

=

12[1+ez2(ax2+by2)]3/2z[ez2(ax2+by2)1]

 

=

ezλ23(ax2+by2).

Hence,

h22

=

[axez2λ23(ax2+by2)2]2+[byez2λ23(ax2+by2)2]2+[ezλ23(ax2+by2)]2

 

=

[ezλ23(ax2+by2)2]2[(axz)2+(byz)2+(ax2+by2)2]

h2

=

[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]1/2.

The three direction-cosines are, then,

γ21=h2(λ2x)

=

axez2λ23(ax2+by2)2[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]1/2

 

=

axz[(axz)2+(byz)2+(ax2+by2)2]1/2

γ22=h2(λ2y)

=

byez2λ23(ax2+by2)2[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]1/2

 

=

byz[(axz)2+(byz)2+(ax2+by2)2]1/2

γ23=h2(λ2z)

=

ezλ23(ax2+by2)[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]1/2

 

=

(ax2+by2)[(axz)2+(byz)2+(ax2+by2)2]1/2.

Direction Cosine Components for T9 Coordinates
n λn hn λnx λny λnz γn1 γn2 γn3
(1) (2) (3) (4) (5) (6) (7) (8) (9)
1 (x2+q2y2+p2z2)1/2 λ13D xλ1 q2yλ1 p2zλ1 (x)3D (q2y)3D (p2z)3D
2A --- 3D(x2+q4y2)1/2 xp2z(x2+q4y2) q2yp2z(x2+q4y2) 1 xp2z3D(x2+q4y2)1/2 q2yp2z3D(x2+q4y2)1/2 (x2+q4y2)3D(x2+q4y2)1/2
2B [1+ez2(ax2+by2)]1/2 [(ax2+by2)2ezλ23]𝔍2B axez2λ23(ax2+by2)2 byez2λ23(ax2+by2)2 ezλ23(ax2+by2) axz𝔍2B byz𝔍2B (ax2+by2)𝔍2B
3 [q2yxq2]1/(q21) (q21)xyλ3(x2+q4y2)1/2 (q2q21)λ3x (q2q21)λ3q2y 0 q2y(x2+q4y2)1/2 x(x2+q4y2)1/2 0

3D

(x2+q4y2+p4z2)1/2,

𝔍02(x2+q4y2)3D2

=

[(xp2z)2+(q2yp2z)2+(x2+q4y2)2],

𝔍2B

[(axz)2+(byz)2+(ax2+by2)2]1/2.


This table titled, "Direction Cosine Components for T9 Coordinates," contains the following information:

  1. The (first) row labeled, n=1, correctly details the scale-factor, h1, and the unit vector expression, e^1=(ı^γ11+ȷ^γ12+k^γ13), that result from the given specification of the λ1 coordinate. By design, the unit vector, e^1, is everywhere normal to the "surface" of the ellipsoid.
  2. The (fourth) row labeled, n=3, correctly details the scale-factor, h3, and the unit vector expression, e^3=(ı^γ31+ȷ^γ32+k^γ33), that result from the given specification of the λ2 coordinate. By design, this unit vector, e^3, has no vertical component — that is, γ33=0 — and, by design, it is everywhere perpendicular to the "surface-normal" unit vector, e^1.
  3. We desire a unit vector, e^2, that is mutually orthogonal to the other two unit vectors; this has been accomplished by examining their cross-product, namely, we have set e^2A=e^3×e^1. Determined in this manner, the expressions for the three direction-cosine components of e^2A have been written in the last three columns of the (second) row labeled, n=2A. While we are confident that the correct specification of e^2 is,

    e^2A=(ı^γ21+ȷ^γ22+k^γ23)

    =

    ı^(xp2z)𝔍0+ȷ^(q2yp2z)𝔍0k^(x2+q4y2)𝔍0,

    as yet (18 February 2021), we have been unable to determine an expression for the coordinate, λ2A(x,y,z), from which all three of these direction-cosine expressions can be simultaneously derived — hence, the dashes in the second column of row 2A. The expression for h2A that has been presented in the third column of row 2A (and framed in pink) is a guess which, when divided into each of the three direction cosines, gives respectively the three guessed partial-derivative expressions shown in columns 4, 5, and 6 of row 2A.

  4. The second column of row 2B contains a guess (framed in yellow) for the coordinate expression, λ2B; this expression contains three unspecified scalar coefficients, a, b and e. The remaining columns of this row contain the three partial derivatives, the associated scale factor, and the three direction cosines that result from this guessed coordinate expression. If we can find values of the three scalar coefficients that give (row 2B) expressions for the three direction cosines that perfectly match the direction cosines written in row 2A, then we will be able to state that λ2B is — at least one form of — our sought-after third coordinate expression.

Yellow-Framed Guess 2B

Referencing the 2B table row, above, we are looking for coefficient values that map,

axzxp2z,

     

byzq2yp2z,

     

(ax2+by2)(x2+q4y2),

and that map,

𝔍2B[(axz)2+(byz)2+(ax2+by2)2]1/2

=

[z2(a2x2+b2y2)+(ax2+by2)2]1/2,

into the expression,

𝔍03D(x2+q4y2)1/2

=

[(x2+q4y2+p4z2)(x2+q4y2)]1/2

 

=

[p4z2(x2+q4y2)+(x2+q4y2)2]1/2

 

=

[x2p4z2+q4y2p4z2+(x2+q4y2)2]1/2.

A portion of these mappings are accomplished by setting a=p2 and b=q2p2, but this pair of specified coefficient values does not satisfy other mappings. Alternatively, a separate subset of mappings — but, again, not all mappings — is satisfied by setting a=1 and b=q4. So the yellow-framed guess 2B does not provide a correct second-coordinate expression.

Blue-Framed Guess 2C

Let's try,

λ2B

=

(cz)2e(ax2+by2)f.


Direction Cosine Components for Additional T9 Coordinate Guesses
n
 
(1)
λn
 
(2)
hn
 
(3)
λn/x
 
(4)
λn/y
 
(5)
λn/z
 
(6)
γn1
 
(7)
γn2
 
(8)
γn3
 
(9)
2A --- --- --- --- --- xp2z𝔍0 q2yp2z𝔍0 (x2+q4y2)𝔍0
2B [1+ez2(ax2+by2)]1/2 [(ax2+by2)2ezλ23]𝔍2B axez2λ23(ax2+by2)2 byez2λ23(ax2+by2)2 ezλ23(ax2+by2) axz𝔍2B byz𝔍2B (ax2+by2)𝔍2B
2C (cz)2e(ax2+by2)f [z(ax2+by2)2eλ2C]𝔍2C 2afxλ2C(ax2+by2) 2bfyλ2C(ax2+by2) 2eλ2Cz (afxze)𝔍2C (bfyze)𝔍2C (ax2+by2)𝔍2C

3D

(x2+q4y2+p4z2)1/2,

𝔍02(x2+q4y2)3D2

=

[(xp2z)2+(q2yp2z)2+(x2+q4y2)2],

𝔍2B

[(axz)2+(byz)2+(ax2+by2)2]1/2,

𝔍2C

[(xz)2(fae)2+(yz)2(fbe)2+(ax2+by2)2]1/2.


Examining just the expression for 𝔍2C, we see that we definitely need: a=1 and b=q4. Also, we need,

afe

=

p2

       

fe

=

p2;

and, separately we need,

fbe

=

q2p2

       

fe

=

p2q2.

These cannot simultaneously be satisfied. So the blue-framed guess 2C does not provide a correct second-coordinate expression. But we are very close! We need one additional scalar coefficient degree of freedom.

Next Thought

More generally, this leads to an expression for the scale-factor of the form,

h22

𝔊(x,y,z)[x2p4z2+q4y2p4z2+(x2+q4y2)2]

λ2x

[𝔊(x,y,z)]1/2[xp2z];

     

λ2y

[𝔊(x,y,z)]1/2[q2yp2z];

     

λ2z

[𝔊(x,y,z)]1/2[x2+q4y2].



Now, if we set, a=p2 and b=q2p2, we have,

(axz)2+(byz)2+(ax2+by2)2

=

(xp2z)2+(q2p2yz)2+(p2x2+q2p2y2)2

 

=

p4[q4y2z2+x2z2+(x2+q2y2)2],

in which case,

Complete Orthogonality Check

If the set of unit vectors is indeed orthogonal, then we must find that,

γmn

=

Mmn,

where the quantity Mmn is the minor of γmn in the determinant, |γmn|. (Note: This last expression is true only for right-handed coordinate systems. If the coordinate system is left-handed, we should find, γmn=Mmn.) More specifically, for any right-handed, orthogonal curvilinear coordinate system we should find:

γ11

=

γ22γ33γ23γ32,

γ12

=

γ23γ31γ21γ33,

γ13

=

γ21γ32γ22γ31,

   

γ21

=

γ32γ13γ33γ12,

γ22

=

γ33γ11γ31γ13,

γ23

=

γ31γ12γ32γ11,

   

γ31

=

γ12γ23γ13γ22,

γ32

=

γ13γ21γ11γ23,

γ33

=

γ11γ22γ12γ21.

and the position vector is,

x

=

ı^x+ȷ^y+k^z

 

=

e^1(γ11x+γ12y+γ13z)+e^2(γ21x+γ22y+γ23z)+e^3(γ31x+γ32y+γ33z).

For (κ1, κ2, κ3)

γ22γ33γ23γ32

=

[q2yp2z3D(x2+q4y2)1/2][0]+[(x2+q4y2)3D(x2+q4y2)1/2][x(x2+q4y2)1/2]=x3D=γ11.       Yes!

γ23γ31γ21γ33

=

[(x2+q4y2)1/23D][q2y(x2+q4y2)1/2][xp2z3D(x2+q4y2)1/2][0]=q2y3D=γ12.       Yes!

γ21γ32γ22γ31

=

[xp2z3D(x2+q4y2)1/2][x(x2+q4y2)1/2]+[q2yp2z3D(x2+q4y2)1/2][q2y(x2+q4y2)1/2]

 

=

[x2p2z3D(x2+q4y2)]+[q4y2p2z3D(x2+q4y2)]=p2z3D=γ13.       Yes!

γ32γ13γ33γ12

=

[x(x2+q4y2)1/2][p2z3D][0][q2y3D]=[xp2z3D(x2+q4y2)1/2]=γ21.       Yes!

γ33γ11γ31γ13

=

[0][x3D]+[q2y(x2+q4y2)1/2][p2z3D]=[q2yp2z3D(x2+q4y2)1/2]=γ22.       Yes!

γ31γ12γ32γ11

=

[q2y(x2+q4y2)1/2][q2y3D][x(x2+q4y2)1/2][x3D]=(x2+q4y2)1/23D=γ23.       Yes!

γ12γ23γ13γ22

=

[q2y3D][(x2+q4y2)1/23D][p2z3D][q2yp2z3D(x2+q4y2)1/2]

 

=

q2y3D2(x2+q4y2)1/2[x2+q4y2+p4z2]=q2y(x2+q4y2)1/2=γ31.       Yes!

γ13γ21γ11γ23

=

[p2z3D][xp2z3D(x2+q4y2)1/2]+[x3D][(x2+q4y2)1/23D]

 

=

x3D2(x2+q4y2)1/2[x2+q4y2+p4z2]=x(x2+q4y2)1/2=γ32.       Yes!

γ11γ22γ12γ21

=

[x3D][q2yp2z3D(x2+q4y2)1/2][q2y3D][xp2z3D(x2+q4y2)1/2]=0=γ33.       Yes!

Given that the prescribed interrelationships between all nine direction cosines are satisfied, we conclude that the (κ1,κ2,κ3) coordinate system is an orthogonal one. Accordingly, the position vector is,

x

=

e^1{x23D+q2y23D+p2z23D}

 

 

+e^2{x2+q2y21p2(x2+q4y2)}p2z3D(x2+q4y2)1/2

 

 

+e^3{xq2y(x2+q4y2)1/2+xy(x2+q4y2)1/2}

 

=

e^1[κ123D]+e^2[x2+q2y21p2(x2+q4y2)]p2z3D(x2+q4y2)1/2+e^3[(1q2)xy(x2+q4y2)1/2]

 

=

e^1[κ123D]+e^2[x2(11p2)q2y2(q21)]p2z3D(x2+q4y2)1/2e^3[(q21)xy(x2+q4y2)1/2].

For (κ1, κ4, κ5)

γ42γ53γ43γ52

=

[xq2yp2z𝒟(1q2y)][p2z3D𝒟(x2q4y2)]+[xq2yp2z𝒟(2p2z)][q2y3D𝒟(p4z2+2x2)]

 

=

x3D𝒟2[(x2q4y2)p4z2+(p4z2+2x2)2q4y2]

 

=

x3D𝒟2[x2p4z2+4x2q4y2+q4y2p4z2]=x3D=γ11.       Yes!

γ43γ51γ41γ53

=

[xq2yp2z𝒟(2p2z)][x3D𝒟(2q4y2+p4z2)][xq2yp2z𝒟(1x)][p2z3D𝒟(x2q4y2)]

 

=

q2y3D𝒟2[2x2(2q4y2+p4z2)p4z2(x2q4y2)]

 

=

q2y3D𝒟2[4x2q4y2+x2p4z2+q4y2p4z2]=q2y3D=γ12.       Yes!

γ41γ52γ42γ51

=

[xq2yp2z𝒟(1x)][q2y3D𝒟(p4z2+2x2)]+[xq2yp2z𝒟(1q2y)][x3D𝒟(2q4y2+p4z2)]

 

=

p2z3D𝒟2[q4y2(p4z2+2x2)+x2(2q4y2+p4z2)]

 

=

p2z3D𝒟2[q4y2p4z2+4x2q4y2+x2p4z2]=p2z3D=γ13.       Yes!

γ52γ13γ53γ12

=

[q2yp2z3D2𝒟(p4z2+2x2)][q2yp2z3D2𝒟(x2q4y2)]

 

=

q2yp2z3D2𝒟[x2+q4y2+p4z2]=q2yp2z𝒟=γ41.       Yes!

γ53γ11γ51γ13

=

[xp2z3D2𝒟(x2q4y2)]+[xp2z3D2𝒟(2q4y2+p4z2)]

 

=

xp2z3D2𝒟[x2+q4y2+p4z2]=xp2z𝒟=γ42.       Yes!

γ51γ12γ52γ11

=

[xq2y3D2𝒟(2q4y2+p4z2)][xq2y3D2𝒟(p4z2+2x2)]

 

=

2xq2y3D2𝒟[x2+q4y2+p4z2]=2xq2y𝒟=γ43.       Yes!

γ12γ43γ13γ42

=

[2xq4y23D𝒟][xp4z23D𝒟]=[x3D𝒟](2q4y2+p4z2)=γ51.       Yes!

γ13γ41γ11γ43

=

[q2yp4z23D𝒟]+[2x2q2y3D𝒟]=[q2y3D𝒟](p4z2+2x2)=γ52.       Yes!

γ11γ42γ12γ41

=

[x2p2z3D𝒟][q4y2p2z3D𝒟]=[p2z3D𝒟](x2q4y2)=γ53.       Yes!

Given that the prescribed interrelationships between all nine direction cosines are satisfied, we conclude that the (κ1,κ4,κ5) coordinate system is an orthogonal one. Accordingly, the position vector is,

x

=

e^1(γ11x+γ12y+γ13z)+e^2(γ41x+γ42y+γ43z)+e^3(γ51x+γ52y+γ53z)

 

=

e^1{x23D+q2y23D+p2z23D}

 

 

+e^2{xq2yp2z𝒟+xyp2z𝒟2xq2yz𝒟}

 

 

+e^3{3Dx2𝒟(2q4y2+p4z2)+3Dq2y2𝒟(p4z2+2x2)+3Dp2z2𝒟(x2q4y2)}

 

=

e^1{κ123D}+e^2{q2p2+p22q2}xyz𝒟

 

 

+e^3{x2(2q4y2+p4z2)+q2y2(p4z2+2x2)+p2z2(x2q4y2)}3D𝒟

 

=

e^1{κ123D}+e^2{q2p2+p22q2}xyz𝒟e^3{2x2q2y2(q21)+x2p2z2(p21)+q2y2p2z2(q2p2)}3D𝒟.

See Also


Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |