Editing 2DStructure/ToroidalCoordinates (section)

===Expression for the Axisymmetric Potential===
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999; hereafter CT99)] derive an expression for the Newtonian gravitational potential in terms of a ''Compact Cylindrical Green's Function'' expansion.  They show &#8212; see, for example, their equation (31) &#8212; that when expressed in terms of cylindrical coordinates, the potential at any meridional location, <math>\varpi = R_*</math> and <math>~Z = Z_*</math>, due to an axisymmetric mass distribution, <math>~\rho(\varpi, Z)</math>, is
<div align="center">
<math>
\Phi(R_*,Z_*) = - \frac{2Gq_0}{R_*^{1/2}}  ,
</math>
</div>
where,
<div align="center">
<math>
q_0 = \int_\Sigma \varpi^{1/2} Q_{-1/2}(\Chi) \rho(\varpi, Z) d\sigma,
</math>
</div>
<math>~d\sigma = d\varpi dZ</math> is a differential area element in the meridional plane, and the dimensionless argument (the modulus) of the special function, <math>~Q_{-1/2}</math>, is,
<div align="center">
<math>
\Chi \equiv \frac{R_*^2 + \varpi^2 + (Z_* - Z)^2}{2R_* \varpi} .
</math>
</div>
Next, following the lead of CT99, we note that according to the Abramowitz &amp; Stegun (1965),
<div align="center">
<math>Q_{-1/2}(\Chi) = \mu K(\mu) \, ,</math>
</div>
where, the function <math>~K(\mu)</math> is the complete elliptic integral of the first kind and, for our particular problem,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mu^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~2(1+\Chi)^{-1}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2\biggl[ 1+\frac{R_*^2 + \varpi^2 + (Z_* - Z)^2}{2R_* \varpi}\biggr]^{-1}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[\frac{4R_*\varpi}{(R_* + \varpi)^2 + (Z_* - Z)^2} \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>

Hence, we can write,
<div align="center">
<math>
q_0 = \int\int \varpi^{1/2} \mu K(\mu) \rho(\varpi, Z) d\varpi dZ \, .
</math>
</div>

As has been explained in [[Appendix/Ramblings/ToroidalCoordinates#Confirmation_Provided_by_Trova.2C_Hur.C3.A9_and_Hersant|an accompanying set of notes]], this is precisely the same expression for the gravitational potential that [http://adsabs.harvard.edu/abs/2012MNRAS.424.2635T A. Trova, J.-M. Hur&eacute; and F. Hersant (2012; MNRAS, 424, 2635)] used in their study of the potential of self-gravitating, axisymmetric discs.  

Our objective, here, is to examine whether or not it might be advantageous to transform this expression to one in which the double integral is performed on a toroidal, rather than a cylindrical, coordinate system.