Friday, 23 March 2007

ca.analysis and odes - spiral of Theodorus

Here's a sketch of a proof that the constant you want exists, and how to find it.



Let
f(n)=arctan(1)+arctan(1/sqrt2)+arctan(1/sqrt3)+ldots+arctan(1/sqrtn).
You want to show that f(n)=sqrtn+C+o(1) for some constant C. (If you're not familiar with the o-notation, think of o(1) as representing some function which goes to 0 as n goes to infinity.)



Then take the power series expansion of arctan(1/sqrtk); this is



()       k1/2frac13k3/2+frac15k5/2+ldots



So summing over 1 to n, we should get
begin{align*}
f(n) = & (1^{-1/2} + 2^{-1/2} + ... + n^{-1/2}) \
- , frac{1}{3} &(1^{-3/2} + 2^{-3/2} + ... + n^{-3/2}) \
+ , frac{1}{5}& (1^{-5/2} + 2^{-5/2} + ... + n^{-5/2}) - ldots
end{align*}
Now, 11/2+21/2+ldots+n1/2 has the asymptotic form
2sqrtn+zeta(1/2)+O(n1/2)
where I cheated a bit and asked Maple, and zeta is the Riemann zeta function. And 1j/2+2j/2+ldots+nj/2 has the asymptotic form
zeta(j/2)O(nj/2+1)
where, if you're not familiar with the O-notation, O(nj/2+1) should be thought of as a function that goes to zero at least as fast as nj/2+1 as n goes to infinity. So, assuming that we can rearrange series however we like,
f(n)=2sqrtn+(zeta(1/2)frac13zeta(3/2)+frac15zeta(5/2)ldots)+o(1).
Since zeta(s) is very close to 1 when s is a large real number, that alternating series should converge. Again cheating and using Maple, I claim it converges to about 2.157782997. This is the constant you call varphi, and what you called K is equal to 2. (An easier way to see that your K is 2 is to note that arctan(1/sqrtn) is about 1/sqrtn, and approximate the sum by an integral.

No comments:

Post a Comment