NOTES ON BERNOULLI NUMBERS AND EULER S SUMMATION FORMULA. B r = [m = 0] r

Similar documents
Existence Of Solutions For Nonlinear Fractional Differential Equation With Integral Boundary Conditions

ERROR ESTIMATES FOR APPROXIMATING THE FOURIER TRANSFORM OF FUNCTIONS OF BOUNDED VARIATION

Week 8 Lecture 3: Problems 49, 50 Fourier analysis Courseware pp (don t look at French very confusing look in the Courseware instead)

Reinforcement Learning

Special Functions. Leon M. Hall. Professor of Mathematics University of Missouri-Rolla. Copyright c 1995 by Leon M. Hall. All rights reserved.

Sequence and Series of Functions

F.Y. Diploma : Sem. II [CE/CR/CS] Applied Mathematics

MATH 104 FINAL SOLUTIONS. 1. (2 points each) Mark each of the following as True or False. No justification is required. y n = x 1 + x x n n

HOMEWORK 6 - INTEGRATION. READING: Read the following parts from the Calculus Biographies that I have given (online supplement of our textbook):

Moment Generating Function

e t dt e t dt = lim e t dt T (1 e T ) = 1

Extremal graph theory II: K t and K t,t

N! AND THE GAMMA FUNCTION

Suggested Solution for Pure Mathematics 2011 By Y.K. Ng (last update: 8/4/2011) Paper I. (b) (c)

UNIVERSITY OF BRISTOL. Examination for the Degrees of B.Sc. and M.Sci. (Level C/4) ANALYSIS 1B, SOLUTIONS MATH (Paper Code MATH-10006)

Review of the Riemann Integral

( a n ) converges or diverges.

LOCUS 1. Definite Integration CONCEPT NOTES. 01. Basic Properties. 02. More Properties. 03. Integration as Limit of a Sum

ONE RANDOM VARIABLE F ( ) [ ] x P X x x x 3

Math 2414 Homework Set 7 Solutions 10 Points

Transient Solution of the M/M/C 1 Queue with Additional C 2 Servers for Longer Queues and Balking

Integration and Differentiation

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

POWER SERIES R. E. SHOWALTER

L-functions and Class Numbers

Supplement: Gauss-Jordan Reduction

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

An Extension of Hermite Polynomials

Review of Sections

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

MA123, Chapter 9: Computing some integrals (pp )

SLOW INCREASING FUNCTIONS AND THEIR APPLICATIONS TO SOME PROBLEMS IN NUMBER THEORY

1 Notes on Little s Law (l = λw)

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

0 otherwise. sin( nx)sin( kx) 0 otherwise. cos( nx) sin( kx) dx 0 for all integers n, k.

{ } { S n } is monotonically decreasing if Sn

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

MATH 507a ASSIGNMENT 4 SOLUTIONS FALL 2018 Prof. Alexander. g (x) dx = g(b) g(0) = g(b),

Comparison between Fourier and Corrected Fourier Series Methods

SUMMATION OF INFINITE SERIES REVISITED

MAS221 Analysis, Semester 2 Exercises

Approximations of Definite Integrals

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods"

Using Linnik's Identity to Approximate the Prime Counting Function with the Logarithmic Integral

Convergence rates of approximate sums of Riemann integrals

SUTCLIFFE S NOTES: CALCULUS 2 SWOKOWSKI S CHAPTER 11

Infinite Series Sequences: terms nth term Listing Terms of a Sequence 2 n recursively defined n+1 Pattern Recognition for Sequences Ex:

Chapter 7 Infinite Series

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

PROGRESSIONS AND SERIES

Convergence rates of approximate sums of Riemann integrals

f(bx) dx = f dx = dx l dx f(0) log b x a + l log b a 2ɛ log b a.

Department of Mathematical and Statistical Sciences University of Alberta

f(t)dt 2δ f(x) f(t)dt 0 and b f(t)dt = 0 gives F (b) = 0. Since F is increasing, this means that

Week 13 Notes: 1) Riemann Sum. Aim: Compute Area Under a Graph. Suppose we want to find out the area of a graph, like the one on the right:

REAL ANALYSIS I HOMEWORK 3. Chapter 1

Functions, Limit, And Continuity

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

Math 6710, Fall 2016 Final Exam Solutions

Sampling Example. ( ) δ ( f 1) (1/2)cos(12πt), T 0 = 1

The limit comparison test

1.3 Continuous Functions and Riemann Sums

4.8 Improper Integrals

LIMITS OF FUNCTIONS (I)

1. (25 points) Use the limit definition of the definite integral and the sum formulas to compute. [1 x + x2

5.1-The Initial-Value Problems For Ordinary Differential Equations

10.5 Test Info. Test may change slightly.

Math 104: Final exam solutions

Mathematical Statistics. 1 Introduction to the materials to be covered in this course

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

Lecture 15 First Properties of the Brownian Motion

ECE-314 Fall 2012 Review Questions

Fermat Numbers in Multinomial Coefficients

Solutions to selected problems from the midterm exam Math 222 Winter 2015

ANALYSIS HW 3. f(x + y) = f(x) + f(y) for all real x, y. Demonstration: Let f be such a function. Since f is smooth, f exists.

We will begin by supplying the proof to (a).

A Note on Random k-sat for Moderately Growing k

Test Info. Test may change slightly.

ON BILATERAL GENERATING FUNCTIONS INVOLVING MODIFIED JACOBI POLYNOMIALS

One of the common descriptions of curvilinear motion uses path variables, which are measurements made along the tangent t and normal n to the path of

Duration Notes 1. To motivate this measure, observe that the duration may also be expressed as. a a T a

[Nachlass] of the Theory of the Arithmetic-Geometric Mean and the Modulus Function.

K3 p K2 p Kp 0 p 2 p 3 p

lecture 16: Introduction to Least Squares Approximation

The Eigen Function of Linear Systems

STK4080/9080 Survival and event history analysis

Certain sufficient conditions on N, p n, q n k summability of orthogonal series

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

INFINITE SEQUENCES AND SERIES

Sequences and Series of Functions

BINOMIAL THEOREM OBJECTIVE PROBLEMS in the expansion of ( 3 +kx ) are equal. Then k =

n 2 + 3n + 1 4n = n2 + 3n + 1 n n 2 = n + 1

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

Lecture 8 April 18, 2018

Linford 1. Kyle Linford. Math 211. Honors Project. Theorems to Analyze: Theorem 2.4 The Limit of a Function Involving a Radical (A4)

SUTCLIFFE S NOTES: CALCULUS 2 SWOKOWSKI S CHAPTER 11

MTH 146 Class 16 Notes

On computing two special cases of Gauss hypergeometric function

Forced Oscillation of Nonlinear Impulsive Hyperbolic Partial Differential Equation with Several Delays

Math 2414 Activity 17 (Due with Final Exam) Determine convergence or divergence of the following alternating series: a 3 5 2n 1 2n 1

Transcription:

NOTES ON BERNOULLI NUMBERS AND EULER S SUMMATION FORMULA MARK WILDON. Beroulli umbers.. Defiiio. We defie he Beroulli umbers B m for m by m ( m + ( B r [m ] r r Beroulli umbers re med fer Joh Beroulli (he mos prolific Beroulli, d he discoverer of he Beroulli effec... Expoeil geerig fucio. If f(z z /! d g(z b z /! he f(zg(z c z /! where he coefficies c re give by c ( r r b r. Thus if we pu β(z B z /! he ( Br [z ]β(z exp z!! r ( r r B r + B! [ ] + B! from which i follows h β(z exp z z + β(z. Therefore ( β(z B z! z exp z. This shows h he rdius of covergece of β(z is π. So lhough i is empig o me subsiuio i log P (e log( e e m m m m e m m m his will o wor, s m is eveully more h π..3. Some vlues. If we dd z/ o boh sides of ( we obi (3 β(z + z z exp z + z/ coh z/, exp z which is eve fucio of z. So B, B / d B if 3 is odd. Some furher vlues re give below. 3 4 5 6 7 8 B / /6 /3 /4 /3 De: November 6, 4.

MARK WILDON Noe h he umerors eed o be ±. For isce B 5/66..4. Coecio wih FLT. Kummer proved h here re o posiive iegrl soluios of x p +y p z p wheever p is odd prime o dividig he umerors of y of he Beroulli umbers B, B 4,..., B p 3. Such primes re sid o be regulr. For isce s B 69/73, he prime 69 is o regulr..5. Esime for B. I [] i is oed h B 85453 38 > 69; ideed he Beroulli umbers re ubouded. The uhors expli h usig Euler s formul for co z (which hs ice elemery proof vi he so-clled Hergloz ric, oe ges z co z z π z β(z + z z coh z z + π + z m ( z for z < π. Hece, comprig coefficies of z, Thus π (! B ( ζ(. (π (π B (! 4 (π. Noe his shows h B /(! is decresig sequece. As hese re coefficies i he coverge geerig fucio β(z, B /(! eds o s eds o..6. Coecio wih uiform disribuio. Le X be disribued uiformly o [, ]. The he mome geerig fucio for X is Ee zx m e z d x (ex β(z. The coefficies of /β(z give us recurrece sisfied by he Beroulli umbers bu s z /β(z ( +! his is jus he defiig recurrece (.

NOTES ON BERNOULLI NUMBERS AND EULER S SUMMATION FORMULA 3.7. Sums of powers. We derive from firs priciples specil cse of he Euler summio formul. Le S m ( m + m +... + ( m m. Cosider he followig geerig fucio where m vries d is fixed: Ŝ(z m S m( zm m!. We hve ( Ŝ(z m z m m! (z m exp z m! exp z m m ( ( exp z. B m+ β(z z! z (m +! zm r Comprig coefficies of z m /m! gives m B r (4 S m ( m! r! m r+ m (m r +! r B r m +.8. Beroulli polyomils. Defie B m (z ( m B z m. m ( m + r m+ r. The firs few Beroulli polyomils re B (z, B (z z, B (z z z +, B 6 3(z z 3 3z +, o so o. Noe h B m( B m d h B m+ ( ( m + B m+ (m + S m ( + B m+ so S m ( B m+( B m+ ( m + which is he obvious defiiio for S m (z. Some useful properies: B m ( ( m B m ( m B + B m [m ] + B m so B m ( B m ( B m uless m i which cse B ( B. For m we hve B m (z ( m B (m z m m ( m B z m mb m (z.. The Euler summio formul From ow o le f : R R be fucio wih s my derivives s eeded. Euler s summio formul is:

4 MARK WILDON Theorem. (Euler. Le m d le < b be iegers. The b m B (5 f( f(xdx +! f ( (x b + R m <b where he remider R m is give by R m ( m+ b B m ({x} f (m (xdx. m!.. Applicio o sums of powers. Se f(x x r d e m > r. The f (m so he remider erm vishes. Puig d b we obi r r+ r+ r + + which grees wih (4. B r+! r r+ B r + ( r + r+.. Firs proof. As ll he erms elescope icely i is sufficie o prove he formul whe d b. This hs he dvge h we c use B m (x rher h B m ({x}. We proceed by iducio o m. The cse m ses h f( f(xdx (f( f( + (x f (xdx. This follows from simple iegrio by prs: (x f (xdx f( For m > we c wrie he righ-hd-side s m f(xdx + B! (f ( ( f ( ( + B m m! (f (m ( f m ( + ( m+ Applyig he iducive hypohesis gives: f( + ( m+ B m (x (m! f (m (x + B m m! (f (m ( f (m ( + ( m+ f(xdx (f( f(. B m (x f (m (xdx m! B m (x f (m (xdx m! Iegrig he fil erm by prs usig he resuls i.8 gives f( + ( m+ B m (x (m! f (m (x + B m m! (f (m ( f (m ( + ( m+ B m(x f (m (x m! + B m (x ( m (m! f (m (xdx. The ouermos wo erms obviously ccel. Ad so do he ier wo, s if m is odd, he s m >, B m.

NOTES ON BERNOULLI NUMBERS AND EULER S SUMMATION FORMULA 5.3. A lerive proof. There used o be secio here climig lerig proof by ig sequece of polyomils covergig uiformly o f o he iervl [, b] d usig uiform covergece o ierchge he iegrl wih he limi of he polyomils: his reduces o he cse where f(x x r. However i is o loger possible o reduce o he cse where d b, d he mi iducive sep i he proof is o subsilly simplified from he previous secio. I is however worh oig h if m > r d we e d b N he he resul comes ou very esily. I his cse Euler s formul sys h <N Usig he ideiy s ( r r N r+ r + + N r+ r + + r+ r r B! r N r ( ( r+, he righ-hd side c be rewrie B r + which is equl o <N r by (4. ( r + N r+.4. Esimes for error erm. We c se Euler s summio formul i he followig form (6 <b f( b f(xdx (f(b f(+ m B (! f ( (x b +R m. I c be show h B m ({x} B m so rough esime for he error erm is give by.5: R m B m (m! b f (m (x dx. If f (m is posiive he his shows h he mgiude of he error erm is mos he mgiude of he fil erm i he sum. Noe h very ofe he remider erm R m will o ed o s he upper limi i he summio, b eds o. However i will ofe hve some o-zero limi. I his cse we hve R m R m ( b B m ({x} f (m (xdx (m! where he righ-hd-side eds o R m ( s b. I is proved i [] p475 h if f (m+ d f (m+4 re posiive for x [, b] he B m+ R m θ m (m +! f m+ (x b where θ m. So he remider erm lies somewhere bewee d wh would hve bee he ex erm i he sum.

6 MARK WILDON 3. Exmples of Euler summio 3.. Firs exmple. We shll emp o use Euler summio o fid S where <. Of course we lredy ow he swer, S ( /(, d i fc Euler summio urs ou o be very ieffecive wy of fidig i! Sill here re some pois of ieres. If f(x x he f (r (x (log r x so by equio (6 we hve ( m x B dx ( + (! (log x + R m ( m log B (log + R m (! for d m. Now somehig lile suprisig hppes: if > e π he log < π d so we c use he expoeil geerig fucio for he Beroulli umbers o obi B (! (log log + d so + ( m+ B (! (log + R m As we ow he swer, he ls equio implies h R m s m. Less rificilly, we c emp o prove his direcly. We fid R m (log m B m ({x} x dx (m! (log m B m (m! ( (log m 4 (π m ( So he error erm R m s m, provided e π < < e π, which holds s we hve lredy ssumed h e π < <. The limiig behviour wih respec o is lso of ieres. Kowig he limi of S llows us o see h lim R m m+ B (log (! This is exmple of he poi mde erlier, h while he remider erm will o usully ed o s eds o, i my well hve some limiig vlue.

NOTES ON BERNOULLI NUMBERS AND EULER S SUMMATION FORMULA 7 3.. Summig squre roos. Le f(x x. Noe h f(xdx ( 3 d h 3 ( f (r (x r! x r. r Euler s summio formul gives ( 3 3 + + + m ( B ( 3 + R m. Seig m we ge ( 3 3 + + + ( + R. We c esime R by usig he boud B ({x} < B 6 : So we hve R 6.8 x 3 4 ( 3 ( 3 + + C + O( for some cos C. I geerl he remr o he previous pge shows h he error is give by ( B m+ R m θ m ( m C(m + O( m m + m where θ m [, ] d C(m does o deped o, This gives ( 3 m ( B 3 + + C(m + 3 + O( m. There is smll simplificio: C(m c be deermied by ig limi wih respec o, so ( C lim ( 3 3 d so does o deped o m. (Exercise 9.7 i [] revels h C ζ( ; i fc he defiiio of ζ(α for α >. 3.3. A esime for P (x. log P (e log( e m m e m m e m. This suggess wys we migh pply Euler summio. (Bu hey ur ou o be equivle.

8 MARK WILDON 3.3.. Te from Kuh fscicle Exercise 5. Le f(x log( e x e mx m. m The we hve log P (e f(. The iegrl of f is give by he Li fucio: x f(udu m e mu m x m e m m m e mx m Li (e Li (e x. The derivives of f re coeced wih he Euleri umbers: le m deoe he umber of permuios i S wih excly m sces (or desces. The we clim Lemm 3.. For we hve f ( (x e x ( ( e x e x. Proof. By iducio o m. If m he he formul is redily verified. Suppose rue for m. The f (m+ (x me x ( m ( e x m+ e x + e x ( ( ( e x e x e x ( + ( e x + e x ( + ( e x + ( ( ( To fiish we pply he ideiy ( m m + ( e x + + (m + m e x e x m. We re ow redy o pply Euler summio. Tig m i (6 we ge f( Li (e Li (e (log( e + log( e + e x ( e x + R. This ideiy c be proved s follows. Le g S hve excly m sces. If ppers i posiio... b... or... he removig gives permuio i S wih m sces. Coversely, give such permuio, here re m desces, d puig bewee y umbers ivolved i desce, or he ed, gives permuio i S wih m sces. Oherwise we hve... b... or... i which cse removig does o chge he umber of sces. Coversely give permuio i S wih m sces, puig bewee he wo umbers ivolved i sce, or he sr, does o chge he umber of sces.

NOTES ON BERNOULLI NUMBERS AND EULER S SUMMATION FORMULA 9 where s he eve derivives re posiive e x ( R e x So leig we ge e e. e log P (e Li (e log( e + O( where O( sds for quiy lwys lyig i [, ]. If we use he ideiy suggesed by Kuh, mely or, equivlely, we ge Li (x + Li ( x ζ( log x log( x Li (e log P (e ζ( ζ( ζ( Li ( e Li ( e + ζ( + log( e 3 log( e + O( log( e + O( ; where O( [, ] log + O( ; where O( [ log, log ]. As Li (y log( y gives Li ( e d log( e log + log( / +... log + O(. I priculr by beig lile more creful wih he O( errors we c ge lower boud for log P (e : π 6 + log log log P (e π 6 + log + log. I fc he righ cos is log( π.989; his follows from log P (e π 6 + log log π + O( where he O( erm is (by he fuciol equio, 4 + log P (e 4π /. Noe h.693 <.989 <.693. 3.3.. Alerive. We migh lso ry o use Euler summio o sum log( e. Perhps surprisigly his urs ou o be equivle o he previous pproch, sice d dx log( e x e x e x e x which is esseilly he fucio we summed erlier. Refereces [] Rold L. Grhm, Dold E. Kuh, d Ore Pshi. Cocree Mhemics. Addiso Wesley, 994.