Tuesday, 17 October 2006

ca.analysis and odes - A two-variable Fourier series and a strange integral

I know how I want to answer this question. I'll write up the easy parts here, and leave the hard part for you :).




First some minor changes. It will be convenient to clear out denominators and work with logleft(4+e2piix+e2piix+e2piiy+e2piiyright). That just changes the constant term of your Fourier series by log2. Next, it is convenient to focus on
int01int01logleft(4+e2piix+e2piix+e2piiy+e2piiyright)e2piimxe2piinydxdy.
A simple linear transformation goes between this and the cosine formulation. Let S=(z,w):|z|=|w|=1. So we are dealing with
frac1(2pii)2intSlogleft(4+z+z1+w+w1right)zm1wn1dzdw.
Dropping out the 4pi2, we want to show the integrand is of the form api+b.



UPDATE: Thanks to fedja for pointing out that I had oversimplified the next paragraph.



Assuming that (m,n)neq(0,0), we can integrate by parts with respect to one of the two variables, let's say z. Once we do that, we will have a quantity of the form
(mboxrationalnumber)cdotintSfrac(zz1)wkzelldwdz4+w+w1+z+z1



So we'd like to show this quantity is of the form a+bpi.



As fedja points out, we need to be careful here. Without the zz1 term, the integral diverges like intintdsdt/(s2+t2) near (1,1).




Whew! Now comes the actual hard part. Let
E:=(z,w)in(mathbbC)2:4+z+z1+w+w1=0.
This is an elliptic curve with four punctures. As Bjorn points out, this is a nodal cubic and can be parameterized as
(z,w)=left(frac1uu(1+u),fracu(u1)1+uright).
We'll come back to this point later.



The 2-form dwdz/(4+z+z1+w+w1) has a simple pole on E. Let omega be the 1-form on E which is the residue of that 2-form.



I think there should be a curve gamma in E such that S is homotopic, in (mathbbC)2setminusE, to a tubular neighborhood of gamma. So
intSfracwkzell(zz1)dwdz4+w+w1+z+z1=intgammaomegawkzell(zz1).



If we substitute in the above parameterization, this will be the integral around a closed loop of some rational function in mathbbQ(u). In particular, we can compute this integral by residues and we will get something of the form a+bpi, as desired.



Actually, it looks to me like we should just get bpi. Maybe the integration by parts doesn't go as well as I hoped?




Obviously, someone should actually work this out explicitly, but I don't think it will be me.

No comments:

Post a Comment