Sources Esope | Cast3M

somme2
C SOMME2    SOURCE    CHAT      05/01/13    03:22:06     5004
      FUNCTION SOMME2 (X,Y,L)
      IMPLICIT INTEGER(I-N)
      REAL*8 SOMME2
************************************************************************
*
*                             S O M M E 2
*                             -----------
*
* FONCTION:
* ---------
*
*     INTEGRATION, PAR LA METHODE DES TRAPEZES, D'UNE FONCTION
*     REPRESENTEE PAR UNE SUITE DE COUPLES DE VALEURS (X,F(X)), DANS
*     L'ORDRE DES "X" CROISSANTS.
*
* PARAMETRES:   (E)=ENTREE   (S)=SORTIE   (+ = CONTENU DANS UN COMMUN)
* -----------
*
*     X       (E)  LISTE DES ABSCISSES.
*     Y       (E)  LISTE DES ORDONNEES.
*     L       (E)  NOMBRE DE POINTS DE DEFINITION DE LA FONCTION.
*     SOMME2  (S)  INTEGRALE.
*
      INTEGER L
      REAL*8 X(*),Y(*)
*
* VARIABLES:
* ----------
*
      REAL*8 XINT
*
* AUTEUR, DATE DE CREATION:
* -------------------------
*
*     PASCAL MANIGOT     22 MARS 1988
*
* LANGAGE:
* --------
*
*     FORTRAN77
*
************************************************************************
*
      L1 = L - 1
      XINT = 0.D0
      DO 20 I=1,L1
         XINT = XINT + (X(I+1)-X(I)) * (Y(I)+Y(I+1)) / 2.D0
  20     CONTINUE
*     END DO
*
      SOMME2 = XINT
*
      END