SOME NEW IDENTITIES INVOLVING π,

Similar documents
Integral Representations and Binomial Coefficients

The Riemann Zeta Function

Zeros of Polynomials

A collection of mathematical formulas involving π

Interesting Series Associated with Central Binomial Coefficients, Catalan Numbers and Harmonic Numbers

arxiv: v1 [math.nt] 28 Apr 2014

Harmonic Number Identities Via Euler s Transform

Created by T. Madas SERIES. Created by T. Madas

In algebra one spends much time finding common denominators and thus simplifying rational expressions. For example:

INFINITE SEQUENCES AND SERIES

ON RUEHR S IDENTITIES

Evaluation of Some Non-trivial Integrals from Finite Products and Sums

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.

arxiv: v3 [math.nt] 24 Dec 2017

Formulas for the Approximation of the Complete Elliptic Integrals

Series with Central Binomial Coefficients, Catalan Numbers, and Harmonic Numbers

-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION

On the Inverse of a Certain Matrix Involving Binomial Coefficients

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

The log-behavior of n p(n) and n p(n)/n

Euler-type formulas. Badih Ghusayni. Department of Mathematics Faculty of Science-1 Lebanese University Hadath, Lebanon

Problem 4: Evaluate ( k ) by negating (actually un-negating) its upper index. Binomial coefficient

3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials

CentralBinomialCoefficients Steven Finch. March 28, A(n +1)=(n +1)A(n), A(0) = 1

Solutions to Homework 1

PROBLEM SET 5 SOLUTIONS 126 = , 37 = , 15 = , 7 = 7 1.

AMS Mathematics Subject Classification : 40A05, 40A99, 42A10. Key words and phrases : Harmonic series, Fourier series. 1.

Continuous Functions

MAT1026 Calculus II Basic Convergence Tests for Series

Curve Sketching Handout #5 Topic Interpretation Rational Functions

6.3 Testing Series With Positive Terms

Chapter 10: Power Series

On the Equivalence of Ramanujan s Partition Identities and a Connection with the Rogers Ramanujan Continued Fraction

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

1 6 = 1 6 = + Factorials and Euler s Gamma function

4x 2. (n+1) x 3 n+1. = lim. 4x 2 n+1 n3 n. n 4x 2 = lim = 3

Complex Analysis Spring 2001 Homework I Solution

Math 113, Calculus II Winter 2007 Final Exam Solutions

Bijective Proofs of Gould s and Rothe s Identities

Simultaneous Generation of Koecher and Almkvist-Granville s Apéry-Like Formulae

Math 132, Fall 2009 Exam 2: Solutions

18.440, March 9, Stirling s formula

On Divisibility concerning Binomial Coefficients

In this section, we show how to use the integral test to decide whether a series

Section 11.8: Power Series

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

Math 155 (Lecture 3)

Series III. Chapter Alternating Series

Math 142, Final Exam. 5/2/11.

Chapter 8. Euler s Gamma function

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Bernoulli numbers and the Euler-Maclaurin summation formula

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Ma 530 Infinite Series I

Q-BINOMIALS AND THE GREATEST COMMON DIVISOR. Keith R. Slavin 8474 SW Chevy Place, Beaverton, Oregon 97008, USA.

Carleton College, Winter 2017 Math 121, Practice Final Prof. Jones. Note: the exam will have a section of true-false questions, like the one below.

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.

SCORE. Exam 2. MA 114 Exam 2 Fall 2017

(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b)

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

Math 113 Exam 3 Practice

(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3

x x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

Sum of cubes: Old proofs suggest new q analogues

Math 128A: Homework 1 Solutions

RAMANUJAN AND PI JONATHAN M. BORWEIN

Chapter 6 Infinite Series

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Ma 530 Introduction to Power Series

Some remarks for codes and lattices over imaginary quadratic

1. Introduction. Notice that we have taken the terms of the harmonic series, and subtracted them from

CERTAIN GENERAL BINOMIAL-FIBONACCI SUMS

e to approximate (using 4

Section 5.5. Infinite Series: The Ratio Test

MATH 10550, EXAM 3 SOLUTIONS

GENERALIZED HARMONIC NUMBER IDENTITIES AND A RELATED MATRIX REPRESENTATION

*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully

Physics 116A Solutions to Homework Set #9 Winter 2012

SCORE. Exam 2. MA 114 Exam 2 Fall 2016

MA131 - Analysis 1. Workbook 9 Series III

Chapter 8. Euler s Gamma function

Infinite Sequences and Series

MATH301 Real Analysis (2008 Fall) Tutorial Note #7. k=1 f k (x) converges pointwise to S(x) on E if and

Roger Apéry's proof that zeta(3) is irrational

The Gamma function Michael Taylor. Abstract. This material is excerpted from 18 and Appendix J of [T].

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

TWO INEQUALITIES ON THE AREAL MAHLER MEASURE

A SOLUTION TO A TANTALIZING PROBLEM. GERT ALMKVIST Institute for Algebraic Meditation, PL 500, S HOOR, Sweden (Submitted November 1984)

COMM 602: Digital Signal Processing

k-generalized FIBONACCI NUMBERS CLOSE TO THE FORM 2 a + 3 b + 5 c 1. Introduction

The integrals in Gradshteyn and Ryzhik. Part 5: Some trigonometric integrals

[ 11 ] z of degree 2 as both degree 2 each. The degree of a polynomial in n variables is the maximum of the degrees of its terms.

COMPUTING THE EULER S CONSTANT: A HISTORICAL OVERVIEW OF ALGORITHMS AND RESULTS

Math 113 Exam 4 Practice

CALCULUS BASIC SUMMER REVIEW

Quadratic Transformations of Hypergeometric Function and Series with Harmonic Numbers

Math 299 Supplement: Real Analysis Nov 2013

NAME: ALGEBRA 350 BLOCK 7. Simplifying Radicals Packet PART 1: ROOTS

Transcription:

SOME NEW IDENTITIES INVOLVING π, HENG HUAT CHAN π AND π. Itroductio The umber π, as we all ow, is defied to be the legth of a circle of diameter. The first few estimates of π were 3 Egypt aroud 9 B.C., 3 8 Babyloia, ad 3 7 < π < 3 Archimedes. 7 Sice the, may mathematicias have devoted their time to the determiatio of the explicit value of π. I this tal, I will preset some recet ew idetities ivolvig π, /π ad /π, with the hope to share with you various mathematical cocepts.. The formulas of Vieta ad Wallis I 65, F. Vieta showed that π = + + + ad i 655, J. Wallis discovered that π = 3 3 5 4 4 5 7 6 6 7 9 8 8. At first glace, these two formulas do ot seem to be related. I 999, over 3 years later, T. Osler [8] succeeded i uitig these two idetities. Osler bega with the simple formula si x = cos x si x. = cos x cos x 4 si x 4 Usig a ot-so-simple formula si θ = θ θ π = = = p cos x cos x cos x p si x p. π θ = π π + θ, π

HENG HUAT CHAN with θ = x p, he deduced that Now, Hece, si x x si x x = p = cos x cos x cos x p p π x p π = cos x = + cos x, cos x = = p π + x p. π + cos x, + cos x p = + + + cos x, } {{ } p radicals + + + p cos x π x p π = }{{} p π + x p. π radicals Settig x = π/ we see that whe p = we have Wallis s formula ad whe p, we have Vieta s formula. 3. The formula of D. Bailey, P. Borwei ad S.Plouffe Aroud 996, usig PSLQ algorithm, D. Bailey, P. Borwei ad S. Plouffe discovered the idetity 4 π = 6 i 8i + 8i + 4 8i + 5. 8i + 6 i= This idetity ca be proved by first observig that x x 8 dx = x +8i dx = / i= Hece, 4 6 i 8i + 8i + 4 8i + 5 = 8i + 6 i= Substitutig y = x, we fid that 6y 6 y 4 y 3 + 4y 4 dy = i= 6 i 8i +. 4 8x 3 4 x 4 8x 5 x 8 dx. 4y y dy 4y 8 y dy = π. y + The BBP formula is very useful i the sese that oe ca determie the -th hexadecimal of π without owig the the first -th hexadecimal of the expasio of π.

SOME NEW IDENTITIES INVOLVING π, π AND π 3 The idea preseted i [3] appears to be useful i derivig -th decimal places of costats provided a formula for those costat exists with powers of i the deomiator of the series. So far, o formula of the type π = T = i i b i has bee foud. However, we all ow that l9/ has a expasio of the above type ad hece the -th decimal places of l9/ ca be foud without owig the -th decimal of l9/. We ed this sectio with the followig observatio: Let Let a = aa + a +. P N,t = 3 5 3 3 A = 45 + 87 3 3! 3 N+t = A. 5 The it appears that the 6N -th to 6N-th decimal places of /π ca be obtaied from the first sixth decimal places of { 6N P N, N }, 9 {x} deotes the fractioal part of x. The above observatio is derived from a series discovered recetly by Cha, W.C. Liaw ad V. Ta [5]. 4. Gosper s formula ad the wor of G. Almvist, C. Krattethaler ad J.Petersso Bill Gosper has a formula for π that taes the form π = 5 6 3. Recetly this formula was used by F. Bellard after the wor of S. Plouffe to compute the -th decimal places of π. The mai tool i their wor is the Chiese Remaider Theorem. Gosper s formula ca be proved as follow: It is well ow that the beta itegral ca be evaluated as Hece x m x dx = 3 = 3 + m!! m + +!. x x dx.

4 HENG HUAT CHAN Therefore, Sice 5 6 3 = x x 5 63 + dx. 5 63 + y = 56y + 97y 3 y 3, we coclude that 5 6 = 8 3 8x 6 56x 5 + 8x 4 97x 3 + 97x 6 x 3 x + 3 dx = π. The above method of proof prompted Almvist, Krattethaler ad Petersso [] to cosider the itegral = m + m p x p x m p dx. They the search for polyomials S of degree d i ad T y of degree d i y so that m + Sy = Lettig y = x p x m p /a, they deduced that S = m a p T y y d+. P x x p x m p dx, a d+ P x is a polyomial of degree md + i x. I order to sum to π oe has to put i restrictios so that x p x m p a has zero at i or + i. Usig LLL-algorithm the authors the costructed may ew series for π. It turs out that Gosper s formula remais to be the simplest of its id. The above experimet led Almvist, Krattethaler ad Petersso to the followig theorem: Theorem 4.. For all, there exist a formula π = S 8 4 4 S is a polyomial i of degree 4 with ratioal coefficiets.

SOME NEW IDENTITIES INVOLVING π, π AND π 5 The proof of this theorem requires very clever methods of computig determiats. A example of such idetity is π = 3 5 7 8 4 S 4 S = 8986 + 3875948 349734 + 47 3 5936 4. is 5. Ramauja s series for /π I [9], Ramauja recorded a total of 7 series for /π. Oe of these series 4 π = 6 + 3! 3 4. Such series coverge rapidly to /π ad ca be used to compute π as well. I fact, the Chudovsys discovered that 643 3/ π = 6!! 3 3! 35949 + 545434 643 3 ad used it to compute over a billio digits of π. The Chudovsys formula exists i the above form because the rig of algebraic itegers of Q 63 is a pricipal ideal domai ad that the correspodig j-ivariat is a algebraic iteger. This is the fastest coverget series to /π of the form C π = α A + Bβ, with α, A, B, β Q. The Borweis were the first to provide proofs to most of Ramauja s series to /π. I both Chudovsys ad Borweis wor, a crucial igrediet is the Clause s formula F a, b; a + b + ; z = 3 F a, b, a + b, a + b + ; a + b; z, mf m a,, a m ; b,, b m ; z := a a m z b b m!. Recetly, Sog Heg, ZhiGuo ad I [4] discovered that Clause s formula is ot eeded for the derivatio of Ramauja s series. Our wor is ispired

6 HENG HUAT CHAN by a series discovered by Taeshi Sato, amely, 5 5. π 4 5 9 + j = j j 3 = j= 5 5 +. Sice o Clause type formula exists for the above series, ew method eeds to be used for the derivatio of series such as 5.. Our ew method yields the followig compaio of 5.: 5. 8 3π = = j= j j j j j + 5. 64 We ed this sectio with two more ew series arisig from this ew method. The first is aother compaio of 5. that I discovered together with H. Verrill after the wor of Almvist ad Zudili []: 3 3 3 + 3 = 3!! 3 4 + = 9 8 3π. The secod oe is oe I discovered recetly with Loo Ko Pig: 3 + = m 3 m m 3 m 6 π 9 8 j i 3 + + G 3, j= i= e π/3 G = e π + e π + + e π + e 4 π + e 3 π + e 6 π + = + 6. 6. The Wilf-Zeilberger algorithm ad Guillera s series for /π A discrete fuctio A, is hypergeometric if A +, A, ad A, + A, are both ratioal fuctios. A pair of fuctios F, ad G, is said to be of WZ after Wilf ad Zeilberger if F ad G are hypergeometric ad F +, F, = G, + G,.

SOME NEW IDENTITIES INVOLVING π, π AND π 7 I this case, Wilf ad Zeilberger showed that there exists a ratioal fuctio C, such that G, = C, F,. The fuctio C, is called the certificate of F, G. Defiig Zeilberger showed that H, = F +, + + G, +, H, = G,. Ehad ad Zeilberger [6] was the first to use the above method to derive a oe page proof of Ramauja s series see before i the previous sectio π = 4 + 3! 3. Motivated by his wor, J. Guillera [7] foud may ew WZ-pairs F, G ad derived ew series ot oly for /π but for /π. Oe of the most elegat formulas of Guillera is 5 8 + 8 + 3 4 = 3 π. This is proved usig Zeilberger s method by exhibitig the WZ-pair, amely, ad F, = 5 3 U, 4 G, = + 84 + 34 + + 3U, U, = cosπ Γ + Γ6 + Γ3 + π 5 Γ + + Γ 3 + Motivated by the elegace of the above series, Guillera discovered empirically usig PSLQ algorithm may series of type C π = C + D + E α H. Most of the series he discovered caot be proved by the WZ-method, which apparetly ca oly be applied whe H is a power of. A example of Guillera s series which remais to be proved is 8 5 π = 5 3 3 6 6! 5 8 3 548 + 693 + 9. Ehad is a computer

8 HENG HUAT CHAN Refereces [] G. Almvist, C. Krattethaler, ad J. Petersso, Some ew formulas for π, Experimet. Math. 3, o. 4, 44 456 [] G. Almvist, W. Zudili, Differetial equatios, mirror maps ad zeta values, preprit. [3] D. Bailey, P. Borwei ad S. Plouffe, O the rapid computatio of various polylogarithmic costats, Math. Comp. 66 997, o. 8, 93 93 [4] H.H. Cha, S.H. Cha ad Z.-G. Liu, Domb s umbers ad Ramauja-Sato type series for /π, Adv. Math. 86 4, o., 396 4. [5] H.H. Cha, W.C. Liaw ad V. Ta Ramauja s class ivariat λ ad a ew class of series for /π. J. Lodo Math. Soc., 64, o., 93-6. [6] S. Ehad, D. Zeilberger, A WZ proof of Ramauja s formula for π, Geometry, aalysis ad mechaics, 7 8, World Sci. Publishig, River Edge, NJ, 994 [7] J. Guillera, Some biomial series obtaied by the WZ-method, Adv. i Appl. Math. 9, o. 4, 599 63. 33F [8] T. Osler, The uio of Vieta s ad Wallis s products for π, Amer. Math. Mothly, 6 999, 774-776. [9] S. Ramauja, Modular equatios ad approximatios to π, Quart. J. Math. Oxford, 45 94, 35-37.