Appendix/Ramblings/T6CoordinatesPt3: Difference between revisions

From jetwiki
Jump to navigation Jump to search
Created page with "__FORCETOC__ <!-- __NOTOC__ will force TOC off --> =Concentric Ellipsoidal (T6) Coordinates= ==Background== Building on our Appendix/Ramblings/DirectionCosines|general intr..."
 
No edit summary
 
Line 1: Line 1:
__FORCETOC__
__FORCETOC__
<!-- __NOTOC__ will force TOC off -->
<!-- __NOTOC__ will force TOC off -->
=Concentric Ellipsoidal (T6) Coordinates=
=Concentric Ellipsoidal (T6) Coordinates (Part 3)=
 
==Background==
Building on our [[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 [[Appendix/Ramblings/EllipticCylinderCoordinates#T5_Coordinates|T5]] (concentric elliptic) coordinate systems, here we explore the creation of a concentric ellipsoidal (T6) coordinate system.  This is motivated by our [[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]].
 
==Orthogonal Coordinates==
 
===Primary (''radial-like'') Coordinate===
 
We start by defining a "radial" coordinate whose values identify various concentric ellipsoidal shells,
<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^2 y^2 + p^2 z^2)^{1 / 2} \, .</math>
  </td>
</tr>
</table>
When <math>~\lambda_1 = a</math>, we obtain the standard definition of an ellipsoidal surface, it being understood that, <math>~q^2 = a^2/b^2</math> and <math>~p^2 = a^2/c^2</math>.  (We will assume that <math>~a > b > c</math>, that is, <math>~p^2 > q^2 > 1</math>.) 
 
A vector, <math>~\bold{\hat{n}}</math>, that is normal to the <math>~\lambda_1</math> = constant surface is given by the gradient of the function,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~F(x, y, z)</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~(x^2 + q^2 y^2 + p^2 z^2)^{1 / 2} - \lambda_1 \, .</math>
  </td>
</tr>
</table>
 
In Cartesian coordinates, this means,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\bold{\hat{n}}(x, y, z)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl( \frac{\partial F}{\partial x} \biggr)
+ \hat\jmath \biggl( \frac{\partial F}{\partial y} \biggr)
+ \hat{k} \biggl( \frac{\partial F}{\partial z} \biggr)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl[ x(x^2 + q^2 y^2 + p^2 z^2)^{- 1 / 2} \biggr]
+ \hat\jmath \biggl[ q^2y(x^2 + q^2 y^2 + p^2 z^2)^{- 1 / 2} \biggr]
+ \hat{k}\biggl[ p^2 z(x^2 + q^2 y^2 + p^2 z^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{x}{\lambda_1} \biggr)
+ \hat\jmath \biggl( \frac{q^2y}{\lambda_1} \biggr)
+ \hat{k}\biggl(\frac{p^2 z}{\lambda_1} \biggr) \, ,
</math>
  </td>
</tr>
</table>
where it is understood that this expression is only to be evaluated at points, <math>~(x, y, z)</math>, that lie on the selected <math>~\lambda_1</math> surface &#8212; that is, at points for which the function, <math>~F(x,y,z) = 0</math>.  The length of this normal vector is given by the expression,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~[ \bold{\hat{n}} \cdot \bold{\hat{n}} ]^{1 / 2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \biggl( \frac{\partial F}{\partial x} \biggr)^2 + \biggl( \frac{\partial F}{\partial y} \biggr)^2 + \biggl( \frac{\partial F}{\partial z} \biggr)^2 \biggr]^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \biggl( \frac{x}{\lambda_1} \biggr)^2
+ \biggl( \frac{q^2y}{\lambda_1} \biggr)^2
+ \biggl(\frac{p^2 z}{\lambda_1} \biggr)^2 \biggr]^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{\lambda_1 \ell_{3D}}
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\ell_{3D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl[ x^2 + q^4y^2 + p^4 z^2 \biggr]^{- 1 / 2} \, .</math>
  </td>
</tr>
</table>
 
It is therefore clear that the ''properly normalized'' normal unit vector that should be associated with any <math>~\lambda_1</math> = constant ellipsoidal surface is,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_1 </math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{ \bold\hat{n} }{ [ \bold{\hat{n}} \cdot \bold{\hat{n}} ]^{1 / 2} }
=
\hat\imath (x \ell_{3D}) + \hat\jmath (q^2y \ell_{3D}) + \hat{k} (p^2 z \ell_{3D}) \, .
</math>
  </td>
</tr>
</table>
From our [[Appendix/Ramblings/DirectionCosines#Scale_Factors|accompanying discussion of direction cosines]], it is clear, as well, that the scale factor associated with the <math>~\lambda_1</math> coordinate is,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~h_1^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\lambda_1^2 \ell_{3D}^2 \, .</math>
  </td>
</tr>
</table>
We can also fill in the top line of our direction-cosines table, namely,
 
<table border="1" cellpadding="8" align="center" width="60%">
<tr>
  <td align="center" colspan="4">
'''Direction Cosines for T6 Coordinates'''
<br />
<math>~\gamma_{ni} = h_n \biggl( \frac{\partial \lambda_n}{\partial x_i}\biggr)</math>
  </td>
</tr>
<tr>
  <td align="center" width="10%"><math>~n</math></td>
  <td align="center" colspan="3"><math>~i = x, y, z</math>
</tr>
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center">&nbsp;<br />
<math>~x\ell_{3D}</math><br />&nbsp;
  <td align="center"><math>~q^2 y \ell_{3D}</math>
  <td align="center"><math>~p^2 z \ell_{3D}</math>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  </td>
</tr>
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  </td>
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  </td>
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  </td>
</tr>
</table>
 
===Other Coordinate Pair in the Tangent Plane===
 
Let's focus on a particular point on the <math>\lambda_1</math> = constant surface, <math>(x_0, y_0, z_0)</math>, that necessarily satisfies the function, <math>F(x_0, y_0, z_0) = 0</math>.  We have already derived the expression for the unit vector that is normal to the ellipsoidal surface at this point, namely,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>\hat{e}_1 </math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>
\hat\imath (x_0 \ell_{3D}) + \hat\jmath (q^2y_0 \ell_{3D}) + \hat\jmath (p^2 z_0 \ell_{3D}) \, ,
</math>
  </td>
</tr>
</table>
where, for this specific point on the surface,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>\ell_{3D}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ x_0^2 + q^4y_0^2 + p^4 z_0^2 \biggr]^{- 1 / 2} \, .</math>
  </td>
</tr>
</table>
 
 
<table border="1" align="center" width="80%" cellpadding="10"><tr><td align="left">
<div align="center">
'''Tangent Plane'''<br />
[See, for example, [http://math.furman.edu/~dcs/courses/math21/ Dan Sloughter's] ([https://www.furman.edu Furman University]) 2001 Calculus III  class lecture notes &#8212; specifically [http://math.furman.edu/~dcs/courses/math21/lectures/l-15.pdf Lecture 15]]
</div>
 
----
 
The two-dimensional plane that is tangent to the <math>\lambda_1</math> = constant surface ''at this point'' is given by the expression,
<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 - x_0) \biggl[ \frac{\partial \lambda_1}{\partial x} \biggr]_0
+ (y - y_0) \biggl[\frac{\partial \lambda_1}{\partial y} \biggr]_0 
+ (z - z_0) \biggl[\frac{\partial \lambda_1}{\partial z} \biggr]_0
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
(x - x_0) \biggl[ \frac{\partial F}{\partial x} \biggr]_0 
+ (y - y_0) \biggl[\frac{\partial F}{\partial y} \biggr]_0 
+ (z - z_0) \biggl[ \frac{\partial F}{\partial z} \biggr]_0
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
(x - x_0) \biggl( \frac{x}{\lambda_1}\biggr)_0 + (y - y_0)\biggl( \frac{q^2 y  }{ \lambda_1 } \biggr)_0 + (z - z_0)\biggl( \frac{ p^2z }{ \lambda_1 } \biggr)_0
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>\Rightarrow~~~
x \biggl( \frac{x}{\lambda_1}\biggr)_0 + y \biggl( \frac{q^2 y  }{ \lambda_1 } \biggr)_0 + z \biggl( \frac{ p^2z }{ \lambda_1 } \biggr)_0
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
x_0 \biggl( \frac{x}{\lambda_1}\biggr)_0 + y_0 \biggl( \frac{q^2 y  }{ \lambda_1 } \biggr)_0 + z_0 \biggl( \frac{ p^2z }{ \lambda_1 } \biggr)_0
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>\Rightarrow~~~
x x_0  + q^2 y y_0  + p^2 z z_0
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
x_0^2  +  q^2 y_0^2 + p^2 z_0^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>\Rightarrow~~~
x x_0  + q^2 y y_0  + p^2 z z_0
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
(\lambda_1^2)_0 \, .
</math>
  </td>
</tr>
</table>
</td></tr></table>
 
Fix the value of <math>\lambda_1</math>.  This means that the relevant ellipsoidal surface is defined by the expression,
<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^2y^2 + p^2z^2 \, .</math>
  </td>
</tr>
</table>
If <math>z = 0</math>, the semi-major axis of the relevant x-y ellipse is <math>\lambda_1</math>, and the square of the semi-minor axis is <math>\lambda_1^2/q^2</math>.  At any other value, <math>z = z_0 < c</math>, the square of the semi-major axis of the relevant x-y ellipse is, <math>~(\lambda_1^2 - p^2z_0^2)</math> and the square of the corresponding semi-minor axis is, <math>(\lambda_1^2 - p^2z_0^2)/q^2</math>.  Now, for any chosen <math>x_0^2 \le (\lambda_1^2 - p^2z_0^2)</math>, the y-coordinate of the point on the <math>~\lambda_1</math> surface is given by the expression,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>y_0^2</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{1}{q^2}\biggl[ \lambda_1^2 - p^2 z_0 -x_0^2 \biggr] \, .</math>
  </td>
</tr>
</table>
The slope of the line that lies in the <math>z = z_0</math> plane and that is tangent to the ellipsoidal surface at <math>(x_0, y_0)</math> is,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>m \equiv \frac{dy}{dx}\biggr|_{z_0}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>- \frac{x_0}{q^2y_0}</math>
  </td>
</tr>
</table>
 
===Speculation1===
 
Building on our experience developing [[User:Tohline/Appendix/Ramblings/T3Integrals#Integrals_of_Motion_in_T3_Coordinates|T3 Coordinates]] and, more recently, [[User:Tohline/Appendix/Ramblings/EllipticCylinderCoordinates#T5_Coordinates|T5 Coordinates]], let's define the two "angles,"
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Zeta</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\sinh^{-1}\biggl(\frac{qy}{x} \biggr)</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp; </td>
  <td align="right">
<math>~\Upsilon</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\sinh^{-1}\biggl(\frac{pz}{x} \biggr) \, ,</math>
  </td>
</tr>
</table>
in which case we can write,
<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(\cosh^2\Zeta + \sinh^2\Upsilon)\, .</math>
  </td>
</tr>
</table>
We speculate that the other two orthogonal coordinates may be defined by the expressions,
<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>~x \biggl[ \sinh\Zeta \biggr]^{1/(1-q^2)}
=
x \biggl[ \frac{qy}{x}\biggr]^{1/(1-q^2)}
=
x \biggl[ \frac{x}{qy}\biggr]^{1/(q^2-1)}
=
\biggl[ \frac{x^{q^2}}{qy}\biggr]^{1/(q^2-1)}
\, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~x \biggl[ \sinh\Upsilon \biggr]^{1/(1-p^2)}
=
x \biggl[ \frac{pz}{x}\biggr]^{1/(1-p^2)}
=
x \biggl[ \frac{x}{pz}\biggr]^{1/(p^2-1)}
=
\biggl[ \frac{x^{p^2}}{pz}\biggr]^{1/(p^2-1)}
\, .</math>
  </td>
</tr>
</table>
 
Some relevant partial derivatives are,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{1}{qy}\biggr]^{1/(q^2-1)} \biggl[ \frac{q^2}{q^2-1} \biggr]x^{1/(q^2-1)}
=
\biggl[ \frac{q^2}{q^2-1} \biggr]\biggl[ \frac{x}{qy}\biggr]^{1/(q^2-1)}
=
\biggl[ \frac{q^2}{q^2-1} \biggr]\frac{\lambda_2}{x}
\, ;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{x^{q^2}}{q}\biggr]^{1/(q^2-1)} \biggl[ \frac{1}{1-q^2} \biggr] y^{q^2/(1-q^2)}
=
- \biggl[ \frac{1}{q^2-1} \biggr] \frac{\lambda_2}{y}
\, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{p^2}{p^2-1} \biggr]\frac{\lambda_3}{x}
\, ;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_3}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \biggl[ \frac{1}{p^2-1} \biggr] \frac{\lambda_3}{z}
\, .</math>
  </td>
</tr>
</table>
And the associated scale factors are,
<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>~
\biggl\{ \biggl[ \biggl( \frac{q^2}{q^2-1} \biggr)\frac{\lambda_2}{x} \biggr]^2 + \biggl[ - \biggl( \frac{1}{q^2-1} \biggr) \frac{\lambda_2}{y} \biggr]^2 \biggr\}^{-1}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl( \frac{q^2}{q^2-1} \biggr)^2 \frac{\lambda_2^2}{x^2}  + \biggl( \frac{1}{q^2-1} \biggr)^2 \frac{\lambda_2^2}{y^2} \biggr\}^{-1}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{x^2 + q^4 y^2 \biggr\}^{-1}
\biggl[ \frac{(q^2 - 1)^2x^2 y^2}{\lambda_2^2} \biggr] \, ;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~h_3^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{x^2 + p^4 z^2 \biggr\}^{-1}
\biggl[ \frac{(p^2 - 1)^2x^2 z^2}{\lambda_3^2} \biggr] \, .
</math>
  </td>
</tr>
</table>
We can now fill in the rest of our direction-cosines table, namely,
 
<table border="1" cellpadding="8" align="center" width="60%">
<tr>
  <td align="center" colspan="4">
'''Direction Cosines for T6 Coordinates'''
<br />
<math>~\gamma_{ni} = h_n \biggl( \frac{\partial \lambda_n}{\partial x_i}\biggr)</math>
  </td>
</tr>
<tr>
  <td align="center" width="10%"><math>~n</math></td>
  <td align="center" colspan="3"><math>~i = x, y, z</math>
</tr>
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center">&nbsp;<br />
<math>~x\ell_{3D}</math><br />&nbsp;
  <td align="center"><math>~q^2 y \ell_{3D}</math>
  <td align="center"><math>~p^2 z \ell_{3D}</math>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center">
<math>~q^2 y \ell_q </math>
  <td align="center">
<math>~-x\ell_q</math>
  <td align="center">
<math>~0</math>
  </td>
</tr>
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center">
<math>~p^2 z \ell_p</math>
  </td>
  <td align="center">
<math>~0</math>
  </td>
  <td align="center">
<math>~-x\ell_p</math>
  </td>
</tr>
</table>
Hence,
<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 \gamma_{21}
+ \hat\jmath \gamma_{22}
+\hat{k} \gamma_{23}
=
\hat\imath (q^2y\ell_q)
- \hat\jmath (x\ell_q) \, ;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\hat{e}_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \gamma_{31}
+ \hat\jmath \gamma_{32}
+\hat{k} \gamma_{33}
=
\hat\imath (p^2z\ell_p)
-\hat{k} (x\ell_p) \, .
</math>
  </td>
</tr>
</table>
 
Check:
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_2 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(q^2y\ell_q)^2
+ (x\ell_q)^2
=
1 \, ;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\hat{e}_3 \cdot \hat{e}_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(p^2z\ell_p)^2
+ (x\ell_p)^2
=
1 \, ;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\hat{e}_2 \cdot \hat{e}_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(q^2y\ell_q)(p^2z\ell_p) \ne 0 \, .
</math>
  </td>
</tr>
</table>
 
===Speculation2===
 
Try,
<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} z^{1/p^2}} \, ,
</math>
  </td>
</tr>
</table>
in which case,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\lambda_2}{x} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{x}{z^{1/p^2}} \biggl(-\frac{1}{q^2}\biggr) y^{-1/q^2 - 1}
=
-\frac{\lambda_2}{q^2 y}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{\lambda_2}{p^2 z}
\, .
</math>
  </td>
</tr>
</table>
The associated scale factor is, then,
<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>~\biggl[
\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
\biggr]^{-1}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[
\biggl( \frac{ \lambda_2}{x} \biggr)^2
+
\biggl( -\frac{\lambda_2}{q^2y} \biggr)^2
+
\biggl( - \frac{\lambda_2}{p^2z} \biggr)^2
\biggr]^{-1}
</math>
  </td>
</tr>
</table>
 
===Speculation3===
 
Try,
<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+p^2 z)^{1 / 2}}{y^{1/q^2} } \, ,
</math>
  </td>
</tr>
</table>
in which case,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{2y^{1/q^2}}\biggl(x + p^2z\biggr)^{- 1 / 2}
=
\frac{\lambda_2}{2(x + p^2z) }
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{\lambda_2}{q^2y}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
 
\, .
</math>
  </td>
</tr>
</table>
 
===Speculation4===
 
====Development====
 
Here we stick with the [[#Primary_.28radial-like.29_Coordinate|primary (radial-like) coordinate as defined above]]; for example,
<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^2 y^2 + p^2 z^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>~\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>
 
<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^2y \ell_{3D}) + \hat{k} (p^2 z \ell_{3D}) \, .
</math>
  </td>
</tr>
</table>
 
<table border="1" width="80%" align="center" cellpadding="10"><tr><td align="left">
Note that, <math>~\hat{e}_1 \cdot \hat{e}_1 = 1</math>, which means that this is, indeed, a properly normalized ''unit'' vector.
</td></tr></table>
 
Then, drawing from our [https://www.phys.lsu.edu/astro/H_Book.current/Appendices/Mathematics/operators.tohline1.pdf earliest discussions of "T1 Coordinates"], we'll try defining the ''second'' coordinate as,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\tan^{-1} u \, ,
</math>
&nbsp; &nbsp; &nbsp; where,
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~u</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>\frac{y^{1/q^2}}{x} \, .</math>
  </td>
</tr>
</table>
The relevant partial derivatives are,
<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{1}{1 + u^2} \biggl[ - \frac{y^{1/q^2}}{x^2} \biggr]
=
- \biggl[ \frac{u}{1 + u^2}\biggr]\frac{1}{x}
=
- \frac{\sin\lambda_3 \cos\lambda_3}{x}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_3}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{1 + u^2} \biggl[ \frac{y^{(1/q^2-1)}}{q^2x} \biggr]
=
\biggl[ \frac{u}{1 + u^2}\biggr]\frac{1}{q^2y}
=
\frac{\sin\lambda_3 \cos\lambda_3}{q^2y} \, ,
</math>
  </td>
</tr>
</table>
which 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>~
\biggl[ \biggl( \frac{\partial \lambda_3}{\partial x} \biggr)^2 + \biggl( \frac{\partial \lambda_3}{\partial y} \biggr)^2 \biggr]^{-1}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{u}{1 + u^2}\biggr]^{-2} \biggl[ \frac{1}{x^2} + \frac{1}{q^4y^2} \biggr]^{-1}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{1 + u^2}{u}\biggr]^{2} \biggl[ \frac{x^2 + q^4y^2}{x^2q^4y^2} \biggr]^{-1}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow~~~h_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{1 + u^2}{u}\biggr]xq^2 y \ell_q = \frac{xq^2 y \ell_q}{\sin\lambda_3 \cos\lambda_3} \, ,
</math>
&nbsp; &nbsp; &nbsp; where,
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\ell_q</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~[x^2 + q^4 y^2]^{-1 / 2} \, .</math>
  </td>
</tr>
</table>
 
The third row of direction cosines can now be filled in to give,
 
<table border="1" cellpadding="8" align="center" width="60%">
<tr>
  <td align="center" colspan="4">
'''Direction Cosines for T6 Coordinates'''
<br />
<math>~\gamma_{ni} = h_n \biggl( \frac{\partial \lambda_n}{\partial x_i}\biggr)</math>
  </td>
</tr>
<tr>
  <td align="center" width="10%"><math>~n</math></td>
  <td align="center" colspan="3"><math>~i = x, y, z</math>
</tr>
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center">&nbsp;<br />
<math>~x\ell_{3D}</math><br />&nbsp;
  <td align="center"><math>~q^2 y \ell_{3D}</math>
  <td align="center"><math>~p^2 z \ell_{3D}</math>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  <td align="center">
&nbsp;<br />
---
<br />&nbsp;
  </td>
</tr>
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center">
<math>~-q^2 y \ell_q</math>
  </td>
  <td align="center">
<math>~x \ell_q</math>
  </td>
  <td align="center">
<math>~0</math>
  </td>
</tr>
</table>
which means that the associated unit vector is,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_3 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\hat\imath (q^2 y \ell_{q}) + \hat\jmath (x \ell_{q})  \, .
</math>
  </td>
</tr>
</table>
 
<table border="1" width="80%" align="center" cellpadding="10"><tr><td align="left">
Note that, <math>~\hat{e}_3 \cdot \hat{e}_3 = 1</math>, which means that this also is a properly normalized ''unit'' vector.  Note, as well, that the dot product between our "first" and "third" unit vectors should be zero if they are indeed orthogonal to each other.  Let's see &hellip;
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_3 \cdot \hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(- q^2y \ell_q)x\ell_{3D} + (x\ell_q) q^2y\ell_{3D} = 0 \, .
</math>
  </td>
</tr>
</table>
Q.E.D.
</td></tr></table>
 
Now, even though we have not yet determined the proper expression for the "second" orthogonal coordinate, <math>~\lambda_2</math>, we should be able to obtain an expression for its associated unit vector from the cross product of the "third" and "first" unit vectors.  Specifically we find,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_2 \equiv \hat{e}_3 \times \hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl[ (e_3)_2 (e_1)_3 - (e_3)_3(e_1)_2 \biggr]
+
\hat\jmath \biggl[ (e_3)_3 (e_1)_1 - (e_3)_1(e_1)_3 \biggr]
+
\hat{k} \biggl[ (e_3)_1 (e_1)_2 - (e_3)_2(e_1)_1 \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat\imath \biggl[ (x \ell_q) (p^2 z \ell_{3D}) - 0 \biggr]
+
\hat\jmath \biggl[ 0 - (-q^2y \ell_q)(p^2z \ell_{3D}) \biggr]
+
\hat{k} \biggl[ (-q^2y \ell_q) (q^2 y \ell_{3D}) - (x\ell_q)(x\ell_{3D}) \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\ell_q \ell_{3D}\biggl[
\hat\imath ( xp^2 z ) + \hat\jmath ( q^2y p^2z ) - \hat{k} ( x^2 + q^4 y^2 )
\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\ell_q \ell_{3D}\biggl[
\hat\imath ( xp^2 z ) + \hat\jmath ( q^2y p^2z ) - \hat{k} \biggl( \frac{1}{\ell_q^2}  \biggr)
\biggr] \, .
</math>
  </td>
</tr>
</table>
 
<table border="1" width="80%" align="center" cellpadding="10"><tr><td align="left">
Note that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_3 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\ell_q^2 \ell_{3D} \biggl[
(- q^2y )x p^2 z 
+ (x) q^2y p^2 z 
\biggr] = 0 \, ;
</math>
  </td>
</tr>
</table>
and,
<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})xp^2z \ell_q \ell_{3D}
+ (q^2y \ell_{3D}) q^2yp^2 z \ell_q \ell_{3D}
- (x^2 + q^4 y^2)\ell_q \ell_{3D} (p^2 z \ell_{3D} )
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \ell_q \ell_{3D}^2 \biggl[
x^2p^2z
+ (q^4y^2 ) p^2 z
- (x^2 + q^4 y^2) (p^2 z )
\biggr] = 0 \, .
</math>
  </td>
</tr>
</table>
We conclude, therefore, that <math>~\hat{e}_2</math> is perpendicular to both of the other unit vectors. <font color="red">'''Hooray!'''</font>
</td></tr></table>
 
 
Filling in the second row of the direction cosines table gives,
 
<table border="1" cellpadding="8" align="center" width="60%">
<tr>
  <td align="center" colspan="4">
'''Direction Cosines for T6 Coordinates'''
<br />
<math>~\gamma_{ni} = h_n \biggl( \frac{\partial \lambda_n}{\partial x_i}\biggr)</math>
  </td>
</tr>
<tr>
  <td align="center" width="10%"><math>~n</math></td>
  <td align="center" colspan="3"><math>~i = x, y, z</math>
</tr>
 
<tr>
  <td align="center"><math>~1</math></td>
  <td align="center">&nbsp;<br />
<math>~x\ell_{3D}</math><br />&nbsp;
  <td align="center"><math>~q^2 y \ell_{3D}</math>
  <td align="center"><math>~p^2 z \ell_{3D}</math>
</tr>
<tr>
  <td align="center"><math>~2</math></td>
  <td align="center">
<math>~x p^2 z\ell_q \ell_{3D}</math>
  <td align="center">
<math>~q^2y p^2 z\ell_q \ell_{3D}</math>
  <td align="center">
<math>~-(x^2 + q^4y^2)\ell_q \ell_{3D} = - \ell_{3D}/\ell_q</math>
  </td>
</tr>
 
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center">
<math>~-q^2 y \ell_q</math>
  </td>
  <td align="center">
<math>~x \ell_q</math>
  </td>
  <td align="center">
<math>~0</math>
  </td>
</tr>
</table>
 
====Analysis====
 
Let's break down each direction cosine into its components.
 
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" colspan="9">'''Direction Cosine Components for T6 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">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center"><math>~\ell_q \ell_{3D} (xp^2z)</math></td>
  <td align="center"><math>~\ell_q \ell_{3D} (q^2 y p^2z) </math></td>
  <td align="center"><math>~- (x^2 + q^4y^2)\ell_q \ell_{3D}</math></td>
</tr>
 
<tr>
  <td align="center"><math>~3</math></td>
  <td align="center"><math>~\tan^{-1}\biggl( \frac{y^{1/q^2}}{x} \biggr)</math></td>
  <td align="center"><math>~\frac{xq^2 y \ell_q}{\sin\lambda_3 \cos\lambda_3}</math></td>
  <td align="center"><math>~-\frac{\sin\lambda_3 \cos\lambda_3}{x}</math></td>
  <td align="center"><math>~+\frac{\sin\lambda_3 \cos\lambda_3}{q^2y}</math></td>
  <td align="center"><math>~0</math></td>
  <td align="center"><math>~-q^2 y \ell_q</math></td>
  <td align="center"><math>~x\ell_q</math></td>
  <td align="center"><math>~0</math></td>
</tr>
</table>
 
Try,
<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>~
\tan^{-1} w \, ,
</math>
&nbsp; &nbsp; &nbsp; where,
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~w</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>\frac{(x^2 + q^2y^2)^{1 / 2}}{z^{1/p^2}}
~~~\Rightarrow~~~\frac{1}{z^{1 / p^2} } = \frac{w}{(x^2 + q^2 y^2)^{1 / 2}}
\, .</math>
  </td>
</tr>
</table>
 
The relevant partial derivatives are,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{1 + w^2} \biggl[ \frac{x}{(x^2 + q^2y^2)^{1 / 2}~z^{1/p^2}} \biggr]
=
\frac{w}{1 + w^2} \biggl[ \frac{x}{(x^2 + q^2y^2)} \biggr]
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{1 + w^2} \biggl[ \frac{q^2y}{(x^2 + q^2y^2)^{1 / 2}~z^{1/p^2}} \biggr]
=
\frac{w}{1 + w^2} \biggl[ \frac{q^2y}{(x^2 + q^2y^2)} \biggr]
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{w}{1 + w^2} \biggl[- \frac{1}{p^2 z} \biggr]
\, ,
</math>
  </td>
</tr>
</table>
 
which means that,
<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>~
\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[ \biggl( \frac{w}{1 + w^2}\biggr)^2 \frac{x^2}{(x^2 + q^2y^2)^2}  \biggr]
+
\biggl[ \biggl( \frac{w}{1 + w^2}\biggr)^2 \frac{q^4 y^2}{(x^2 + q^2y^2)^2}  \biggr]
+
\biggl[ \biggl( \frac{w}{1 + w^2}\biggr)^2 \frac{1}{p^4 z^2} \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{w}{1 + w^2}\biggr)^2
\biggl[ \frac{(x^2 + q^4y^2)(p^4 z^2) + (x^2 + q^2y^2)^2}{(x^2 + q^2y^2)^2~p^4 z^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{1 + w^2}{w}\biggr) \biggl\{ \frac{(x^2 + q^2y^2)~p^2 z}{ \mathcal{D}}  \biggr\}
\, ,
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\mathcal{D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~[(x^2 + q^4y^2)(p^4 z^2) + (x^2 + q^2y^2)^2]^{1 / 2} \, .</math>
  </td>
</tr>
</table>
Hence, the trio of associated direction cosines are,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\gamma_{21} = h_2 \biggl( \frac{\partial \lambda_2}{\partial x} \biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{1 + w^2}{w}\biggr) \biggl\{ \frac{(x^2 + q^2y^2)~p^2 z}{ \mathcal{D}}  \biggr\}\frac{w}{1 + w^2} \biggl[ \frac{x}{(x^2 + q^2y^2)} \biggr]
=
\biggl\{ \frac{x~p^2 z}{ \mathcal{D}}  \biggr\} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\gamma_{22} = h_2 \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{1 + w^2}{w}\biggr) \biggl\{ \frac{(x^2 + q^2y^2)~p^2 z}{ \mathcal{D}}  \biggr\} \frac{w}{1 + w^2} \biggl[ \frac{q^2y}{(x^2 + q^2y^2)} \biggr]
=
\biggl\{ \frac{q^2 y~p^2 z}{ \mathcal{D}}  \biggr\} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\gamma_{23} = h_2 \biggl( \frac{\partial \lambda_2}{\partial z} \biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{1 + w^2}{w}\biggr) \biggl\{ \frac{(x^2 + q^2y^2)~p^2 z}{ \mathcal{D}}  \biggr\}\frac{w}{1 + w^2} \biggl[- \frac{1}{p^2 z} \biggr]
=
\biggl\{- \frac{(x^2 + q^2y^2)}{ \mathcal{D}}  \biggr\} \, .
</math>
  </td>
</tr>
</table>
 
<font color="red">'''VERY close!'''</font>
 
Let's examine the function, <math>~\mathcal{D}^2</math>.
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{1}{\ell_{3D}^2 \ell_d^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(x^2 + q^4 y^2)(x^2 + q^4 y + p^4 z)
=
(x^2 + q^4 y^2)p^4 z + (x^2 + q^4 y^2)^2 \, .
</math>
  </td>
</tr>
</table>
 
===Eureka (NOT!)===
 
 
Try,
<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>~
\tan^{-1} w \, ,
</math>
&nbsp; &nbsp; &nbsp; where,
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~w</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>\frac{(x^2 + q^2y^2)^{1 / 2}}{p^2 z}
~~~\Rightarrow~~~\frac{1}{p^2 z } = \frac{w}{(x^2 + q^2 y^2)^{1 / 2}}
\, .</math>
  </td>
</tr>
</table>
 
The relevant partial derivatives are,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{1 + w^2} \biggl[ \frac{x}{(x^2 + q^2y^2)^{1 / 2}~p^2 z} \biggr]
=
\frac{w}{1 + w^2} \biggl[ \frac{x}{(x^2 + q^2y^2)} \biggr]
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{1 + w^2} \biggl[ \frac{q^2y}{(x^2 + q^2y^2)^{1 / 2}~p^2z} \biggr]
=
\frac{w}{1 + w^2} \biggl[ \frac{q^2y}{(x^2 + q^2y^2)} \biggr]
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{1 + w^2} \biggl[- \frac{(x^2 + q^2y^2)^{1 / 2}}{~p^2z^2} \biggr]
=
\frac{w}{1 + w^2} \biggl[- \frac{1}{z} \biggr]
\, ,
</math>
  </td>
</tr>
</table>
 
which means that,
<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>~
\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{w}{1 + w^2} \biggr]^2 \biggl\{
\biggl[  \frac{x}{(x^2 + q^2y^2)} \biggr]^2
+
\biggl[ \frac{q^2y}{(x^2 + q^2y^2)}  \biggr]^2
+
\biggl[ - \frac{1}{z} \biggr]^2
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{w}{1 + w^2} \biggr]^2 \biggl\{
\frac{x^2 + q^4y^2}{(x^2 + q^2y^2)^2}
+
\frac{1}{z^2}
\biggr\}
</math>
  </td>
</tr>
</table>
 
===Speculation5===
 
====Spherical Coordinates====
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~r\cos\theta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~z \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~r\sin\theta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(x^2 + y^2)^{1 / 2} \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\tan\varphi</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{y}{x} \, .</math>
  </td>
</tr>
</table>
 
====Use &lambda;<sub>1</sub> Instead of r ====
 
Here, as above, we define,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\lambda_1^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~x^2 + q^2 y^2 + p^2 z^2 </math>
  </td>
</tr>
</table>
 
Using this expression to eliminate "x" (in favor of &lambda;<sub>1</sub>) in each of the three spherical-coordinate definitions, we obtain,
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~r^2 \equiv x^2 + y^2 + z^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\lambda_1^2 - y^2(q^2-1) - z^2(p^2-1) \, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\tan^2\theta \equiv \frac{x^2 + y^2}{z^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{z^2}\biggl[
\lambda_1^2 -y^2(q^2-1) -p^2z^2
\biggr] \, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{1}{\tan^2\varphi} \equiv \frac{x^2}{y^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\lambda_1^2 - p^2z^2}{y^2} - q^2 \, .
</math>
  </td>
</tr>
</table>
 
After a bit of additional algebraic manipulation, we find that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>\frac{z^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ (1+\tan^2\varphi)}{\mathcal{D}^2} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{y^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[\frac{ \mathcal{D}^2 \tan^2\varphi -  p^2 \tan^2\varphi (1+\tan^2\varphi)}{(1+q^2\tan^2\varphi) \mathcal{D}^2}  \biggr]
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{x^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1 - q^2\biggl(\frac{y^2}{\lambda_1^2} \biggr) - p^2\biggl(\frac{z^2}{\lambda_1^2}\biggr)
\, ,
</math>
  </td>
</tr>
</table>
 
where,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\mathcal{D}^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\biggl[
(1 + q^2\tan^2\varphi)(p^2 + \tan^2\theta) - p^2(q^2-1)\tan^2\varphi
\biggr] \, .
</math>
  </td>
</tr>
</table>
 
<table border="1" cellpadding="8" align="center" width="80%"><tr><td align="left">
As a check, let's set <math>~q^2 = p^2 = 1</math>, which should reduce to the normal spherical coordinate system.
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\lambda_1^2</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
r^2 \, ,
</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp; </td>
  <td align="right">
<math>~\mathcal{D}^2</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
\biggl[
(1 + \tan^2\varphi)(1 + \tan^2\theta)
\biggr] \, .
</math>
  </td>
</tr>
</table>
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{z^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
\frac{1}{1+\tan^2\theta} = \cos^2\theta = \frac{z^2}{r^2} \, ;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{y^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
\biggl[\frac{(1 + \tan^2\varphi)(1 + \tan^2\theta) \tan^2\varphi -  \tan^2\varphi (1+\tan^2\varphi)}{(1+\tan^2\varphi) (1 + \tan^2\varphi)(1 + \tan^2\theta)}  \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\tan^2\varphi}{(1 + \tan^2\varphi)}
\biggl[\frac{\tan^2\theta }{ (1 + \tan^2\theta)}  \biggr]
= \sin^2\theta \sin^2\varphi
= \frac{y^2}{r^2} \,;
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{x^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
1 - \biggl(\frac{y^2}{\lambda_1^2} \biggr) - \biggl(\frac{z^2}{\lambda_1^2}\biggr)
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
1 - \sin^2\theta \sin^2\varphi - \cos^2\theta
=
- \sin^2\theta \sin^2\varphi + \sin^2\theta
=
\sin^2\theta \cos^2\varphi = \frac{x^2}{r^2}
\, .
</math>
  </td>
</tr>
</table>
</td></tr></table>
 
====Relationship To T3 Coordinates====
 
If we set, <math>~q = 1</math>, but continue to assume that <math>~p > 1</math>, we expect to see a representation that resembles our previously discussed, [[User:Tohline/Appendix/Ramblings/T3Integrals#Integrals_of_Motion_in_T3_Coordinates|T3 Coordinates]].  Note, for example, that the new "radial" coordinate is,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\lambda_1^2</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
(\varpi^2 + p^2z^2) \, ,
</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp; </td>
  <td align="right">
<math>~\mathcal{D}^2</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
\biggl[
(1 + \tan^2\varphi)(p^2 + \tan^2\theta)
\biggr]  \, .
</math>
  </td>
</tr>
</table>
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{p^2z^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~\frac{ p^2}{(p^2 + \tan^2\theta)} = \frac{ 1}{(1 + p^{-2} \tan^2\theta)}\, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\varpi^2}{\lambda_1^2} = \frac{x^2}{\lambda_1^2} + \frac{y^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~1 -  p^2 \biggl( \frac{z^2}{\lambda_1^2}\biggr)
=
\biggl[1 -  \frac{ 1}{(1 + p^{-2}\tan^2\theta)} \biggr] \, .</math>
  </td>
</tr>
</table>
 
We also see that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\varpi^2}{p^2z^2}</math>
  </td>
  <td align="center">
<math>~\rightarrow</math>
  </td>
  <td align="left">
<math>~
(1 + p^{-2}\tan^2\theta)\biggl[1 -  \frac{ 1}{(1 + p^{-2}\tan^2\theta)} \biggr]
=
p^{-2}\tan^2\theta \, .
</math>
  </td>
</tr>
</table>
 
 
====Again Consider Full 3D Ellipsoid====
Let's try to replace everywhere, <math>~[\varpi/(pz)]^2 = p^{-2}\tan^2\theta</math> with <math>~\lambda_2</math>.  This gives,
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\mathcal{D}^2}{p^2}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\biggl[
(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi
\biggr] \, .
</math>
  </td>
</tr>
</table>
 
which means that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>\frac{p^2 z^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ (1+\tan^2\varphi)}{[(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi]}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ (1+\tan^2\varphi)/\tan^2\varphi}{[q^2 \lambda_2 (1 + q^2\tan^2\varphi)/(q^2\tan^2\varphi) - (q^2-1)]}
= \frac{1/\sin^2\varphi}{[q^2\lambda_2 Q^2 - (q^2-1) ]}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{q^2y^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ q^2 \tan^2\varphi }{(1+q^2\tan^2\varphi) } 
-
\frac{ q^2 \tan^2\varphi (1+\tan^2\varphi)}{(1+q^2\tan^2\varphi) [(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi] }
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ q^2 \tan^2\varphi }{(1+q^2\tan^2\varphi) } \biggl\{1 
-
\frac{ (1+\tan^2\varphi)}{ [(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi] }
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ q^2 \tan^2\varphi }{(1+q^2\tan^2\varphi) } \biggl\{1 
-
\frac{ (1+\tan^2\varphi)/\tan^2\varphi}{ [q^2\lambda_2(1 + q^2\tan^2\varphi)/(q^2\tan^2\varphi) - (q^2-1) ] }
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ 1}{Q^2 } \biggl\{1 
-
\frac{ 1/\sin^2\varphi}{ [q^2\lambda_2 Q^2 - (q^2-1) ] }
\biggr\}
= \frac{1}{Q^2[~~]} \biggl[ [~~] - \frac{1}{\sin^2\varphi} \biggr]
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{x^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1
- \frac{ q^2 \tan^2\varphi }{(1+q^2\tan^2\varphi) } \biggl\{1 
-
\frac{ (1+\tan^2\varphi)}{ [(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi] }
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
- \frac{ (1+\tan^2\varphi)}{[(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi]}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1
- \frac{ q^2 \tan^2\varphi }{(1+q^2\tan^2\varphi) }
- \biggl\{\frac{ (1+\tan^2\varphi)}{ [(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi] }
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1
- \frac{ 1 }{ Q^2 }
- \biggl\{\frac{ 1/\sin^2\varphi}{ [q^2\lambda_2 Q^2 - (q^2-1) ] }
\biggr\}
=
\frac{1}{Q^2 [~~] } \biggl\{
Q^2[~~] - [~~] - \frac{Q^2}{\sin^2\varphi}
\biggr\}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{x^2 + q^2y^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~1
- \biggl[1 + \frac{ q^2 \tan^2\varphi }{(1+q^2\tan^2\varphi) }\biggr] \biggl\{
\frac{ (1+\tan^2\varphi)}{ [(1 + q^2\tan^2\varphi)\lambda_2 - (q^2-1)\tan^2\varphi] }
\biggr\} \, .
</math>
  </td>
</tr>
</table>
 
Now, notice that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{q^2 y^2 Q^2}{\lambda_1^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~1 - \frac{1}{[~~]\sin^2\varphi} \, ,</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{x^2}{\lambda_1^2} + \frac{1}{Q^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~1 - \frac{1}{[~~]\sin^2\varphi} \, .</math>
  </td>
</tr>
</table>
Hence,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{x^2}{\lambda_1^2} + \frac{1}{Q^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{q^2 y^2 Q^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>~Q^4 \biggl( \frac{q^2 y^2}{\lambda_1^2} \biggr) - Q^2 \biggl( \frac{x^2}{\lambda_1^2} \biggr) - 1</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~Q^4  - Q^2 \biggl( \frac{x^2}{q^2 y^2} \biggr) - \biggl( \frac{\lambda_1^2}{q^2 y^2} \biggr) \, ,</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~Q^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{1 + q^2\tan^2\varphi}{q^2\tan^2\varphi} \, .</math>
  </td>
</tr>
</table>
Solving the quadratic equation, we have,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~Q^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{2} \biggl\{ \biggl( \frac{x^2}{q^2 y^2} \biggr) \pm \biggl[ \biggl( \frac{x^2}{q^2 y^2} \biggr)^2 + 4\biggl( \frac{\lambda_1^2}{q^2 y^2} \biggr) \biggr]^{1 / 2}  \biggr\}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{x^2}{2q^2 y^2} \biggr) \biggl\{ 1  \pm \biggl[ 1 + 4\biggl( \frac{\lambda_1^2 q^2 y^2}{x^4} \biggr) \biggr]^{1 / 2}  \biggr\} \, .</math>
  </td>
</tr>
</table>
 
<table border="1" align="center" width="80%" cellpadding="10"><tr><td align="left">
<div align="center">'''Tentative Summary'''</div>
<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^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^2 + y^2)^{1 / 2}}{pz}
\, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\lambda_3 = Q^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{x^2}{2q^2 y^2} \biggr) \biggl\{ 1  \pm \biggl[ 1 + 4\biggl( \frac{\lambda_1^2 q^2 y^2}{x^4} \biggr) \biggr]^{1 / 2}  \biggr\} \, .</math>
  </td>
</tr>
</table>
 
</td></tr></table>
 
====Partial Derivatives &amp; Scale Factors====
=====First Coordinate=====
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_1}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{x}{\lambda_1} \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_1}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{q^2 y}{\lambda_1} \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_1}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{p^2 z}{\lambda_1} \, .</math>
  </td>
</tr>
</table>
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~h_1^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{\partial \lambda_1}{\partial x} \biggr)^2
+ \biggl( \frac{\partial \lambda_1}{\partial y} \biggr)^2
+ \biggl( \frac{\partial \lambda_1}{\partial z} \biggr)^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{x}{\lambda_1} \biggr)^2
+ \biggl( \frac{q^2 y}{\lambda_1} \biggr)^2
+ \biggl( \frac{p^2 z}{\lambda_1} \biggr)^2
\, .</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ h_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left" colspan="3">
<math>~
\lambda_1 \ell_{3D} \, ,
</math>
  </td>
</tr>
</table>
where,
<div align="center"><math>\ell_{3D} \equiv (x^2 + q^4y^2 + p^4z^2)^{-1 / 2} \, .</math></div>
As a result, the associated unit vector is,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath} h_1 \biggl( \frac{\partial \lambda_1}{\partial x} \biggr)
+ \hat{\jmath} h_1 \biggl( \frac{\partial \lambda_1}{\partial y} \biggr)
+ \hat{k} h_1 \biggl( \frac{\partial \lambda_1}{\partial z} \biggr)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </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^2 z \ell_{3D}
\, .
</math>
  </td>
</tr>
</table>
 
Notice that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(x^2 + q^4 y^2 + p^4 z^2) \ell_{3D}^2 = 1 \, .
</math>
  </td>
</tr>
</table>
 
 
=====Second Coordinate (1<sup>st</sup> Try)=====
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{pz} \biggl[ \frac{x}{(x^2 + y^2)^{1 / 2}} \biggr] \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{pz} \biggl[ \frac{y}{(x^2 + y^2)^{1 / 2}} \biggr]
\, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \frac{(x^2 + y^2)^{1 / 2}}{pz^2}
\, .</math>
  </td>
</tr>
</table>
 
<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>~
\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{1}{pz} \biggl[ \frac{x}{(x^2 + y^2)^{1 / 2}} \biggr] \biggr\}^2
+ \biggl\{ \frac{1}{pz} \biggl[ \frac{y}{(x^2 + y^2)^{1 / 2}} \biggr] \biggr\}^2
+ \biggl\{ \frac{(x^2 + y^2)^{1 / 2}}{pz^2} \biggr\}^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{x^2}{(x^2 + y^2)p^2 z^2} \biggr] \biggr\}
+ \biggl\{ \biggl[ \frac{y^2}{(x^2 + y^2)p^2 z^2} \biggr] \biggr\}
+ \biggl\{ \frac{(x^2 + y^2)}{p^2 z^4} \biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{p^2 z^2}
+ \frac{(x^2 + y^2)}{p^2 z^4}
=
\frac{(x^2 + y^2 + z^2)}{p^2 z^4}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow~~~h_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{p z^2}{r }
</math>
  </td>
</tr>
</table>
As a result, the associated unit vector is,
<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} h_2 \biggl( \frac{\partial \lambda_2}{\partial x} \biggr)
+ \hat{\jmath} h_2 \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)
+ \hat{k} h_2 \biggl( \frac{\partial \lambda_2}{\partial z} \biggr)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath}  \biggl[ \frac{xz}{r(x^2 + y^2)^{1 / 2}} \biggr]
+ \hat{\jmath}  \biggl[ \frac{yz}{r(x^2 + y^2)^{1 / 2}} \biggr]
- \hat{k} \biggl[ \frac{(x^2 + y^2)^{1 / 2}}{r} \biggr] \, .
</math>
  </td>
</table>
 
Notice that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_2 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x^2 z^2}{r^2(x^2 + y^2)} \biggr]
+ \biggl[ \frac{y^2 z^2}{r^2(x^2 + y^2)} \biggr]
+ \biggl[ \frac{(x^2 + y^2)}{r^2} \biggr]
= 1 \, .
</math>
  </td>
</tr>
</table>
 
Let's check to see if this "second" unit vector is orthogonal to the "first."
<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{xz}{r(x^2 + y^2)^{1 / 2}} \biggr]
+ q^2 y\ell_{3D} \biggl[ \frac{yz}{r(x^2 + y^2)^{1 / 2}} \biggr]
- p^2 z \ell_{3D} \biggl[ \frac{(x^2 + y^2)^{1 / 2}}{r} \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\ell_{3D} \biggl\{ \biggl[ \frac{x^2z}{r(x^2 + y^2)^{1 / 2}} \biggr]
+ \biggl[ \frac{q^2 y^2 z}{r(x^2 + y^2)^{1 / 2}} \biggr]
-  \biggl[ \frac{p^2 z(x^2 + y^2)^{1 / 2}}{r} \biggr]
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{z\ell_{3D}}{r (x^2 + y^2)^{1 / 2}} \biggl\{ \biggl[ x^2\biggr]
+ \biggl[ q^2 y^2 \biggr]
-  \biggl[ p^2 (x^2 + y^2) \biggr]
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\ne</math>
  </td>
  <td align="left">
<math>~
0 \ .
</math>
  </td>
</tr>
</table>
 
 
=====Second Coordinate (2<sup>nd</sup> Try)=====
 
Let's try,
<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>~
\biggl[\frac{(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2}}{pz}  \biggr]
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow~~~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{x}{pz(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2} } = \frac{x}{p^2 z^2 \lambda_2}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{q^2 y}{pz(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2} } = \frac{q^2y}{p^2 z^2 \lambda_2}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\mathfrak{f}\cdot p^2z}{pz(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2} }
-
\frac{(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2}}{pz^2} 
= \frac{1}{p^2z^2 \lambda_2 } \biggl( \mathfrak{f}\cdot p^2z \biggr) - \frac{\lambda_2 }{z}
= \frac{1}{p^2z^2 \lambda_2 } \biggl( \mathfrak{f}\cdot p^2z  - p^2z \lambda_2^2 \biggr)
\, .
</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>~
\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{x}{p^2 z^2 \lambda_2} \biggr]^2
+ \biggl[ \frac{q^2y}{p^2 z^2 \lambda_2} \biggr]^2
+ \biggl[ \frac{ \mathfrak{f} }{z \lambda_2 }  - \frac{\lambda_2 }{z}\biggr]^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x^2 + q^4 y^2}{p^4 z^4 \lambda_2^2} \biggr]
+ \biggl[ \frac{1}{z\lambda_2}\biggl( \mathfrak{f} - \lambda_2^2 \biggr) \biggr]^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{1}{p^4 z^4 \lambda_2^2}
\biggl[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^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>~
\frac{p^2 z^2 \lambda_2}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \, .
</math>
  </td>
</tr>
</table>
So, the associated unit vector is,
<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} h_2 \biggl( \frac{\partial \lambda_2}{\partial x} \biggr)
+ \hat{\jmath} h_2 \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)
+ \hat{k} h_2 \biggl( \frac{\partial \lambda_2}{\partial z} \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}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ \hat{\jmath}  \biggl\{ \frac{q^2 y}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ \hat{k} \biggl\{ \frac{p^2z(\mathfrak{f}-\lambda_2^2)}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\} \, .
</math>
  </td>
</table>
 
Checking orthogonality &hellip;
<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{x}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ q^2 y\ell_{3D} \biggl\{ \frac{q^2 y}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ p^2 z \ell_{3D} \biggl\{ \frac{p^2z(\mathfrak{f}-\lambda_2^2)}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^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^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} }
\biggl\{ x^2 + q^4y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^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^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} }
\biggl\{ x^2 + q^4y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)\biggr\} \, .
</math>
  </td>
</tr>
</table>
 
If <math>~\mathfrak{f} = 0</math>, we have &hellip;
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~p^2 z (\mathfrak{f} - \lambda_2^2) </math>
  </td>
  <td align="center">
&nbsp; &nbsp; <math>~~~\rightarrow ~~~ </math>&nbsp; &nbsp;
  </td>
  <td align="left">
<math>~
\biggl[- p^2 z \lambda_2^2\biggr]_{\mathfrak{f} = 0}
=
- p^2 z\biggl[\frac{(x^2 + q^2y^2 + \cancelto{0}{\mathfrak{f} \cdot }p^2 z^2 )^{1 / 2}}{pz}  \biggr]^2
=
- \frac{(x^2 + q^2y^2  )}{z} \, ,</math>
  </td>
</tr>
</table>
which, in turn, means &hellip;
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~[ x^2 + q^4 y^2 + p^4 z^2 (\cancelto{0}{\mathfrak{f}} - \lambda_2^2)^2 ]^{1 / 2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[ x^2 + q^4 y^2 + p^4 z^2 \lambda_2^4 ]^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ x^2 + q^4 y^2 + p^4 z^2 \biggl[\frac{(x^2 + q^2y^2 + \cancelto{0}{\mathfrak{f}\cdot} p^2 z^2)^{1 / 2}}{pz}  \biggr]^4 \biggr\}^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ x^2 + q^4 y^2 + \biggl[\frac{(x^2 + q^2y^2 )^{2}}{z^2}  \biggr] \biggr\}^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(x^2 + q^4 y^2)^{1 / 2} \biggl[ 1 + \frac{(x^2 + q^2y^2 )}{z^2} \biggr]^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{(x^2 + q^4 y^2)^{1 / 2}}{z} \biggl[ z^2 + (x^2 + q^2y^2 ) \biggr]^{1 / 2} \, ,
</math>
  </td>
</tr>
</table>
and,
<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>~\frac{\ell_{3D}}{ [ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} }
\biggl\{ x^2 + q^4y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)\biggr\} \, .
</math>
  </td>
</tr>
</table>
 
===Speculation6===
 
====Determine &lambda;<sub>2</sub>====
This is very similar to the [[#Speculation2|above, Speculation2]].
Try,
<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}}{ z^{2/p^2}} \, ,
</math>
  </td>
</tr>
</table>
in which case,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\lambda_2}{x} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{x}{z^{2/p^2}} \biggl(\frac{1}{q^2}\biggr) y^{1/q^2 - 1}
=
\frac{\lambda_2}{q^2 y}
\, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{2\lambda_2}{p^2 z}
\, .
</math>
  </td>
</tr>
</table>
The associated scale factor is, then,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~h_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[
\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
\biggr]^{-1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[
\biggl( \frac{ \lambda_2}{x} \biggr)^2
+
\biggl( \frac{\lambda_2}{q^2y} \biggr)^2
+
\biggl( - \frac{2\lambda_2}{p^2z} \biggr)^2
\biggr]^{-1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{\lambda_2}\biggl[
\frac{ 1}{x^2}
+
\frac{1}{q^4y^2}
+
\frac{4}{p^4z^2}
\biggr]^{-1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{\lambda_2}\biggl[
\frac{ (q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2}{x^2 q^4 y^2 p^4 z^2}
\biggr]^{-1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{\lambda_2}\biggl[
\frac{x q^2 y p^2 z}{ \mathcal{D}}
\biggr] \, .
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" 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>
 
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} h_2 \biggl( \frac{\partial \lambda_2}{\partial x} \biggr)
+ \hat{\jmath} h_2 \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)
+ \hat{k} h_2 \biggl( \frac{\partial \lambda_2}{\partial z} \biggr)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{x q^2 y p^2 z}{\mathcal{D}} \biggl\{
\hat{\imath} \biggl( \frac{1}{x} \biggr)
+ \hat{\jmath}  \biggl( \frac{1}{q^2 y} \biggr)
+ \hat{k} \biggl( -\frac{2}{p^2 z} \biggr)
\biggr\} \ .
</math>
  </td>
</tr>
</table>
 
Recalling that the unit vector associated with the "first" coordinate is,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_1 </math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\hat\imath (x \ell_{3D}) + \hat\jmath (q^2y \ell_{3D}) + \hat{k} (p^2 z \ell_{3D}) \, ,
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\ell_{3D}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2 + p^4 z^2 )^{- 1 / 2} \, ,</math>
  </td>
</tr>
</table>
let's check to see whether the "second" unit vector is orthogonal to the "first."
<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>~
\frac{(x q^2 y p^2 z) \ell_{3D}}{\mathcal{D}} \biggl[
1 + 1 - 2
\biggr] = 0 \, .
</math>
  </td>
</tr>
</table>
 
<font color="red">'''Hooray!'''</font>
 
====Direction Cosines for <i>Third</i> Unit Vector====
Now, what is the unit vector, <math>~\hat{e}_3</math>, that is simultaneously orthogonal to both these "first" and the "second" unit vectors?
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_3 \equiv \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_{2y} )( e_{1z}) ) \biggl]
+ \hat\jmath \biggl[ ( e_{1z} )( e_{2x}) - ( e_{2z} )( e_{1x}) ) \biggl]
+ \hat{k} \biggl[ ( e_{1x} )( e_{2y}) - ( e_{2x} )( e_{1y}) ) \biggl]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{(x q^2 y p^2 z) \ell_{3D}}{\mathcal{D}}
\biggl\{
\hat\imath \biggl[ \biggl( -\frac{2 q^2y}{p^2 z} \biggr) - \biggl( \frac{p^2z}{q^2y} \biggr)  \biggl]
+ \hat\jmath \biggl[ \biggl( \frac{p^2z}{x} \biggr) - \biggl(-\frac{2x}{p^2z} \biggr)  \biggl]
+ \hat{k} \biggl[ \biggl( \frac{x}{q^2y} \biggr) - \biggl( \frac{q^2y}{x} \biggr)  \biggl]
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{(x q^2 y p^2 z) \ell_{3D}}{\mathcal{D}}
\biggl\{
-\hat\imath \biggl[ \frac{2 q^4y^2 + p^4z^2}{q^2 y p^2 z} \biggl]
+ \hat\jmath \biggl[ \frac{p^4z^2 + 2x^2}{xp^2 z}  \biggl]
+ \hat{k} \biggl[ \frac{x^2 - q^4y^2}{x q^2y} \biggl]
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}}{\mathcal{D}}
\biggl\{
-\hat\imath \biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]
+ \hat\jmath \biggl[ q^2 y(p^4z^2 + 2x^2 )  \biggl]
+ \hat{k} \biggl[ p^2z( x^2 - q^4y^2 ) \biggl]
\biggr\} \, .
</math>
  </td>
</tr>
</table>
 
Is this a valid unit vector?  First, note that &hellip;
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)
(x^2 + q^4y^2 + p^4 z^2 )
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(x^2 q^4 y^2 p^4 z^2 + x^4 p^4 z^2 + 4x^4q^4y^2)
+ (q^8 y^4 p^4 z^2 + x^2 q^4y^2 p^4 z^2 + 4x^2q^8y^4)
+(q^4 y^2 p^8 z^4 + x^2 p^8 z^4 + 4x^2q^4y^2 p^4 z^2)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
6x^2 q^4 y^2 p^4 z^2 + x^4(p^4 z^2 + 4 q^4y^2)
+ q^8 y^4(p^4 z^2  + 4x^2)
+p^8z^4(x^2 + q^4 y^2 )\, .
</math>
  </td>
</tr>
</table>
 
<span id="Eureka">Then we have,</span>
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2}\hat{e}_3 \cdot \hat{e}_3</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>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x^2(4 q^8y^4 + 4q^4y^2p^4z^2 + p^8z^4 )
+
q^4 y^2(p^8z^4 + 4x^2p^4z^2 + 4x^4 )
+
p^4z^2( x^4 - 2x^2q^4 y^2 + q^8y^4 )
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
4 x^2 q^8y^4 + 4x^2 q^4y^2p^4z^2 + x^2 p^8z^4
+
q^4 y^2p^8z^4 + 4x^2q^4 y^2p^4z^2 + 4x^4q^4 y^2
+
x^4p^4z^2 - 2x^2q^4 y^2p^4z^2 + q^8y^4p^4z^2 
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
6x^2 q^4y^2p^4z^2
+  p^8z^4 (x^2  +q^4 y^2)
+  x^4(4q^4 y^2 + p^4z^2)
+ q^8 y^4(4 x^2  + p^4z^2 )
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2} \, ,
</math>
  </td>
</tr>
</table>
which means that, <math>~\hat{e}_3\cdot \hat{e}_3 = 1</math>.  &nbsp; &nbsp;<font color="red">'''Hooray! Again (11/11/2020)!'''</font>
 
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" colspan="9">'''Direction Cosine Components for T6 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}}{ z^{2/p^2}}</math></td>
  <td align="center"><math>~\frac{1}{\lambda_2}\biggl[\frac{x q^2 y p^2 z}{ \mathcal{D}}\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>~-\frac{2\lambda_2}{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>~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{\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>~\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>
 
<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>
  </td>
</tr>
</table>
 
 
Let's double-check whether this "third" unit vector is orthogonal to both the "first" and the "second" unit vectors.
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}^2}{\mathcal{D}}
\biggl\{
-x \biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]
+ q^2 y  \biggl[ q^2 y(p^4z^2 + 2x^2 )  \biggl]
+ p^2 z  \biggl[ p^2z( x^2 - q^4y^2 ) \biggl]
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}^2}{\mathcal{D}}
\biggl\{
- (2 x^2q^4y^2 + x^2p^4z^2 )
+ (q^4 y^2 p^4z^2 + 2x^2 q^4 y^2)
+ ( x^2p^4z^2  - q^4y^2 p^4z^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>
 
and,
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_2 \cdot \hat{e}_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}}{\mathcal{D}} \cdot \frac{x q^2 y p^2 z}{\mathcal{D}}
\biggl\{
- \biggl[ (2 q^4y^2 + p^4z^2 ) \biggl]
+  \biggl[ (p^4z^2 + 2x^2 )  \biggl]
- \biggl[ 2( x^2 - q^4y^2 ) \biggl]
\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>
 
<font color="red">'''Q. E. D.'''</font>
 
====Search for <i>Third</i> Coordinate Expression====
Let's try &hellip;
<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>~\mathcal{D}^n \ell_{3D}^m </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)^{n / 2} (x^2 + q^4y^2 + p^4 z^2 )^{- m / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{\partial \lambda_3}{\partial x_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\ell_{3D}^m \biggl[ \frac{n}{2} \biggl(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2 \biggr)^{n / 2 - 1} \biggr] \frac{\partial}{\partial x_i} \biggl[(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+ \mathcal{D}^n \biggl[ - \frac{m}{2} (x^2 + q^4y^2 + p^4 z^2 )^{- m / 2 - 1} \biggr] \frac{\partial}{\partial x_i} \biggl[(x^2 + q^4y^2 + p^4 z^2 )\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\mathcal{D}^n \ell_{3D}^m \biggl[ \frac{n}{2} \biggl(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2 \biggr)^{- 1} \biggr] \frac{\partial}{\partial x_i} \biggl[(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+ \mathcal{D}^n \ell_{3D}^m \biggl[ - \frac{m}{2} (x^2 + q^4y^2 + p^4 z^2 )^{- 1} \biggr] \frac{\partial}{\partial x_i} \biggl[(x^2 + q^4y^2 + p^4 z^2 )\biggr] \, .
</math>
  </td>
</tr>
</table>
 
Hence,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{1}{\mathcal{D}^n \ell_{3D}^m} \cdot \frac{\partial \lambda_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{n}{2\mathcal{D}^2}\frac{\partial}{\partial x} \biggl[q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2 \biggr]
-
\frac{m \ell_{3D}^2}{2} \frac{\partial}{\partial x} \biggl[ x^2 + q^4y^2 + p^4 z^2 \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x \biggl\{
\frac{n}{\mathcal{D}^2}\biggl[p^4 z^2 + 4q^4y^2 \biggr]
-
m \ell_{3D}^2
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x \biggl\{
\frac{n (p^4 z^2 + 4q^4y^2)}{ ( q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2 ) }
-
\frac{m}{ ( x^2 + q^4y^2 + p^4 z^2 ) }
\biggr\}
</math>
  </td>
</tr>
</table>
 
This is overly cluttered!  Let's try, instead &hellip;
 
 
<table border="1" cellpadding="8" align="center" width="80%"><tr><td align="left">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~A \equiv \ell_{3D}^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2 + p^4 z^2 ) \, ,</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp; </td>
  <td align="right">
<math>~B \equiv \mathcal{D}^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)</math>
  </td>
</tr>
</table>
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{\partial A}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2x \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial A}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2q^4 y \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial A}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 2p^4 z\, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~ \frac{\partial B}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2x( 4q^4y^2 + p^4 z^2 ) \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial B}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2q^4 y (p^4 z^2 + 4x^2) \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{\partial B}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 2p^4 z(q^4 y^2 + x^2)\, .</math>
  </td>
</tr>
</table>
 
</td></tr></table>
 
 
Now, let's assume that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{A}{B} \biggr)^{1 / 2} \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow~~~ \frac{ \partial \lambda_3}{\partial x_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{2 (AB)^{1 / 2}} \cdot \frac{\partial A}{\partial x_i}
-
\frac{A^{1 / 2}}{2 B^{3 / 2}} \cdot \frac{\partial B}{\partial x_i}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\lambda_3}{2AB}
\biggl[
B \cdot \frac{\partial A}{\partial x_i}
-
A \cdot \frac{\partial B}{\partial x_i}
\biggr]
\, .
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl[ \frac{2AB}{\lambda_3} \biggr] \frac{\partial \lambda_3}{\partial x_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2) \cdot \frac{\partial A}{\partial x_i}
-
(x^2 + q^4y^2 + p^4 z^2 ) \cdot \frac{\partial B}{\partial x_i} \, .
</math>
  </td>
</tr>
</table>
 
<table border="1" cellpadding="8" align="center" width="80%"><tr><td align="left">
Looking ahead &hellip;
<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{\lambda_3}{2AB} \biggl[
B \cdot \frac{\partial A}{\partial x}
-
A \cdot \frac{\partial B}{\partial x}
\biggr] \biggr\}^2
+
\biggl\{ \frac{\lambda_3}{2AB} \biggl[
B \cdot \frac{\partial A}{\partial y}
-
A \cdot \frac{\partial B}{\partial y}
\biggr] \biggr\}^2
+
\biggl\{ \frac{\lambda_3}{2AB} \biggl[
B \cdot \frac{\partial A}{\partial z}
-
A \cdot \frac{\partial B}{\partial z}
\biggr] \biggr\}^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \biggl[\frac{2AB}{\lambda_3}  \biggr]^2 h_3^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[
B \cdot \frac{\partial A}{\partial x}
-
A \cdot \frac{\partial B}{\partial x}
\biggr]^2
+
\biggl[
B \cdot \frac{\partial A}{\partial y}
-
A \cdot \frac{\partial B}{\partial y}
\biggr]^2
+
\biggl[
B \cdot \frac{\partial A}{\partial z}
-
A \cdot \frac{\partial B}{\partial z}
\biggr]^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \biggl[\frac{\lambda_3}{2AB}  \biggr] h_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{\biggl[
B \cdot \frac{\partial A}{\partial x}
-
A \cdot \frac{\partial B}{\partial x}
\biggr]^2
+
\biggl[
B \cdot \frac{\partial A}{\partial y}
-
A \cdot \frac{\partial B}{\partial y}
\biggr]^2
+
\biggl[
B \cdot \frac{\partial A}{\partial z}
-
A \cdot \frac{\partial B}{\partial z}
\biggr]^2
\biggr\}^{-1 / 2}
</math>
  </td>
</tr>
</table>
Then, for example,
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\gamma_{31} \equiv h_3 \biggl(\frac{\partial \lambda_3}{\partial x} \biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[
B \cdot \frac{\partial A}{\partial x}
-
A \cdot \frac{\partial B}{\partial x}
\biggr]
\biggl\{\biggl[
B \cdot \frac{\partial A}{\partial x}
-
A \cdot \frac{\partial B}{\partial x}
\biggr]^2
+
\biggl[
B \cdot \frac{\partial A}{\partial y}
-
A \cdot \frac{\partial B}{\partial y}
\biggr]^2
+
\biggl[
B \cdot \frac{\partial A}{\partial z}
-
A \cdot \frac{\partial B}{\partial z}
\biggr]^2
\biggr\}^{-1 / 2}
</math>
  </td>
</tr>
</table>
 
</td></tr></table>
 
As a result, we have,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl[ \frac{2AB}{\lambda_3} \biggr] \frac{\partial \lambda_3}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2x \biggl[
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2) 
-
(x^2 + q^4y^2 + p^4 z^2 ) ( 4q^4y^2 + p^4 z^2 ) \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2x \biggl[
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2) 
-
(4x^2q^4y^2 + x^2 p^4 z^2 + 4q^8y^4 + q^4y^2 p^4z^2 + 4q^4y^2 p^4 z^2 + p^8z^4) \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2x \biggl[ 
-
(4q^8y^4  + 4q^4y^2 p^4 z^2 + p^8z^4) \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-2x (2q^4y^2  + p^4z^2)^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-8x \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl[ AB \biggr] \frac{\partial \ln\lambda_3}{\partial \ln{x}}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-\biggl[ 2x \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)\biggr]^2 \, ;
</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl[ \frac{2AB}{\lambda_3} \biggr] \frac{\partial \lambda_3}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2q^4y\biggl[
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2) 
-
(x^2 + q^4y^2 + p^4 z^2 ) (p^4 z^2 + 4x^2)
\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2q^4y\biggl[
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2) 
-
( x^2p^4z^2 + 4x^4 + q^4y^2p^4z^2 + 4x^2q^4y^2 + p^8z^4 + 4x^2p^4z^2)
\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-2q^4y( 4x^4 + p^8z^4 + 4x^2p^4z^2)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-2q^4y( 2x^2 + p^4z^2 )^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl[ AB \biggr] \frac{\partial \ln\lambda_3}{\partial \ln{y} }</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \biggl[2q^2y\biggl( x^2 + \frac{p^4z^2}{2} \biggr) \biggr]^2 \, ;
</math>
  </td>
</tr>
</table>
 
and,
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl[ \frac{2AB}{\lambda_3} \biggr] \frac{\partial \lambda_3}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2p^4 z \biggl[
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2) 
-
(x^2 + q^4y^2 + p^4 z^2 ) (q^4 y^2 + x^2)
\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2p^4 z \biggl[
(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2) 
-
(x^2q^4y^2 + x^4 + q^8y^4 + x^2q^4y^2 + q^4y^2p^4z^2 + x^2p^4z^2)
\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2p^4 z \biggl[
( 2x^2q^4y^2) 
-
( x^4 + q^8y^4 )
\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-2p^4 z \biggl[
x^4 + q^8y^4
- 2x^2q^4y^2 
\biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-2p^4 z (x^2 - q^4y^2 )^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl[ AB \biggr] \frac{\partial \ln \lambda_3}{\partial \ln{z}}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-4 \biggl[ \biggl( \frac{p^4 z^2}{4} \biggr) (x^2 - q^4y^2 )^2 \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-\biggl[ 2\biggl( \frac{p^2 z}{2} \biggr) (x^2 - q^4y^2 ) \biggr]^2
\, .
</math>
  </td>
</tr>
</table>
 
<font color="red">'''Wow! &nbsp; Really close!''' (13 November 2020)</font>
 
 
Just for fun, let's see what we get for <math>~h_3</math>.  It is given by the expression,
<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{\partial \lambda_3}{\partial x} \biggr)^2
+\biggl( \frac{\partial \lambda_3}{\partial y} \biggr)^2
+\biggl( \frac{\partial \lambda_3}{\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_3}{ABx} \biggl[ 2x \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)\biggr]^2 \biggr\}^2
+\biggl\{ \frac{\lambda_3}{ABy} \biggl[2q^2y\biggl( x^2 + \frac{p^4z^2}{2} \biggr) \biggr]^2 \biggr\}^2
+\biggl\{ \frac{\lambda_3}{ABz} \biggl[ 2\biggl( \frac{p^2 z}{2} \biggr) (x^2 - q^4y^2 ) \biggr]^2 \biggr\}^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow~~~ \biggl[ \frac{AB}{\lambda_3}\biggr]^2 h_3^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \frac{1}{x^2} \biggl[ 2x \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)\biggr]^4 \biggr\}
+\biggl\{ \frac{1}{y^2} \biggl[2q^2y\biggl( x^2 + \frac{p^4z^2}{2} \biggr) \biggr]^4 \biggr\}
+\biggl\{ \frac{1}{z^2} \biggl[ 2\biggl( \frac{p^2 z}{2} \biggr) (x^2 - q^4y^2 ) \biggr]^4 \biggr\}
</math>
  </td>
</tr>
</table>
 
====Fiddle Around====
 
Let &hellip;
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\mathcal{L}_x</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
- \biggl[
B \cdot \frac{\partial A}{\partial x}
-
A \cdot \frac{\partial B}{\partial x}
\biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~8x \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2}{x} \biggl[ 2x \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)\biggr]^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~8x~\mathfrak{F}_x(y,z)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\mathcal{L}_y</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
- \biggl[
B \cdot \frac{\partial A}{\partial y}
-
A \cdot \frac{\partial B}{\partial y}
\biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
8q^4y\biggl( x^2 + \frac{p^4z^2}{2} \biggr)^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{2}{y} \biggl[ 2q^2y\biggl( x^2 + \frac{p^4z^2}{2} \biggr)\biggr]^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~8y~\mathfrak{F}_y(x,z)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\mathcal{L}_z</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
-\biggl[
B \cdot \frac{\partial A}{\partial z}
-
A \cdot \frac{\partial B}{\partial z}
\biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2p^4 z \biggl(x^2 - q^4y^2 \biggr)^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{2}{z} \biggl[p^2 z \biggl(x^2 - q^4y^2 \biggr)\biggr]^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~8z~\mathfrak{F}_z(x,y)
</math>
  </td>
</tr>
</table>
With this shorthand in place, we can write,
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\hat{e}_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}}{\mathcal{D}}
\biggl\{
-\hat\imath \biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]
+ \hat\jmath \biggl[ q^2 y(p^4z^2 + 2x^2 )  \biggl]
+ \hat{k} \biggl[ p^2z( x^2 - q^4y^2 ) \biggl]
\biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{(AB)^{1 / 2}}
\biggl\{
-\hat\imath \biggl[ \frac{x \mathcal{L}_x}{2} \biggl]^{1 / 2}
+ \hat\jmath \biggl[ \frac{y \mathcal{L}_y}{2} \biggl]^{1 / 2}
+ \hat{k} \biggl[ \frac{z \mathcal{L}_z}{2} \biggl]^{1 / 2}
\biggr\}
\, .
</math>
  </td>
</tr>
</table>
 
We therefore also recognize that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~h_3 \biggl(\frac{\partial \lambda_3}{\partial x}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{1}{(AB)^{1 / 2}} \biggl[ \frac{x \mathcal{L}_x}{2} \biggl]^{1 / 2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\frac{1}{(AB)^{1 / 2}} \biggl[ 4x^2 ~\mathfrak{F}_x \biggl]^{1 / 2} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~h_3 \biggl(\frac{\partial \lambda_3}{\partial y}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{(AB)^{1 / 2}} \biggl[ \frac{y \mathcal{L}_y}{2} \biggl]^{1 / 2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{(AB)^{1 / 2}} \biggl[ 4y^2~\mathfrak{F}_y \biggl]^{1 / 2} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~h_3 \biggl(\frac{\partial \lambda_3}{\partial z}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{(AB)^{1 / 2}} \biggl[ \frac{z \mathcal{L}_z}{2} \biggl]^{1 / 2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{(AB)^{1 / 2}} \biggl[ 4z^2~\mathfrak{F}_z \biggl]^{1 / 2} \, .
</math>
  </td>
</tr>
</table>
 
Now, if &#8212; and it is a BIG "if" &#8212; <math>~h_3 = h_0(AB)^{-1 / 2}</math>, then we have,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~h_0 \biggl(\frac{\partial \lambda_3}{\partial x}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-\biggl[ 4x^2 ~\mathfrak{F}_x \biggl]^{1 / 2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-2x \biggl[ \mathfrak{F}_x \biggl]^{1 / 2} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~h_0 \biggl(\frac{\partial \lambda_3}{\partial y}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ 4y^2~\mathfrak{F}_y \biggl]^{1 / 2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2y \biggl[ \mathfrak{F}_y \biggl]^{1 / 2} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~h_0 \biggl(\frac{\partial \lambda_3}{\partial z}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ 4z^2~\mathfrak{F}_z \biggl]^{1 / 2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2z\biggl[ \mathfrak{F}_z \biggl]^{1 / 2} \, ,
</math>
  </td>
</tr>
</table>
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Rightarrow~~~ h_0 \lambda_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-x^2 \biggl[ \mathfrak{F}_x \biggl]^{1 / 2}
+ y^2 \biggl[ \mathfrak{F}_y \biggl]^{1 / 2}
+ z^2 \biggl[ \mathfrak{F}_z \biggl]^{1 / 2} \, .
</math>
  </td>
</tr>
</table>
But if this is the correct expression for <math>~\lambda_3</math> and its three partial derivatives, then it must be true 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>~
\biggl(\frac{\partial \lambda_3}{\partial x}\biggr)^2
+
\biggl(\frac{\partial \lambda_3}{\partial y}\biggr)^2
+
\biggl(\frac{\partial \lambda_3}{\partial z}\biggr)^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \biggl( \frac{h_3}{h_0}\biggr)^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
4x^2 \biggl[ \mathfrak{F}_x \biggl]
+
4y^2 \biggl[ \mathfrak{F}_y \biggl]
+
4z^2\biggl[ \mathfrak{F}_z \biggl]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
4x^2 \biggl[ \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)^2 \biggl]
+
4y^2 \biggl[ q^4\biggl( x^2 + \frac{p^4z^2}{2} \biggr)^2\biggl]
+
4z^2\biggl[ \frac{p^4}{4}\biggl(x^2 - q^4y^2 \biggr)^2 \biggl]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x^2 (2q^4y^2  + p^4z^2 )^2
+
q^4 y^2( 2x^2 + p^4z^2 )^2
+
p^4 z^2 (x^2 - q^4y^2 )^2
</math>
  </td>
</tr>
</table>
 
Well &hellip; the right-hand side of this expression is identical to the right-hand side of the [[#Eureka|above expression]], where we showed that it equals <math>~(\ell_{3D}/\mathcal{D})^{-2}</math>.  That is to say, we are now showing that,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\biggl( \frac{h_3}{h_0}\biggr)^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2} = [AB]</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{h_3}{h_0}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(AB)^{-1 / 2} \, .</math>
  </td>
</tr>
</table>
And this is ''precisely'' what, just a few lines above, we hypothesized the functional expression for <math>~h_3</math> ought to be.  <font color="red">'''EUREKA!'''</font>
 
====Summary====
 
In summary, then &hellip;
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\lambda_3</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
-x^2 \biggl[ \mathfrak{F}_x \biggl]^{1 / 2}
+ y^2 \biggl[ \mathfrak{F}_y \biggl]^{1 / 2}
+ z^2 \biggl[ \mathfrak{F}_z \biggl]^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-x^2 \biggl[\biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)^2 \biggl]^{1 / 2}
+ y^2 \biggl[ q^4\biggl( x^2 + \frac{p^4z^2}{2} \biggr)^2\biggl]^{1 / 2}
+ z^2 \biggl[ \frac{p^4}{4} \biggl(x^2 - q^4y^2 \biggr)^2 \biggl]^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-x^2 \biggl(q^4y^2  + \frac{p^4z^2}{2} \biggr)
+ q^2y^2 \biggl( x^2 + \frac{p^4z^2}{2} \biggr)
+ \frac{p^2 z^2}{2}  \biggl(x^2 - q^4y^2 \biggr) \, ,
</math>
  </td>
</tr>
</table>
and,
 
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~h_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(AB)^{-1 / 2} </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[
(x^2 + q^4y^2 + p^4 z^2 )(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)
\biggr]^{-1 / 2} \, .</math>
  </td>
</tr>
</table>
 
No!  Once again this does not work.  The direction cosines &#8212; and, hence, the components of the <math>~\hat{e}_3</math> unit vector &#8212; are not correct!
 
 
===Speculation7===
 
 
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" colspan="9">'''Direction Cosine Components for T6 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}}{ z^{2/p^2}}</math></td>
  <td align="center"><math>~\frac{1}{\lambda_2}\biggl[\frac{x q^2 y p^2 z}{ \mathcal{D}}\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>~-\frac{2\lambda_2}{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>~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{\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>~\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>
 
<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>
  </td>
</tr>
</table>
 
 
On my white-board I have shown that, if
<div align="center">
<math>~\lambda_3 \equiv \ell_{3D} \mathcal{D} \, ,</math>
</div>
then everything will work out as long as,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\mathcal{L}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{\mathcal{D}}{\ell_{3D}} \biggr)^2 \frac{1}{\ell_{3D}^4} \, ,
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\mathcal{L}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
x^2 (2q^4 y^2 + p^4z^2 )^4
+
q^8 y^2 (2x^2 + p^4 z^2)^4
+
p^8z^2( x^2 - q^4y^2)^4
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x^2 (4q^8 y^4 + 4q^4y^2p^4z^2 + p^8z^4 )^2
+
q^8 y^2 (4x^4 + 4x^2p^4z^2 + p^8 z^4)^2
+
p^8z^2( x^4 - 2x^2q^4y^2 + q^8y^4)^2
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x^2
[16q^{16}y^8 + 16q^{12}y^6p^4z^2 + 4q^8y^4p^8z^4 + 16q^{12}y^6p^4z^2 + 16q^8y^4p^8z^4 + 4q^4y^2p^{12}z^6 + 4q^8y^4p^8z^4 + 4q^4y^2 p^{12}z^6 + p^{16}z^8]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+
q^8 y^2
[16x^8 + 16x^6p^4z^2 + 4x^4p^8z^4 + 16x^6p^4z^2 + 16x^4p^8z^4 + 4x^2p^{12}z^6 + 4x^4p^8z^4 + 4x^2p^{12}z^6 + p^{16}z^8]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+
p^8z^2
[x^8 - 2x^6q^4y^2 + x^4q^8y^4 - 2x^6q^4y^2 + 4x^4q^8y^4 - 2x^2q^{12}y^6 + x^4q^8y^4 - 2x^2q^{12}y^6 + q^{16}y^8]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x^2
[16q^{16}y^8 + 32q^{12}y^6p^4z^2 + 24q^8y^4p^8z^4 + 8q^4y^2p^{12}z^6 + p^{16}z^8]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+
q^8 y^2
[16x^8 + 32x^6p^4z^2 + 24x^4p^8z^4  + 8x^2p^{12}z^6 + p^{16}z^8]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+
p^8z^2
[x^8 - 4x^6q^4y^2 + 6x^4q^8y^4 - 4x^2q^{12}y^6 + q^{16}y^8]
</math>
  </td>
</tr>
</table>
 
Let's check this out.
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\mathrm{RHS}~\equiv \biggl(\frac{\mathcal{D}}{\ell_{3D}} \biggr)^2 \frac{1}{\ell_{3D}^4}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)(x^2 + q^4y^2 + p^4 z^2 )^3</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(x^2 + q^4y^2 + p^4 z^2 )(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)[x^4 + 2x^2q^4y^2 + 2x^2p^4z^2 +q^8y^4 + 2q^4y^2p^4z^2 + p^8z^4]</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[(6x^2q^4 y^2 p^4 z^2 + x^4 p^4 z^2 + 4x^4q^4y^2)
+
(q^8 y^4 p^4 z^2 + 4x^2q^8 y^4)
+
(q^4 y^2 p^8 z^4 + x^2 p^8 z^4 )]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~\times [x^4 + 2x^2q^4y^2 + 2x^2p^4z^2 +q^8y^4 + 2q^4y^2p^4z^2 + p^8z^4]</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[6x^2q^4 y^2 p^4 z^2 + x^4(p^4 z^2 + 4q^4y^2)
+
q^8 y^4(p^4 z^2 + 4x^2)
+
p^8 z^4 (q^4 y^2 + x^2 )]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~\times [x^4 + 2x^2q^4y^2 + 2x^2p^4z^2 +q^8y^4 + 2q^4y^2p^4z^2 + p^8z^4]</math>
  </td>
</tr>
</table>


==Best Thus Far==
==Best Thus Far==

Latest revision as of 18:19, 23 July 2021

Concentric Ellipsoidal (T6) Coordinates (Part 3)

Best Thus Far

Part A

Direction Cosine Components for T6 Coordinates
n λn hn λnx λny λnz γn1 γn2 γn3
1 (x2+q2y2+p2z2)1/2 λ13D xλ1 q2yλ1 p2zλ1 (x)3D (q2y)3D (p2z)3D
2 xy1/q2z2/p2 1λ2[xq2yp2z𝒟] λ2x λ2q2y 2λ2p2z xq2yp2z𝒟(1x) xq2yp2z𝒟(1q2y) xq2yp2z𝒟(2p2z)
3 --- --- --- --- --- 3D𝒟[x(2q4y2+p4z2)] 3D𝒟[q2y(p4z2+2x2)] 3D𝒟[p2z(x2q4y2)]

3D

(x2+q4y2+p4z2)1/2,

𝒟

(q4y2p4z2+x2p4z2+4x2q4y2)1/2.

 

A3D2

=

(x2+q4y2+p4z2),

      and,      

B𝒟2

=

(q4y2p4z2+x2p4z2+4x2q4y2)

Ax

=

2x,

Ay

=

2q4y,

Az

=

2p4z;

Bx

=

2x(4q4y2+p4z2),

By

=

2q4y(p4z2+4x2),

Bz

=

2p4z(q4y2+x2).

Try …

λ3

(BA)m/2=(D3D)m

[ABmλ3]λ3xi

12{ABxiBAxi}

[ABm]lnλ3lnxi

xi2{ABxiBAxi}.

In this case we find,

x2{}x

=

x2(2q4y2+p4z2)2,

y2{}y

=

q4y2(2x2+p4z2)2,

z2{}z

=

p4z2(x2q4y2)2.

The scale factor is, then,

h32

=

i=13(λ3xi)2

 

=

i=13{[mλ32AB][ABxiBAxi]}2

 

=

[mλ3AB]2{[x(2q4y2+p4z2)2]2+[q4y(2x2+p4z2)2]2+[p4z(x2q4y2)2]2}

h3

=

[ABmλ3]{[x(2q4y2+p4z2)2]2+[q4y(2x2+p4z2)2]2+[p4z(x2q4y2)2]2}1/2.

Part B (25 February 2021)

Now, from above, we know that,

(𝒟3D)2=AB

=

[x(2q4y2+p4z2)]2+[q2y(p4z2+2x2)]2+[p2z(x2q4y2)]2.

Example:
(q2,p2)=(2,5)     and     (x,y,z)=(0.7,0.23,0.1)λ1=1.0

[x(2q4y2+p4z2)]2 [q2y(p4z2+2x2)]2 [p2z(x2q4y2)]2 (𝒟3D)2
2.14037 1.39187 0.04623 3.57847

As an aside, note that,

AB

=

[ABm]lnλ3lnx+[ABm]lnλ3lny+[ABm]lnλ3lnz

m

=

lnλ3lnx+lnλ3lny+lnλ3lnz.

We realize that this ratio of lengths may also be written in the form,

(𝒟3D)2

=

6x2q4y2p4z2+x4(4q4y2+p4z2)+q8y4(4x2+p4z2)+p8z4(x2+q4y2).

Same Example, but Different Expression:
(q2,p2)=(2,5)     and     (x,y,z)=(0.7,0.23,0.1)λ1=1.0

6x2q4y2p4z2 x4(4q4y2+p4z2) q8y4(4x2+p4z2) p8z4(x2+q4y2) (𝒟3D)2
0.67620 0.94359 1.87054 0.08813 3.57847


Let's try …

λ5x

=

x(2q4y2+p4z2),

λ5y

=

q2y(p4z2+2x2),

λ5z

=

p2z(x2q4y2).

This means that the relevant scale factor is,

h52

=

[x(2q4y2+p4z2)]2+[q2y(p4z2+2x2)]2+[p2z(x2q4y2)]2=(𝒟3D)2

h5

=

(3D𝒟),

and the three associated direction cosines are,

γ51=h5(λ5x)

=

x(2q4y2+p4z2)(3D𝒟),

γ52=h5(λ5y)

=

q2y(p4z2+2x2)(3D𝒟),

γ53=h5(λ5z)

=

p2z(x2q4y2)(3D𝒟).

These direction cosines exactly match what is required in order to ensure that the coordinate, λ5, is everywhere orthogonal to both λ1 and λ4. GREAT! The resulting summary table is, therefore:

Direction Cosine Components for T10 Coordinates
n λn hn λnx λny λnz γn1 γn2 γn3
1 (x2+q2y2+p2z2)1/2 λ13D xλ1 q2yλ1 p2zλ1 (x)3D (q2y)3D (p2z)3D
4 xy1/q2z2/p2 1λ2[xq2yp2z𝒟] λ2x λ2q2y 2λ2p2z xq2yp2z𝒟(1x) xq2yp2z𝒟(1q2y) xq2yp2z𝒟(2p2z)
5 --- 3D𝒟 x(2q4y2+p4z2) q2y(p4z2+2x2) p2z(x2q4y2) 3D𝒟[x(2q4y2+p4z2)] 3D𝒟[q2y(p4z2+2x2)] 3D𝒟[p2z(x2q4y2)]

3D

(x2+q4y2+p4z2)1/2,

𝒟

(q4y2p4z2+x2p4z2+4x2q4y2)1/2.


Try …

λ5

=

x2q4y2q2+yq2p4zq4p2+xp4zp2

 

=

y2q2x2q4+yq2p4zq4p2+zp4xp2

 

=

1x2q4+p2zq4p2{[xp2]y2q2[zq4p2]+[x2q4+p2]yq2p4+[x2q4]zp4+q4p2}

 

=

1x2q4+p2zq4p2{F5}.

This gives,

λ5x

=

2q4x(y2q2x2q4)p4x(zp2xp2)

 

=

1x2q4+p2+1[2q4y2q2xp2+p4x2q4zp2].

Or, given that,

x2q4+p2+1

=

xzq4p2{F5λ5},

we can also write,

λ5x

=

zq4p2x{λ5F5}[2q4y2q2xp2+p4x2q4zp2]

Similarly,

λ5y

=

2q2y(y2q2x2q4)+q2p4y(yq2p4zq4p2)

 

=

1x2q4zq4p2[2q2y(y2q2zq4p2)+q2p4y(yq2p4x2q4)]

 

=

xp2y{λ5F5}[2q2(y2q2zq4p2)+q2p4(yq2p4x2q4)];

λ5z

=

q4p2z(yq2p4zq4p2)+p4z(zp4xp2)

 

=

1xp2zq4p2[p4z(zp4+q4p2)q4p2z(xp2yq2p4)]

 

=

x2q4z{λ5F5}[p4(zp4+q4p2)q4p2(xp2yq2p4)]

Understanding the Volume Element

Let's see if the expression for the volume element makes sense; that is, does

(h1h4h5)dλ1dλ4dλ5

=

dxdydz?

First, let's make sure that we understand how to relate the components of the Cartesian line element with the components of our T10 coordinates.

Line Element

MF53 claim that the following relation gives the various expressions for the scale factors; we will go ahead and incorporate the expectation that, since our coordinate system is orthogonal, the off-diagonal elements are zero.

ds2=dx2+dy2+dz2

=

i=1,4,5hi2dλi2.

Let's see. The first term on the RHS is,

h12dλ12

=

h12[(λ1x)dx+(λ1y)dy+(λ1z)dz]2

 

=

h12[(λ1x)2dx2+(λ1y)2dy2+(λ1z)2dz2

 

 

+2(λ1x)(λ1y)dxdy+2(λ1x)(λ1z)dxdz+2(λ1x)(λ1y)dydz];

the other two terms assume easily deduced, similar forms. When put together and after regrouping terms, we can write,

i=1,4,5hi2dλi2

=

[h12(λ1x)2+h42(λ4x)2+h52(λ5x)2]dx2

 

 

+[h12(λ1y)2+h42(λ4y)2+h52(λ5y)2]dy2

 

 

+[h12(λ1z)2+h42(λ4z)2+h52(λ5z)2]dz2.

Given that this summation should also equal the square of the Cartesian line element, (dx2+dy2+dz2), we conclude that the three terms enclosed inside each of the pair of brackets must sum to unity. Specifically, from the coefficient of dx2, we can write,

h12(λ1x)2

=

1h42(λ4x)2h52(λ5x)2.

Using this relation to replace h12 in each of the other two bracketed expressions, we find for the coefficients of dy2 and dz2, respectively,

[1h42(λ4x)2h52(λ5x)2](λ1y)2

=

[1h42(λ4y)2h52(λ5y)2](λ1x)2;

[1h42(λ4x)2h52(λ5x)2](λ1z)2

=

[1h42(λ4z)2h52(λ5z)2](λ1x)2.

We can use the first of these two expressions to solve for h42 in terms of h52, namely,

(λ1y)2h42(λ4x)2(λ1y)2h52(λ5x)2(λ1y)2

=

(λ1x)2h42(λ4y)2(λ1x)2h52(λ5y)2(λ1x)2

h42[(λ4y)2(λ1x)2(λ4x)2(λ1y)2]

=

(λ1x)2(λ1y)2+h52(λ5x)2(λ1y)2h52(λ5y)2(λ1x)2

Analogously, the second of these two expressions gives,

h42[(λ4z)2(λ1x)2(λ4x)2(λ1z)2]

=

(λ1x)2(λ1z)2+h52(λ5x)2(λ1z)2h52(λ5z)2(λ1x)2

Eliminating h4 between the two gives the desired overall expression for h5, namely,

0

=

[(λ1x)2(λ1z)2+h52(λ5x)2(λ1z)2h52(λ5z)2(λ1x)2][(λ4y)2(λ1x)2(λ4x)2(λ1y)2]

 

 

[(λ1x)2(λ1y)2+h52(λ5x)2(λ1y)2h52(λ5y)2(λ1x)2][(λ4z)2(λ1x)2(λ4x)2(λ1z)2]

 

=

h52{[(λ5x)2(λ1y)2(λ5y)2(λ1x)2][(λ4z)2(λ1x)2(λ4x)2(λ1z)2]

 

 

[(λ5x)2(λ1z)2(λ5z)2(λ1x)2][(λ4y)2(λ1x)2(λ4x)2(λ1y)2]}

 

 

[(λ1x)2(λ1z)2][(λ4y)2(λ1x)2(λ4x)2(λ1y)2]

 

 

+[(λ1x)2(λ1y)2][(λ4z)2(λ1x)2(λ4x)2(λ1z)2]

 

=

h52h14h42h52{[γ512γ122γ522γ112][γ432γ112γ412γ132][γ512γ132γ532γ112][γ422γ112γ412γ122]}

 

 

+1h42h12{(λ1x)2[γ432γ112γ412γ132](λ1y)2[γ432γ112γ412γ132](λ1x)2[γ422γ112γ412γ122]+(λ1z)2[γ422γ112γ412γ122]}

 

=

h52h14h42h52{[γ512γ122γ522γ112][γ43γ11+γ41γ13][γ43γ11γ41γ13][γ512γ132γ532γ112][γ42γ11+γ41γ12][γ42γ11γ41γ12]}

 

 

+1h42h12{(λ1x)2[γ43γ11+γ41γ13][γ43γ11γ41γ13](λ1y)2[γ43γ11+γ41γ13][γ43γ11γ41γ13](λ1x)2[γ42γ11+γ41γ12][γ42γ11γ41γ12]+(λ1z)2[γ42γ11+γ41γ12][γ42γ11γ41γ12]}

 

=

1h14h42{[γ51γ12+γ52γ11]γ43[γ43γ11+γ41γ13]γ52+[γ51γ13+γ53γ11]γ42[γ42γ11+γ41γ12]γ53}

 

 

+1h14h42{γ122[γ43γ11+γ41γ13]γ52γ112[γ43γ11+γ41γ13]γ52γ112[γ42γ11+γ41γ12]γ53+γ132[γ42γ11+γ41γ12]γ53}

 

=

1h14h42{[(γ51γ12γ52γ11)γ43+γ122γ112](γ43γ11+γ41γ13)γ52+[(γ51γ13+γ53γ11)γ42γ112+γ132](γ42γ11+γ41γ12)γ53}

… Not sure this is headed anywhere useful!

Volume Element

(h1h4h5)dλ1dλ4dλ5

=

(h1h4h5)[(λ1x)dx+(λ1y)dy+(λ1z)dz][(λ4x)dx+(λ4y)dy+(λ4z)dz][(λ5x)dx+(λ5y)dy+(λ5z)dz]

 

=

[(γ11)dx+(γ12)dy+(γ13)dz][(γ41)dx+(γ42)dy+(γ43)dz][(γ51)dx+(γ52)dy+(γ53)dz]

COLLADA

Here we try to use the 3D-visualization capabilities of COLLADA to test whether or not the three coordinates associated with the T6 Coordinate system are indeed orthogonal to one another. We begin by making a copy of the Inertial17.dae text file, which we obtain from an accompanying discussion. When viewed with the Mac's Preview application, this group of COLLADA-based instructions displays a purple ellipsoid with axis ratios, (b/a, c/a) = (0.41, 0.385). This means that we are dealing with an ellipsoid for which,

qab

=

2.44

      and,      

pac

=

2.60.

λ1

(x2+q2y2+p2z2)1/2;

λ2

=

xy1/q2z2/p2;

3D

[x2+q4y2+p4z2]1/2;

𝒟

(q4y2p4z2+x2p4z2+4x2q4y2)1/2.

First Trial

First Trial
(specified variable values have bgcolor="pink")
x y z λ1 3D 𝒟
0.5 0.35493 0.00000 1 0.46052 2.11310

Unit Vectors

e^1

=

ı^(x03D)+ȷ^(q2y03D)+k^(p2z03D)

 

=

ı^(0.23026)+ȷ^(0.97313);

e^2

=

x0q2y0p2z0𝒟{ı^(1x0)+ȷ^(1q2y0)+k^(2p2z0)}

 

=

k^(2x0q2y0𝒟)

 

=

k^(1) ;

e^3

=

3D𝒟{ı^[x0(2q4y02+p4z02)]+ȷ^[q2y0(p4z02+2x02)]+k^[p2z0(x02q4y02)]}

 

=

2q2x0y03D𝒟{ı^(q2y0)+ȷ^(x0)}=(1)3D{ı^(q2y0)+ȷ^(x0)}

 

=

ı^(0.97313)+ȷ^(0.23026).

Tangent Plane

From our above derivation, the plane that is tangent to the ellipsoid's surface at (x0,y0,z0) is given by the expression,

xx0+q2yy0+p2zz0

=

(λ12)0.

For this First Trial, we have (for all values of z, given that z0=0) …

(0.5)x+(2.11310)y

=

1

y

=

(10.5x)2.11310.

So let's plot a segment of the tangent plane whose four corners are given by the coordinates,

Corner x y z
A x_0 - 0.25 = +0.25 0.41408 -0.25
B x_0 + 0.25 = +0.75 0.29577 -0.25
C x_0 - 0.25 = +0.25 0.41408 +0.25
D x_0 + 0.25 = +0.75 0.29577 +0.25

Now, in order to give some thickness to this tangent-plane, let's adjust the four corner locations by a distance of ±0.1 in the e^1 direction.

Eight Corners of Tangent Plane

Corner 1: Shift surface-point location (x0,y0,z0) by (+Δe1) in the e^1 direction, by (+Δe2) in the e^2 direction, and by by (+Δe3) in the e^3 direction. This gives …

x1

=

x0+(Δe1)0.23026(Δe2)0.97313

Second Trial

Second Trial(q=2.44,p=2.60)
[specified variable values have bgcolor="pink"]
x_0 y_0 z_0 λ1 3D 𝒟
0.5 0.35493 0.00000 1 0.46052 2.11310

Generic Unit Vector Expressions

Let's adopt the notation,

e^i

=

ı^[eix]+ȷ^[eiy]+k^[eiz]

      for,       i=1,3.

Then, for the T6 Coordinate system, we have,

e1x

x03D;

     

e1y

q2y03D;

     

e1z

p2z03D;

e2x

q2y0p2z0𝒟;

     

e2y

x0p2z0𝒟;

     

e2z

2x0q2y0𝒟;

e3x

x0(2q4y02+p4z02)3D𝒟;

     

e3y

q2y0(p4z02+2x02)3D𝒟;

     

e3z

p2z0(x02q4y02)3D𝒟.


Second Trial
  x y z
e1 0.23026 0.97313 0.0
e2 0.0 0.0 -1.0
e3 - 0.97313 0.23026 0.0


What are the coordinates of the eight corners of a thin tangent-plane? Let's say that we want the plane to extend …

  • From (Δ1) to (+Δ1) in the e^1 direction … here we set Δ1=0.05;
  • From (Δ2) to (+Δ2) in the e^2 direction … here we set Δ2=0.25;
  • From (Δ3) to (+Δ3) in the e^3 direction … here we set Δ3=0.25.

Δx

=

Δ1e1x+Δ2e2x+Δ3e3x=0.23177;

Δy

=

Δ1e1y+Δ2e2y+Δ3e3y=+0.10622;

Δz

=

Δ1e1z+Δ2e2z+Δ3e3z=0.25000.

Tangent Plane Schematic
vertex x y z   x y z
0 x0|Δx| y0|Δy| z0|Δz| 0.26823 0.24871 -0.25
1 x0|Δx| y0+|Δy| z0|Δz| 0.26823 0.46115 -0.25
2 x0|Δx| y0|Δy| z0+|Δz| 0.26823 0.24871 0.25
3 x0|Δx| y0+|Δy| z0+|Δz| 0.26823 0.46115 0.25
4 x0+|Δx| y0|Δy| z0|Δz| 0.73177 0.24871 -0.25
5 x0+|Δx| y0+|Δy| z0|Δz| 0.73177 0.46115 -0.25
6 x0+|Δx| y0|Δy| z0+|Δz| 0.73177 0.24871 0.25
7 x0+|Δx| y0+|Δy| z0+|Δz| 0.73177 0.46115 0.25

In the figure on the left, vertex 0 is hidden from view behind the (yellow) solid rectangle.

Third Trial

GoodPlane01

Third Trial(q=2.44,p=2.60)
[specified variable values have bgcolor="pink"]
x_0 y_0 z_0 λ1 3D 𝒟
0.8 0.24600 0.00000 1 0.59959 2.34146

Again, for the T6 Coordinate system, we have,

e1x

x03D;

     

e1y

q2y03D;

     

e1z

p2z03D;

e2x

q2y0p2z0𝒟;

     

e2y

x0p2z0𝒟;

     

e2z

2x0q2y0𝒟;

e3x

x0(2q4y02+p4z02)3D𝒟;

     

e3y

q2y0(p4z02+2x02)3D𝒟;

     

e3z

p2z0(x02q4y02)3D𝒟.


Third Trial
  x y z ΔTP
e1 0.47967 0.87745 0.0 0.02
e2 0.0 0.0 -1.0 0.25
e3 - 0.87753 0.47952 0.0 0.25

In constructing the Tangent-Plane (TP) for a 3D COLLADA display, we first move from the point that is on the surface of the ellipsoid, x0=(x0,y0,z0)=(0.8,0.246,0.0), to

vertex
"m"
Pm Components
xm=ı^Pm ym=ȷ^Pm zm=k^Pm
0 x0Δ1e^1Δ2e^2Δ3e^3

x0Δ1e1xΔ2e2xΔ3e3x

y0Δ1e1yΔ2e2yΔ3e3y

z0Δ1e1zΔ2e2zΔ3e3z

   

0.8 - 0.02 (0.47952) - 0.25 (-0.87752) = 1.00979

0.24590 - 0.02 (0.87753) - 0.25 (0.47952) = 0.10847

- 0.25 (-1.0) = + 0.25

1 x0Δ1e^1+Δ2e^2Δ3e^3

x0Δ1e1x+Δ2e2xΔ3e3x

y0Δ1e1y+Δ2e2yΔ3e3y

z0Δ1e1z+Δ2e2zΔ3e3z

   

0.8 - 0.02 (0.47952) - 0.25 (-0.87752) = 1.00979

0.24590 - 0.02 (0.87753) - 0.25 (0.47952) = 0.10847

+ 0.25 (-1.0) = - 0.25

2 x0Δ1e^1Δ2e^2+Δ3e^3

x0Δ1e1xΔ2e2x+Δ3e3x

y0Δ1e1yΔ2e2y+Δ3e3y

z0Δ1e1zΔ2e2z+Δ3e3z

   

0.8 - 0.02 (0.47952) + 0.25 (-0.87752) = 0.56307

0.24590 - 0.02 (0.87753) + 0.25 (0.47952) = 0.34823

- 0.25 (-1.0) = + 0.25

3 x0Δ1e^1+Δ2e^2+Δ3e^3

x0Δ1e1x+Δ2e2x+Δ3e3x

y0Δ1e1y+Δ2e2y+Δ3e3y

z0Δ1e1z+Δ2e2z+Δ3e3z

   

0.8 - 0.02 (0.47952) + 0.25 (-0.87752) = 0.57103

0.24590 - 0.02 (0.87753) + 0.25 (0.47952) = 0.34823

+ 0.25 (-1.0) = - 0.25

4  

0.8 + 0.02 (0.47952) - 0.25 (-0.87752) = 1.0290

0.24590 + 0.02 (0.87753) - 0.25 (0.47952) = 0.1436

- 0.25 (-1.0) = + 0.25

5  

0.8 + 0.02 (0.47952) - 0.25 (-0.87752) = 1.0290

0.24590 + 0.02 (0.87753) - 0.25 (0.47952) = 0.1436

+ 0.25 (-1.0) = - 0.25

6  

0.8 + 0.02 (0.47952) + 0.25 (-0.87752) = 0.59021

0.24590 + 0.02 (0.87753) + 0.25 (0.47952) = 0.38333

- 0.25 (-1.0) = + 0.25

7  

0.8 + 0.02 (0.47952) + 0.25 (-0.87752) = 0.59021

0.24590 + 0.02 (0.87753) + 0.25 (0.47952) = 0.38333

+ 0.25 (-1.0) = - 0.25


Tangent Plane Schematic Vertex Locations via Excel
x0=0.8,z0=0.0,y0=0.246,λ1=1.0
Tangent Plane Schematic
Δ1=0.02,Δ2=0.25,Δ3=0.25

GoodPlane02

x0=0.075,z0=0.0,y0=0.4089,λ1=1.0
Tangent Plane Schematic
Δ1=0.02,Δ2=0.25,Δ3=0.25

GoodPlane03

x0=0.25,z0=0.20,y0=0.33501,λ1=1.0
Tangent Plane Schematic Tangent Plane Schematic
Δ1=0.02,Δ2=0.25,Δ3=0.25 Δ1=0.02,Δ2=0.10,Δ3=0.25

CAPTION:   The image on the right differs from the image on the left in only one way — Δ2 = 0.1 instead of 0.25. It illustrates more clearly that the e^3 (longest) coordinate axis is not parallel to the z-axis when z00.

GoodPlane04

x0=0.25,z0=1/3,y0=0.1777,λ1=1.0
Tangent Plane Schematic
Δ1=0.02,Δ2=0.10,Δ3=0.25

Further Exploration

Let's set: x0=0.25,y0=0.33501,z0=0.2λ1=1.00000,λ2=0.33521.


qab

=

2.43972

      and,      

pac

=

2.5974.

λ1

(x2+q2y2+p2z2)1/2;

λ2

=

xy1/q2z2/p2;

3D

[x2+q4y2+p4z2]1/2;

𝒟

(q4y2p4z2+x2p4z2+4x2q4y2)1/2.

Next, let's examine the curve that results from varying z while λ1 and λ2 are held fixed. From the expression for λ2 we appreciate that,

x

=

λ2z2/p2y1/q2;

and from the expression for λ1 we have,

x2

=

λ12q2y2p2z2.

Hence, the relationship between y and z is,

[λ2z2/p2y1/q2]2

=

λ12q2y2p2z2

λ22z4/p2

=

y2/q2[λ12q2y2p2z2].


Alternatively,

y

=

[λ2z2/p2x]q2.

Hence, the relationship between x and z is,

x2

=

λ12p2z2q2[λ2z2/p2x]2q2

x2(q2+1)

=

x2q2[λ12p2z2]q2[λ2z2/p2]2q2

x2

=

{x2q2[λ12p2z2]q2[λ2z2/p2]2q2}1/(q2+1)

Here are some example values …

λ1=1    and,     λ2=0.33521
z 1st Solution   2nd Solution lambda_3 coordinate
y1 x1 y2 x2
0.01 0.407825695 0.0995168 - -
0.03 0.40481851 0.138 - -
0.04 0.40309223 0.1503934 - -
0.08 0.393779065 0.1854283 - -
0.12 0.37990705 0.2103761 - -
0.16 0.36067787 0.23111 1.04123×10-4 0.9095546
0.2 0.33500747 0.2500033 2.23778×10-4 0.85448
0.22 0.31923525 0.2592611 3.36653 ×10-4 0.82065
0.24 0.30106924 0.2686685 5.2327 ×10-4 0.78192
0.26 0.2799962 0.2784963 8.53243 ×10-4 0.73752
0.28 0.25521147 0.2891526 1.491545 ×10-3 0.68634
0.3 0.22530908 0.3013752 2.89262 ×10-3 0.62671
0.32 0.1873233 0.3168808 6.6223 ×10-3 0.55579
0.34 0.13149897 0.3423994 2.09221 ×10-2 0.46637
0.343 0.1191543 0.3490285 0.026458 0.4496
0.344 0.1145 0.3517 0.02880 0.4435
0.345 0.1093972 0.354688 0.03155965 0.4371186
0.3485 0.0847372 0.3713588 0.0480478 0.4085204

See Also


Tiled Menu

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