A fixed-point iteration method for the number Pi with arbitrary odd order of convergence based on the sine function

Abstract

In this paper, we present a fixed point method for high-precision computation of number $\ddot{I}€$ based on the sine function. Let $P\in \mathbb{N}$. We define the function: \[ S\left(x\right) =x+\sum_{k=1}^{P}\left(\prod_{\ell=1}^{k-1}\frac {2\,\ell-1}{2\,\ell}\right)\frac{\sin\left(x\right)^{2\,k-1}}{2\,k-1}\,. \] For every initial value $x_0$ sufficiently close to $\ddot{I}€$, the sequence \[x_{n+1}=x_n+S\left(x_{n}\right)\;;\,n=0,1,\ldots\] is converging to $\ddot{I}€$ with order of convergence exactly $\left(2P+1\right)$. The computational tests we performed demonstrate the efficiency of the method. \[\] \[\textbf{Zusammenfassung}\] In dieser Abhandlung stellen wir ein Fixpunktverfahren zur Berechnung der Kreiszahl $\ddot{I}€$ auf Basis der sinus Funktion vor. Es sei $P\in \mathbb{N}$. Wir definieren die Funktion: \[ S\left(x\right) =x+\sum_{k=1}^{P}\left(\prod_{\ell=1}^{k-1}\frac {2\,\ell-1}{2\,\ell}\right)\frac{\sin\left(x\right)^{2\,k-1}}{2\,k-1}\;. \] Für jeden Startwert $x_0$ hinreichend nahe bei $\ddot{I}€$ konvergiert die Folge \[x_{n+1}=x_n+S\left(x_{n}\right)\;;\,n=0,1,\ldots\] gegen $\ddot{I}€$ mit Konvergenzordnung genau $\left(2P+1\right)$. Anhand von praktischen Berechnungen zeigen wir die Effizienz des Verfahrens. \[\text{Deutsche Version ab Seite 19}\]

bilingual: English (18 pages) and German (17 pages)