Letzter Kommentar: vor 12 Jahren3 Kommentare2 Personen sind an der Diskussion beteiligt
Hallo,
hab das Beispiel mal in c++ und Fortran programmiert. Bekomme die gleichen Nullstellen, aber in anderer Reihenfolge.
Meine Tabelle lautet:
It.-Nr
p
q
r
0
1.0000000000000000+0.0000000000000000*i
0.40000000000000002 +0.90000000000000002*i
-0.65000000000000002+0.71999999999999997*i
1
1.3607734793516519 +2.0222302921553132*i
-1.3982133295376751-0.69356635962504332*i
3.0374398501860242-1.3286639325302698*i
2
0.98096328371966823 +1.3474626910848719*i
-0.33525193260127306 -0.64406860772816377*i
2.3542886488816048-0.70339408335670806*i
3
0.31718054925650652+0.93649454851955727*i
0.49001572078718691-0.96614107903072544*i
2.1928037299563061+2.96465305111680566E-002*i
4
0.20901563897345732 +1.5727420147652911*i
4.12060386626913466E-002-1.5275192097633474*i
2.7497783223638517-4.52228050019438110E-002*i
5
0.21297050700971853+1.3948274731404162*i
0.18467846583682393-1.3845653821841168*i
2.6023510271534578-1.02620909562992635E-002*i
6
0.20653075193800677+1.3748787427714850*i
0.20600107336130224 -1.3746529207714702*i
2.5874681747006911-2.25822000014998941E-004*i
Hab ich was falsch verstanden, kann nummerisches Rauschen zum "Austausch" der Nullstellen führen, oder ist die Tabelle auf der Wiki-Seite falsch?
Mein Quelltext:
DK-Fraktal Was genau passiert ist, musst Du Bo Jacoby auf der englischen Seite fragen, ich habe die Tabelle aus dem englischen Artikel übernommen. An Deinem Code ist nichts zu meckern. Die nachprüfbaren Fakten sind auf Deiner Seite. Macht man das Experiment, die Konvergenz bei variablen Starttripeln farbig darzustellen, ergibt sich nebenstehende Grafik. Das Gitter ist ganzzahlig. Der Punkt z=1, um den es hier geht, sieht dabei sehr stabil aus, sicher entfernt von den numerisch instabilen fraktalen Rändern. Wenn Du Deine Tabelle in den Artikel übernimmst, bitte noch sinnvoll auf 4 oder 6 Nachkommastellen runden, ist übersichtlicher.--LutzL (Diskussion) Diskussion:Weierstra%C3%9F-(Durand-Kerner)-Verfahren#c-LutzL-2012-06-08T17:08:00.000Z-129.69.55.84-2012-06-08T14:14:00.000Z11Beantworten
Auf der englischen Diskussionsseite ist der Code zur Tabelle zu finden. Die Zahlen stammen von einer sequentiellen Implementierung a la Gauß-Seidel-Iteration, bei der die neu gefundene Approximation einer Nullstelle sofort in die Berechnung der darauffolgenden Nullstellen eingesetzt wird. Ändert man also den Code so ab, dass das Hilfsarray gestrichen und in Zuweisungen durch das Nullstellenarray ersetzt wird, dann ergeben sich die Werte aus der alten Tabelle.--LutzL (Diskussion) Diskussion:Weierstra%C3%9F-(Durand-Kerner)-Verfahren#c-LutzL-2012-06-13T11:57:00.000Z-129.69.55.84-2012-06-08T14:14:00.000Z11Beantworten