Esope source files

   1 : C INT6P2    SOURCE    CHAT      05/01/13    00:39:32     5004           
   2 :       SUBROUTINE INT6P2(P, RES)
   3 : C
   4 : C     Procedure d'integration de polynomes de 2 variables sur le
   5 : C         domaine (-1., 1.)x(-1., 1.)
   6 : C     Entree : P(30) : premier polynome : les coefficients sont
   7 : C                       ranges comme suit :
   8 : C               indice : 1 2  3    4    5    6      7       8    9
   9 : C               terme :  1 T1 T2 T1*T2 T1^2 T2^2 T1*T2^2 T1^2*T2 T1^3
  10 : C               indice : 10   11      12        13      14   15  
  11 : C               terme :  T2^3 T1*T2^3 T1^2*T2^2 T1^3*T2 T1^4 T2^4
  12 : C               indice : 16      17        18        19      20   21
  13 : C               terme :  T1*T2^4 T1^2*T2^3 T1^3*T2^2 T1^4*T2 T1^5 T2^5
  14 : C               indice : 22      23        24        25        26
  15 : C               terme :  T1*T2^5 T1^2*T2^4 T1^3*T2^3 T1^4*T2^2 T1^5*T2
  16 : C               indice : 27        28        29        30
  17 : C               terme :  T1^2*T2^5 T1^3*T2^4 T1^4*T2^3 T1^5*T2^2
  18 : C     Sortie : RES : REAL*8 : resultat de l'integration
  19 : C
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8(A-H,O-Z)
  22 :       DIMENSION P(30), C1(30)
  23 :       DATA C1/0.5, 0.166666666666667D0, 0.166666666666667D0,
  24 :      1  0.041666666666667D0, 0.083333333333333D0, 0.083333333333333D0,
  25 :      2  0.016666666666667D0, 0.016666666666667D0, 0.05D0, 0.05D0,
  26 :      3  0.008333333333333D0, 0.005555555555556D0, 0.008333333333333D0,
  27 :      4  0.033333333333333D0, 0.033333333333333D0, 0.004761904761904D0,
  28 :      5  0.002380952380952D0, 0.002380952380952D0, 0.004761904761904D0,
  29 :      6  0.023809523809524D0, 0.023809523809524D0, 0.002976190476D0,
  30 :      7  0.001190476190D0, 0.0008928571429D0, 0.001190476190D0,
  31 :      8  0.002976190476D0, 0.0006613756614D0, 0.0003968253968D0,
  32 :      9  0.0003968253968D0, 0.0006613756614D0/ 
  33 :       RES = 0.D0
  34 :       DO 30 I = 1, 30
  35 :   30     RES = RES + C1(I)*P(I)
  36 :       RETURN
  37 :       END
  38 : 
  39 :  
  40 :