THESE - Université Ferhat Abbas de Sétif

CALL dqdag(FEIM,Xa,Xb,ERRABS,ERRREL,IRULE,SOMI1,ERREST(
YB=1.071046d0*RSZ
Xa=YA*WL
Xb=YB*WL
CALL dqdag(FEI,Xa,Xb,ERRABS,ERRREL,IRULE,SOMI2,ERREST(
YC=1.773407d0*RSZ
Xa=YB*WL
Xb=YC*WL
CALL dqdag(FEIM,Xa,Xb,ERRABS,ERRREL,IRULE,SOMI3,ERREST(
Xa=YC*WL
Xb=4.0d0*WL*DFLOAT(2*LL+1)
CALL dqdag(FEI,Xa,Xb,ERRABS,ERRREL,IRULE,SOMI4,ERREST(
EI=SOMI1+SOMI2+SOMI3+SOMI4
FI=(ER*ER+EI*EI)
RETURN
end
c calcul de la fonction a integrer en cos
FUNCTION FER(XX)
IMPLICIT DOUBLE PRECISION (a-h,o-z)
c REAL*8 DBSJ0,ER,LAGO2,FER
DOUBLE PRECISION LAGO2
EXTERNAL DBSJ0,LAGO2,RJB
COMMON/PAM1/WL,Z,alambd,ft,DRL,LL
COMMON/PAM3/R,RSZ,W0
COMMON/PHASE/GOUY
RAU=XX
c XL est la variable R normalise pour le poly de Laguerre
XL=2.0d0*RAU*RAU/WL/WL
pi=4.d0*datan(1.d0)
C R est la variable radiale dans le plan z
c RAU est la variable radiale dans le plan de l"EOD
c calcul avec lentille
AGRC=pi*RAU*RAU*(DRL+(1.0d0/Z)-(1.0d0/ft))/alambd
c calcul sans lentille
c AGRC=pi*RAU*RAU*(DRL+(1.0d0/Z))/alambd
c CVAR=DCOS(GOUY*DFLOAT(2*LL-1)+AGRC(
CVAR=DCOS(AGRC)
!!!!IMPORTANT : l'introduction du déphasage de Gouy perture totalement
!!le calcul. En particulier on ne retrouve pas qu'une gaussienne qui
cc traverse une lentille est une gaussienne transformée mais une
cc gaussienne tout de même. De plus la phase de Gouy ne dépendant pas
cc de la coordonnée radiale n'apporte rien au calcul de diffraction.
ANG=2.d0*pi*R*RAU/alambd/Z
CDV=(W0/WL)*2.0d0*pi/alambd/Z
95