Construction Of Binomial Sums For π And Polylogarithmic Constants Inspired By BBP Formulas
|
|
- Miles Jason Daniels
- 5 years ago
- Views:
Transcription
1 Applied Mathematics E-Notes, 77, 7-46 c ISSN 67-5 Available free at mirror sites of amen/ Construction Of Binomial Sums For π And Polylogarithmic Constants Inspired By BBP Formulas Boris Gourévitch, Jesús Guillera Goyanes Received 8 October 6 Abstract We present new sums involving binomial coefficients for π and various logarithms and polylogarithms constants. These sums are a generalization of BBP formulas first introduced by D. Bailey, P. Borwein and S. Plouffe in 995. In this paper, we describe how to find and prove such sums using the Beta function at integer and rational arguments. Introduction In 995, David Bailey, Peter Borwein and Simon Plouffe discovered a new seminal formula for π 4 π 6 n 8n + 8n +4 8n +5. 8n +6 n This amazing formula allows extracting the n-th binary digit of π without knowledge of previous digits []. Such sums are called BBP Bailey-Borwein-Plouffe formulas. Many new BBP formulas for π, G Catalan s constant or ζ Riemann Zeta function were subsequently discovered by Adamchik [], Bellard [], Broadhurst [4], Lupas [5] and Huvent [6]. The main tool for finding BBP formulas is an algorithm designed to detect linear combinations on the set of natural numbers N like PSLQ from its use of a partial-sum-of-squares PS vector and lower-diagonal-orthogonal LQ matrix factorization or LLL Lenstra Lenstra Lovász [7], [8]. To then prove such formulas, an integral is evaluated giving rise to logarithms or polylogarithms Li p z n zn, n p [9]. In October, Almkvist and colleagues presented new sums for π involving binomial coefficients such as π n Mathematics Subject Classifications: 4C 7 rue Montera, 75 Paris, boris@pi4.net Av. Cesareo Alierta, esc. izda 4oVA, Zaragoza 5n 6 n n n. 7
2 8 Binomial Sums for π and Polylogarithmic Constants They proved this formula using the so-called Beta method which is based on use of Beta function Br, s xr x s for positive integer values r and s []. In this paper, using an adapted Beta method, we show how to prove some more general sums of the form T q n mn n b mn + + b mn b k mn + k + + b m mn +m where T is a constant often related to π and the natural logarithm lnz. These sums involve binomial coefficients and they are similar to BBP formulas. We therefore call them BBP binomial formulas. One example of such formulas is π 687 n 7n n n n + 6 7n + 7n + 5 7n n n +6 4 that we discovered with LLL algorithm. An explicit proof of this formula is detailed in this paper. Then we show how to generate an infinite set of similar formulas for π and the natural logarithm lnz. We also present outlines for proofs of formulas showing that Beta method is not restrained to integer parameters. Proofs of BBP Binomial Formulas Consider the Eq.. Proof of such formula is reduced to the following proposition: PROPOSITION The Beta Method. There exists a polynomial P X Q m [X] polynomial of degree m with coefficients in the set of rational numbers Q such that qpx T 5 q x p x m p PROOF. The Beta function Br, s is Br, s ΓrΓs Γr + s 6 where Γz tz e t dt is the Gamma function. With parameters m, p, n N, we use Γn n! to obtain mn + mn x x m. 7
3 B. Gourévitch and J. Guillera Goyanes 9 Summing Eq. 7 on n with mn+ mn on the left side, we obtain n q n mn mn + q n n x x m n x p x m p n q q. 8 q x p x m p We can evaluate this integral by partial fraction decomposition. More generally, to obtain the elements in Eq. we write mn + k x +k x m mn+k +k mn k mn + mn +...mn + k. 9 For a given k, the partial fraction decomposition of k mn+mn+...mn+k in n gives a k,j j,..,k such that x k x x m mn ak, mn + + a k, mn a k,k mn + k with a k,j Q. For values k ranging from to m, we linearly combine a k,j j,..,k in order to find the m values b k of Eq.. This is equivalent to solving the system a, a,... a m, b a,... a m, V b a m,k b m Vector V contains the coefficients of the polynomial P X V t. m, such that q n mn n n when we use Eq.. b mn + + b mn P x x x m q n qpx q x p x m p X... X m b m mn +m of degree
4 4 Binomial Sums for π and Polylogarithmic Constants REMARK. Almkvist et al. [] use the relation n n n + xn x n in order to prove the formula in Eq.. In general, they consider formulas of the form Sn π q n n where S Q m [X] and q N. However, such forms of sums lead to the intricate evaluation of fx n mn +Sn x p x n. m p q fx is indeed an integral whose denominator is a power of the one considered in Eq. 8. Our method uses on the left side of Eq. 7, which produces formulas of the general form of mn+ mn Eq.. Proof of Eq. 4 We write Eq. 4 as mn U n π Using Proposition, we are looking for the polynomial P x such that x n x 5n U n P x P x n n n x x 5 P x x x +x 5 x 4 +x + x +. Because n 4 x x + π, it is likely that P x x 5 x 4 +x + x +. Let us verify this conjecture. To obtain the equivalent linear system in, and therefore vector V, we calculate all coefficients a k,i i columns of left matrix in Eq. with help from the computation of the integral in Eq.. For example, for k nd column of matrix, we have x.x n x 5n B5n +, n + n + 7n 7n + 7n + n 5 7n 7 7n n + n Consequently, a, 5 7 and a, 7. The whole linear system is V
5 B. Gourévitch and J. Guillera Goyanes 4 Its solution is V. Finally, we have P X V t X... X 5 which is the expected polynomial. +X +X X 4 + X 5 4 Construction and Prediction of BBP Binomial Formulas qpx Using proposition, we choose a suitable P x that simplifies. Knowing P then gives rise to a BBP formula by using Eq., after combination and q x m p x p simplification. As an example, we find formulas for π choosing P X such that due to +x qpx q x m p x p +x π 4. To find m and p, we have to know when + X divides q X m p X p i.e. when i and i are roots of q X m p X p. Case m p in Eq. is excluded since it only provides classical BBP formulas without binomial coefficients. Let us choose p,m, then i i i + i. Now there exists a formula for π with n n, p and q. Moreover, i and i are fourth roots of, so we can multiply x x by x 4k and keep roots i and i. Thus there exists formulas for π involving C+4kn n, k N. More generally, if KX q X m p X p, { } {i and i are roots of K, p < m} q ±i m p i p,p<m { ± i p m+p q±i m,p<m} { + i p m+p qi m,p<m} by complex conjugate { p } q and p π 4 m + pπ + m π [π],p<m with modulus and argument, where [x] designates modulo x { p q, p < m and p m [8] } { p k, k N, k q, p < m and 7p m [8] }. This idea also leads to other constants using the root ln the polynomial x x +π 9 x+x. For non-alternating sums b n mn+ + b mn bm mn+m with equivalent integral representation q n mn, or +x qpx q x p x m p, we obtain the following table for
6 4 Binomial Sums for π and Polylogarithmic Constants p<m: Polynomial Constant Condition given by roots Existence of a sum for p<m +x π i m p i p q p k, k q, 7p m [8] +x ln m p p q m p +k, q p x + x π p p+m +i q q,mp +6k +x ln z +x ln z+ z +z x ln z For alternating sums n bm with equiva- mn+m lent integral representation m p p q m p +k, q k p z m p z p q m p +k, q z k z p m p p z + z q m p +k, q z p m q n mn b mn+ + b mn z p+k qpx q+x p x m p, we get the following table: Polynomial Constant Condition given by roots Existence of a sum for p<m +x π i m p i p q p k, k q, 7p m + 4 [8] +x ln m p+ p q m p +k +,q p x + x π p+ +i p+m q q,mp ++6k +x ln z +x ln z+ z +z x ln z m p p q m p +k +,q k p z m p z p q m p +k +,q z k z p m p p z + z q m p +k +,q z p z p+k 5 Generalized Beta Function Method Albeit more intricate, some similar formulas can also be found using fractional parameters in Beta function. Considering Beta function and a b, c Q, we have d x + a b x mn+ Γ + a c d b + Γ mn + c d + Γ p + mn + a b + c d + a c b + d mn+ Γ a b Γ c d a b + d c Γ a p+mn+ b + c d Γx + n where x n xx+ x+n is Pochhammer symbol. Γ a Γx b is often known for small integers values of a et b. All following formulas stem from summing this relation on n for particular values of a b, c d. We now emphasize three particular cases of interest.
7 B. Gourévitch and J. Guillera Goyanes 4 5. Case : a b, c d We have q x q x p x m n x n x x mn q n q n x x mn q n Γ + Γmn + Γ p + mn + n + n!mn +!mn!. 4 q n mn +!! n EXAMPLE. Case p,m. We have π n xx + x x +x x n n x x +x x x x n n x n x n n x x +x x n x x n 5 So now using Eq. 4, we have x x +x x n x n x n [,,, ] x n+[ 5,,, ] T x n Γ n + [ 7, 5,, ] T Γn + [,,, ] Γ 4n + [ 9, 7, 5, ] T 6n!4n!n!4n 885 n!8n! 5 8n n n +5 + in simplified notations. Finally with Eq. 5 and Eq. 6, we have 6n 9 π n n n 885 8n n n n 4n 4n n 5 8n n Case : a + c Z positive and negative integers b d This case allows us to obtain special values of Γ a b Γ c d equal to π multiplied by an algebraic number. In some cases, π is absent and we obtain amazing formulas for algebraic numbers.. 6
8 44 Binomial Sums for π and Polylogarithmic Constants EXAMPLE. The following relations are based on Beta values for a b 4, c d or a b, c d or a b 5, c d n n 7 n n 4π n n 7 n n 8π which gives the formula 8 n x n 4 x n π x n 5 x n 4π n n n 54 n n +. n x n x n x n x n 5. Case : a + c Q\N rationals but not naturals b d We obtain Pochhammer symbols in sums then factors Γ a b Γ c d Γ a b + c in the outcome. d EXAMPLE. Considering x n 4 x n Γ n + 4 Γ n + Γ n x n+ 4 x n Γ n Γ n + Γ n and the following integral n x n xn 4 x n 4 π x xx + we obtain finally n n 4 n 8 n n n n 4 n 8n + + Γ 8n K π where Kk π/ dθ is the complete elliptic integral of the first kind, thus k sin θ making the connection with singular values of the elliptic integral. 6 Higher Order Constants One point of interest is to generalize this method to higher order constants. A constant is of order p if it is equal to a linear combination in N of p-th order polylogarithms. For example, π or ln are of order whereas π ln or G are of order and ζ or π
9 B. Gourévitch and J. Guillera Goyanes 45 are of order. A BBP binomial sum involves a constant of order p, as for BBP sums, if the polynomial in n in the denominator is of degree p. For example, we have n n n + ln + ζ. 7 n + n n n However, the straightforward application of integral representation does not completely help since, for order for instance, we have q n mn mn + n q dy. q y m x m x m p The evaluation is not often easy, there is no known proof of Eq. 7 to our knowledge. Nevertheless, there is evidence for existence of BBP binomial formulas for higher order constants. Acknowledgements. The authors would like to thank Typhenn Brichieri-Colombi, Professor Sui Sun Cheng and the anonymous reviewer for their helpful comments and corrections. References [] D. Bailey, P. Borwein and S. Plouffe, On The Rapid Computation of Various Polylogarithmic Constants, Math. Comp., 66997, 9 9. [] V. Adamchik, π : A -Year Search Changes Direction, Education and Research, 5996, 9 [] F. Bellard, Computation of the n th digit of π in any base in O n, preprint 997, n.ps [4] D. J. Broadhurst, Polylogarithmic ladders, hypergeometric series and the ten millionth digits of ζ and ζ 5, preprint 998, [5] A. Lupas, Some BBP functions, preprint, Sibiu University, Romania. [6] G. Huvent, Formules BBP, Seminar ANO, IRMA, Lille University, huvent.pdf [7] A. K. Lenstra, H. W. Jr. Lenstra, L. Lovász, Factoring polynomials with rational coefficients, Mathematische Annalen 6498, [8] Pari-GP, software pour PC, MacOS et Unix, documentation at [9] L. Lewin, Polylogarithms and Associated Functions, North-Holland, New York, 98.
10 46 Binomial Sums for π and Polylogarithmic Constants [] G. Almkvist, C. Krattenthaler and J. Petersson, Some new formulas for π, preprint, Matematiska Institutionen, Lunds Universitet, Sweden. and Institut fur Mathematik der Universitat Wien, Austria,
Formulae for some classical constants
Formulae for some classical constants Alexandru Lupaş (to appear in Proceedings of ROGER-000 The goal of this paper is to present formulas for Apéry Constant, Archimede s Constant, Logarithm Constant,
More informationA gentle introduction to PSLQ
A gentle introduction to PSLQ Armin Straub Email: math@arminstraub.com April 0, 200 Abstract This is work in progress. Please let me know about any comments and suggestions. What PSLQ is about PSLQ is
More informationInteger Relation Detection and Lattice Reduction David H. Bailey 1. To appear in Computing in Science and Engineering
Integer Relation Detection and Lattice Reduction David H. Bailey To appear in Computing in Science and Engineering. Introduction Let x =(x,x,,x n ) be a vector of real or complex numbers. x is said to
More informationThe Miraculous Bailey-Borwein-Plouffe Pi Algorithm
Overview: 10/1/95 The Miraculous Bailey-Borwein-Plouffe Pi Algorithm Steven Finch, Research and Development Team, MathSoft, Inc. David Bailey, Peter Borwein and Simon Plouffe have recently computed the
More informationBinary BBP-Formulae for Logarithms and Generalized Gaussian-Mersenne Primes
3 47 6 3 Journal of Integer Sequences, Vol. 6 003), Article 03.3.7 Binary BBP-Formulae for Logarithms and Generalized Gaussian-Mersenne Primes Marc Chamberland Department of Mathematics and Computer Science
More informationInteger Relation Methods : An Introduction
Integer Relation Methods : An Introduction Special Session on SCIENTIFIC COMPUTING: July 9 th 2009 Jonathan Borwein, FRSC www.carma.newcastle.edu.au/~jb616 Laureate Professor University of Newcastle, NSW
More informationExperimental mathematics and integration
Experimental mathematics and integration David H. Bailey http://www.davidhbailey.com Lawrence Berkeley National Laboratory (retired) Computer Science Department, University of California, Davis October
More informationFormulas for Odd Zeta Values and Powers of π
3 47 6 3 Journal of Integer Sequences, Vol. 4 (0), Article..5 Formulas for Odd Zeta Values and Powers of π Marc Chamberland and Patrick Lopatto Department of Mathematics and Statistics Grinnell College
More informationBackground and Definitions...2. Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and Functions...
Legendre Polynomials and Functions Reading Problems Outline Background and Definitions...2 Definitions...3 Theory...4 Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and
More informationπ in terms of φ via the Machin s Route
in terms of φ via the Machin s Route Hei-Chi Chan Mathematical Science Program, University of Illinois at Springfield Springfield, IL 62703-507 email: chan.hei-chi@uis.edu Abstract In this paper, we prove
More informationParametric Euler Sum Identities
Parametric Euler Sum Identities David Borwein, Jonathan M. Borwein, and David M. Bradley September 23, 2004 Introduction A somewhat unlikely-looking identity is n n nn x m m x n n 2 n x, valid for all
More informationExtracting Hexadecimal Digits of π
Extracting Hexadecimal Digits of π William Malone William Malone is currently a senior Mathematics major at Ball State University. Upon graduation he plans to pursue graduate studies starting in the fall
More informationBINARY BBP-FORMULAE FOR LOGARITHMS AND GENERALIZED GAUSSIAN-MERSENNE PRIMES. Marc Chamberland
1 BINARY BBP-FORMULAE FOR LOGARITHMS AND GENERALIZED GAUSSIAN-MERSENNE PRIMES Marc Chamberland Department of Mathematics and Computer Science Grinnell College Grinnell, IA, 5011 E-mail: chamberl@math.grinnell.edu
More informationExperimental Determination of Apéry-like Identities for ζ(2n +2)
Experimental Determination of Apéry-lie Identities for ζ(n + David H. Bailey, Jonathan M. Borwein, and David M. Bradley CONTENTS. Introduction. Discovering Theorem. 3. Proof of Theorem.. An Identity for
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More informationClosed-form Second Solution to the Confluent Hypergeometric Difference Equation in the Degenerate Case
International Journal of Difference Equations ISS 973-669, Volume 11, umber 2, pp. 23 214 (216) http://campus.mst.edu/ijde Closed-form Second Solution to the Confluent Hypergeometric Difference Equation
More informationQuadratic Transformations of Hypergeometric Function and Series with Harmonic Numbers
Quadratic Transformations of Hypergeometric Function and Series with Harmonic Numbers Martin Nicholson In this brief note, we show how to apply Kummer s and other quadratic transformation formulas for
More informationAlgorithms for Experimental Mathematics I David H Bailey Lawrence Berkeley National Lab
Algorithms for Experimental Mathematics I David H Bailey Lawrence Berkeley National Lab All truths are easy to understand once they are discovered; the point is to discover them. Galileo Galilei Algorithms
More informationThe Value of the Zeta Function at an Odd Argument
International Journal of Mathematics and Computer Science, 4(009), no., 0 The Value of the Zeta Function at an Odd Argument M CS Badih Ghusayni Department of Mathematics Faculty of Science- Lebanese University
More informationSimplifying Rational Expressions and Functions
Department of Mathematics Grossmont College October 15, 2012 Recall: The Number Types Definition The set of whole numbers, ={0, 1, 2, 3, 4,...} is the set of natural numbers unioned with zero, written
More informationComputer discovery of new mathematical facts and formulas
Computer discovery of new mathematical facts and formulas David H. Bailey http://www.davidhbailey.com Lawrence Berkeley National Laboratory (retired) University of California, Davis, Department of Computer
More informationSums and Products. a i = a 1. i=1. a i = a i a n. n 1
Sums and Products -27-209 In this section, I ll review the notation for sums and products Addition and multiplication are binary operations: They operate on two numbers at a time If you want to add or
More informationSolutions to Homework 2
Solutions to Homewor Due Tuesday, July 6,. Chapter. Problem solution. If the series for ln+z and ln z both converge, +z then we can find the series for ln z by term-by-term subtraction of the two series:
More informationMathematics 324 Riemann Zeta Function August 5, 2005
Mathematics 324 Riemann Zeta Function August 5, 25 In this note we give an introduction to the Riemann zeta function, which connects the ideas of real analysis with the arithmetic of the integers. Define
More informationFinding and Excluding b-ary Machin-Type Individual Digit Formulae
DRAFT: Canad. J. Math. July 6, 2004 14:27 File: borwein3298 pp.1 29 Page 1 Sheet 1 of 29 Canad. J. Math. Vol. XX (Y), ZZZZ pp. 1 29 Finding and Excluding b-ary Machin-Type Individual Digit Formulae Jonathan
More informationThe Generating Functions for Pochhammer
The Generating Functions for Pochhammer Symbol { }, n N Aleksandar Petoević University of Novi Sad Teacher Training Faculty, Department of Mathematics Podgorička 4, 25000 Sombor SERBIA and MONTENEGRO Email
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationAnalytical expressions of 3 and 4-loop sunrise Feynman integrals and 4-dimensional lattice integrals
3 September 8 Analytical expressions of 3 and 4-loop sunrise Feynman integrals and 4-dimensional lattice integrals arxiv:83.7v3 [hep-ph] Mar 8 S. Laporta Museo Storico della Fisica e Centro Studi e Ricerche
More informationHYPERGEOMETRIC BERNOULLI POLYNOMIALS AND APPELL SEQUENCES
HYPERGEOMETRIC BERNOULLI POLYNOMIALS AND APPELL SEQUENCES ABDUL HASSEN AND HIEU D. NGUYEN Abstract. There are two analytic approaches to Bernoulli polynomials B n(x): either by way of the generating function
More informationMULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng
Indian J. pure appl. Math., 39(: 37-55, April 008 c Printed in India. MULTISECTION METHOD AND FURTHER FORMULAE FOR π De-Yin Zheng Department of Mathematics, Hangzhou Normal University, Hangzhou 30036,
More informationrama.tex; 21/03/2011; 0:37; p.1
rama.tex; /03/0; 0:37; p. Multiple Gamma Function and Its Application to Computation of Series and Products V. S. Adamchik Department of Computer Science, Carnegie Mellon University, Pittsburgh, USA Abstract.
More informationThe Riemann and Hurwitz zeta functions, Apery s constant and new rational series representations involving ζ(2k)
The Riemann and Hurwitz zeta functions, Apery s constant and new rational series representations involving ζ(k) Cezar Lupu 1 1 Department of Mathematics University of Pittsburgh Pittsburgh, PA, USA Algebra,
More informationSummation Techniques, Padé Approximants, and Continued Fractions
Chapter 5 Summation Techniques, Padé Approximants, and Continued Fractions 5. Accelerated Convergence Conditionally convergent series, such as 2 + 3 4 + 5 6... = ( ) n+ = ln2, (5.) n converge very slowly.
More informationON THE RAPID COMPUTATION OF VARIOUS POLYLOGARITHMIC CONSTANTS
MATHEMATICS OF COMPUTATION Volume 66 Number 218 April 1997 Pages 903 913 S 0025-5718(97)00856-9 ON THE RAPID COMPUTATION OF VARIOUS POLYLOGARITHMIC CONSTANTS DAVID BAILEY PETER BORWEIN AND SIMON PLOUFFE
More information1. Introduction Interest in this project began with curiosity about the Laplace transform of the Digamma function, e as ψ(s + 1)ds,
ON THE LAPLACE TRANSFORM OF THE PSI FUNCTION M. LAWRENCE GLASSER AND DANTE MANNA Abstract. Guided by numerical experimentation, we have been able to prove that Z 8 / x x + ln dx = γ + ln) [cosx)] and to
More information27 Wyner Math 2 Spring 2019
27 Wyner Math 2 Spring 2019 CHAPTER SIX: POLYNOMIALS Review January 25 Test February 8 Thorough understanding and fluency of the concepts and methods in this chapter is a cornerstone to success in the
More informationAnalytic Aspects of the Riemann Zeta and Multiple Zeta Values
Analytic Aspects of the Riemann Zeta and Multiple Zeta Values Cezar Lupu Department of Mathematics University of Pittsburgh Pittsburgh, PA, USA PhD Thesis Overview, October 26th, 207, Pittsburgh, PA Outline
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationRamanujan-Type Series Related to Clausen s Product
Ramanujan-Type Series Related to Clausen s Product arxiv:.9v [math.nt] 7 Jan John M. Campbell York University Toronto, ON Canada maxwell8@yorku.ca Abstract Infinite series are evaluated through the manipulation
More informationLaura Chihara* and Dennis Stanton**
ZEROS OF GENERALIZED KRAWTCHOUK POLYNOMIALS Laura Chihara* and Dennis Stanton** Abstract. The zeros of generalized Krawtchouk polynomials are studied. Some interlacing theorems for the zeros are given.
More informationIntegral approximations to π with nonnegative integrands
Integral approximations to π with nonnegative integrands S.K. Lucas Department of Mathematics and Statistics James Madison University Harrisonburg VA 2287 Email: lucassk@jmu.edu May 27 One of the more
More informationMathematical Induction Assignments
1 Mathematical Induction Assignments Prove the Following using Principle of Mathematical induction 1) Prove that for any positive integer number n, n 3 + 2 n is divisible by 3 2) Prove that 1 3 + 2 3 +
More informationFor many years, researchers have dreamt INTEGER RELATION DETECTION T HEME ARTICLE. the Top
the Top T HEME ARTICLE INTEGER RELATION DETECTION Practical algorithms for integer relation detection have become a staple in the emerging discipline of experimental mathematics using modern computer technology
More informationFRACTIONAL HYPERGEOMETRIC ZETA FUNCTIONS
FRACTIONAL HYPERGEOMETRIC ZETA FUNCTIONS HUNDUMA LEGESSE GELETA, ABDULKADIR HASSEN Both authors would like to dedicate this in fond memory of Marvin Knopp. Knop was the most humble and exemplary teacher
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationINTEGRATION WORKSHOP 2003 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 23 COMPLEX ANALYSIS EXERCISES DOUGLAS ULMER 1. Meromorphic functions on the Riemann sphere It s often useful to allow functions to take the value. This exercise outlines one way to
More informationON THE EQUATION X n 1 = BZ n. 1. Introduction
ON THE EQUATION X n 1 = BZ n PREDA MIHĂILESCU Abstract. We consider the Diophantine equation X n 1 = BZ n ; here B Z is understood as a parameter. We give new conditions for existence of solutions to this
More informationFastTrack - MA109. Exponents and Review of Polynomials
FastTrack - MA109 Exponents and Review of Polynomials Katherine Paullin, Ph.D. Lecturer, Department of Mathematics University of Kentucky katherine.paullin@uky.edu Monday, August 15, 2016 1 / 25 REEF Question
More informationf(g(x)) g (x) dx = f(u) du.
1. Techniques of Integration Section 8-IT 1.1. Basic integration formulas. Integration is more difficult than derivation. The derivative of every rational function or trigonometric function is another
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (I)
Laplace s Equation in Cylindrical Coordinates and Bessel s Equation I) 1 Solution by separation of variables Laplace s equation is a key equation in Mathematical Physics. Several phenomena involving scalar
More informationON THE RAPID COMPUTATION OF VARIOUS POLYLOGARITHMIC CONSTANTS. David Bailey, Peter Borwein 1 and Simon Plouffe
ON THE RAPID COMPUTATION OF VARIOUS POLYLOGARITHMIC CONSTANTS David Bailey, Peter Borwein 1 and Simon Plouffe Abstract. We give algorithms for the computation of the d-th digit of certain transcendental
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationORTHOGONALITY AND IRRATIONALITY. This version was printed on November 3, 2013
ORTHOGONALITY AND IRRATIONALITY PETER BORWEIN, TAMÁS ERDÉLYI, AND JOHN ZHANG This version was printed on November 3, 23 Abstract. The fairly recent proofs of the irrationality of log2,π 2,ζ(3) and various
More informationCongruences for algebraic sequences
Congruences for algebraic sequences Eric Rowland 1 Reem Yassawi 2 1 Université du Québec à Montréal 2 Trent University 2013 September 27 Eric Rowland (UQAM) Congruences for algebraic sequences 2013 September
More informationRelations and Functions (for Math 026 review)
Section 3.1 Relations and Functions (for Math 026 review) Objective 1: Understanding the s of Relations and Functions Relation A relation is a correspondence between two sets A and B such that each element
More informationA Note about the Pochhammer Symbol
Mathematica Moravica Vol. 12-1 (2008), 37 42 A Note about the Pochhammer Symbol Aleksandar Petoević Abstract. In this paper we give elementary proofs of the generating functions for the Pochhammer symbol
More informationMultiple Zeta Values of Even Arguments
Michael E. Hoffman U. S. Naval Academy Seminar arithmetische Geometrie und Zahlentheorie Universität Hamburg 13 June 2012 1 2 3 4 5 6 Values The multiple zeta values (MZVs) are defined by ζ(i 1,..., i
More informationExpectations on Fractal Sets
Expectations on Fractal Sets David H. Bailey http://www.davidhbailey.com Lawrence Berkeley Natl. Lab. (retired) Computer Science Dept., University of California, Davis Co-authors: Jonathan M. Borwein (CARMA,
More informationMahler measure and special values of L-functions
Mahler measure and special values of L-functions Matilde N. Laĺın University of Alberta mlalin@math.ualberta.ca http://www.math.ualberta.ca/~mlalin October 24, 2008 Matilde N. Laĺın (U of A) Mahler measure
More informationNumeration and Computer Arithmetic Some Examples
Numeration and Computer Arithmetic 1/31 Numeration and Computer Arithmetic Some Examples JC Bajard LIRMM, CNRS UM2 161 rue Ada, 34392 Montpellier cedex 5, France April 27 Numeration and Computer Arithmetic
More informationComplex Numbers: Definition: A complex number is a number of the form: z = a + bi where a, b are real numbers and i is a symbol with the property: i
Complex Numbers: Definition: A complex number is a number of the form: z = a + bi where a, b are real numbers and i is a symbol with the property: i 2 = 1 Sometimes we like to think of i = 1 We can treat
More informationDifference Equations for Multiple Charlier and Meixner Polynomials 1
Difference Equations for Multiple Charlier and Meixner Polynomials 1 WALTER VAN ASSCHE Department of Mathematics Katholieke Universiteit Leuven B-3001 Leuven, Belgium E-mail: walter@wis.kuleuven.ac.be
More informationLecture 3. Frits Beukers. Arithmetic of values of E- and G-function. Lecture 3 E- and G-functions 1 / 20
Lecture 3 Frits Beukers Arithmetic of values of E- and G-function Lecture 3 E- and G-functions 1 / 20 G-functions, definition Definition An analytic function f (z) given by a powerseries a k z k k=0 with
More information17 Advancement Operator Equations
November 14, 2017 17 Advancement Operator Equations William T. Trotter trotter@math.gatech.edu Review of Recurrence Equations (1) Problem Let r(n) denote the number of regions determined by n lines that
More informationA P-ADIC SUPERCONGRUENCE CONJECTURE OF VAN HAMME
A P-ADIC SUPERCONGRUENCE CONJECTURE OF VAN HAMME ERIC MORTENSON Abstract. In this paper we prove a p-adic supercongruence conjecture of van Hamme by placing it in the context of the Beuers-lie supercongruences
More informationSome algebraic number theory and the reciprocity map
Some algebraic number theory and the reciprocity map Ervin Thiagalingam September 28, 2015 Motivation In Weinstein s paper, the main problem is to find a rule (reciprocity law) for when an irreducible
More informationOn A Central Binomial Series Related to ζ(4).
On A Central Binomial Series Related to ζ. Vivek Kaushik November 8 arxiv:8.67v [math.ca] 5 Nov 8 Abstract n n n 7π 3 by We prove a classical binomial coefficient series identity n evaluating a double
More informationA PROBLEM ON THE CONJECTURE CONCERNING THE DISTRIBUTION OF GENERALIZED FERMAT PRIME NUMBERS (A NEW METHOD FOR THE SEARCH FOR LARGE PRIMES)
A PROBLEM ON THE CONJECTURE CONCERNING THE DISTRIBUTION OF GENERALIZED FERMAT PRIME NUMBERS A NEW METHOD FOR THE SEARCH FOR LARGE PRIMES) YVES GALLOT Abstract Is it possible to improve the convergence
More informationMath221: HW# 7 solutions
Math22: HW# 7 solutions Andy Royston November 7, 25.3.3 let x = e u. Then ln x = u, x2 = e 2u, and dx = e 2u du. Furthermore, when x =, u, and when x =, u =. Hence x 2 ln x) 3 dx = e 2u u 3 e u du) = e
More informationPrelim Examination 2010 / 2011 (Assessing Units 1 & 2) MATHEMATICS. Advanced Higher Grade. Time allowed - 2 hours
Prelim Examination 00 / 0 (Assessing Units & ) MATHEMATICS Advanced Higher Grade Time allowed - hours Read Carefully. Calculators may be used in this paper.. Candidates should answer all questions. Full
More informationMathematics 136 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 19 and 21, 2016
Mathematics 36 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 9 and 2, 206 Every rational function (quotient of polynomials) can be written as a polynomial
More informationPrecise error estimate of the Brent-McMillan algorithm for the computation of Euler s constant
Precise error estimate of the Brent-McMillan algorithm for the computation of Euler s constant Jean-Pierre Demailly Université de Grenoble Alpes, Institut Fourier Abstract Brent and McMillan introduced
More informationOn the power-free parts of consecutive integers
ACTA ARITHMETICA XC4 (1999) On the power-free parts of consecutive integers by B M M de Weger (Krimpen aan den IJssel) and C E van de Woestijne (Leiden) 1 Introduction and main results Considering the
More informationAS PURE MATHS REVISION NOTES
AS PURE MATHS REVISION NOTES 1 SURDS A root such as 3 that cannot be written exactly as a fraction is IRRATIONAL An expression that involves irrational roots is in SURD FORM e.g. 2 3 3 + 2 and 3-2 are
More informationarxiv: v3 [math.nt] 9 May 2011
ON HARMONIC SUMS AND ALTERNATING EULER SUMS arxiv:02.592v3 [math.nt] 9 May 20 ZHONG-HUA LI Department of Mathematics, Tongji University, No. 239 Siping Road, Shanghai 200092, China Graduate School of Mathematical
More informationA New Shuffle Convolution for Multiple Zeta Values
January 19, 2004 A New Shuffle Convolution for Multiple Zeta Values Ae Ja Yee 1 yee@math.psu.edu The Pennsylvania State University, Department of Mathematics, University Park, PA 16802 1 Introduction As
More informationAbstract. 2. We construct several transcendental numbers.
Abstract. We prove Liouville s Theorem for the order of approximation by rationals of real algebraic numbers. 2. We construct several transcendental numbers. 3. We define Poissonian Behaviour, and study
More informationSign changes of Fourier coefficients of cusp forms supported on prime power indices
Sign changes of Fourier coefficients of cusp forms supported on prime power indices Winfried Kohnen Mathematisches Institut Universität Heidelberg D-69120 Heidelberg, Germany E-mail: winfried@mathi.uni-heidelberg.de
More informationITALIAN JOURNAL OF PURE AND APPLIED MATHEMATICS N ( ) 451
ITALIAN JOURNAL OF PURE AND APPLIED MATHEMATICS N. 39 2018 451 464 451 ON THE k-normal ELEMENTS AND POLYNOMIALS OVER FINITE FIELDS Mahmood Alizadeh Department of Mathematics Ahvaz Branch Islamic Azad University
More informationAMICABLE PAIRS AND ALIQUOT CYCLES FOR ELLIPTIC CURVES OVER NUMBER FIELDS
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 46, Number 6, 2016 AMICABLE PAIRS AND ALIQUOT CYCLES FOR ELLIPTIC CURVES OVER NUMBER FIELDS JIM BROWN, DAVID HERAS, KEVIN JAMES, RODNEY KEATON AND ANDREW QIAN
More informationq-binomial coefficient congruences
CARMA Analysis and Number Theory Seminar University of Newcastle August 23, 2011 Tulane University, New Orleans 1 / 35 Our first -analogs The natural number n has the -analog: [n] = n 1 1 = 1 + +... n
More informationSeries Solutions of Linear Differential Equations
Differential Equations Massoud Malek Series Solutions of Linear Differential Equations In this chapter we shall solve some second-order linear differential equation about an initial point using The Taylor
More informationThe Mertens conjecture revisited
The Mertens conjecture revisited Herman te Riele, CWI, Amsterdam (joint work with Tadej Kotnik, Univ. of Ljubljana, Slovenia) 7th Algorithmic Number Theory Symposium Technische Universität Berlin, 23 28
More informationGaussian vectors and central limit theorem
Gaussian vectors and central limit theorem Samy Tindel Purdue University Probability Theory 2 - MA 539 Samy T. Gaussian vectors & CLT Probability Theory 1 / 86 Outline 1 Real Gaussian random variables
More informationarxiv: v1 [math.nt] 25 Nov 2009
Transformations of Jesus Guillera s formulas for 1/. Gert Almvist arxiv:0911.4849v1 [math.nt] 2 Nov 2009 Introduction. Jesus Guillera has found nine formulas of Ramanuan type for 1/ (see [] [6]). The first
More informationExplicit solution of a class of quartic Thue equations
ACTA ARITHMETICA LXIV.3 (1993) Explicit solution of a class of quartic Thue equations by Nikos Tzanakis (Iraklion) 1. Introduction. In this paper we deal with the efficient solution of a certain interesting
More informationPUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes
PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES Notes. x n+ = ax n has the general solution x n = x a n. 2. x n+ = x n + b has the general solution x n = x + (n )b. 3. x n+ = ax n + b (with a ) can be
More informationAn iteration procedure for a two-term Machin-like formula for pi with small Lehmer s measure
arxiv:706.08835v3 [math.gm] 6 Jul 07 An iteration procedure for a two-term Machin-like formula for pi with small Lehmer s measure S. M. Abrarov and B. M. Quine July 6, 07 Abstract In this paper we present
More informationMULTI-RESTRAINED STIRLING NUMBERS
MULTI-RESTRAINED STIRLING NUMBERS JI YOUNG CHOI DEPARTMENT OF MATHEMATICS SHIPPENSBURG UNIVERSITY SHIPPENSBURG, PA 17257, U.S.A. Abstract. Given positive integers n, k, and m, the (n, k)-th m- restrained
More informationANALOGUES OF THE TRIPLE PRODUCT IDENTITY, LEBESGUE S IDENTITY AND EULER S PENTAGONAL NUMBER THEOREM
q-hypergeometric PROOFS OF POLYNOMIAL ANALOGUES OF THE TRIPLE PRODUCT IDENTITY, LEBESGUE S IDENTITY AND EULER S PENTAGONAL NUMBER THEOREM S OLE WARNAAR Abstract We present alternative, q-hypergeometric
More informationarxiv: v1 [math.co] 18 Oct 2016
ON SOME CLASS OF SUMS ANDREI K. SVININ Abstract. We consider some class of the sums which naturally include the sums of powers of integers. We suggest a number of conjectures concerning a representation
More informationCongruences for combinatorial sequences
Congruences for combinatorial sequences Eric Rowland Reem Yassawi 2014 February 12 Eric Rowland Congruences for combinatorial sequences 2014 February 12 1 / 36 Outline 1 Algebraic sequences 2 Automatic
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More information(1) = 0 = 0,
GALLERY OF WALKS ON THE SQUARE LATTICE BY A TURING PLOTTER FOR BINARY SEQUENCES RICHARD J. MATHAR Abstract. We illustrate infinite walks on the square lattice by interpretation of binary sequences of zeros
More informationZero estimates for polynomials inspired by Hilbert s seventh problem
Radboud University Faculty of Science Institute for Mathematics, Astrophysics and Particle Physics Zero estimates for polynomials inspired by Hilbert s seventh problem Author: Janet Flikkema s4457242 Supervisor:
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationODD MAGNE OGREID AND PER OSLAND
DESY 97-45 hep-th/9868 December 997 SUMMING ONE- AND TWO-DIMENSIONAL SERIES RELATED TO THE EULER SERIES arxiv:hep-th/9868v 6 Jan 998 ODD MAGNE OGREID AND PER OSLAND Abstract We present results for some
More informationDeterministic Polynomial Time Equivalence of Computing the RSA Secret Key and Factoring
Deterministic Polynomial Time Equivalence of Computing the RSA Secret Key and Factoring Jean-Sébastien Coron and Alexander May Gemplus Card International 34 rue Guynemer, 92447 Issy-les-Moulineaux, France
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
More information