Sources Esope | Cast3M

intfr2
C INTFR2    SOURCE    CB215821  17/11/30    21:16:29     9639           
      SUBROUTINE INTFR2(CC,SS,XF,XD)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
C =====================================================================
C SOUS-PROGRAMME FORTRAN APPELE PAR ELFGR1 POUR CALCULER LES INTEGRALES
C
C                  / XF
C   CC=1/SQRT(2PI)/   COS(T)/SQRT(T) DT
C                / XD
C
C                  / XF
C   SS=1/SQRT(2PI)/   SIN(T)/SQRT(T) DT
C                / XD
C
C         AVEC      1 &lt; XF ET XD &lt; 10
C
C   PROGRAMMEUR : GUILBAUD
C   CREATION    : 13/01/87
C =====================================================================
C
-INC CCREEL
      DEUPI=2.D0*XPI
      DX2=XPI/2000.D0
C
C  INTEGRATION PAR LA METHODE DES TRAPEZES
C
      KMAX=INT((XF-XD)/DX2)
      IF(KMAX.EQ.0) KMAX=1
      DX=(XF-XD)/(KMAX+1)
      SSX=SIN(XD)*0.5D0/SQRT(XD)
      SCX=COS(XD)*0.5D0/SQRT(XD)
      XX=XD
      DO 20 K=1,KMAX
      XX=XX+DX
      SSX=SSX+SIN(XX)/SQRT(XX)
      SCX=SCX+COS(XX)/SQRT(XX)
  20  CONTINUE
      SS=(SSX+0.5D0*SIN(XF)/SQRT(XF))*DX/SQRT(DEUPI)
      CC=(SCX+0.5D0*COS(XF)/SQRT(XF))*DX/SQRT(DEUPI)
      RETURN
      END