Sources Esope | Cast3M

intqp2
C INTQP2    SOURCE    CHAT      05/01/13    00:41:39     5004
      SUBROUTINE INTQP2(P, RES)
C
C     Procedure d'integration de polynomes de 2 variables sur le
C         domaine (-1., 1.)x(-1., 1.)
C     Entree : P(30) : premier polynome : les coefficients sont
C                       ranges comme suit :
C               indice : 1 2  3    4    5    6      7       8    9
C               terme :  1 T1 T2 T1*T2 T1^2 T2^2 T1*T2^2 T1^2*T2 T1^3
C               indice : 10   11      12        13      14   15
C               terme :  T2^3 T1*T2^3 T1^2*T2^2 T1^3*T2 T1^4 T2^4
C               indice : 16      17        18        19      20   21
C               terme :  T1*T2^4 T1^2*T2^3 T1^3*T2^2 T1^4*T2 T1^5 T2^5
C               indice : 22      23        24        25        26
C               terme :  T1*T2^5 T1^2*T2^4 T1^3*T2^3 T1^4*T2^2 T1^5*T2
C               indice : 27        28        29        30
C               terme :  T1^2*T2^5 T1^3*T2^4 T1^4*T2^3 T1^5*T2^2
C     Sortie : RES : REAL*8 : resultat de l'integration
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION P(30), C1(30)
      DATA C1/4.D0, 0.D0, 0.D0, 0.D0, 1.333333333333333D0,
     1        1.333333333333333D0, 0.D0, 0.D0, 0.D0, 0.D0, 0.D0,
     2        0.444444444444444D0, 0.D0, 0.8D0, 0.8D0, 0.D0, 0.D0,
     3        0.D0, 0.D0, 0.D0, 0.D0, 0.D0, 0.266666666666667D0, 0.D0,
     4        0.266666666666667D0, 0.D0, 0.D0, 0.D0, 0.D0, 0.D0/
      RES = 0.D0
      DO 30 I = 1, 30
  30     RES = RES + C1(I)*P(I)
      RETURN
      END