Sources Esope | Cast3M

intega
C INTEGA    SOURCE    CHAT      05/01/13    00:40:08     5004
      SUBROUTINE INTEGA(G,F2,A,B,C,D,N,D2,CO)
C----------------------------------------------------------------------
C
C   CALCUL DU FACTEUR RELATIF A DEUX SEGMENTS EN 3D
C   UNE INTEGRATION ANALYTIQUE
C
C----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION G(11,10),A(3),B(3),C(3),D(3)
      F2=0
      N1=(N+1)/2
      DO 10 I=1,N1
      M=2*I-1
      YAMB1=(G(M,N)+1)/2
      CALL PROGA(F1,G,YAMB1,A,B,C,D,D2,CO)
      F2=F2+F1*G(M+1,N)
          IF(G(M,N).NE.0)THEN
          YAMB1=(1-G(M,N))/2
          CALL PROGA(F1,G,YAMB1,A,B,C,D,D2,CO)
          F2=F2+F1*G(M+1,N)
          ENDIF
 10   CONTINUE
      F2=F2/2
      RETURN
      END