Saturday, 16 March 2013

calculus - Two sums for pi




I was looking for a simple way to evaluate the integral 0sinxxdx ( a belated look at this question). There are symmetries to be exploited, for one thing. So I had an idea, but the idea depends, optimistically, on two expressions for π that I cannot prove, and my question is whether anyone sees reasonable proofs of either expression (1, 2 below).



Put



0sinxxdx=π0k=k=0{sinxx+2kπsinxx+2kπ+π}dx



The intuition behind this step is that as we make our way around the unit circle from 0 to π, for each angle α we have an angle α+π for which the value of sinxxwill be smaller and negative. No need to integrate beyond π, and for angles greater than 2π we are simply dividing up the original integral's range by adding multiples of π to the value of x.



WLOG we can form Riemann sums by dividing the range of integration into an even number of equi-angle subranges dxn. For example, if n=6, we have




π0k=k=0sinx{1x+2kπ1x+2kπ+π}dx



{k=06n=1sin(nπ6)(1π)(1n6+2k1n6+2k+1)}{π6}.



Taking this as true, we note that pairs of summands of the Riemann sums symmetric about π2 sum to 1, and the value of the summand at π2 is 1/2. Again, taking this as true, we have that



f(xn)=12+1n22+sinπf(π)=n12,


and
f(xn)dxn=(n12)πn.




So for the approximation above for n=6, we expect a value of about 5π12 for high k.



In words, we are only counting n - 1 summands, since the last summand at π is 0. Since we divided the interval into n subintervals, we are stuck with the sum above and the conditional limit (conditioned at least on the truth of the two limits below):



0sinxxdx=limn=1n12πn=π2.



The problem is this: we need to prove, to begin with, that



k=0(1π/2+2πk1π/2+2kπ+π)=12




or



1π11/2+2k11/2+2kπ+1=12,



that is,



(1) k=043+16k+16k2=π2.



We must also prove for 0<j<n that




(2) sin(jπn)(k=0(1jπn+2πk1jπn+2kπ+π)+k=0(1(nj)πn+2πk1(nj)πn+2kπ+π))=1.



or , factoring 1π out of the sum and mutiplying:



sin(jπn)(k=0(1jn+2k1jn+2k+1)+k=0(1(nj)n+2k1(nj)n+2k+1))=π.



(1) might be a special case of something at Wolfram's site, but I didn't spot it. (2) looks messy, but maybe induction?



If this is otherwise correct, I think one could work back from the relations for π and establish the value of the integral. It is important not to forget that k starts at 0.




Edit: the following is equivalent to (2) and maybe easier to scan. For $0

(2)k=02(sinπr)(1+4k+4k2r+r2)(1+2k+r)(1+2kr)(2+2kr)(2k+r)=π



This was derived assuming r rational, but I don't think it matters.


Answer



Nk=01α+2k=Ψ(N+1+α/2)Ψ(α/2)2



where Ψ(t)=ln(t)+O(1/t) as t, so
k=0(11/2+2k13/2+2k)=Ψ(3/4)Ψ(1/4)2=π2

However, I think the evaluation of Ψ(3/4)Ψ(1/4) may be just as hard as your integral.



EDIT: It turns out that Ψ has the integral form
Ψ(t)=0(es1(s+1)t)dss


so that
Ψ(3/4)Ψ(1/4)=0((s+1)3/4(s+1)1/4)dss=4arctan((s+1)1/4)|0=π


No comments:

Post a Comment

real analysis - How to find limhrightarrow0fracsin(ha)h

How to find limh0sin(ha)h without lhopital rule? I know when I use lhopital I easy get $$ \lim_{h\rightarrow 0}...