C ELPDK1    SOURCE    CHAT      05/01/12    23:36:51     5004
      SUBROUTINE ELPDK1(XLP0
     &                  ,XD,XNU
     &                  ,NTRAP2,CS1,CB)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-B,D-H,O-Z)
      IMPLICIT COMPLEX*16(C)
************************************************************************
*
*    PLAQUES PAR EQUATION INTEGRALE:
*
*       ON INTEGRE PAR LA METHODE DES TRAPEZES EN STATIQUE
*       ET EN DYNAMIQUE
*    op (w-ws)
************************************************************************
*
*
          DIMENSION P0(2)
          DIMENSION QI(2)
          DIMENSION XOP(8)
          DIMENSION COPD(8)
*
          DIMENSION XN (2)
*
          DIMENSION XT (2)
*
*
          DIMENSION CS1 (8)
*
* 1-   TAILLE DES TRAPEZES ( DEMI SEGMENT SEULEMENT)
*
        X0 = XLP0/(4*NTRAP2)
        X1 = XLP0 /(2*NTRAP2)
*
* 2-   POSITIONS DE REFERENCE
*
        P0(1) = 0D0
        P0(2) = 0D0
        XN(1) = 1D0
        XN(2) = 0D0
        XT(1) = 0D0
        XT(2) = 1D0
        QI(1) = 0D0
        QI(2) = 0D0
*
* 2-   INTEGRATION
*
         DO 10 J=1,8
             CS1(J)  = 0D0
             COPD(J) = 0D0
10       CONTINUE
         DO 100 I=1,NTRAP2
           QI(2) = X0 + (I-1)*X1

           CALL ELPDKP(P0,XN,QI,XN,XT,XD,XNU,COPD,CB)

           DO 110 J = 1,8
             CS1(J) = CS1(J) + 2*( X1 * COPD(J))
110        CONTINUE
100      CONTINUE
*
       RETURN
      END