Esope source files

   1 : C INTEGA    SOURCE    CHAT      05/01/13    00:40:08     5004           
   2 :       SUBROUTINE INTEGA(G,F2,A,B,C,D,N,D2,CO)
   3 : C----------------------------------------------------------------------         
   4 : C                                                                               
   5 : C   CALCUL DU FACTEUR RELATIF A DEUX SEGMENTS EN 3D
   6 : C   UNE INTEGRATION ANALYTIQUE                                                  
   7 : C
   8 : C----------------------------------------------------------------------         
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8 (A-H,O-Z)
  11 :       DIMENSION G(11,10),A(3),B(3),C(3),D(3)
  12 :       F2=0
  13 :       N1=(N+1)/2
  14 :       DO 10 I=1,N1
  15 :       M=2*I-1
  16 :       YAMB1=(G(M,N)+1)/2
  16 :       CALL PROGA(F1,G,YAMB1,A,B,C,D,D2,CO)
  17 :       F2=F2+F1*G(M+1,N)
  18 :           IF(G(M,N).NE.0)THEN
  19 :           YAMB1=(1-G(M,N))/2
  19 :           CALL PROGA(F1,G,YAMB1,A,B,C,D,D2,CO)
  20 :           F2=F2+F1*G(M+1,N)
  21 :           ENDIF
  22 :  10   CONTINUE
  23 :       F2=F2/2
  24 :       RETURN
  25 :       END
  26 :