Serie de fourier double

[jgorsse]

Je dois develloper un chargement p(x,y) en serie de fourier pour resoudre un probleme de mecanique des deformations.
p(x,y) est defini sur :
-a/2-b/2p(x,y)=p0 sur un disque de centre (0,0) et de rayon R et 0 ailleurs.
De part les symetries on a:
p(x,y)= sum(sum(pmn*cos(m*Pi*x/a)*cos(n*Pi*y/b),n=0..inf),m=0..inf)
avec
pmn=4/(a*b)*int(int(p(x,y)*cos(m*Pi*x/a)*cos(n*Pi*y/b),y=-b/2..b/2)x=-a/2..a/2).
Mon probleme se situe dans le calcul des pmn :mur: si vous savez comment faire... :help:

a+