Post by Marko van Hausereine Tastenkombination habe ich bis jetzt nicht gefunden, kann noch wer helfen ?
===========
Vorschlag: ein Algorithmus zu finden sqrt[p]{A} , A>0 ,
p natuerliche Zahl, p >= 2 .
Es sei H:(0,infty)--->(0,infty) mit
==================================
(p-1)*X^p +(p+1)*A
H(X)= X*-------------------- .
(p+1)*X^p +(p-1)*A
================================
Weiter 0< y_0 < sqrt[p]{A} und x_0 := (p-1+A)/p .
Die iterative Folgen (x_n) , (y_n) wobei
==================================================
(1) x_{k+1}=H(x_k) , y_{k+1}=H(y_k) , k=0,1,...,n,...,
==================================================
die folgende EIgenschaften haben :
0 < y_0< ...< y_{n} < y_{n+1}< ...< sqrt[p]{A}<...<
<x_{n+1}<x_n<...< x_0 .
----------------------------------------------------------
Hinweise:
1) fuer (A-x^p)(H(x)-x)>0 ,
2) H'(x) >= 0 fuer x>0 ,
3) Man weiss dass die Geometrische-mittel(:= Geom) einer Zahlen , z.B.
(*) a_1=a_2=...=a_{p-1}=1 , a_p=A ,
ist kleiner als die Arithmetische-mittel (:= Arithm) .
Fuer die Zahlen aus (*)
Geom=sqrt[p]{A} =< Arithm = (p-1+A)/p .
------------------------------------------------------------
SATZ. Es gibt eine positive Konstante r=r(p) so dass
|y_{n+1}-y_n| =< r*|y_n-y_{n-1}|^3 .
Genauer (y_n) wird gegen sqrt[p]{A} convergiert mit Ordnung 3.
BEWEIS/HINWEIS: wenn L:=sqrt[p]{A} dann
H(L)=L , H'(L)=H"(L)=0 , H^{(3)}(L)=(p^2-1)/(2L^2) .
Weiter verwenden Sie die Taylor'sche Entwicklung, z.B.
y_{n+1}-L =H(y_n) - H(L)=(y_n-L)H'(L)+(y_n-L)^2H"(L)/2!+
+ (y_n-L)^3*H^{(3)}(c)/3! mit
c zwischen y_n und L , oder
y_{n+1}-y_n=H(y_n)-H(y_{n-1})=......
Gruss, Alex.
=========