Esope source files

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