1 : C INTQP2 SOURCE CHAT 05/01/13 00:41:39 5004 2 : SUBROUTINE INTQP2(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/4.D0, 0.D0, 0.D0, 0.D0, 1.333333333333333D0, 24 : 1 1.333333333333333D0, 0.D0, 0.D0, 0.D0, 0.D0, 0.D0, 25 : 2 0.444444444444444D0, 0.D0, 0.8D0, 0.8D0, 0.D0, 0.D0, 26 : 3 0.D0, 0.D0, 0.D0, 0.D0, 0.D0, 0.266666666666667D0, 0.D0, 27 : 4 0.266666666666667D0, 0.D0, 0.D0, 0.D0, 0.D0, 0.D0/ 28 : RES = 0.D0 29 : DO 30 I = 1, 30 30 : 30 RES = RES + C1(I)*P(I) 31 : RETURN 32 : END 33 : 34 : 35 :
© Cast3M 2003 - All rights reserved.
Disclaimer