Esope source files

   1 : C INTFR2    SOURCE    CHAT      05/01/13    00:41:05     5004           
   2 :       SUBROUTINE INTFR2(CC,SS,XF,XD)                                            
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C                                                                               
   6 : C =====================================================================         
   7 : C SOUS-PROGRAMME FORTRAN APPELE PAR ELFGR1 POUR CALCULER LES INTEGRALES         
   8 : C                                                                               
   9 : C                  / XF                                                         
  10 : C   CC=1/SQRT(2PI)/   COS(T)/SQRT(T) DT                                         
  11 : C                / XD                                                           
  12 : C                                                                               
  13 : C                  / XF                                                         
  14 : C   SS=1/SQRT(2PI)/   SIN(T)/SQRT(T) DT                                         
  15 : C                / XD                                                           
  16 : C                                                                               
  17 : C         AVEC      1 < XF ET XD < 10                                           
  18 : C                                                                               
  19 : C   PROGRAMMEUR : GUILBAUD                                                      
  20 : C   CREATION    : 13/01/87                                                      
  21 : C =====================================================================         
  22 : C                                                                               
  23 : -INC CCREEL                                                                     
  24 :       DEUPI=2.D0*XPI                                                            
  25 :       DX2=XPI/2000.D0                                                           
  26 : C                                                                               
  27 : C  INTEGRATION PAR LA METHODE DES TRAPEZES                                      
  28 : C                                                                               
  29 :       KMAX=IDINT((XF-XD)/DX2)                                                   
  30 :       IF(KMAX.EQ.0) KMAX=1                                                      
  31 :       DX=(XF-XD)/(KMAX+1)                                                       
  32 :       SSX=SIN(XD)*0.5D0/SQRT(XD)                                                
  33 :       SCX=COS(XD)*0.5D0/SQRT(XD)                                                
  34 :       XX=XD                                                                     
  35 :       DO 20 K=1,KMAX                                                            
  36 :       XX=XX+DX                                                                  
  37 :       SSX=SSX+SIN(XX)/SQRT(XX)                                                  
  38 :       SCX=SCX+COS(XX)/SQRT(XX)                                                  
  39 :   20  CONTINUE                                                                  
  40 :       SS=(SSX+0.5D0*SIN(XF)/SQRT(XF))*DX/SQRT(DEUPI)                            
  41 :       CC=(SCX+0.5D0*COS(XF)/SQRT(XF))*DX/SQRT(DEUPI)                            
  42 :       RETURN                                                                    
  43 :       END                                                                       
  44 :  
  45 :  
  46 :