20130320, 01:05  #1 
Dec 2003
Hopefully Near M48
2×3×293 Posts 
Can Euler's Basel Solution be Saved?
I'm thinking of giving a talk to an undergrad math club, and I had the idea of presenting Euler's solution to the Bessel Problem (1 + 1/4 + 1/9 + ... = pi^2/6) as the topic. While there are some very grave gaps, it a deliciously elegant solution, and it would be great if there were a way to repair those gaps. My (very rough draft of a) plan is to present Euler's solution as I read it, explain what the gaps are, and then try to fix them. While I know how to repair some of the problems, there are some I don't, which is why I'm posting this.
Anyway, here's Euler's solution: 1) Lemma 1: For all x, sin(x) = x  x^3/3! + x^5/5!  ... Proof: Taken for granted. 2) Lemma 2: Let f be a complex polynomial with f(0) = 1 and distinct roots r_1, ... r_n. Then f(x) = (1  x/r_1)(1 x/r_2)...(1x/r_n) Proof: Since f has n distinct roots, deg(f) = n. We have specified the values of f at n+1 distinct inputs; it is easy to check that there is at most one f of degree n which can satisfy those conditions. The given f satisfies those conditions. 3) Define sinc(x) = sin(x)/x when x =/= 0 and 1 when x = 0. Then sinc(x)= 1 x^2/3! + x^4/5!  ... 4) Clearly, sinc(0) = 1 and the roots of sinc(x) are just the roots of sin(x) with 0 removed (distinct because the roots of sin(x) are distinct), namely +/npi for n = 0, 1, 2, etc. 5) Regard sinc(x) as an (infinite) polynomial. By Lemma 2, sinc(x) = (1x/pi)(1+x/pi)(1x/2pi)(1+x/2pi)... 6) So (1x/pi)(1+x/pi)(1x/2pi)(1+x/2pi)... = 1  x^2/3! + x^4/5!  ... 7) (1x^2/pi^2)(1x^2/4pi^2)(1  x^2/9pi^2)... = 1  x^2/3! + x^4/5!  ... 8) Multiplying out the lefthand side and collecting like terms, we get 1 + (1/pi^2  1/4pi^2  1/9pi^2  ...)x^2 + O(x^4) = 1  x^2/3! + x^4/5!  ... 9) Equate x^2 coefficients to get 1/pi^2  1/4pi^2  1/pi^2 = ... = 1/3! 10) Multiply both sides by pi^2 to get the desired result.  Problems with the solution and my ideas for fixing them  4) We need to look for all complex zeroes of the function, not just the wellknown ones on the real line. But this is easily fixed by using the formula sin(z) = e^(iz)  e^(iz)/2i, setting the numerator equal to zero, splitting z into real and imaginary parts, and solving. It is quickly seen that Im(z) = 0, which brings us back to the wellknown case. 5) This is the biggest problem. sinc(x) is of course an entire function, not a polynomial; so we can't legally just apply the lemma. To make matters much worse, the conditions in Lemma 2 aren't even enough to uniquely determine an entire function. Replacing sinc(x) with e^(h(x))sinc(x) for any entire function h with h(0) = 0 yields another function that satisfies exactly the same conditions; but is a very different function and of course will have a very different power series. What I need is an analogue of Lemma 2 for entire functions that allows me to uniquely pin down the function using its roots, its value at the origin, and some extra data that needs to be easy to get for sinc(x). 7) Warning: Messing with an infinite product in such a nontrivial way is in general not ok. But in this case, it is ok assuming we have patched up the proof up till Line 6 because we know that the infinite product is convergent. In a convergent infinite product, multiplication is associative because the new sequence of partial products is just a subsequence of the original sequence of partial products. 8) Multiplying out an infinite product like that seems kind of fishy to me; can't think of a theorem to justify it off the top of my head. But I would think it's something that can be justified using standard undergrad analysis techniques. Am I right?  So in summary, the main thing I'm missing is a way to justify step 5. A justification for step 8 would also be nice. Can it be done? Thanks 
20130320, 09:06  #2 
"Nathan"
Jul 2008
Maryland, USA
5×223 Posts 
Nitpickoftheday: It's actually the Basel problem (after Euler's hometown of Basel, Switzerland). There are, of course, Bessel functions, and knowing Euler, he may well have been involved with them too, but they play no part in this problem. Anyway...
Take a look at this paper/talk on the same subject. It discusses a lot of the complex analysis behind this problem, and specifically, states that the issue in your step (5) is resolved via the Weierstrass Factorization Theorem, which is a stronger result that comes from the Fundamental Theorem of Algebra. I also looked at Wikipedia's article on the Basel problem, and it explains that your step (8) is legal because of Newton's identities. The latter article does not resemble anything that I have seen in an undergraduate curriculum, so you might want to just make mention of the buzzword "Newton's identities" in explaining why step (8) is possible. It's worth mentioning to the "kids" that the upshot is that, as with many problems in the earliest days of calculus and the study of the infinite and infinitesimal, mathematicians were often able to discover these beautiful and unexpected results, but were unable to put them on any sort of rigorous foundation for at least another 100+ years when analysis would be formally established. Nonetheless, having results such as this one would also prove useful as a check on the new rigorous methods that were being implemented. Another interesting point to mention: Why doesn't this work for the odd powers? What happens if we try to use the cosine series instead of the sine series? Zeta(3) is irrational [cf. Apery, 1978]...but that's about all we know... Obviously, this result also gives us an approximation for pi (which was part of a talk that I gave at a student session at MathFest in 2004): the quantity sqrt(6 * sum(1/n^2)) should approximate pi for large enough n. How large? How fast is the convergence? (Hint: it's not too bad!) You might also recommend that they check out William Dunham's books Journey through Genius and Euler: Master of Us All. Both of these books provide nice sketches of Euler's work on this problem (he didn't stop with Zeta(2); rather, he calculated Zeta(n) for even n as high as 18, long before pocket calculators), as well as discussing how he played fast and loose with infinite series. Good luck! This is a really neat topic and a fine introduction to Euler's amazing and prodigious mathematics. 
20130320, 12:57  #3 
Dec 2003
Hopefully Near M48
2×3×293 Posts 
Thanks for the links. I already read Journey through Genius in high school; that's where I learned of this solution!
Going to have to think about the Weierstrass Factorization and it's application to this case; but as for the other concern, I still don't quite see how to prove the validity of (8) from Newton's identities. I mean, I'm comfortable with applying Newton's identities to get the coefficients of x^2, x^4, x^6 etc. if we had a finite product. The problem here is the product being infinite. Seems to me it boils down to three steps: 8.1) The sequence of partial products (1x^2/pi^2), (1x^2/pi^2)(1x^2/4pi^2), (1x^2/pi^2)(1x^2/4pi^2)(1x^2/9pi^2), etc. converges to sinc(x). [Pointwise this is ok, but is that going to be enough?] 8.2) Use Newton's identities to find the x^2 (and x^4, x^6, etc. if I feel like it) coefficients of each of these partial products. 8.3) Argue that if {f_n} is a sequence of entire functions converging (is pointwise enough?) to an entire function f, then the x^2 (and higher) coefficients of the Taylor expansions of f_n converge to the x^2 coefficient of the Taylor expansion of f. [Come to think of it, I'm pretty sure I heard that in the real setting, even f_n > f uniformly does not imply f_n' > f'. In the complex setting, err... Not sure...] Thanks Last fiddled with by jinydu on 20130320 at 13:10 
20130320, 16:03  #4 
"Lucan"
Dec 2006
England
2×3×13×83 Posts 
Bessel

20130320, 21:15  #5  
Basketry That Evening!
"Bunslow the Bold"
Jun 2011
40<A<43 89<O<88
3·29·83 Posts 
I can't say much on this particular solution (I haven't had analysis yet), but I can say that another method I've seen to calculate this series (and which generalizes well to larger even n, I believe) is a combination of Fourier analysis/linear algebra/Parseval's identity (the latter is essentially the Pythagorean theorem for an arbitrary inner product space), as well as being quite accessible to undergrads with little analysis.
The following were questions on a lin. alg. quiz I took in high school: Quote:
For problem 4), we apply Parseval's identity with our orthonormal basis being . We also note that (where ), or rewriting to put it in the form of Parseval's identity, . Thus . Meanwhile, , so we have , or . The only squishy part (AFAICT) is proving Parseval's identity (esp. in the infinite dimensional case). We can also see why we get only even powers  because of the squaring of the fourier coefficients. It ought to be easy (though I haven't really thought about it) to use f(x)=x^2, x^3, ... to get . (The only part that I haven't confirmed is that the Fourier coefficients take the form I think they take for the higher order f's.) Last fiddled with by Dubslow on 20130320 at 21:20 Reason: (esp. in the infinite dimensional case) 

Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Intermediate files not saved when service is stopped  mPat  Software  0  20170326 10:42 
GHz days saved update  Chuck  GPU to 72  5  20120222 06:22 
Euler (6,2,5) details.  Death  Math  10  20110803 13:49 
Two worker threads, only one is saved  Unregistered  Information & Answers  7  20110331 04:57 
r***** saved file  bayanne  Software  5  20061108 19:01 