Self-replication and Borwein-like algorithms

Jesús Guillera

DOI: https://doi.org/10.48550/arXiv.1702.05378

2017-02-21

Abstract:Using a self-replicating method, we generalize with a free parameter some Borwein algorithms for the number $\pi$. This generalization includes values of the Gamma function like $\Gamma(1/3)$, $\Gamma(1/4)$ and of course $\Gamma(1/2)=\sqrt{\pi}$. In addition, we give new rapid algorithms for the perimeter of an ellipse.

Number Theory

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is how to use the self - replicating method to generalize some Borwein algorithms for calculating the value of pi ($\pi$) and other special functions, such as the gamma function $\Gamma(1/3)$, $\Gamma(1/4)$ and $\Gamma(1/2)=\sqrt{\pi}$. In addition, the paper also proposes a new fast algorithm for calculating the perimeter of an ellipse.

Specifically:

1. **Generalizing Borwein algorithms**:
- The author generalizes some classical Borwein algorithms to a broader form that can calculate special function values by introducing a free parameter.
- These generalizations include the calculation of the gamma function $\Gamma(1/3)$, $\Gamma(1/4)$ and $\Gamma(1/2)=\sqrt{\pi}$.

2. **Proposing new algorithms**:
- The paper proposes new quadratic, cubic and quartic iterative algorithms, which can quickly converge to $\pi$ and other special function values.
- At the same time, the author also gives a new fast algorithm for calculating the perimeter of an ellipse. These algorithms are based on the self - replicating method rather than the traditional arithmetic - geometric mean (AGM) method.

### Specific formulas and derivations

- **Landen transformation**:
\[
\sum_{k = 0}^{\infty}\frac{\left(\frac{1}{2}\right)_k^2}{(1)_k^2}x^{2k}=(1 + t)\sum_{k = 0}^{\infty}\frac{\left(\frac{1}{2}\right)_k^2}{(1)_k^2}t^{2k},\quad t=\frac{1-\sqrt{1 - x^2}}{1+\sqrt{1 - x^2}}
\]
This is a quadratic transformation because there is a relationship of $t\sim4^{-1}\cdot x^2$ between $t$ and $x$.

- **Cubic transformation**:
\[
\sum_{k = 0}^{\infty}\frac{\left(\frac{1}{3}\right)_k\left(\frac{2}{3}\right)_k}{(1)_k^2}x^{3k}=(1 + 2t)\sum_{k = 0}^{\infty}\frac{\left(\frac{1}{3}\right)_k\left(\frac{2}{3}\right)_k}{(1)_k^2}t^{3k},\quad t=\frac{1-\sqrt[3]{1 - x^3}}{1+2\sqrt[3]{1 - x^3}}
\]

- **Quartic transformation**:
\[
\sum_{k = 0}^{\infty}\frac{\left(\frac{1}{2}\right)_k^2}{(1)_k^2}x^{4k}=(1 + t)^2\sum_{k = 0}^{\infty}\frac{\left(\frac{1}{2}\right)_k^2}{(1)_k^2}t^{4k},\quad t=\frac{1-\sqrt[4]{1 - x^4}}{1+\sqrt[4]{1 - x^4}}
\]

### Algorithm for the perimeter of an ellipse

For the perimeter $P(a, b)$ of an ellipse, where $a$ and $b$ are the semi - major axis and semi - minor axis of the ellipse respectively, its perimeter can be represented by the following elliptic integral:
\[
P(a, b)=4a\int_{0}^{\frac{\pi}{2}}\sqrt{1 - x^2\cos^2\phi}\,d\phi,\quad x = \sqrt{1-\frac{b^2}{a^2}}
\]
Through the self - replicating method, the author derives quadratic and quartic algorithms for calculating the perimeter of an ellipse.

### Summary

This paper generalizes Borwein algorithms through the self - replicating method and proposes new fast algorithms for calculating $\pi$ and other special function values. At the same time, it also provides...