[soon]

Zitat von Struktron

... neulich einen Hinweis:
"in 2004 hans de vries posted a terse, elegant and (still, to date)
accurate formula for the fine structure constant:

http://tinyurl.com/devriesconst

a formula in python2.7 would be as follows:

from math import pi,e

a = 0.007 # start off arbitrarily close to alpha
for x in range(1, 15):
t = 0.0
g = 0.0
for i in range(70):
g = g + pow(a, i)/pow(2*pi, t)
t = t + i
a = pow(g, 2)/(pow(e, (pow(pi, 2)/2)))
return 1/a

[update: when running this with 30 loops and using python
BigFloat (libmpfr) at 150 decimal places of precision
the following values are obtained:
137.03599909582961049584812391855684195959053299
0.007297352568653858214986442807070539674568087313 6
thus demonstrating that the algorithm is still within the
margin of error for at least the CODATA 2010 value of alpha]".

Wenn ich ein Spaßvogel wäre, würde ich auch soetwas versuchen:

Ich nehme eine Naturkonstante, bastel eine Iteration, iteriere, bis ich einen lustigen Wert bekomme, bezeichne diesen als mysteriösen Startwert und iteriere rückwärts. Wenn dann die Naturkonstante nicht wieder herauskommt, erzähle ich etwas von 150 Nachkommastellen, die ich brauche, so daß niemand den Quatsch nachprüfen mag.

- Ist nur eine Vermutung, da ich auf andere Werte kam, nachdem ich den Quelltext in Delphi eingetippt habe.

(procedure TForm1.Button1Click(Sender: TObject);
var
a,t,g,pi,e : extended;
x,i : integer;
begin
a := 0.007;
pi := 3.141592653589793238;
e := 2.718281828459045235;

for x := 1 to 15 do
begin
t := 0;
g := 0;

for i := 0 to 70 do
begin
g := g + power(a,i) / power(2*pi,t);
t := t + 1;
a := power(g,2) / (power(e,(power(pi,2)/2)));
end;
a := 1/a;

memo1.lines.add(floattostrF(a,fffixed,100,100));
end;

Ergebnis:
138,727507278841451000
138,727507278841451000
138,727507278841451000
138,727507278841451000
...


LG soon