Borwein's Algorithm (Others)

Jonathan and Peter Borwein devised various algorithms to calculate the value of π. The most prominent and oft-used one is explained under Borwein's algorithm. Other algorithms found by them include the following:
  1. Cubical convergence, 1991:
    • Start out by setting
      a_0 = \frac{1}{3}
      s_0 = \frac{\sqrt{3} - 1}{2}
    • Then iterate
      r_{k+1} = \frac{3}{1 + 2(1-s_k^3)^{1/3}}
      s_{k+1} = \frac{r_{k+1} - 1}{2}
      a_{k+1} = r_{k+1} a_k - 3^k(r_{k+1}^2-1)
    Then ak converges cubically against 1/π; that is, each iteration approximately triples the number of correct digits.
  2. Quartical convergence, 1984:
    • Start out by setting
      a_0 = \sqrt{2}
      b_0 = 0
      p_0 = 2 + \sqrt{2}
    • Then iterate
      a_{n+1} = \frac{\sqrt{a_n} + 1/\sqrt{a_n}}{2}
      b_{n+1} = \frac{\sqrt{a_n} (1 + b_n)}{a_n + b_n}
      p_{n+1} = \frac{p_n b_{n+1} (1 + a_{n+1})}{1 + b_{n+1}}
    Then pk converges quartically against π; that is, each iteration approximately quadruples the number of correct digits. The algorithm is not self-correcting; each iteration must be performed with the desired number of correct digits of π.
  3. Quintical convergence:
    • Start out by setting
      a_0 = \frac{1}{2}
      s_0 = 5(\sqrt{5} - 2)
    • Then iterate
      x_{n+1} = \frac{5}{s_n} - 1
      y_{n+1} = (x_{n+1} - 1)^2 + 7
      z_{n+1} = \left(\frac{1}{2} x_{n+1}\left(y_{n+1} + \sqrt{y_{n+1}^2 - 4x_{n+1}^3}\right)\right)^{1/5}
      s_{m+1} = \frac{25}{(z_{n+1} + x_{n+1}/z_{n+1} + 1)^2 s_n}
      a_{n+1} = s_n^2 a_n - 5^n\left(\frac{s_n^2 - 5}{2} + \sqrt{s_n(s_n^2 - 2s_n + 5)}\right)
    Then ak converges quintically against 1/π (that is, each iteration approximately quintuples the number of correct digits), and the following condition holds:
    0 < a_n - \frac{1}{\pi} < 16\cdot 5^n\cdot e^{-5^n}\pi
  4. Nonical convergence:
    • Start out by setting
      a_0 = \frac{1}{3}
      r_0 = \frac{\sqrt{3} - 1}{2}
      s_0 = (1 - r_0^3)^{1/3}
    • Then iterate
      t_{n+1} = 1 + 2r_n
      u_{n+1} = (9r_n (1 + r_n + r_n^2))^{1/3}
      v_{n+1} = t_{n+1}^2 + t_{n+1}u_{n+1} + u_{n+1}^2
      w_{n+1} = \frac{27 (1 + s_n + s_n^2)}{v_{n+1}}
      a_{n+1} = w_{n+1}a_n + 3^{2n-1}(1-w_{n+1})
      s_{n+1} = \frac{(1 - r_n)^3}{(t_{n+1} + 2u_{n+1})v_{n+1}}
      r_{n+1} = (1 - s_{n+1}^3)^{1/3}
    Then ak converges nonically against 1/π; that is, each iteration approximately multiplies the number of correct digits by nine.

 

<< PreviousWord BrowserNext >>
canton of zug
john thomas
urban trad
aswini
parjanya
prithivi mata
misiones province
puntland
romsdal
johrei
dormancy
bartolomeo eustachi
eustachian tube
list of politics by country articles
gabriele falloppio
giovanni filippo ingrassias
jean fernel
ruble
east china normal university
pollination
muslims in britain
insecurity
tacit knowledge
rational unified process
phantom zone
serissa
cryoelectronics
list of korea related topics
eurovision song contest 2004
claudia beni
spherical earth
david bronstein
world geodetic system
chiton
peruvian apple cactus
history of geodesy
nach dir, herr, verlanget mich
figure of the earth
royal flying corps
list of greek phrases
rockwell kent
delaware and raritan canal
george of podebrady
paul guldin