[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: Spheres vs Toroids



Original poster: "Godfrey Loudner" <ggreen-at-gwtc-dot-net> 

Hello Antonio

The messy part of the exact formula is
f(x) = Sum[Q(n-1/2,x)/P(n-1/2,x), {n, 1, infinity}]
where x = (D-d)/d.
Set f(x,k) = Sum[Q(n-1/2,x)/P(n-1/2,x), {n, 1, k}], i.e.,
the sum of the first k terms. By playing around with integral
inequalities, one can show the below.

absolutevalue[f(x)-f(x,k)] < e whenever x > 1 and

k > ln[Pi(2x)^(1/2)/e(x-1)(x^2-1)^(1/4)]/ln[x].

This can be slightly improved if one admits having
Q(-1/2,x) into the estimate. The problem with my
estimates is that they may over shoot the number of
terms necessary to achieve a given accuracy. I don't
know how optimal the estimates are.

Let's take your D = 90cm and d = 30cm, then x = 2.
Now with e = 1/10^m with accuracy to m places in
mind.

place value accuracy     number of terms
          1                 6 or more
          2                 9 or more
          3                13 or more
          4                16 or more
          5                19 or more
          6                23 or more

 >From my notes, Mathematica gives 40.56477807 pf

However with x very close to 1 from above, my estimate
indicates that an astronomical number of terms might be
required. If D = 2d, we have a torus with no hole and x = 1.
I have tried to compute the limit of the formula as
x approaches 1 from above, but without success so far.

Godfrey Loudner







Example:
Major diameter: 90 cm
Minor diameter: 30 cm
I get: Exact capacitance: 40.5648 pF (all the digits correct)
The document lists:       40.46   pF

There is also a table in the documentation of the program Mandk,
but that table appears to use an approximate formula, the same
listed in Gary's paper.

Antonio Carlos M. de Queiroz